A numerical comparison between the approximation and the estimation of the reversed hazard rate function for the inverted Kumaraswamy distribution

The unknown shape parameter of the inverse Kumaraswamy distribution and some other life-time properties such as the reversed hazard rate function (RHRF) deriving based on complete data, and imposed another shape parameter is known were introduced. Three methods for estimating the reversed hazard rate function were used. Maximum likelihood estimation (MLE) method was used for the non-Bayesian estimator, as well as using the Bayesian estimators with two informative priors (Gamma and Exponential) under symmetric (squared error loss (SEL) function) and asymmetric (entropy loss (EL) function) to estimate the reversed hazard rate function, which was approximated using a numerical method (Boubaker polynomials (BbP) method). The different estimators were compared experimentally by adopting the Monte Carlo simulation study, also comparing the estimators and the approximation for the reversed hazard rate function depending on the integrated mean squared error (IMSE). The program (MATLAB 2015) was used to obtain the numerical calculations.