Bivariate Step-Stress Accelerated Life Tests for the Kavya–Manoharan Exponentiated Weibull Model under Progressive Censoring with Applications

: A new three-parameter survival model is proposed using the Kavya–Manoharan (KM) transformation family and the exponentiated Weibull (EW) distribution. The shapes of the pdf for the new model can be asymmetric and symmetric shapes, such as unimodal, decreasing, right-skewed and symmetric. In addition, the shapes of the hrf for the suggested model can be increasing, decreasing, constant and J-shaped. Statistical properties are obtained: quantile function, mode, moments, incomplete moments, residual life time, reversed residual life time, probability weighted moments, order statistics and entropy. We discuss the maximum likelihood estimation for the model. The relevance and ﬂexibility of the model are demonstrated using two real datasets. The distribution is very ﬂexible, and it outperforms many known distributions, such as the three-parameter exponentiated Weibull, the modiﬁed Weibull model, the Kavya–Manoharan Weibull, the extended Weibull, the odd Weibull inverse Topp–Leone and the extended odd Weibull inverse Nadarajah–Haghigh model. A bivariate step-stress accelerated life test based on progressive type-I censoring (PTIC) using the model is presented. This pattern is noticed when a particular number of lifetime test units are routinely eliminated from the test at the conclusion of each post-test period of time. Minimizing the asymptotic variance of the MLE of the log of the scale parameter at design stress under PTIC yields an expression for the ideal test plan under PTIC.


Introduction
In order to enhance survival data modeling, statisticians and applied researchers have been increasingly interested in developing flexible lifetime models. As a consequence, significant progress has been achieved in generalizing and applying a number of wellknown lifetime models. The Weibull distribution, introduced by Weibull [1], is one of the most used distributions for modeling lifetimes. The Weibull has been utilized effectively as a purely empirical model in many applications due to its flexible form and ability to predict a wide variety of failure rates. The Weibull model may be theoretically developed as a type of extreme value distribution that governs the time until occurrence of the "weakest link" among several competing failure processes. This might account for its performance in applications such as capacitor, ball bearing, relay and material strength failures. The hazard rate function (hrf) of the Weibull distribution can only be increasing, decreasing or constant. As a result, it is not useful for modeling human mortality or machine lifetimes with a non-monotonic hrf. Many variants that generalize the Weibull distribution have emerged. Mudholkar et al. [2] were the first to extend the Weibull distribution with a bathtub-shaped hazard rate, and this extended Weibull has been successfully applied by Xie et al. [3]. Other extensions and modified forms of the Weibull with greater flexibility have been proposed; for example, the additive Weibull distribution with a bathtub-shaped hrf by Xie and Lai [4], the EW distribution by Nadarajah et al. [5], the transmuted additive Weibull by Elbatal et al. [6], the Kumaraswamy transmuted exponentiated modified Weibull by Al-Babtain et al. [7], the Burr X EW model by Khalil et al. [8], the Marshall-Olkin powergeneralized Weibull distribution by Afify et al. [9], truncated Cauchy power Weibull-G by Alotaibi et al. [10], a new version of Weighted Weibull distribution by Alahmadi et al. [11], exponentiated power generalized Weibull power series family by Aldahlan et al. [12], exponentiated truncated inverse Weibull-G by Almarashi et al. [13], odd inverse power generalized Weibull by Al-Moisheer et al. [14], extended inverse Weibull distribution by Alkarni et al. [15], four-parameter Weibull distribution by Abouelmagd et al. [16], exponentiated Weibull Weibull distribution by Hassan and Elgarhy [17] among others.
Mudholkar and Srivastava [18] presented the three-parameter EW distribution as an extension of the Weibull family, which is a particularly flexible class of probability distribution functions. Mudholkar et al. [2] demonstrated the use of the EW distribution in reliability and survival tests.
The cumulative function (cdf) and density function (pdf) of the EW model are given by and g EW (x; α, β, θ) = βθα β x β−1 e −(αx) β 1 − e −(αx) β θ−1 , where β and θ are two shape parameters and α is the scale parameter. For β = 1, it represents the exponentiated exponential (EE) distribution, for β = 2, it represents the exponentiated Rayleigh (ER) distribution, for β = θ = 1, it represents the exponential (E) distribution, for β = 2, θ = 1, it represents the Rayleigh (R) distribution [19] and for θ = 1, it represents the Weibull distribution. Thus, EW is a generalization of the EE distribution, as well as the Weibull distribution. The EW distribution also has a very good physical interpretation. We have three reasons for using the EW distribution. First, it expands the popular Weibull and EE models. Second, the hazard rate function of this distribution has several forms, one of which is the bathtub form. As a result, it could be a helpful distribution for analyzing biological and mortality data. Third, if a parallel system has n components and their lifetimes are independently and identically distributed as EW, then the system lifetime is also EW.
Additional parameters give greater flexibility, but they also increase the complexity of the estimation. To counter this, ref. [20] proposed the Dinesh-Umesh-Sanjay (DUS) transformation to obtain new parsimonious classes of distributions. This is as follows. If G(x) is the baseline cdf, the DUS transformation generates a new cdf F(x) expressed as The merit of using this transformation is that the resulting distribution is parameterparsimonious because no extra parameters are added. In this way, Maurya et al. [21] proposed a new class of distributions that includes many flexible hazard rates. They explored using the DUS transformation by using the exponentiated cdf, introducing the generalized DUS (GDUS) transformation. Kavya and Manoharan [22] proposed a generalized lifetime model based on the DUS transformation, with the cdf of the GDUS transformation given by where α > 0. The associated pdf is given by where G(x; ζ) is the baseline distribution in the GDUS family distribution. This approach will always create a parsimonious distribution because it is a transformation rather than a generalization, so no additional parameters beyond those in the baseline distribution are introduced.
Recently, Kavya and Manoharan [23] introduced a new transformation, the KM transformation family of distributions. The cdf and pdf are, respectively, and The hrf is provided via Using a given baseline distribution, this family generates new lifetime models or distributions.
Kavya and Manoharan [23] used the exponential and Weibull distributions as baseline distributions because they are widely used in reliability theory and survival analysis.
Censored data arise in real-world testing trials when investigations, including the lifetime of test units, must be terminated before complete observation. Censoring is a common and unavoidable routine activity for a variety of reasons, including time restrictions and cost savings. Censorship in its different forms has been carefully examined, with types I and II censorship being the most prevalent. In comparison to traditional censorship designs, a generalized censorship style known as progressive censored schemes has recently attracted substantial attention in the literature due to its efficient use of available resources. One of these progressive censored schemes is the PTIC. This pattern is noticed when a particular number of lifetime test units are routinely eliminated from the test at the conclusion of each post-test period of time. It has the capacity to determine the termination time realistically and there is additional design freedom by allowing test units to be terminated during non-terminal time periods, as studied by Balakrishnan et al. [24]. Bayesian and classical inference for an odd Lindley Burr XII model are proposed in the study by Korkmaz et al. [25] and Topp-Leone NH distribution is proposed in the study by Yousof et al. [26].
Suppose that a life testing experiment has n units. Presume that X 1 , X 2 , . . . , X n represent the lifetime of all n units taken from a population. Assume that x (1) < x (2) < . . . < x (n) denotes the respective ordered lifetime recorded from the life test. At the end of the preset period of censoring Tqi, Ri items are excluded from the surviving items in accordance with qi th Qs, i = 1, 2, . . . , m, where m signifies the number of testing phases, T q i > T q i−1 and n = r + ∑ m i=1 R i . The quantities T q i should always be determined ahead of time:

1.
Predicated on the experimenter's past knowledge and skills with the objects under consideration, as in the study by Balasooriya and lwa [27];

2.
Or the quantiles of the lifetimes distribution, qith, which are possibly computed by using the provided expression Within those cases, Ri, Tqi and n are fixed and constant, whereas li represents the quantity of surviving entities at a particular moment in time, and Tqi and r = ∑ i = 1 m li are random variables. The likelihood function is represented as is the observable lifetime of the ith OS [28]. This censorship mechanism is depicted in Figure 1; see Balakrishnan and Cramer [29]. The type-I and type-II conventional censoring methods are generalized by these progressive censoring schemes, which also allow for partial withdrawals from the ongoing experiment. The fact that they assume a constant review of the life test experiment is a significant problem with these progressive censoring schemes, as well as with all other classic censoring schemes. Complete samples, as well as the type-I censoring scheme, can indeed be viewed as special examples of this censoring technique.
The PTIC approach for the Weibull distribution was proposed in the study by Balakrishnan and Cramer [29]. The maximum likelihood estimates (MLEs) and Bayesian estimates for the unknown parameters of the generalized inverse exponential distribution under PTIC were obtained in the study by Mahmoud et al. [30]. Abo-Kasem et al. [31] discussed inferential survival analysis for inverted NH distribution under adaptive progressive hybrid censoring. For the PTIC, there are two publications that are closely linked. The first was the MLEs and asymptotic confidence estimates for the parameters of the extended inverse exponential model, predicated on the premise that there are two kinds of failures, as in the study by Mahmoud et al. [32]. Ahmad et al. [33] discussed inference on reliability estimation for a multi-component stress-strength model under power Lomax distribution. Almetwally et al. [34] obtained the best optimal plan of multi-stress-strength reliability Bayesian and non-Bayesian methods for the alpha power exponential model using progressive first failure. The second was the MLEs and Bayesian estimates for the unknown parameters of the extended inverse exponential model, as in the study by Mahmoud et al. [35]. Algarni et al. [36] investigate the statistical inference of the inverse Weibull model under PTIC. Elbatal et al. [37] investigate the Bayesian and non-Bayesian estimation of the Nadarajah-Haghighi distribution under PTIC.
Turning now to accelerated life tests (ALTs), these are a means of gathering more information in less time than would otherwise be possible by subjecting items to more stress than in usual operating conditions. Such testing can save a significant amount of time and money. The step-stress accelerated life test (SSALT) is one type of ALT; see Hakamipour [38]. During the test, the experimenter progressively increases the stress levels at pre-specified time intervals, typically starting with a stress level that is slightly above normal condition. The test continues until either the entire sample of items fails or the time limit for the duration of the test is reached and censoring occurs.
When stress(es) are applied items in only two phases, this is referred to as a simple SSALT. Meeker [39], as well as Nelson [40], are good resources for interested readers. Many different types of SSALT censorship have been used in the literature. El-Sherpieny et al. [41] introduced PTIIC for bivariate distributions based on ALT. The bulk of these studies concentrate on types I and II censoring. When items are eliminated before they fail on test, the cost of the test is lowered since these deleted samples can be utilized elsewhere or in other tests. This is referred to as progressive censorship. Most tests employ just one accelerating stress variable. It is recommended to employ many stress factors since using only one variable may result in insufficient failure data. Increasing the number of stress components would result in a better knowledge of the simultaneous effects of the stress variables, as well as more failure data.
Han [42] addressed the best design of a basic SSALT for an exponential distribution using PTIC. Li and Fard [43] investigated SSALT for two stress variables using Weibull failure times and type-I censoring. To obtain optimal hold durations, Ling et al. [44] built an SSALT for two stress variables in a type-I hybrid censoring scheme. "Bivariate SSALT" refers to an SSALT with two stress components; for more details of this scheme, see Figure 2. The study by [45] obtained the best design for a bivariate SSALT model with type-II censoring for the Gompertz distribution. Hakamipour [38] discussed the optimum design for a bivariate SSALT with a generalized exponential distribution under PTIC. Khan and Chandra [46] discussed bivariate SSALT estimation and the optimum design under PTIC. Alotaibi et al. [47] introduced the optimal design for a bivariate SSALT for alpha power exponential distribution based on PTIC samples. The test is repeated until T is attained, at which time, n 2 units fail this stage. At time T, all of the remaining surviving units c 3 = N − n 1 − c 1 − n 2 − c 2 − n 3 are removed from the test. For more explanation, see Figure 2. In the article under consideration, our primary focus lies in introducing a new lifetime model called the KMEW model as a new three-parameter lifetime model based on the KM transformation family, EW distribution. The following arguments give enough motivation to study the proposed model. We specify them as follows: (i) the new suggested distribution is very flexible and contains some distributions as sub-models; (ii) the shapes of the pdf for the new model can be asymmetric (unimodal, decreasing and right-skewed) and symmetric. In addition, the shapes of the hrf for the suggested model can be increasing, decreasing, constant and J-shaped; (iii) the new suggested model has a closed form for the quantile function and this makes the calculation of some properties, such as skewness and kurtosis, very easy, and generating random numbers from the new suggested model also becomes easy; (iv) some statistical and mathematical properties of the new suggested model are explored; (v) a maximum likelihood method of estimation is produced to estimate the parameters of the FB model; (vi) the KMEW is a good alternative to several lifetime distributions, such as the three-parameter exponentiated Weibull, the modified Weibull model, the Kavya-Manoharan 9 Weibull, the extended Weibull, the odd Weibull inverse Topp-Leone and the extended odd Weibull 10 inverse Nadarajah-Haghigh model for modeling skewed data in applications; (vii) we propose a bivariate SSALTs under PTIC using the KMEW model. For our proposed bivariate SSALT under PTIC, the best test plan is obtained by minimizing the asymptotic variance (AV) of the MLEs of the log of the scale parameter. This paper is organized as follows. In Section 2, a new lifetime model using an exponentiated Weibull distribution as the baseline distribution in the KM transformation family and its different distributions are presented. In Section 3, we demonstrate the statistical features of the KMEW model. The maximum likelihood inference for model parameters is discussed in Section 4. Its application to three real datasets is discussed in Section 5. The bivariate SSALT under the PTIC model is introduced in Section 6.

Kavya-Manoharan Exponentiated Weibull Distribution
In this section, we construct a new flexible distribution called the Kavya-Manoharan transformation exponential Weibull (KMEW) distribution by inserting (1) into (3) to obtain and the corresponding pdf is The survival function and hrf for the KMEW distribution are and The cumulative hazard rate function is given by Figure 3 shows the forms of the KMEW pdf (7) and hrf in (9) using various specific parameter settings. It demonstrates that the pdf of the KMEW distribution can be unimodal and right-skewed, with very heavy tails.    Figure 3. Forms of the hrf include increasing, decreasing and reversed-J form, as shown in the right of Figure 3. The KMEW distribution is a very flexible model that provides different distributions when its parameters are changed. It contains the following special models: 1.
Using the following generalized binomial we can obtain a useful linear representation for the pdf in Equation (7) as follows:

Statistical Measures
In this section, we study the statistical measures of the KMEW distribution, including quantile function, mode, moments, probability weighted moments, order statistics and entropy.

Quantile Function
Quantile functions are used in theoretical aspects, statistical applications and Monte Carlo methods. The p th quantile function of KMEW distribution is given by where p ∈ (0, 1). In addition, the quantile function allows us to define several shape measures, such as skewness and kurtosis. Furthermore, the analysis of the variability of the skewness and kurtosis on the shape parameters α, β and θ can be investigated based on quantile measures. The median can be obtained by putting p = 0.5 in (11) as

Mode
The mode of the KMEW (α, β, θ) can be obtained as a solution of the following non- It is not possible to obtain the explicit solution in the general case. It has to be obtained numerically.

Moments
Statistical moments of different orders are important in defining uncertainty characteristics of the distributions. Using the expansion (10), the r th moment of X is attained as Setting y = (j + 1)(αx) β , after some algebra, the r th moment is provided via (12) Table A1 shows the numerical values of the E(X), E(X 2 ), E(X 3 ), E(X 4 ), mode, Var (X), skewness (SK) and coefficient of variation (CV) of the KMEW distribution.
The s th incomplete moment of KMEW distribution can be written as where γ(s, t) = t 0 x s−1 e −x dx is the lower incomplete gamma function. The conditional moments KMEW distribution is defined by E(X s | X > t). It can be written as and Γ(s, t) = ∞ t x s−1 e −x dx is the upper incomplete gamma function. The momentgenerating function of the KMEW distribution can be written as The r th moment of the residual lifetime is provided via Adding (10) and the binomial expansion of (x − t) r to the previous calculation yields Similarly, the r th moment of the reversed residual life may be determined utilizing method Adding (10) and the binomial expansion of (x − t) r to the previous calculation yields .
The mean past (inactivity) time is denoted by

Probability Weighted Moments
The probability weighted moments (PWMs) of random variable X whose Cdf is F(x), say ξ r,s , is provided via Combining Equations (6) and (7), we may compute Thus, the (r, s) th PWMs of X can indeed be expressed simply by

Order Statistics
Moments of order statistics (OS) (see [48]) are very useful in quality control and reliability when the prediction of the hazard of failure of future items depends on the times of few early failures. Assume that X (1:n) ≤ · · · ≤ X (n:n) denotes the OS of a random sample X 1 , . . . , X n from a continuous population with cdf F(x) and pdf f (x); then, the pdf of k th OS, denoted by X k:n , is provided with The pdf of the k th OS for KMEW distribution can be written as where and f Ew (x; α, β, θ(j + 1)) is an EW pdf with parameters α, β and θ(j + 1).Therefore, the pdf of the KMEW OS can be expressed as a linear combination of the EW densities. Thus, according to (19), we can find some mathematical properties of X k:n from those properties of EW. The r th moment of X k:n is provided via where Z r θ(j+1) has EW distribution with parameters α, β and θ(j + 1).

Entropy
The entropy of a random variable X is a measure of variation of uncertainty and has been used in many fields, such as physics, engineering and economics. The Rényi entropy [49] is characterized by From (7), the Rényi entropy of X can be obtained as follows: using power series, and, after some simplifications, we obtain

Maximum Likelihood Estimation
In this section, we use maximum likelihood (see [50]) to estimate the parameters of the KMEW distribution. Let x 1 , . . . , x n be a random sample of size n from the KMEW distribution given by (7). The log-likelihood function (log-LLF) of the KMEW distribution is provided via In order to obtain the MLEs of parameters α, β and θ, we compute the first partial derivatives of (20) with respect to α, β and θ, and, equating them to 0, we derive the accompanying three equations as: and Then, the MLEs of the parameters α, β and θ are obtained by solving a non-linear system of equations ∂L ∂α = 0, ∂L ∂β = 0 and ∂L ∂θ = 0. The MLEs can be evaluated numerically using statistical software. For interval estimation on the model parameters, we require the observed information matrix

Application of Real Data
Three genuine data sets from engineering and medical science are used to demonstrate the importance and adaptability of the proposed KMEW model. We utilize the "maxLik" package to compute likelihood estimates in the R package using the Newton-Raphson (NR) algorithms; see Henningsen and Toomet [51]. The first dataset is the waiting times (in minutes) of 100 bank clients until the service is provided. It was originally used by [52]. The data are as follows: 0. 8 [54], we conclude that the our distribution is better than the beta-flexible Weibull distribution and many other distributions.The third dataset is remission periods (in months) of a random sample of 128 bladder cancer patients [55]. In this example, we compare the KMEW distribution with some competing distributions, such as the EW [56], modified Weibull (MW) [57], Kavya-Manoharan Weibull (KMW) [23], generalized Weibull (GW) [58], extended Weibull (ExW) [59], odd Weibull inverse Topp-Leone (OWITL) [60] and extended odd Weibull inverse Nadarajah-Haghighi (EOWINH) [61] distributions. What follows is the application of the proposed distribution to the three datasets and an observation of the findings in Tables 1-3, which show standard errors (SEs), Akaike information criterion (AINC), Bayesian information criterion (BINC), Cramér-von-Mises (CVM) [62,63], Kolmogorov-Smirnov (KS) [64,65] with its p-value (PV) and Anderson-Darling (AD) [66,67]. We find that our proposed distribution is the best model because it has the lowest values for all measurements except PV. Although the PV of KMEW is equal to the GW distribution, KMEW is better than GW because the values of AINC, BINC, CVM and AD for KMEW distribution are smaller than the values of GW. The best model is the one that has the smallest values of AINC, BINC, CVM, AD and KS, and the largest value of PV.
For waiting times data, we discussed Figures 4-8. For time between failures of secondary reactor pumps data, we discussed Figures 9-13. For bladder cancer data, we discussed Figures 14-18. Figures 5, 10 and 15 support this result, demonstrating that the estimated pdf is close to the probability histogram and the calculated cdf is close to the actual cdf. We sketched the log-likelihood for every pair of parameters, while fixing the other parameters, as shown in Figures 8, 13 and 18. We plot cdf and pdf of alternative models in Figures 6, 7, 11, 12, 16, and 17 for the three real data sets.. We find that the root of the parameter pair is a single point, which assures that the roots are unique. Figure 4, 9 and 14 show TTT plots and the estimated hazard for the three real data sets.

SSALT Based on PTIC
In this section, the bivariate SSALT under PTIC is described. In addition, the MLE of model parameters is discussed.

Model Assumptions
The bivariate SSALT under PTIC is as follows: when employing the bivariate SSALT, each stress variable (SV) has two levels. Let H s be the stress level (SL) s of variable l, where l = 1, 2 and s = 0, 1, 2. Normal operational situations are represented by the SLs H 10 and H 20 . Allow all n units of the experiment to begin at the 1 st step with SLs (K 11 , K 21 ) for a time τ 1 , during which, n 1 failures will be recorded. At time τ 1 , the 1 st SV is raised from H 11 to H 12 , c 1 units are randomly removed from the surviving N − n 1 units and the 1 st SV is raised from H 11 to H 12 . The 2 nd phase is repeated until the predefined time τ 2 is calculated at time τ 2 from the remaining N − n 1 − c 1 − n 2 units. At the end of the 2 nd step, the other SV is raised from H 21 to H 22 . The test is repeated until T is attained, at which time, n 2 units fail this stage. At time T, all of the remaining surviving units c 3 = N − n 1 − c 1 − n 2 − c 2 − n 3 are removed from the test. For more of an explanation, see Figure 2.
In the 1 st phase, the life of test units is determined by a KMEW distribution with cdf in Equation (6).
The scale parameter α i at test step i has a log-linear function of SLs for i = 1, 2, 3.
Step 3: ln(α 3 where B 0 , B 1 and B 2 are unknown parameters that differ depending on the product and the test method. It is considered that the two stresses have no relationship. The cumulative exposure model is also taken into account. The remaining life in this model is solely determined by the present cumulative failure probability and the current SL, regardless of how the probability is calculated ( [40]).
For all SLs, the shape parameter β remains constant. The cdf of the lifetimes of test units for the bivariate SSALT and cumulative exposure models is then where i = 1, 2, 3. The pdf of bivariate SSALT for this can be written as and

Likelihood Function of Bivariate SSALT model
In a bivariate SSALT, assume that x ij , i = 1, 2, 3, j = 1, 2, . . . , n i denote the observations derived from a PTIC sample with random deletions. The number of items eliminated from the test at any one moment is distributed binomially, and each unit is excluded with the same probability p. In other words, If C i is independent of x ij for all i, then the joint log-LLF of the bivariate SSALT model under the PTIC sample is as follows: where where F i (x ij ) and f i (x ij ) will be substituted for F i (x ij ) and f i (x ij ) in (24)- (27). Finally, the log-LLF is constructed by inserting L 1 (θ). Maximizing the logarithm of the likelihood function instead of the likelihood function itself is frequently easier. The log-LLFs of KMEW based on the bivariate SSALT model under the PTIC sample is Table 4 shows the SLs for lieftime data after bivariate SSALT and PTIC. When τ 1 = 8, τ 2 = 13, T = 20, Table 5 shows the ML estimators of the model. We may deduce from the data in Table 5 that the model's efficiency improves when the chance of binomial elimination increases and the AINC and BINC values decrease. Table 6 shows the SLs for the time between failures of secondary reactor pumps data after bivariate SSALT and PTIC. When τ 1 = 0.6, τ 2 = 1, T = 4.5, Table 7 shows the ML estimates of the model. We may deduce from the results in Table 7 that the model's efficiency improves when the chance of binomial elimination increases and the AINC and BINC values decrease.    The SLs for bladder cancer data under bivariate SSALT and PTIC are shown in Table 8. The ML estimators of the model when τ 1 = 6, τ 2 = 10, T = 20 are shown in Table 9. From the results in Table 9, we conclude that the efficiency of the model increases when the probability of binomial removal increases, where the AINC and BINC values are smaller. In addition, we conclude that the model with p = 0.5 is better than another model. The bladder cancer data set performs quite well, as the five roots of the parameters are global maximums, as seen in Figure 19.

Concluding Remarks
In this study, we explored a new three-parameter model that is called the Kavya-Manoharan exponentiated Weibull model. Its statistical and mathematical features (moments, incomplete moments, mean deviations, Bonferroni and Lorenz curves, entropy and order statistics) were derived. The maximum likelihood estimation of the parameters was discussed. The relevance and flexibility of the KMEW model were demonstrated using two real datasets. A bivariate SSALT based on PTIC for the KMEW model was introduced. Minimizing the asymptotic variance of the MLE of the log of the scale parameter at design stress under PTIC data yields an expression of an ideal test plan under PTIC. In Tables 6 and 8, we compared the different schemes based on different values of binomial removals. From these results, we note that, if the value of binomial removals increases, then the efficiency of this model increases. For future directions, the new suggested model is very interesting to study. Many authors can estimate the parameters of the suggested model under different censored schemes. This work requires a remarkable increase in investigation, which we will delegate to future research.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: