Classical and Bayesian Estimation of the Inverse Weibull Distribution: Using Progressive Type-I Censoring Scheme

+e challenge of estimating the parameters for the inverse Weibull (IW) distribution employing progressive censoring TypeI (PCTI) will be addressed in this study using Bayesian and non-Bayesian procedures. To address the issue of censoring time selection, qauntiles from the IW lifetime distribution will be implemented as censoring time points for PCTI. Focusing on the censoring schemes, maximum likelihood estimators (MLEs) and asymptotic confidence intervals (ACI) for unknown parameters are constructed. Under the squared error (SEr) loss function, Bayes estimates (BEs) and concomitant maximum posterior density credible interval estimations are also produced. +e BEs are assessed using two methods: Lindley’s approximation (LiA) technique and the Metropolis-Hasting (MH) algorithm utilizing Markov Chain Monte Carlo (MCMC). +e theoretical implications of MLEs and BEs for specified schemes of PCTI samples are shown via a simulation study to compare the performance of the different suggested estimators. Finally, application of two real data sets will be employed.


Introduction
Keller and Kamath [1] were the ones to propose the IW model as a sustainable idea for describing the deterioration of structural devices in diesel engines. e IW distribution gives an excellent match to various real data sets, according to [2]. In the perspective of a mechanical system's load-strength relationship, Calabria and Pulcini [3] gave an essential explanation of this distribution. e IW distribution, which was created to explain failures of structural devices influenced by degradation phenomena, plays a critical part in reliability engineering and lifetime testing. It has been looked into from a variety of angles. On the basis of the PC type-II data set, Musleh and Helu [4] used both conventional and Bayesian estimation techniques to estimate parameters from the IW distribution. Singh et al. [5] evaluated simulated hazards of several estimators, with a focus on the Bayesian approach. De Gusmao et al. [6] and Elbatal and Muhammed [7] focused their efforts on its comprehensive versions research, which included both generalized and exponentiated generalized IW distributions. e distribution of IW has been studied from many angles. Khan et al. [8] visually and quantitatively depicted several aspects of this distribution, including mean (M), variance (V), kurtosis, and skewness. Erto [9] calculated the IW distribution using a fresh prior distribution, taking into account the shape parameter's range and the predicted value of a quantile (Q) of the sampling distribution. Sultan et al. [10] go into great depth on how to include two IW distributions into a hybrid model. e IW distribution's density function (pdf ), cumulative function (cdf ), and Q function (Qf ) are shown as follows: (2) Note that when c � 1 and c � 2 the IW model reduces to inverted exponential (IE) inverted Rayleigh (IR) models.
Censored data arises in real-life testing trials when the experiments, which include the lifetime of test units, must be stopped before acquiring complete observation. For a variety of reasons, including time constraints and cost minimization, the censoring method is frequent and inescapable in practice. Various types of censorship have been explored in depth, with Type-I and Type-II censoring being the most prevalent. In comparison to classic censoring designs, a generalized type of censoring known as PC schemes has recently garnered substantial attention in the works as a result of its efficient use of available resources. PCTI is one of these PC schemes. When a certain number of lifetime test units are continually eliminated from the test at the conclusion of each of the post periods of time, this pattern is seen [11].
Assume there are n units in a life testing experiment. Assume that X 1 , X 2 , . . . , X n indicate the lifetime of all these n units drawn from a population. e equivalent ordered lifetimes recorded from the life test are denoted by x (1) < x (2) < · · · < x (n) . Eventually R i items are omitted from the surviving items at the predetermined period of censoring T q i in accordance with q th i Qs, i � 1, 2, . . . , m, where m denotes the number of testing stages, T q i > T q i− 1 and n � r + m i�1 R i . e values T q i must be established beforehand: (1) According on the experimenter's prior knowledge and expertise with the items on test [12].
(2) e Qs of lifetimes distribution, q i th, which may be constructed using the given formula In these situations, R i , T q i , and n are preset and constant, whereas l i is the number of surviving objects at a given point in time T q i and r � m i�1 l i are random variables. e LL function is indicated by where x (i) is the observed lifetime of the i th order statistic [13]. Figure 1 describes this scheme of censoring [14].
One can observe that complete samples and also Type-I censoring scheme can be considered as special cases of this scheme of censoring.
Cohen [13] introduced PCTI scheme for the Weibull distribution. Mahmoud et al. [15] derived the MLLEs and the BEs for the parameters of the generalized IE model under PCTI. ere are two closely related papers for the PCTI. e first one was the MLEs and ACI estimates for the unknown parameters of the generalized IE model under the idea that there are two types of failures [16]. e MLLEs and BEs for the unknown parameters of the generalized IE model [17] were the second. e purpose of this paper is to look at the PCTI scheme when the lifetimes have their own IW model. We use two distinct techniques to drive the MLEs and BEs and derive the ACI of these different parameters: MCMC and Lindley Approximation. We look at a simulation outcome and a real data set to see how the various models perform in practice. e following is how the remaining of the article is structured: e MLLE and confidence intervals are discussed in Section 2. In Section 3, the Metropolis-Hasting (MH) algorithm and LiA are used to explore the Bayesian estimation technique, fully accrediting the gamma distribution as a prior distribution for unknown parameters. A simulated outcome and a real data set are utilized to demonstrate the theoretical conclusions in Section 5. Finally, there are some final observations and a summary.

Estimation Using Method of Maximum Likelihood
According on the PCTI method, for the unknown parameters of the IW distribution, the MLLE technique of estimate is examined in this section. is is how the PCTI system can be put into practice:    (1) , x (2) , . . . , x (r) ) that indicate the reported lifetime of r units under that same censoring procedure.
By applying equation (1) of IW distribution in equation (4) of LL function under PCTI, the connected LL function of c and λ given the PCTI data, x, may be interpreted as Take logarithm of ℓ(c, λ) to obtain log-LL L as First partial derivatives of log-LL function L in terms of c and λ are computed as follows: where A j � T Focused on the empirical distribution of the MLLE of the parameters, CIs for the unknown parameters c and λ will be computed. It is established from the empirical distribution of the MLLE of the parameters that where N 2 (·) is bivariate normal distribution and I(·) is the the Fisher information matrix which is defined in equation (12). Considering specific regularity constraints, the twosided 100(1 − α)%, 0 < α < 1, ACIs for the unknown parameters c and λ can be obtained as c ± Z α/2 ��� σ 11 √ and λ ± Z α/2 ��� σ 22 √ , where σ 11 and σ 22 are the asymptotic Vs of the MLEs of c and λ, respectively; here Z α/2 is the upper α/2 th percentile of the standard normal distribution.

Bayesian Estimation
In this part, we will look at how to use Bayesian estimation to estimate the unknown parameters of an IW distribution using a PCTI method. e SEr loss function will be used for Bayesian parameter estimation. It is possible to use separate gamma priors for both parameters of the IW distribution c and λ with pdfs e hyperparameters a 1 , b 1 , a 2 , b 2 are used to represent past knowledge of the unknown parameters in this situation. e joint prior (JP) for c and λ is as follows: Hyperparameter elicitation: the informative priors (IPs) will be used to elicit the hyperparameters. e above IPs will indeed be deduced from the MLEs for (c, λ) by equating the M and V of (c j , λ j ) with both the M and V of the regarded priors (Gamma priors), where j � 1, 2, . . . , k and k is the number of observations from the IW distribution that are Advances in Civil Engineering available. us, on equating M and V of (c j , λ j ) with the M and V of gamma priors, we acquire (29) e calculated hyperparameters may now be expressed as after solving the preceding two equations According to the observed data, the appropriate posterior density (PD) x � (x (1) , x (2) , . . . , x (r) ) can indeed be expressed as e PD function is denoted by the symbol where As a result, the PD may be rewritten as follows: Under the SEr, the Bayes estimator of any function, such as g(c, λ), is provided by Consequently, equation (24) cannot be calculated for general g(c, λ). As a result, we recommend using the most widely used approximate BEs of c and λ MCMC.
where g(c, λ) is function of c and λ only and L(c, λ) is the log-LL given in (5) and ρ(c, λ) is the log-JP distribution. Using the approach developed by [19], the ratio L(X) can be expressed as where Partial derivatives of L in equations (7)- (11) are also where e log-JP is given as us, the partial derivatives of log-JP distribution are Under the SEr function, the BE of c is already provided as By substituting equations (7)-(11), (18)-(33) in (26), the estimates of c and λ can be written as c Lindley and λ Lindley .

Metropolis-Hasting Algorithm.
We need to specify IW model and beginning values for the unknown parameters c and λ to run the MH method for the IW distribution. We explore a bivariate normal distribution for the proposal distribution, that is, q((c′, λ ′ )|(c, λ)) � N 2 ((c, λ), S c,λ ). We may get negative observations, which are undesirable, if Sc, λ represents the variance-covariance matrix. We use the MLE for c and λ to determine the starting values, that is, is the Fisher information matrix. It is worth mentioning that the choice of S c,λ is critical in the MH Algorithm 1, since the acceptance rate is determined by it. e following steps are followed used by the MH method to draw a sample from the PD given by equation (24) supplied in the following manner.
Eventually, using the PD's random samples of size M, part of the initial samples can indeed be eliminated (burnin), and the remaining samples can be used to produce BEs. Extra precisely equation (24) can be estimated as where l B is the total number of burn-in samples.

Highest Posterior Density (HDP).
We use the samples generated from the proposed MH method in the preceding subsection to create HPD intervals for the unknown parameters c and λ of the IW distribution under PCTI. Consider the following scenario: c (δ) and λ (δ) are the δ th Q of c and λ, respectively, that is, where 0 < δ < 1 and Π(·) is the posterior distribution function of c and λ. It is worth noting that, for specific c * and λ * , a simulation accurate estimator of π((c, λ)|x) may be calculated as Here I (c,λ)≤(c * ,λ * ) is the indicator function. en the appropriate estimate is computed as where Let us now compute a 100(1 − δ)% HPD credible interval for c and λ represents the biggest integer that is less than or equal to a. en select HPD j * from among all HPD j ′ s with the shortest width.

Simulation Study.
In this part, we use a Monte Carlo simulation study to evaluate the performance of estimation approaches, namely, MLLE and Bayesian estimation using MCMC and Lindley's approximation, for the IW distribution using a PCTI scheme. We create 1000 data sets from the IW distribution for the MLEs under the next assumptions: (1) Two initial values are IW(c � 2, λ � 1.5) and IW(c � 1, λ � 2) (2) Sample sizes are n � 25, 50, 100 (3) Number of stages of PCTI is m � 3, 4, 5 (4) Censoring time T q j is calculated on the basis for the identified q th j Qs from IW(2, 1.5) and IW(1, 2): (a) At m � 3 and q j � (10%, 40%, 70%) (b) At m � 4 and q j � (10%, 30%, 50%, 70%) (c) At m � 5 and q j � (10%, 25%, 40%, 55%, 70%) (5) Omitted items R j are proposed based on a fraction of sample size n where the number of removed items in each stage of censoring can be computed from us, the proposed schemes of removing items can be and r is the number of failure items and ⌈⌉ is the ceiling function. It is indicated that scheme I represents Type-I censoring scheme where R m � n − r.
We construct MLEs and related 95% asymptotic CIs premised on the data that is generated. When constructing MLEs, the initial estimate values are assumed to be the same as the real parameter values.
Under the informative prior (IP), we calculate BEs using the MH algorithm for Bayesian estimation. us, we have the following. To compute the desired estimations, the aforementioned IP values are entered in. We use the MLEs as starting guess values, as well as the related V-C matrix S θ of (ln(λ), ln(α)) when using the MH algorithm. Finally, we eliminated 2000

First
Step Put initial value of θ as θ (0) � (c, λ). Second Step For i � 1, 2, . . . , M the necessary phases should be repeated: 2.1: Set θ � θ (i− 1) . 2.2: Create a new value for the candidate parameter δ from N 2 (ln θ, S θ ). 2.3: Set θ′ � exp(δ). 2.4: Compute β � (π(θ′|x)/π(θ′x)), where π(·) is the PD. 2.5: Construct a sample u from the uniform U(0, 1) model 2.6: Accept or deny the new request according to θ′ ALGORITHM 1: Algorithm of MCMC. 6 Advances in Civil Engineering burn-in samples from the total of 10,000 samples generated by the PD and then used [20] approach to get BEs and HPD interval estimations. Tables 1-6 show the average estimates for both techniques. In addition, the first row displays average estimates (AVEs) and interval estimates (IEs), while the second row displays related mean square errors (MSErs) and average lengths (ALs) with coverage probabilities (CPrs). It can be seen from the table of values that, depending on MSErs, larger values of n result in better estimates. It is also worth noting that MLEs outperform Lindley BEs and that the performance of BEs produced using MCMC outperforms Bayes estimates obtained via LiA. It is also worth noting that MCMC's ALs and related CPrs for HPD intervals outperform LiA. Furthermore, when the units are eliminated at an early stage in scheme II and scheme III, MSErs and ALs of related interval estimations are typically smaller.

Real Data Application.
Two real data set are investigated for illustration and also to examine the statistical performance of the MLEs and BEs for the IW distribution under different PCTI censoring schemes.

Data Set I.
e accompanying basic data set corresponds to an unfiltered data set. e data collection includes 34 observations of vinyl chloride data from [21], which indicates cleanup gradient ground-water monitoring wells in mg/L.
We begin by determining if the IW distribution is appropriate for studying this data set. To assess the quality of fit, we provide the MLEs of the parameters as well as the value of the Kolmogorov-Smirnov (KS) test statistic. e estimated KS and p-value for the IW distribution are 0.116 1 and 0.749 6, respectively, where c � 0.8431 and λ � 0.7121, which indicate that this distribution can be considered as an adequate model for the given data set.
From the raw data, one may construct, for example, three PCTI censored samples with m � 3, 4, 5 stages at time censoring T q j matching to the chosen q th j , Qs, where j � 1, . . . , m.
e sequences of deleted elements are presented for three fractions of eliminating f � 0%, 25%, 50% of total sample size (n � 34). Table 7 describes these patterns of removal process.
Here, (1 * 4 , 0) indicates, for example, that the censoring scheme used is (1, 1, 1, 1, 0)\enleadertwodots In Table 8, the MLEs of the parameters c and λ, as well as their related standard error and ACI at suggested schemes for PCTI samples in the provided real data set, have been computed.
Additionally, BEs were computed using the MH algorithm under the noninformative prior, i.e., a1 � b1 � a2 � b2 � 0. It is said while using the MH method to generate samples from the posterior distribution, starting values of (c, λ) are regarded as (c (0) , λ (0) ) � (c, λ) and c, λ are the MLEs of the parameters c, λ. Afterwards, 2000 burnin samples were eliminated from the overall 10 000 samples generated by the PD, and BEs and HPD intervals were produced. BEs based on LiA have been calculated based on the MLEs. e convergence of MCMC estimate is in the case of PCTI scheme III for the data set I where the percentage of removal is f � 50%. e BEs utilizing MCMC are converged through three subgraphs, as shown in Figure 2: scatter plot, histogram, and cumulative mean of the 10 000 estimates.

Data Set II.
A real data set of the carbonation depth of pier of a reinforced concrete girder bridge was analyzed under progressive Type-I censoring scheme [22] First, we check whether the IW distribution is suitable for analyzing this data set. Also, we provide the MLEs of the parameters as well as the value of the Kolmogorov-Smirnov (KS) test statistic. e estimated KS and p-value for the IW distribution are 0.207 4 and 0.195 7, respectively, where c � 5.2000 and λ � 157.8003, which indicate that this distribution can be considered as an adequate model for the given data set.
Two different PCTI censored samples are assumed for the given data set with m � 3, 5 stages at time censoring T q j matching to the chosen q th j , Qs, where j � 1, . . . , m. Table 9 describes the patterns of removing elements for three fractions of eliminating f � 0%, 25%, 50% of total sample size (n � 27).
In Table 10, the MLEs of the parameters c and λ, as well as their related ACI estimations and ACI at suggested schemes for PCTI samples in the provided real data set, have been computed.
As in data set I, BEs were computed using the MH algorithm under the noninformative prior. e starting values of (c, λ) are regarded as (c (0) , λ (0) ) � (c, λ). Finally, 2000 burn-in samples were eliminated from the overall 10 000 samples generated by the PD, and BEs and HPD intervals were produced. BEs based on LiA have been calculated based on the MLEs as a prior estimates. Figure 3 illustrates the convergence of MCMC estimate in the case of PCTI scheme II for the real data set II where the percentage of removal is f � 50%. e BEs utilizing MCMC are converged through three subgraphs, as shown in scatter plot, histogram, and cumulative mean of the 10 000 estimates.

Summary and Conclusion
We investigated the topic of IW distribution estimation and prediction under PCTI from both classical and Bayesian perspectives in this work. For the unknown parameters of the IW distribution, we calculated maximum likelihood estimates and associated asymptotic confidence intervals. en, using informative priors, we produced Bayes estimates and the related HPD interval estimates. In addition, when an informative prior is taken into account, a discussion of how to pick the values of hyperparameters in Bayesian estimation is examined based on historical samples. e simulation           12 Advances in Civil Engineering outcomes demonstrate that MLEs informative Bayes estimates using Lindley approximation perform better than both MLEs and informative prior using Lindley approximation and that estimates under informative prior using MCMC perform better than both MLEs and informative prior using Lindley approximation. We used Bayesian estimation with the squared error loss function for future work, but other loss functions can also be used. In addition, the current approach may be extended to the construction of an optimum progressive censoring, as well as alternative   14 Advances in Civil Engineering censoring methods. Neutrosophic statistics can be an extended work on area of progressive censoring schemes under the assumed distribution and PCTTI.

Data Availability
Interested parties can reach out to the author in order to receive a numerical data set used to perform the research described in the study.

Conflicts of Interest
e authors declare that there are no conflicts of interest.