Open Journal of Applied Sciences
Vol.09 No.08(2019), Article ID:94721,9 pages
10.4236/ojapps.2019.98054
Bayesian Estimations with Fuzzy Data to Estimation Inverse Rayleigh Scale Parameter
Shurooq Ahmed Kareem AL-Sultany
Department of Mathematics, College of Science, Mustansiriyah University, Baghdad, Iraq
Copyright © 2019 by author(s) and Scientific Research Publishing Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
Received: July 28, 2019; Accepted: August 26, 2019; Published: August 29, 2019
ABSTRACT
In this paper, Bayesian computational method is used to estimate inverse Rayleigh Scale parameter with fuzzy data. Based on imprecision data, the Bayes estimates cannot be obtained in explicit form. Therefore, we provide Tierney and Kadane’s approximation to compute the Bayes estimates of the scale parameter under Square error and Precautionary loss function using Non-informative Jefferys Prior. Also, we provide compared numerically through Monte-Carlo simulation study to obtained estimates of the scale parameter in terms of mean squared error values.
Keywords:
Inverse Rayleigh Distribution, Imprecision Data, Modified Newton Method, Tierney and Kadane’s Approximation
1. Introduction
The Rayleigh distribution (RD) is originated from two parameter Weibull distribution and it is an appropriate model for life-testing. It can be shown by transformation of random variable that if the random variable X has Rayleigh distribution,
Then the random variable has an inverse Rayleigh distribution (IRD)
[1]. The Inverse Rayleigh distribution (IRD) has been introduced by Trayer (1964) [2]. The distribution of life times of several types of experimental units can be approximated by the IRD [3]. The IRD plays an important role in many applications, including life test and reliability studies [4]. A random variable Y is said to have a one-parameter (IRD) if it has the following (PDF),
(1)
and (CDF), is given by:
(2)
where is the scale parameter.
2. Maximum Likelihood Estimators (MLE)
Given be an (i.i.d.) random vector of a random sample of size m from (IRD), the complete-data likelihood function is:
(3)
Now if is not observed precisely. Then, we can compute its probability by using Zadeh’s definition of an imprecision event [5]. The observed-data log-likelihood function can then be obtained as,
(4)
where is the Borel measurable membership function.
Now, by take the natural logarithm for the likelihood function and differentiating with respect to and then equating to zero we get:
(5)
Since, the (MLE) of is the solution of Equation (5), so, we used the modified Newton’s Method to determine the MLE of the parameter .
Where, at iteration
(6)
and
(7)
3. Bayes Estimator
In this section, we describe Bayesian method to estimate the parameter . In Bayesian opinion the parameter itself is considered as a random variable from a given probability distribution whose variability can be described by the prior distribution.
Assume that the prior distribution of the unknown scale parameter of IRD defined as using Jeffery’s prior information , which is given by [2] :
where
, a is a constant,
where
Now, the posterior density function of given imprecision data is:
(8)
In this study we consider non-informative prior density for based on square error and precautionary loss function as the following:
3.1. Bayes Estimator Based on Square Error Loss Function
Bayes estimation of any function of the scale parameter say , based on a squared error loss function, may be written as,
(9)
3.2. Bayes Estimator Based on Precautionary Loss Function
Precautionary loss function was proposed by Norstrom (1996) [6] , as follows:
,
where is an estimate of .
Bayes estimation of any function of the scale parameter say , based on a precautionary error loss function, may be written as,
(10)
Note that, Bayes estimator in (9) and (10) cannot be simplified in to a closed form. Therefore, we consider Tierney and Kadane’s approximation form to obtain Bayes estimator of of IRD.
4. Tierney and Kadane’s Approximation Form
Tierney and Kadane (1986) [7] proposed an alternative method for the evaluation of the ratio of integrals of the form (9) and (10).
Setting
(11)
(12)
Now, set
(13)
And
(14)
(15)
(16)
Now, the Equation (15) and Equation (16) can be written as
(17)
(18)
where, : is the minus the inverses of the second derivative of or at depending on what loss function have been used. : is the minus the inverses of the second derivative of at . And maximize and as well as maximize
Now, the function is given by,
(19)
where,
(20)
and that maximize , can be obtained by solving the following equation,
(21)
It is clear there is no explicit solution to Equation (21). Therefore, modified Newton method is applied to solve the required equation.
(22)
where
(23)
then,
Now, following the same argument with
4.1. Tierney and Kadane’s Approximation of λ Based on Square Error Loss Function (TKS)
Set , Equation (13) will be,
(24)
where k is a constant as in (20).
Now, that maximize in (24) can be obtained by solving following equation iteratively as in
(25)
where
And
(26)
Now, Bayes estimate of of IRD based on square error loss function, denoted by , can be obtained from Equation (17), where all the H and elements are evaluated in and respectively.
4.2. Bayes Estimate of λ Based on Precautionary Loss Function (TKP)
Set, , Equation (14) will be,
(27)
where k is a constant as in (20).
Now, that maximize in (27) can be obtained by solving following equation
iteratively as in
(28)
where
And
(29)
Now, Bayes estimate of of IRD based on prec. loss function, denoted by , can be obtained from Equation (18), where all the H and elements are evaluated in and respectively.
5. Simulation Study
In trying to illustrate and compare the methods as described above, a Monte-Carlo simulation study was perform to generate an (i.i.d) random samples, say , according to IRD through the adoption of inverse transformation method with size n = 10, 30 and 90 to take care of small, medium and large data sets. The scale parameter λ = 0.3, 0.5, 1, 1.5, 2. Then, each observation of was made Imprecision based on an appropriate selected membership function among four membership functions in the Imprecision Information System as the following Figure 1.
The simulation program has been written by using MATLAB (R2010b) program. The results of Monte-Carlo simulation have been summarized in Table 1.
The initial values required for proceeding modified Newton-Raphson method chosen to be the symmetrical rank regression estimators. The comparisons between the parameter estimates were based on values from MSE where [8] :
Figure 1. Imprecision information system.
Table 1. MSE values for estimates of the scale parameter ( ) of IRD with different cases.
n: sample size; : maximum likelihood estimate of λ by newton-raphson; : Bayes estimate of λ of IRD based on square error loss function; : Bayes estimate of λ of IRD based on prec. loss function.
(30)
: is the estimate of respectively at the jth run.
L: is the number of sample replicated chosen to be (500).
6. Conclusions and Recommendations
The most important conclusions of Monte-Carlo simulation results are:
Tierney and Kadane’s approximation based on square error loss function (TKS) estimate introduced the best perform compared with the different estimates for all sample sizes and for all cases except and , where Bayes Estimate based on Precautionary loss function (TKP) is the best.
Based on this, we recommend,
1) Using the TKS estimate to compute estimates of the scale parameter of IRD for all sample sizes and with cases , and .
2) Using the TKP estimate to compute estimates of the scale parameter of IRD for all sample sizes and with the cases and .
3) For further study, we suggest such type of work can be done by using other informative priors for the parameter of the IRD and also the parameter can be estimated by other methods.
4) Research can be applied to real data and demonstrate the importance of this distribution in practice.
Conflicts of Interest
The authors declare no conflicts of interest regarding the publication of this paper.
Cite this paper
AL-Sultany, S.A.K. (2019) Bayesian Estimations with Fuzzy Data to Estimation Inverse Rayleigh Scale Parameter. Open Journal of Applied Sciences, 9, 673-681. https://doi.org/10.4236/ojapps.2019.98054
References
- 1. Rao, G.S. and Mbwambo, S. (2019) Exponentiated Inverse Rayleigh Distribution and an Application to Coating Weights of Iron Sheets Data. Journal of Probability and Statistics, 2019, Article ID: 7519429. https://doi.org/10.1155/2019/7519429
- 2. Trayer, V.N. (1964) Doklady Acad, Nau, Belorus, U.S.S.R.
- 3. Rasheed, H.A., Ismail, S.Z. and Jabir, A.G. (2015) Acomparison of the Classical Estimators with the Bayes Estimators of One Parameter Inverse Rayleigh Distribution, International Journal of Advanced Research, 3, 738-749.
- 4. Rasheed, H.A. and Aref, R.K.H. (2016) Bayesian Approach in Estimation of Scale Parameter of Inverse Rayleigh Distribution. Mathematics and Statistics Journal, 2, 8-13.
- 5. Khoolenjani, N.B. and Shahsanaei, F. (2016) Estimating the Parameter of Exponential Distribution under Type-II Censoring from Fuzzy Data. Journal of Statistical Theory and Applications, 15, 181-195. https://doi.org/10.2991/jsta.2016.15.2.8
- 6. Norstrom, J.G. (1996) The Use of Precautionary Loss Function in Risk Analysis. IEEE Transactions on Reliability, 45, 400-403. https://doi.org/10.1109/24.536992
- 7. Tierney, L. and Kadane, J.B. (1986) Accurate Approximations for Posterior Moments and Marginal Densities. Journal of the American Statistical Association, 81, 82-86. https://doi.org/10.1080/01621459.1986.10478240
- 8. Pak, A., Parham, G.A. and Saraj, M. (2013) Inference for the Weibull Distribution Based on Fuzzy Data. Revista Colombiana de Estadistica, 36, 339-358.