Abstract

In this study, the estimation of parameters of a three-parameter flexible reduced logarithmic-inverse Lomax (FRL-IL) distribution based on progressive type-II right censored sample is studied. These methods include maximum likelihood estimations (MLEs) and Bayesian estimators. Approximate confidence intervals (ACIs) for the reliability and hazard functions are estimated based on the asymptotic distribution of maximum likelihood estimates (MLEs). In addition, two bootstrap CIs are also proposed. Bayesian estimates are obtained for symmetric and asymmetric loss functions such as squared error loss (SEL) and linear-exponential (LINEX) loss functions. The Gibbs within Metropolis-Hasting sampler procedure is applied using the Markov Chain Monte Carlo (MCMC) technique to get the Bayes estimates of the unknown parameters and their credible intervals (CRIs). Finally, a real-life dataset that represents a group of patients with bladder cancer is considered an application of the proposed methods.

1. Introduction

In life testing and reliability experiments, experiments must be often terminated before all units on the test have failed. For example, units may accidentally break in an industrial experiment, individuals may drop out of the study in a clinical trial, or the study may have to be terminated for lack of funds. For many reasons, the removal of units before failure is very often a procedure due to limitations of time and cost associated with the experiment. In such cases, in a lifetime test, it is impossible for experimenters to obtain the complete sample information in a short period of time. It leads us to the area of censoring, namely, some surviving units are removed from the experiment following a specific censoring scheme, and the final sample obtained is the censored sample, see Klein and Moeschberger [1]. Reducing the total test time and the associated test is one of the major reasons for censoring. In the life-testing experiment, if all the items under the test are observed until failure, the life test is said to be a complete life test. The available data in most practical situations are not complete. Therefore, censoring is very important in lifetime data analysis, either due to test time or budget constraints. Well-known censoring schemes are divided into two, type-I and type-II, censors. And then, a progressive censoring scheme is further proposed. This, in turn, allows the experimenters to remove the remaining units at different stages during the test, which is more flexible and realistic. The progressive type-II censoring scheme can be demonstrated as follows: the tester assigns independent and identical units on the life test. Suppose there are units to be tested at the beginning of the experiment, the lifetime test is terminated at the failure time of the mth unit. Once the first failure occurs, time is registered, and units are randomly withdrawn from the remaining survival units. Consequently, time will be registered and units are randomly withdrawn from the remaining survival units at the occurrence of the second failure. This experiment stops at the failure, which is determined in advance, at time , and (Figure 1).

There are a lot of authors who have studied progressive type-II censoring with different breakdowns of failure time, including Mann et al. [2], Lawless [3], Meeker and Escobar [4], and Balakrishnan [5] studied the properties of progressively censored order statistics and provided an overview of various developments in inferential procedures based on progressive type-I and type-II right censored samples and identified some interesting possible problems for further research. Balakrishnan and Cramer [6] provide a comprehensive literature review on progressive censoring, details on this progressive censoring scheme, and its different applications. Tse et al. [7] indicated that, for example, the number of patients who drop out of a clinical test at each stage is random and cannot be predetermined. In some reliability trials, the experimenter may decide that it is inappropriate or dangerous to perform the test on some of the units tested even though these units did not fail. In such instances, the pattern of removal is random with each failure. For further information on the progressively censored samples, see Balakrishnan and Sandhu [8] and Balakrishnan and Aggarwala [9]. Aggarwala and Balakrishnan [10] have discussed inference for the case of progressive type II censored when lifetime distributions are Weibull, log normal, and exponential. Balakrishnan and Sandhu [11], and Aggarwala and Balakrishnan [12] have developed an algorithm to simulate general progressive type-II censored samples from uniform or any other continuous distributions. An interesting real application of progressive type II censored has been carried out by Montanari et al [13], and Eryilmaz and Bairamov [3] and also has been studying the estimation of parameters from different lifetime distributions based on progressive type-II censoring by many authors, including Balakrishnan and Kannan [14], Ali Mousa and Jaheen [15, 16], and Ali Mousa and Al-sagheer [17]. Recently, Salemi et al. [18] studied A-optimal and D-optimal censoring plans in progressive type-II right censored order statistics, Qin et al. [19] propose a new test statistic based on spacings to test whether the general progressive type-II censored samples are from an exponential distribution, Maiti and Kayal [20] studied estimation of the unknown parameters, reliability, and hazard functions of a log-logistic distribution under the progressive type-II censored sample, also, a study of the generalized Rayleigh distribution was presented when a controlled stepwise sample of the second type was available by Maiti and Kayal [21], and parametric (classical and the Bayes) point estimation procedures have been proposed to estimate the reliability characteristics, viz, reliability function and mean time to failure for XGD based on the complete observations, using four different estimation methods proposed by Saha and Yadav [22].

Bladder cancer is one of the most prevalent malignant tumors and one of the most expensive to treat per patient. Bladder cancer is chronically under-recognized as a public health concern and scientifically underfunded, despite its high prevalence, morbidity, death, and related management costs. In the last 30 years, 5-year survival rates for patients with prostate and kidney cancer have dramatically improved, but progress in bladder cancer has slowed a renewed interest from the clinical and research communities. Improved diagnostics, therapies, and health services for patients with bladder cancer will emerge as awareness and funding for the disease growth. For this reason, since our study will be on patients with bladder cancer, which requires us to reduce the number of experimental patients and the time of the lifetime experiment of these patients, we will apply the progressive type-II censoring scheme, which satisfies the requirements of obtaining good estimation by reducing the experimental time of patients and keeping some experimental patients from failure.

In this study, the three-parameter flexible reduced logarithmic-inverse Lomax (FRL-IL) distribution proposed by Buzaridah et al. [23] is considered as a good model for the failure times of the bladder cancer patient components. The random variable X has a (FRL-IL) flexible reduced logarithmic-inverse Lomax distribution if its probability density function (PDF) and cumulative distribution function (CDF) are given bywhile the reliability and hazard function of FRL-IL, respectively, are given by

Buzaridah et al. [23] discussed the statistical properties for the FRL-IL, including the behavior of the function of probability density, the reliability or survival function, the function of the hazard rate, the function of the reversed hazard rate, and the residual (reversed) life. The moments and the moments generating function, quantile function, and skewness and kurtosis are obtained. In addition, they investigated the maximum likelihood estimators of the unknown parameters and their asymptotic confidence intervals based on complete data.

The rest of this study is organized as follows: In Section 2, the MLEs of , , , and are formulated. Asymptotic confidence intervals based on the maximum likelihood estimates are presented in Section 3. The confidence intervals for the unknown parameters, reliability function, and hazard function are constructed using two parametric bootstrap procedures in Section 4. Bayes estimates for the aforementioned parameters and functions are also obtained for different loss functions such as squared error loss (SEL) and LINEX loss functions in Section 5. The real dataset has been analyzed in Section 6. A simulation study is performed in Section 7 to assess the quality of the different estimators developed in this study. Finally, we conclude the study in Section 8.

2. Maximum Likelihood Estimators

In this section, we study the estimation problem of the FRL-IL parameters under progressive type-II censored samples using an estimation method called the maximum likelihood estimators (MLEs).

Let denote the progressive type-II censored sample of size from a sample of with scheme , drawn from the FRL-IL distribution. The likelihood function based on a progressive type-II censored sample from FRL-IL is given bywhere is a constant defined as .

and are cdf and pdf defined in equations (1) and (2), respectively.

Substituting cdf and pdf into equation (4), then, we can write the log-likelihood function of , and based on progressive type-II censored sample as follows:

The log-likelihood function without the constant term can be written as

Calculating the first partial derivatives of with respect to , , and and equating each to zero, we get the likelihood equations as

Since equations (7)–(9) do not have closed-form solutions, the Newton-Raphson iteration method is used to obtain the estimates. Ahmed [24] and Mahmoud et al. [25] describe the steps in which the algorithm works in detail.

The algorithm is described as follows:(1)Use the method of moments or any other methods to estimate the parameters , , and as starting point of iteration, denote the estimates as (, , ), and set .(2)Calculate and the observed Fisher information matrix , given in the next section.(3)Update as(4)Set and then go back to Step 1.(5)Continue the iterative steps until is smaller than a threshold value. The final estimates of are the MLE of the parameters denoted as .

Moreover, the invariance property of MLEs, and the MLEs of and can be obtained after replacing , , and by , , and as

3. Asymptotic Confidence Intervals

The asymptotic variances and covariances of the MLEs, , and are given by the entries of the inverse of the Fisher information matrix , where and . Unfortunately, exact closed forms for the expectations listed above are difficult to come by. Therefore, the observed Fisher information matrix , which is obtained by dropping the expectation operator , will be used to construct confidence intervals (CIs) for the parameters. As a result, the observed information matrix is given by

Therefore, the approximate asymptotic variance-covariance matrix for the MLEs is obtained by inverting the observed information matrix ; therefore,

It is well known that under some regularity conditions (see Lawless [26]), is approximately distributed as multivariate normal with mean and covariance matrix . Then, the two-sided confidence intervals of , and can be given bywhere is the percentile of the standard normal distribution with right-tail probability .

To construct the asymptotic confidence interval of the reliability and hazard functions, which are functions in the parameters , and , it is necessary to compute their variances. The delta method is used to compute the approximate estimates of the variance of and (for more details about the delta method, see Greene [27], Ahmed [24], and El-Sagheer [28]). According to this method, the variance of and can be approximated, respectively, bywhere and are, respectively, the gradient of and with respect to , and .

Then, the two-sided confidence intervals of and can be given by

4. Bootstrap Confidence Intervals

A parametric bootstrap interval, as opposed to a point estimate, substantially provides more information about the population value of the quantity of interest. In addition, it is clear that confidence intervals based on asymptotic results perform poorly for a small sample size. Two parametric bootstrap algorithms are offered to calculate the bootstrap confidence intervals of and . The first one is the percentile bootstrap (Boot-p) confidence interval, proposed by Efron [29]. The second one is the bootstrap-t (Boot-t) confidence interval, which is presented in Hall [30]. Boot-t was created using a studentized “pivot” and requires an estimator of the variance of the MLE of and .

4.1. Parametric Boot-p

(1)Based on the original data , , and are obtain by maximizing equations (7)–(9).(2)Based on the prespecified progressive censoring scheme , a type-II progressive censoring sample from the FRL-IL distribution with parameters , and are generated, using the algorithm described in Balakrishnan and Sandhu [11].(3)The MLEs are obtained based on the bootstrap sample, and this bootstrap estimate is denoted by (in our case could) be and .(4)Steps (2) and (3) are repeated boot times, and is obtained, where , .(5) is arranged in ascending orders, and is obtained.

Let be the cumulative distribution function of . is defined for given . The approximate bootstrap-p CI of is given by

4.2. Parametric Boot-t

(1)–(3) The same as the parametric Boot-p.(4) Based on the asymptotic variance-covariance matrix (19) and delta method (21), respectively, the variance-covariance matrix and the approximate estimates of the variance and are computed.(5) The statistic is computed, which is defined as(6) Steps 2–5 are repeated NBoot times, and is obtained.(7) is sorted in ascending orders, and the ordered sequences are obtained.

Let be the cumulative distribution function of for given . is defined.

Then, the approximate bootstrap-t CI of is given by

5. Bayes Estimation

In this section, we obtain Bayesian estimates that deal with the parameters as random and uncertainties in the parameters are described by a joint prior distribution, which was developed before the collected failure data. The Bayesian approach is highly useful in reliability analysis because it may incorporate previous knowledge into the analysis. This is important because one of the main challenges associated with reliability analysis is the limited availability of data. Bayesian estimates of the unknown parameters , , and and some lifetime parameters and against the SEL and LINEX loss functions are developed. It is assumed here that the parameters , , and are independent and follow the gamma prior distributions,where all the hyperparameters and , are assumed to be non-negative and known.

The posterior distribution of the parameters , and denoted by can be obtained by combining the likelihood function in equation (4) with the priors in equation (23), and it can be written as

A commonly used loss function is the SEL, which is a symmetrical loss function that assigns equal losses to overestimation and underestimation. If is the parameter to be estimated by an estimator , then the square error loss function is defined as

Therefore, the Bayes estimate of any function of , , and , say under the SEL function, can be obtained aswhere

Variance [31] considered a LINEX (linear exponential) loss function for a parameter is given by

This loss function is suitable for situations where overestimation is more costly than its underestimation. Zellner [32] discussed Bayesian estimation and prediction using LINEX loss. Hence, under LINEX loss function in equation (30), the Bayes estimate of a function iswherewhere the calculation of the multiple integrals in equations (30) and (33) cannot be analytically solved. Thus, the MCMC technique is used to generate samples from the joint posterior density function in equation (27). In order to be able to implement the MCMC technique, we consider the Gibbs within the Metropolis-Hasting sampler procedure. The M-H and Gibbs sampling are two useful MCMC methods that have been widely used in statistics. These methods recently have a high degree of discussions. (for details, see Ahmed [33], Panahi and Asadi [34], and Abushal [35]).

From equation (27), the joint posterior distribution can be written as follows:

The posterior conditionals for , , and are

But, the conditional posteriors of , and in equations (35)–(37) do not present standard forms, so Gibbs sampling is not a straightforward option. Therefore, the use of the Metropolis-Hasting sampler is required for the implementation of MCMC technique. Within Gibbs sampling, the Metropolis-Hasting algorithm is as follows:(1)Start with an initial guess .(2)Set .(3)Using the following M-H algorithm, generate , and from , and with the normal proposal distributionswhere , and can be obtained from the main diagonal in inverse Fisher information matrix in equation (12).(4)Generate a proposal from from , and from .(i)Evaluate the acceptance probabilities(ii)Generate a , , and from a uniform distribution.(iii)If , accept the proposal and set , else set .(iv)If , accept the proposal and set , else set .(v)If , accept the proposal and set , else set .(5)Compute the reliability function and hazard function as(6)Set .(7)Repeat steps (3)–(6) times and obtain and .(8)To compute the CRs of , , , and as , then the CRIs of is

The first simulated variants are discarded in order to ensure convergence, and the affection of initial value selection is removed. Then, the selected samples are , for sufficiently large .

Based on SEL function, the approximate Bayes estimates of are given byThe approximate Bayes estimates for , under LINEX loss function, from equation (33), are

6. Applications

In this section, the proposed estimation methods are applied to the failure data for a group of 30 patients with bladder cancer. Bladder cancer is where a growth of abnormal tissue, known as a tumor, develops in the bladder lining. The tumor may migrate into the bladder muscle in some cases (see Lee and Wang [36]). The Kolmogorov-Smirnov (K-S) distance between the empirical distribution of failure data and the CDF of the FRL-IL distribution is 0.111778 with -value that equals 0.0755062. Hence, the FRL-IL distribution fits well to the given data. A progressive type-II censored sample of effective size was randomly selected from 30 failure observations with a progressive censored scheme .

The MLEs of parameters, reliability, and hazard functions based on progressive type-II data presented in Table 1 are obtained to be and at and are displayed in Table 2. Using the algorithms described in Section 4 of the bootstrap methods, the mean of Boot-p (Bp) and Boot-t (Bt) samples of the lifetime parameters is also displayed in Table 2. The Bayes estimates relative to both SEL and LINEX functions are computed, for different values of the shape parameter of LINEX loss function for the parameters , , and and the reliability and hazard functions at , and are also displayed in Table 2.

The ACIs and CRIs for the parameters , , and are computed. Moreover, the delta method is used to obtain the ACIs and CRIs for the reliability and hazard functions, and the results are displayed in Table 3. Also, the 95% bootstrap (Boot-p and Boot-t) confidence intervals (CIs) are displayed in Table 3.

It is known that the LINEX loss function becomes symmetric for tending to zero. From Table 3, the results of SEL and LINEX loss functions, at , are the same, indicating that the recommended approaches are accurate.

7. Simulation Study

In order to compare the estimators of parameters and some lifetime parameters’ reliability function and hazard function of the FRL-IL distribution, a simulation study was performed utilizing 1,000 progressive type-II censored samples for each simulation. To generate progressive type-II censored samples from the FRL-IL distribution, with initial values , and , the comparison of the different approaches of the estimators of the results and , at , has been considered in their mean square error (MSE), which is computed for and , where is the number of simulated samples. Another criterion is used to compare the CIs obtained by using asymptotic distributions of the MLEs and CRIs. The comparison between them is made in terms of the average confidence interval lengths (ACLs) and coverage probability (CP). A confidence interval’s CP is the percentage of time the interval contains the starting value of interest. The following progressive schemes are studied in this study:Scheme I: for .Scheme II: for and .Scheme III: for .

The results of estimate parameters and their MSE are shown in Tables 48, and the results of ACL and CP of CIs are shown in Tables 9 and 10.

From the results, the following conclusions can be drawn:(1)It is observed that from Tables 48, as sample size increases, the MSEs decrease and Bayes estimates have the smallest MSEs for , , and . Hence, Bayes estimates perform better than the MLE methods in all cases considered.(2)Bayes estimates under LINEX with provide better estimates in the sense of having smaller MSEs.(3)For fixed values of the sample n and effective sizes , scheme I performs better than schemes II and III in the sense of having smaller MSEs.(4)From Tables 9 and 10 for various sample sizes, observed failures, and schemes, it can be seen that the CRIs provide more accurate results than the ACIs.

8. Conclusions

The aim of this study is to develop various methods for estimating and deriving confidence intervals for parameters of the FRL-IL distribution under progressive type-II samples. The different estimation methods proposed in this study were used to calculate the time period (in months) to recovery for some patients with bladder cancer. Using asymptotic distributions and parametric bootstrap methods, the MLEs of the unknown parameters are obtained and different confidence intervals proposed. On the other hand, the estimation results were used to calculate the failure risk of these components. A simulation study was conducted to evaluate the quality of the proposed estimations, and this study showed that Bayesian methods have good performance in all different cases. In addition, scheme I performs better than schemes II and III in the sense of having smaller MSEs. Since the distribution of FRL-IL is well suited to the applied data, this study can effectively contribute to the study of failure analysis of patients with bladder cancer.

In future works, we are going to do an estimation of the parameters of the flexible reduced logarithmic-inverse Lomax distribution under adaptive progressive type-II censored data, and also, we will make the comparison to all types of censored algorithms that we applied.

Data Availability

The bladder cancer data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.