ISRN Applied Mathematics

Volume 2012 (2012), Article ID 250393, 18 pages

http://dx.doi.org/10.5402/2012/250393

## Bayesian Estimation for Burr Distribution Type III Based on Trimmed Samples

^{1}Department of Mathematical Statistics, Institute of Statistical Studies and Research, Cairo University, 5 Dr. Ahmed Zewail Street, Dokki, Giza 12613, Egypt^{2}Faculty of Science, King Abdulaziz University, 21589 Jeddah, Saudi Arabia

Received 27 June 2012; Accepted 2 September 2012

Academic Editors: S.-W. Chyuan, S. W. Gong, D. Kuhl, G. Kyriacou, X. Wen, and K.-V. Yuen

Copyright © 2012 A. M. Abd-Elfattah and A. H. Alharbey. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Trimmed samples are widely employed in several areas of statistical practice, especially when some sample values at either or both extremes might have been contaminated. The problem of estimating the parameters of Burr distribution type III based on a trimmed samples and prior information will be considered. In this paper, we study the estimation of unknown parameters based on doubly censored type II. The problem discussed using maximum likelihood method and Bayesian approach to estimate the shape parameters of Burr type III distribution. The numerical illustration requires solving nonlinear equations, therefore, MathCAD (2001) statistical package used to asses these effects numerically.

#### 1. Introduction

Burr [1] introduced a family of twelve cumulative distribution functions for modeling lifetime data. The two important members of the family are Burr types III and XII. The two important distributions, Burr type III and Burr type XII, are interrelated through simple transformation. Burr type III distribution allows for a wider region for the skewness and kurtosis plane, which covers several distributions including the log-logistic, and the Weibull and Burr type XII distributions.

However, outliers may occur in the data set. Trimmed samples are widely employed in several areas of statistical practice, especially when some sample values at either or both extremes might have been contaminated. The problem of estimating the parameters of Burr distribution type III based on a trimmed sample and prior information is considered.

Many authors discuss different methods of estimation for Burr type XII distribution. Al-Hussaini and Jaheen [2] and Al-Hussaini et al. [3] used different techniques for obtaining Bayes estimates of the shape parameters and , reliability and failure rate functions based on type II censored samples. Al-Hussaini et al. [4] obtained the maximum likelihood, uniformly minimum variance unbiased, Bayes and empirical Bayes estimators for the parameter and reliability function when is known. Wingo [5] developing the theory for the ML point estimation of the parameters of the Burr distribution when Type II singly censored sample is at hand. Ali-Mousa and Jaheen [6] obtained interval estimates of the parameter and reliability function when is known, using a Bayesian approach based on Type II censored data. Ali-Mousa [7] obtained empirical Bayes estimation of the parameter and the reliability function based on accelerated Type II censored data. Gupta et al. [8] discusses analysis of failure time data by Burr distribution. Wang et al. [9] obtained the maximum likelihood estimation of the Burr XII distributions parameter with censored and uncensored data. Hossain and Nath [10] deal with unweighted least squares estimation of the parameters. They compared the results with the maximum likelihood and maximum product of spacing methods. Ali-Mousa and Jaheen [11] obtained the maximum likelihood and Bayes estimates for two parameters , and the reliability function of the Burr Type XII distribution based on progressive type II censored samples.

The objective of this paper is to obtain the estimators of the unknown shape parameters of Burr type III based on doubly censored type II. The problem discussed using maximum likelihood method and Bayesian approach to estimate the shape parameters of Burr type III distribution requires solving nonlinear equations, therefore, numerical study is carried out to asses these effects using MathCAD (2001) statistical package.

#### 2. Burr Type III Distribution

The Burr family of distributions has, in recent years, assumed an important position in the field of life testing because of its uses to fit business failure data, also it includes the well-known exponential and Weibull distributions as special cases. Burr [1] has suggested this family of distributions by solving the following differential equation: where , , in the range over which the solution is being satisfied. By using different forms of and Pearson systems, Burr obtained twelve distribution functions which listed in Burr [1].

Burr family of distributions includes twelve types of cumulative distribution functions, which yield a variety of density shapes. Many standard theoretical distributions, including the Weibull, exponential logistic, generalized logistic, Gompertz, normal, extreme value, and uniform distributions are special cases or limiting cases of Burr family of distributions. The simple closed form of these distributions has been applied to studies in simulation. Burr [1] developed the family of Burr distributions as an outgrowth into methods for fitting cumulative frequency functions rather than probability density functions to frequency data. Some of the forms of Burr distributions are related by simple transformation For example, the Burr type III distribution can be obtained from Burr type II distribution by replacing with , see Johnson et al. [12]. Similarly, Burr type XII distribution can be derived from Burr type III distribution by replacing with , see Burr and Cislak [13]. The Burr type I family is more commonly known as the uniform distribution. The Burr III, IV, V, IX, and XII families have a variety of density shapes. Types III and XII are the simplest functionally and therefore, the two distributions are the most desirable for statistical modeling. The Burr type III distribution is much richer than Burr type XII distribution.

The twelve distributions of Burr are listed in Burr [1] and Johnson et al. [12].

*Properties of Type III Burr *

The distribution of Burr type III is
where the parameters and are the shape parameters of the distribution. Its density function is
The th moment
where is the standard beta function.

The expectation of the distribution is obtained as follows.

If in (2.4) we have: where , .

The variance

The mode

The median

#### 3. Bayesian and Non-Bayesian Estimation Methods

In this section we will obtain the estimation of Burr parameters based on trimmed samples using Bayesian and non Bayesian methods of estimation.

##### 3.1. Non-Bayesian Estimation for Burr Parameters Based on Trimmed Sample

In this section, the estimation problem of Burr type III with two parameters () under type II double censored data (trimmed samples) are obtained. Some data may not be observed, a known number of observation in an ordered sample are missing at both ends in failure censored experiments, the observations are the smallest and the largest are random then the data collected will be and the likelihood function in double censored type II takes the following form: It is usually easier to use maximize the natural logarithm of the likelihood function rather than the likelihood function itself. Therefore, the logarithm of the likelihood function is where . The maximum likelihood estimators of are the solutions of the system of equations obtained by letting the first partial derivatives of the total log likelihood with respect to , be zero. The systems of equations are as follows: From (3.3) the maximum likelihood estimator of is expressed by From (3.4) the maximum likelihood estimator of is expressed by Since the closed form solution to nonlinear equations (3.5) and (3.6) is very hard to obtain, Newton-Raphson method is applied for solving the nonlinear equations simultaneously to obtain .

###### 3.1.1. Asymptotic Variances Covariance Matrix

The asymptotic variances covariance matrix of the parameters ( and ) is obtained by inverting the Fisher information matrix where or . Hence the approximate variance-covariance matrix is given by The elements of matrix are the negative of second derivatives of the natural logarithm of likelihood function defined in (3.2).

The elements of matrix are given as follows: where denotes . The maximum likelihood estimators , have asymptotic variance-covariance matrix defined by inverting the Fisher information matrix.

##### 3.2. Numerical Illustration

In estimation problem, it is required to study the properties of the derived expressions for the estimators theoretically or analytically. Sometimes it seems very difficult to study the properties of the estimators theoretically because of the complicated formula of the estimators. Consequently, a simulation study will be set up, treating separately the sampling distribution of the estimators. Simulation studies have been performed using MathCAD (2001) for illustrating the theoretical results of the estimation problem. The simulation procedures will described below.

*Step 1. *Generate a random sample of size 10, 20, 40, 60, 80, and 100 from Burr type III distribution. The generation of Burr type III distribution is very simple if has a uniform random number, then follows a Burr type III distribution. In double censored type II, the true parameters selected values are , , , and .

*Step 2. *Choose censored failure .

*Step 3. *For each sample and for the four sets of parameters distribution were estimated under doubly censored type II.

*Step 4. *Newten-Raphson method was used for solving the nonlinear equations for and defined in (3.4) and (3.6).

*Step 5. *Equations (3.9) and (3.10) were used to obtain the variance covariance matrix of ().

Results are tabulated in Table 1.

From Table 1, we note that the standard deviation decreases when is increasing. Also, we note that, when and increase the standard deviation of the estimators decrease. Similarly, the maximum likelihood estimator of the parameters has the same behaviors when the sample size becomes large and the properties of two parameters and at (1.2 and 1.5), respectively is better than the other values.

##### 3.3. Bayesian Analysis for Burr Distribution Type III Based on Trimmed Samples

The Bayesian approach allows both sample and prior information to be incorporated into analysis, which will improve the quality of the inferences. In this section, Bayesian estimators and posterior variance of the shapes parameters are obtained in the case of double censored type II. The prior distribution could be “none informative” such as a flat, uniform distribution, assuming equal probability for the parameter value within a realistic range, and when the prior could be more informative such as a normal distribution or other possible distribution that represents an initial assessment of what is known about the parameter before the collection of data and the analysis. Moreover, a numerical examples are given for illustration study.

###### 3.3.1. Bayesian Estimation for Burr Distribution Type III Parameters Based on Trimmed Sample in Case of Noninformative Prior

The likelihood function takes the form in (3.1). Assumed that the parameters , have independent prior distribution and let the noninformative prior (NIP) for and are, respectively, given by: Consequently, the joint (NIP) will be defined as follows: The joint posterior density functions of , under double censored sample type II will be where is the normalized constant defined as follows: Now, the marginal posterior of one parameter is obtained by integrating the joint posterior distribution with respect to the other parameter, hence the marginal posterior probability density function of will be Similarly integrating the joint posterior (3.14) with respect to . Then, the marginal posterior of is

It is well known that under a squared error loss function, the Bayes estimator of the parameter will be its posterior expectation. To obtain the posterior mean and posterior variance, a numerical integration is required. Then, the posterior mean and posterior variance of the shape parameters are expressed as follows: Equations (3.18) are very difficult to be solved exactly so that an iterative procedure is needed to solve these equations numerically. Using the statistical package, MathCAD (2001), the posterior mean and posterior variance of the shapes parameter will be obtained.

The posterior mean and the posterior variance of , are obtained numerically as described below.

Repeat Steps 1 and 2 in Section 3.2.

*Step 3. *Solve the nonlinear equations to obtain posterior variance of the shape parameter using (3.18).

*Step 4. *The posterior mean and the posterior variance of the estimators for the shape parameters for all sample size and for sets of parameter values were obtained.

Numerical results are summarized in Table 2.

###### 3.3.2. Bayes Estimator under LINEX Loss Function

Under the assumption that minimal loss function occurs , the LINEX loss function for can be expressed as defined by Zellner [14], the Bayes estimator of under LINEX loss function is

Now, in (3.19) put , then the Bayes estimate of parameter relative to the LINEX loss function is Set in (3.19), then the Bayes estimate of parameter relative to the LINEX loss function is The equations (3.20) to (3.21) are very difficult to obtain their solutions exactly, an iterative procedure is needed to solve these equations numerically using MathCAD (2001) statistical package to obtain posterior variance of shapes parameters . The numerical procedures will describe as follow.

Repeat Steps 1 and 2 in Section 3.2.

*Step 3. *Solve the nonlinear equations to obtain the posterior variance of the shape parameter in (3.20) and (3.21).

*Step 4. *The posterior mean and the posterior variance of the estimators for the shape parameter for all sample size and for sets of parameters were obtained.

Numerical results are summarized in Table 3.

###### 3.3.3. Bayes Estimator under General Entropy (GE) Loss Function

The Bayes estimator under GE loss function is Now, in (3.22) put , then the Bayes estimate of parameter relative to the GE loss function is In (3.22), set , then the Bayes estimate of parameter relative to the GE loss function is

The equations (3.23) to (3.24) cannot have exact solutions, an iterative procedure is needed to solve these equations numerically using MathCAD (2001) statistical package to obtain posterior variance of shapes parameter. The numerical procedures will be described as follow.

Repeat Steps 1 and 2 in Section 3.2.

*Step 3. * Solving the nonlinear Bayesian for posterior variance of the shape parameter in (3.23) and (3.24).

*Step 4. *The posterior mean and the posterior variance of the estimators for the shape parameter for all sample size and for sets of parameters were obtained.

Numerical results are summarized in Table 4.

##### 3.4. Bayesian Estimation for Burr Distribution Type III Parameters Based on Trimmed Sample in Case of Informative Prior (IP)

Assume that the parameters , have independent prior distribution and let the informative prior (IP) for the parameter is given by the gamma density as follow: and the prior density function of is the exponential density as follow: Then the joint prior density of and , given by The joint posterior density function of , using double censored sample type II will be where is the normalized constant equal to Now, the marginal posterior of one parameter is obtained by integrating the joint posterior distribution with respect to the other parameter, hence the marginal posterior probability density function of is Similarly integrating the joint posterior (3.28) with respect to the marginal posterior of is

Hence, under a squared error loss function, the Bayes estimator of the parameter will be its posterior expectation. To obtain the posterior mean and posterior variance a numerical integration is required. Then, the posterior mean and posterior variance of the shape parameters are expressed as follows: Nonlinear equations (3.32) are very difficult to be solved exactly an iterative procedure is needed to solve these equations numerically. By using the statistical package, MathCAD (2001), the posterior mean and posterior variance of the shape parameters will be obtained.

The posterior mean and the posterior variance of , are obtained numerically in the following procedures as follow.

Repeat Steps 1 and 2 in Section 3.2.

*Step 3. *Solve the nonlinear Bayesian for posterior variance of the shape parameter in (3.32).

*Step 4. *The posterior mean and the posterior variance of the estimators for the shape parameter for all sample size and for sets of parameters were obtained.

Numerical results are summarized in Table 5.

It is noted that the posterior mean decreases when is increasing. Similarly the posterior variance of the parameters has the same behaviors when the sample size becomes large.

#### Acknowledgment

This project was founded by the deanship of scientific research (DSR), King Abdulaziz University, Jeddah under Grant no. (124/130/1432). The authors, therefore, acknowledge with thanks DSR technical and financial support.

#### References

- I. W. Burr, “Cumulative frequency functions,”
*Annals of Mathematical Statistics*, vol. 13, pp. 215–232, 1942. View at Google Scholar - E. K. Al-Hussaini and Z. F. Jaheen, “Bayesian estimation of the parameters, reliability and failure rate functions of the Burr type XII failure model,”
*Journal of Statistical Computation and Simulation*, vol. 41, no. 1-2, pp. 31–40, 1992. View at Publisher · View at Google Scholar - E. K. AL-Hussani, M. A. M. Ali-Mousa, and Z. E. Jaheen, “Approximate Bayes estimators applied to tile Burr model,”
*Communications in Statistics*, vol. 23, no. 1, pp. 99–121, 1994. View at Google Scholar - E. K. Al-Hussaini, M. A. M. Ali-Mousa, and Z. F. Jaheen, “Estimation under the Burr type failure model based on censored data: a comparative study,”
*Test*, vol. 1, no. 1, pp. 47–60, 1992. View at Publisher · View at Google Scholar - D. R. Wingo, “Maximum likelihood estimation of Burr XII distribution parameters under type II censoring,”
*Microelectronics Reliability*, vol. 33, no. 9, pp. 1251–1257, 1993. View at Google Scholar - M. A. M. Ali Mousa and Z. F. Jaheen, “Interval estimation under the Burr type XII failure model based on censored data,”
*Journal of the Egyptian Mathematical Society*, vol. 2, pp. 67–73, 1994. View at Google Scholar - M. A. M. Ali-Mousa, “Empirical Bayes estimators for the Burr Type XII accelerated life testing model based on Type-2 censored data,”
*The Journal of Statistical Computation and Simulation*, vol. 52, pp. 95–103, 1995. View at Google Scholar - P. L. Gupta, R. C. Gupta, and S. J. Lvin, “Analysis of failure time data by Burr distribution,”
*Communications in Statistics*, vol. 25, no. 9, pp. 2013–2024, 1996. View at Publisher · View at Google Scholar - F. K. Wang, J. B. Keats, and W. J. Zimmer, “Maximum likelihood estimation of the burr XII parameters with censored and uncensored data,”
*Microelectronics Reliability*, vol. 36, no. 3, pp. 359–362, 1996. View at Publisher · View at Google Scholar · View at Scopus - A. M. Hossain and S. K. Nath, “Estimation of parameters in the presence of outliers for a Burr XII distribution,”
*Communications in Statistics*, vol. 26, no. 3, pp. 637–652, 1997. View at Google Scholar · View at Scopus - M. A. M. Ali-Mousa and Z. F. Jaheen, “Statistical inference for the Burr model based on progressively censored data,”
*Computers & Mathematics with Applications*, vol. 43, no. 10-11, pp. 1441–1449, 2002. View at Publisher · View at Google Scholar - N. L. Johnson, S. Kotz, and N. Balakrishnan,
*Continuous Univariate Distributions*, vol. 2, John Wiley & Sons, New York, NY, USA, 2nd edition, 1995. - I. W. Burr and P. J. Cislak, “On a general system of distributions: I. Its curve-shape characteristics. II The sample median,”
*Journal of the American Statistical Association*, vol. 63, pp. 627–635, 1968. View at Google Scholar - A. Zellner, “Bayesian estimation and prediction using asymmetric loss functions,”
*Journal of the American Statistical Association*, vol. 81, no. 394, pp. 446–451, 1986. View at Google Scholar