Parametric inference of Akash distribution for Type-II censoring with analyzing of relief times of patients

Abstract: In this paper, the problem of estimating the parameter of Akash distribution applied when the lifetime of the product follow Type-II censoring. The maximum likelihood estimators (MLE) are studied for estimating the unknown parameter and reliability characteristics. Approximate confidence interval for the parameter is derived under the s-normal approach to the asymptotic distribution of MLE. The Bayesian inference procedures have been developed under the usual error loss function through Lindleys technique and Metropolis-Hastings algorithm. The highest posterior density interval is developed by using Metropolis-Hastings algorithm. Finally, the performances of the different methods have been compared through a Monte Carlo simulation study. The application to set of real data is also analyzed using proposed methods.


Introduction
Firstly, Akash distribution, a continuous one parameter distribution which has been introduced by Shanker [1]. Akash distribution has many real life application in studies pertaining to medical, engineering and survival analysis. Akash distribution has been proved useful in making statistical inferences for such situations. They studied the shapes of the density moments, distribution of order statistics, Renyi entropy measure, skewness and kurtosis. In particular, the author have obtained reliability characteristics and stress-strength reliability of the proposed model. They have computed the MLE, method of moment estimate for the unknown parameter and studied their behavior numerically. Two real lifetime data sets are analysed from medical science.
Akash distribution is described respectively by the probability density function (PDF) and cumulative distribution function (CDF), take the following forms: This distribution can be easily expressed as a combination of exponential (β) and gamma (3, β) with their mixing proportions β 2 β 2 +2 and 2 β 2 +2 respectively. Consequently, the hazard rate, reliability and mean residual life functions take the following forms , t > 0. (1.5) In last few years, many authors have studied inference methods for Akash distribution. Using one parameter family of distribution such as Akash, exponential and Lindley distributions have been innovated and considered the concept of modeling of lifetime data by Shanker et al. [2]. They have also obtained some inferential problems. Many real life data sets are used to reflected its exibility over the exponential distribution. Shanker and Shukla [3] studied two-parameter Akash distribution and calculated its different statistical properties, estimation problem and application to this distribution. Shanker et al. [4] has introduced a generalized Akash distribution and studied their statistical aspects in this paper. The method of moments and the maximum likelihood method for estimating its parameters have been obtained. Also for real data, the author have fitted this distribution and compareed with other several distribution. Abushal [5] studied the classical and Bayesian property of the unknown parameters and reliability characteristic of Shanker distribution. However, two real data sets are studied and the applicability of Shanker distribution have been presented. Owing to this we can treat Akash distribution as an alternative lifetime model in reliability analysis.
In the literature, Several censoring schemes have been discussed. Even though, Type-I and Type-II censoring schemes are most popular censoring. Consider a life test where n independent units taken from a Akash distribution are placed under observation and failure time of each unit is recorded. Suppose that the test is terminated when r th , (1 ≤ r ≤ n, r is prefixed) unit fails. These observed failure times, say (X 1 , X 2 , . . . , X r ); X 1 ≤ X 2 ≤ · · · ≤ X r , is a Type-II censored sample of size r. In this censoring scheme n−r units remain unobserved and survive beyond the time of termination. In Type-II censoring the time of termination is a random variable. The likelihood function based on X 1 , X 2 , . . . , X r is given by (see, Cohen [6] for detail) The number of failure times r is fixed in Type-II censoring whereas the random observation is the endpoint X (r) . One may refer to the book by Lawless [7] for extensive literature and applications of this censoring.
The current paper is concerned with estimating the unknown parameter β, the reliability function R(t), the hazard function h(t) and the mean residual life function m(t) of an Akash distribution grounded on Type-II censored sample. Section 2, displays the estimates of the unknown parameter using maximum likelihood estimate. Approximate confidence intervals (CIs) are also derived. Section 3, considered the Bayesian estimates and Highest posterior density (HPD) interval are discussed. Furthermore, a Simulation study is presented using Monte Carlo simulations in section 4 and section 6, presents a real life data to explain all the methods discussed. At last, the conclusion shown in section 7.

Maximum likelihood estimators
Assume that n independent observed values taken of Akash distribution as presented in (1.1) are put on a test. Using the Type-II censoring, we obtained the ordered r failures. If the ordered r failures are X (1) , X (2) , . . . , X (r) then the likelihood function of β, under Type-II censored data drawn of an Akash distribution, is obtained as follows: We consider the logarithm of likelihood function (2.1) as log L(β) ∝ 3 r log β − n log(β 2 + 2) − β s + (n − r) log(β 2 x 2 (r) + 2β x (r) + β 2 + 2).
MLEs of β is a solution of Eq (2.2) accomplished by addressing the first partial derivatives of the total log-likelihood to be zero. So, we consider the equation as follows, The closed form solutions to the nonlinear Eq (2.3) is difficult to reach and a numerical method must be applied to solve these simultaneous equation for obtaining the MLEβ of β. Iterative techniques like Newton-Raphson and Broydan used to obtain the desired estimate. In particular, the nonlinear equations are solved using nleqslv package of the R statistical software. We refer to Pradhan and Kundu [8] for a detailed discussion for the initial guess and convergence of the iterative process. We have taken true parameter values as our initial guess.
Finally, We have obtained the MLE of reliability characteristic as follows: The interval estimation for the parameter that gets in the observed Fisher's information I(β) for the likelihood is obtain by Frequently, following the asymptotic variance of MLE for Akash distribution parameter is computed using observed Fisher's information i.e., The sampling distribution of (β−β) √

Var(β)
can be approximating by a standard normal distribution. That is, the approximate of a two sided 100(1 − η)% confidence intervals for β given byβ ± z η 2 Var(β), where 0 < η < 1 and z η 2 is the η 2 th percentile of standard normal distribution. The coverage probability using simulation is presented by

The Bayesian estimation
In this section, we have derived Bayes estimates for parameter β and reliability characteristics of an Akash distribution. We consider the Bayesian inference under squared error loss function which is also called the quadratic loss function. This is symmetrical loss function and it has been considered equal important. This loss function may be mathematically expressed as: Here d(µ) denotes some parametric function of µ. thus, the Bayesian estimate; denoted byd(µ) can be determined by the posterior mean of d(µ).
Suppose X (1) , X (2) , . . . , X (r) be a Type-II censored order statistics of a random sample of size n, of Akash distribution. Any conjugate and noninformative prior distribution can used if unknown parameter does not exit. Kundu and Pradhan [9] discussed in detail that Gamma distribution can accommodate variety of shapes depending upon parameter values. Thus the family of gamma distributions is highly flexible in nature and can be considered as suitable prior β. It is assumed that β has Gamma (a, b) prior i.e., Based on the prior, the posterior distribution of β is derived as Hence, x = (x (1) , x (2) , . . . , x (r) ). Therefore, the Bayes estimate of β under the square error loss function is computed as, Next, the Bayes estimate of reliability characteristic can be addressed as Generally speaking, the ratio of the above integrals of all Bayes estimates of the parameter and reliability characteristic is not to be achieved in a closed forms. Therefore, the next section focuses on employing two famouse approximation techniques.

Lindley's approximation
It may be noted that the above Bayes estimates takes a ratio of two integrals can't be reduced to a closed form. Lindley [10] introduced an alternative method to approximate these integrals into a finite values. Consequently, suppose I(x), indicates the posterior exception of β and reliability characteristic gathering with the posterior distribution, that is expressed in the following form: here g(β) is any function of β. and ρ(β) = log π(β) and l(β) is the log-likelihood. Lindley's approximation schemas, I(x) becomes here g ββ represents the second derivative of g(β) with respect to β andĝ ββ denotes the expression obtained for β =β. All other terms coming concerning I(x) can be addressed aŝ The corresponding Bayes estimate of β, we have g(β) = β, g β = 1, g ββ = 0,β l =β + ρ βσββ + 0.5σ 2 ββl βββ .
Consequently, we use the square error loss function, to compute the Bayes estimates for reliability characteristic in a similar manner.
We adopt the Metropolis-Hastings (MH) algorithm and derive Bayes estimates in the next part. Focuses on Hastings [11] and Metropolis et al. [12].

Metropolis-Hastings algorithm
Therefore, in the following algorithm, we employ the well-known Metropolis-Hastings MH algorithm gathering normal proposal distribution to generate samples from this distributions.
Step 1: Let initial value of the parameter to be β 0 and set l=1.
Step 4: We obtain a sample α with the help of uniform distribution with parameter 0 and 1.
Step 6: Repeat steps (2-5) B times and we get adequate number of replicates.
We are able to obtain sample from the posterior distribution of β by this way. Let us assume that B represent the total number of generated sample and B 0 represent the initial burn-in sample.
The Bayes estimate of β and R(t) against SELF are computed by: We get the Bayes estimates of h(t) and m(t) under square error loss function in a similar manner. By the idea of Chen and Shao [13], we compute HPD interval of β using the sample generated from the above algorithm. Order β 1 , β 2 , ..., β B , then consider the 100(1 − η)% symmetric credible interval is computed as (β 1 , β (1−η)B+1 )....(β Bη , β B ). Here, δ represent the largest integer less than or equal to δ. The shortest one in among all such credible intervals is the HPD interval.

Simulation study
The aim of the simulation is to discuss the properties of the derived estimators. Consequently, we use Monte Carlo simulation study to compare the different estimators discussed in section 2 and section 3. The true value of parameter can be taken such as β = 1.5 and samples are generated from Akash distribution under Type-II censoring. We carry out the simulation procedures for different combinations of (n,r) and Monte-Carlo simulation study of 10000 samples. For n = 30 the r is 20, 24, 28, 30 and for n = 50 the r is 40, 44, 48, 50. We calculated behavior of the MLE and Bayes estimates for different combinations of n and r. The hyper-parameter assigned for the non prior density is a = 0 and b = 0 and prior density is a = 3 and b = 2. we set Q = 10000 replications to compute the Bayes estimates and HPD interval utilizing MH algorithm.
The performances of the different estimators are compared regarding their average estimates and their mean squared errors(MSEs) ofβ,β l ,β m , based on 10000 simulated samples of Akash distribution are written in Tables 1 and 2. Thus, based on the simulated sampes the coverage probabilities and average lengths of CI and HPD interval Considering 95% of the true coverage probability obtained for comparison purposes. In Table 9, the coverage probabilities and the corresponding average lengths are computed. By the Tables 1-9, the following conclusions obtained based on calculated values of the average estimates, estimated MSE, coverage probability and average length.
From the simulation of different estimates established in Tables 1-8, it can be seen that the MLE and Bayes estimates of unknown parameters and the reliability characteristics of Akash distribution are very good in terms of minimum MSEs of serval values for n and r. The Bayes estimates using informative prior (IP) and non informative prior (NIP) are dscussed. The Bayes estimates using the MH algorithm for gamma informative prior (IP) are preferable as they include prior information than MLE. In most cases, it is noticed that the Bayes estimates using the MH algorithm have performance better than Lindley estimates on the basis of minimum MSEs. For fixed n, the MSE of all estimates become a little bit small by increasing r. That is, when sample size n is increase then expected MSE values of all estimate decrease.   Table 3. Results of simulation study of average estimates and MSEs values for several values of r and T when n = 30.  Table 4. Results of simulation study of average estimates and MSEs values for several values of r and T when n = 50.  Table 5. Results of simulation study of average estimates and MSEs values of proposed estimates of h(t) for several values of r and T when n = 30.  Table 6. Results of simulation study of average estimates and MSEs values for several values of r and T when n = 50.    In Table 9, we have calculated several 95% interval estimates of β for several values of n and r. The asymptotic intervals, informative and noninformative HPD intervals of the unknown parameter have been derived. By Table 9 we can say that informative HPD intervals have a better results compared to the other intervals in terms of the average length (AL) of these intervals. In addition, the asymptotic intervals compete good when compared with noninformative HPD intervals. Also observed that the ALs of all intervals tend to be decreased as n increases. For all intervals, the Coverage probabilities (CP) is reasonably vary about the nominal level.
Therefore, the Bayes estimate of Akash distribution of the parameter and reliability characteristics with MH algorithm are recommended.

Relief times of patients data
This section focuses on a data set from Gross and Clark [14]. The observed data of the lifetimes data relating to relief times (in minutes) of patients receiving an analgesic are 1. 4 This data is related to 20 patents and also recently considered by Shanker et al. [2]. In this paper the author fitted this real data to Akash distribution and obtain that this data is best fit for Akash distribution. They also computed mathematical and inference problems for this distribution. Based on several combinations of (n, r), we considered different Type-II censoring schemes and obtained all estimates. We Obtained MLEs and Bayes estimates with SELF using Lindley's approximation and MH algorithm of β and reliability characteristic. Because there is no prior infromation so we considered the non-infromative prior (NIP), i.e., a=b=0. In Table 10, the MLEs and Bayes estimates of β and reliability characteristic are computed. In Table 11, based on true coverage probability of 95%, we derived the approximate CI and non informative HPD intervals. We found that all the estimates the parameter β and reliability characteristic are close to each other in Table 10. By the Table 11, we can see that among all the intervals, HPD interval performs better in respect of length.

Conclusions
This paper discussed the classical and Bayesian inferential methods for the Type-II censored data from Akash distribution. Both MLE and Bayes estimates are provided. Also we have derived the CIs and HPD interval of the parameter and reliability characteristic. we calculated numerical MLE and Bayes estimators of parameter β and reliability characteristic and compared these estimates in terms of their MSE by using Monte Carlo simulation. Although the consequences have been derived under the Type-II censoring scheme, similar methods can be applied to the rest of censoring schemes. We believe that more work can be done in these contexts. Therefore, additional research is needed.
Remarks: We can use Tierney-Kadane's [15] approximation rather than Lindley's method to obtain Bayes estimates with minor mean squared errors, see Abdel-Hamid and AL-Hussaini [16].