Bayesian Estimation for Burr Distribution Type III Based on Trimmed Samples

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.


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 ISRN Applied Mathematics
Bayes estimates of the shape parameters c and k, 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 k and reliability function when c 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 k and reliability function when c is known, using a Bayesian approach based on Type II censored data.Ali-Mousa 7 obtained empirical Bayes estimation of the parameter k 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 c, k and the reliability function r t 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.

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 wellknown exponential and Weibull distributions as special cases.Burr 1 has suggested this family of distributions by solving the following differential equation: where g x, y > 0, 0 ≤ y ≤ 1, x in the range over which the solution is being satisfied.By using different forms of g x, y 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 X with ln X , see Johnson et al. 12 .Similarly, Burr type XII distribution can be derived from Burr type III distribution by replacing X with 1/X , 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 c, k
The distribution of Burr c, k type III is where the parameters c > 0 and k > 0 are the shape parameters of the distribution.Its density function is The rth moment where β a, b is the standard beta function.
The expectation of the distribution is obtained as follows.
If r 1 in 2.4 we have: where The mode The median

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.

Non-Bayesian Estimation for Burr Parameters Based on Trimmed Sample
In this section, the estimation problem of Burr type III with two parameters c, k 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 r and the largest r s are random then the data collected will be x r 1 ≤ x r 2 ≤ • • • ≤ x n−s and the likelihood function in double censored type II takes the following form:

3.1
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 w n!/ r−1 !n−s ! .The maximum likelihood estimators c, k of c, k are the solutions of the system of equations obtained by letting the first partial derivatives of the total log likelihood with respect to c, k be zero.The systems of equations are as follows:

3.4
From 3.3 the maximum likelihood estimator of k is expressed by From 3.4 the maximum likelihood estimator of c is expressed by

3.6
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 c, k.

Asymptotic Variances Covariance Matrix
The asymptotic variances covariance matrix of the parameters c and k is obtained by inverting the Fisher information matrix where θ 1 , θ 2 c or k.Hence the approximate variance-covariance matrix is given by

3.8
The elements of matrix I are the negative of second derivatives of the natural logarithm of likelihood function defined in 3.2 .

ISRN Applied Mathematics
The elements of matrix I are given as follows:

3.11
where . The maximum likelihood estimators c, k have asymptotic variance-covariance matrix defined by inverting the Fisher information matrix.

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 2. Choose censored failure r 0.6, s 6 .
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 c and k defined in 3.4 and 3.6 .
Step 5. Equations 3.9 and 3.10 were used to obtain the variance covariance matrix of c, k .
Results are tabulated in Table 1.
From Table 1, we note that the standard deviation decreases when n is increasing.Also, we note that, when c and k 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 c and k at 1.2 and 1.5 , respectively is better than the other values.

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.

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 c, k have independent prior distribution and let the noninformative prior NIP for c and k are, respectively, given by:

3.12
Consequently, the joint NIP will be defined as follows: The joint posterior density functions of c, k under double censored sample type II will be where λ is the normalized constant defined as follows:

3.15
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 c will be

3.16
Similarly integrating the joint posterior 3.14 with respect to k.Then, the marginal posterior of k is

3.17
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 c, k are expressed as follows: dc dk.

3.18
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 c, k will be obtained.The posterior mean and the posterior variance of c, k 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 c, k using 3.18 .
Step 4. The posterior mean and the posterior variance of the estimators for the shape parameters c, k for all sample size and for sets of parameter values were obtained.Numerical results are summarized in Table 2.

Bayes Estimator under LINEX Loss Function
Under the assumption that minimal loss function occurs θ θ, the LINEX loss function for θ θ c, k can be expressed as defined by Zellner 14 , the Bayes estimator θ of θ under LINEX loss function is

3.21
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 c, k .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 c, k in 3.20 and 3.21 .
Step 4. The posterior mean and the posterior variance of the estimators for the shape parameter c, k for all sample size and for sets of parameters were obtained.
Numerical results are summarized in Table 3.

Bayes Estimator under General Entropy (GE) Loss Function
The Bayes estimator θ under GE loss function is

3.22
Now, in 3.22 put θ c, then the Bayes estimate c of parameter c relative to the GE loss function is

3.23
In 3.22 , set θ k, then the Bayes estimate k of parameter k 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.
Step 3. Solving the nonlinear Bayesian for posterior variance of the shape parameter c, k in 3.23 and 3.24 .
Step 4. The posterior mean and the posterior variance of the estimators for the shape parameter c, k for all sample size and for sets of parameters were obtained.
Numerical results are summarized in Table 4.

Bayesian Estimation for Burr Distribution Type III Parameters Based on Trimmed Sample in Case of Informative Prior (IP)
Assume that the parameters c, k have independent prior distribution and let the informative prior IP for the parameter k is given by the gamma density as follow: and the prior density function of c is the exponential density as follow: Then the joint prior density of c and k, given by

3.27
The joint posterior density function of c, k using double censored sample type II will be where E is the normalized constant equal to

3.29
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 c is

ISRN Applied Mathematics
Similarly integrating the joint posterior 3.28 with respect to k the marginal posterior of k is

3.31
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 c, k are expressed as follows:

3.32
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 c, k will be obtained.
The posterior mean and the posterior variance of c, k 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 c, k in 3.32 .
Step 4. The posterior mean and the posterior variance of the estimators for the shape parameter c, k for all sample size and for sets of parameters were obtained.
Numerical results are summarized in Table 5.

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

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 U has a uniform 0, 1 random number, then Y U −1/k − 1 −1/c follows a Burr type III distribution.In double censored type II, the true parameters selected values are c 1.2, k 1.5 , c 1.2, k 2 , c 2.5, k 1.5 , and c 2.5, k 2 .

Table 1 :
The maximum likelihood estimator, the standard deviation, and covariance of the Burr distribution with two parameter under doubly censored sample when a c 1.2, k 1.5 , b c 1.2, k 2 , c c 2.5, k 1.5 , d c 2.5, k 2 .

Table 2 :
The posterior mean and posterior variance of the Burr distribution with two parameter distribution under doubly censored sample when c 1.2, k 1.5 in case of noninformative prior.
Now, in 3.19 put θ c, then the Bayes estimate c of parameter c relative to the LINEX loss function is c − 1 a ln e −ac g 1 c | k, x dc.3.20 Set θ k in 3.19 , then the Bayes estimate k of parameter k relative to the LINEX loss function is k − 1 a ln e −ak g 2 k | c, x dk.

Table 3 :
The expectation mean and variance of the Burr distribution with two parameter distribution under doubly censored sample when c 1.2, k 1.5 under LINEX loss function.

Table 4 :
The expectation mean and variance of the Burr distribution with two parameter distribution under doubly censored sample when c 1.2, k 1.5 general entropy loss function.

Table 5 :
The posterior mean and posterior variance of the Burr distribution with two parameter under doubly censored sample when c 1.2, k 1.5 in case of informative prior.