Next Article in Journal
Quaternion Quantum Mechanics II: Resolving the Problems of Gravity and Imaginary Numbers
Next Article in Special Issue
Symmetrical and Asymmetrical Distributions in Statistics and Data Science
Previous Article in Journal
Exact Boundary Correction Methods for Multivariate Kernel Density Estimation
Previous Article in Special Issue
The Reliability of Stored Water behind Dams Using the Multi-Component Stress-Strength System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reliability Analysis of Kavya Manoharan Kumaraswamy Distribution under Generalized Progressive Hybrid Data

1
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
2
Department of Statistics, Faculty of Business Administration, Delta University for Science and Technology, Gamasa 11152, Egypt
3
Department of Mathematical Statistics, Faculty of Graduate Studies for Statistical Research, Cairo University, Cairo 12613, Egypt
4
Department of Statistics, Al-Azhar University, Cairo 11751, Egypt
*
Author to whom correspondence should be addressed.
Submission received: 25 July 2023 / Revised: 19 August 2023 / Accepted: 22 August 2023 / Published: 30 August 2023

Abstract

:
Generalized progressive hybrid censoring approaches have been developed to reduce test time and cost. This paper investigates the difficulties associated with estimating the unobserved model parameters and the reliability time functions of the Kavya Manoharan Kumaraswamy (KMKu) distribution based on generalized type-II progressive hybrid censoring using classical and Bayesian estimation techniques. The frequentist estimators’ normal approximations are also used to construct the appropriate estimated confidence intervals for the unknown parameter model. Under symmetrical squared error loss, independent gamma conjugate priors are used to produce the Bayesian estimators. The Bayesian estimators and associated highest posterior density intervals cannot be derived analytically since the joint likelihood function is provided in a complicated form. However, they may be evaluated using Monte Carlo Markov chain (MCMC) techniques. Out of all the censoring choices, the best one is selected using four optimality criteria.

1. Introduction

The progressive type-II censoring (PCS-T2) method is the most popular scheme in reliability and survival analysis. Compared with the traditional type-II censoring method, it is better. Progressive censoring is advantageous in a variety of real-world applications, including business, medical research, and therapeutic settings. Up until the test’s conclusion, it permits the removal of any remaining experimental units. Assume that n units are used in a life test and that it is not desirable to record every failure because of financial and time constraints. Consequently, only a portion of unit failures are seen. A sample like this is known as a censored sample. Assume that one of the units was accidentally damaged after the test started but before they all burned out. This unit needs to be taken out of the life test if the experiment is still going on. In this situation, a framework for analyzing this kind of data is provided by the progressive censoring scheme. A few examples of primary references are [1,2].
PCS-T2 has drawn a lot of attention in the literature as a very flexible censoring system (see [3] for further details). When testing n independent units at a time T = 0 , the failure number to be noticed s and the progressive censored samples, R _ = R 1 , R 2 ,   , R s   , where n = i = 1 s R i + s , are specified. When the initial failure is seen (suppose that Y 1 : s : n ), the other surviving units n 1 are chosen at random, and R 1 of those units is disqualified from the test. Similarly, at the moment of the second failure (suppose that Y 2 : s : n ), R 2   of n R 1 2 are selected at random and deleted from the test, and so on. At the time of the s t h failure (suppose that Y s : s : n ), every survivor unit still present R s = n s j = 1 s 1 R j is removed from the experiment.
Whenever the test units are particularly reliable, the major drawback of this censoring is that it could take longer to finish the progressively type-II hybrid censored samples (PHCS-T2). The authors of [4] proposed a progressive type-I hybrid censored strategy (PHCS-T1) as a remedy for this issue. This method combines PCS-T2 with conventional type-I censoring. Under PHCS-T1, the trial period is stopped at T , maximum likelihood estimators (MLEs) were not always available due to the fact that relatively a few failures might occur before time T in PHCS-T1. To resolve this issue, [5] presented the PHCS-T2 scheme. At T * = m a x Y s : s : n ,   T , the experiment comes to an end under PHCS-T2. It can take some time until such s t h failures are really observed, despite the fact that PHCS-T2 promises a fixed number of failures.
It could take a while to gather the needed failures, even though the PHCS-T2 ensures an effective number of observable failures. Thus, [6] devised the generalized progressive type-II hybrid censoring (GPHC-T2). Assume that the thresholds T i , i = 1 , and 2, as well as the integer s , are preassigned in such a way that 0 < T 1 < T 2 < and 1 < s < n . c 1 and c 2 represent the overall number of failures up to periods T 1   and   T 2 , respectively. Then, at Y 1 : s : n , R 1 of n 1 are arbitrarily excluded from the test, followed by R 2 of n R 1 2 , and so on.
The experiment is over, and all remaining units are deleted at T * = m a x T 1 , m i n Y s : s : n ,   T 2 . If Y s : n < T 1 , failures are observed without any further withdrawals up until time T 1 (Case-I); if T 1 < Y s : s : n < T 2 , the test is terminated at time Y s : s : n (Case-II); or, if not, the test is terminated at time T 2 (Case-III). Keep in mind that the GPHCS-T2 modifies the PHCS-T2 by guaranteeing that the test is completed at the scheduled time T 2 . T 2 demonstrates the longest period of time the researcher is willing to let the experiment continue. As a result, one of the following three data types will be visible to the experimenter:
Y _ , R _ = C a s e   I ; Y 1 : n , R 1 , , Y s 1 : s : n , R s 1 , Y s : s : n , 0 , , Y c 1 : n , 0 C a s e   I I ; Y 1 : s : n , R 1 , , Y c 1 : n , R c 1 , Y s 1 : s : n , R s 1 , , Y s : s : n , R s C a s e   I I I ; Y 1 : s : n , R 1 , , Y c 1 : n , R c 1 , Y c 2 1 : n , R c 2 1 , , Y c 2 : n , R c 2 .
Figure 1 indicates the cases of generalized type-II progressive hybrid sample as follows:
Assume that in a distribution with a cumulative distribution function (cdf) F(.), and probability density function (pdf) f(.), the variables Y _ and R _ represent the respective lifetimes. As a result, the GPHCS-T2 likelihood function is expressed as follows:
L φ θ , β y _ = C φ j = 1 D φ f Y j : s : n ; θ , β 1 F Y j : s : n ; θ , β R j ψ φ T τ ; θ , β ,
where τ = 1,2 ,   φ = 1,2 , 3 ,   stand in for Case-I, Case-II, and Case-III, respectively, and ψ φ . is a combination form of dependability functions. Table 1 displays the GPHCS-T2 notations from Equation (1). Many censoring techniques can also be inferred as particular examples from Equation (1), including
  • With T 1 setting to 0, use PHCS-T1.
  • T 2 . by setting PHCS-T2.
  • You may do hybrid type-I censoring by setting T 1 0 ,   R j = 0 ,   j = 1,2 , , s 1 ,   R s = n s .
  • T 2 ,   R j = 0 ,   j = 1,2 , , s 1 , R s = n s can be used to do hybrid type-II censoring.
  • To do type-I censoring, set T 1 = 0 ,   s = 1 ,   R j = 0 ,   j = 1 ,   2 , ,   s 1 ,   R s = n s .
  • A type-II censored sample is produced by setting T 1 = 0 ,   T 2 ,   s = 1 ,   R j = 0 ,   j = 1 ,   2 , ,   s 1 ,   R s = n s .
On the basis of GPHCS-T2, more studies have been conducted. For instance, Ref. [7] investigated the prediction issue of forthcoming Burr-XII distribution failure rates. The authors of [8] created the Weibull distribution with little data with an objective Bayesian analysis. The authors of [9] addressed the competing risks from exponential data, and [10] more recently examined both the point and interval estimations of the Burr-XII parameters. Last but not least, [11] addressed the Fréchet distribution’s optimality under generalized censoring schemes. In this paper, the KMKu model under generalized censoring samples is studied. Where the KMKu model was initially proposed by [11]. Also, they found that the Kumaraswamy model’s and KMKu shape forms in the pdf for different parameter values are comparable. It may be asymmetric, unimodal, increasing, or decreasing. In addition, the bathtub, U-shape, J-shape, or increasing shapes of the hazard rate function (hrf) for the KMKu model are all possible. But suppose that Y is the lifespan random variable of a test item adheres to the KMKu distribution, denoted by the notation K M K u ( θ , β ) , where θ > 0 , β > 0 are the shape parameters. Therefore, it is supplied by its pdf, cdf, reliability function (RF), R ( . ) , and hrf, all represented by the letters f ( . ) , F ( . ) , and h ( . ) accordingly:
f y ; θ , β = θ β y θ 1 e 1 1 y θ β 1 e 1 y θ β , 0 < y < 1 ; θ β > 0 ,
F y ; θ , β = e e 1 1 e 1 e 1 y θ β ,
R y ; θ , β = 1 e e 1 1 e 1 e 1 y θ β ,
and
h y ; θ , β = θ β y θ 1 e 1 1 y θ β 1 1 e 1 y θ β .
Although the KMKu model has a lot of flexibility because of its different shapes of hrf and pdf, to our knowledge, no studies have yet been done under censorship. Particularly, the generalized type-II progressively hybrid censoring scheme has not produced any data for the new KMKu lifetime model’s survival traits and model parameters. To fill this gap, the following are the objectives of this study: Firstly, the probability inference for any function of the unknown KMKu parameters, such as R(t) or h(t), is derived. The second objective is to derive independent gamma priors from the squared error (SE) loss and produce Bayes estimates for the same unknown parameters, employing the provided estimation procedures, such as classical and Bayesian approaches. The unknown parameters of the KMKu distribution are discovered using the approximation confidence intervals (ACIs) and highest posterior density (HPD) interval estimators. The acquired estimates are computed using the R programming language’s “maxLik” and “coda” packages because the theoretical findings of θ   and β obtained by the suggested estimation techniques cannot be represented in closed form. [12,13] offered these packages. Using four optimality criteria, the ultimate aim is to develop the most efficient progressively censored sample technique. The effectiveness of the different estimators is investigated using a Monte Carlo simulation with the entire sample size, which can be combined in a variety of ways, effective sample size, threshold timings, and progressively censored samples. We compare the average confidence lengths (ACLs), mean relative absolute biases (MRABs), and simulated root mean squared errors (RMSEs) of the derived estimators. The optimal censoring tactic should be chosen after evaluating how effectively the given techniques will function in practice. The remaining portions of this study are structured as follows: The maximum likelihood, Bayes inferences, and reliability functions of the unknown parameters are presented in Section 2 and Section 3, respectively. The credible and asymptotic intervals are built into Section 4. Section 5 goes into depth about the results of the Monte Carlo simulation. The optimal methods for progressive censoring are discussed in Section 6. Two actual data sets are indicated in Section 7. Finally, the conclusion and discussion are given in Section 8.

2. Likelihood Estimation

Assume that the representation of a GPHCS-T2 sample of size   c 2 taken from K M K u ( θ , β ) is Y = ( ( Y 1 : s : n ,   R 1 ) , , ( Y c 1 : n ,   R c 1 ) , , ( Y c 2 : n ,   R c 2 ) ) . The probability function of GPHCS-T2 may be represented by substituting y j for y j : s : n in Equation (1) and adding Equations (2) and (3); for more information, see [14].
L φ θ , β Y _ j = 1 D φ θ β y j θ 1 e 1 1 y j θ β 1 e 1 y j θ β 1 e e 1 + e 1 e 1 y θ β R i ψ φ T τ ; θ , β ,
where
ψ 1 T 1 ; θ , β = 1 e e 1 1 e 1 e 1 T 1 θ β R c 1 + 1 * ,   ψ 2 T τ ; θ , β = 1 and ψ 3 T 2 ; θ , β = 1 e e 1 1 e 1 e 1 T 2 θ β R c 2 + 1 * .
The proper log-likelihood function for Equation (6) is l φ . l n   L φ . as follows:
l φ θ , β Y _ D φ l n θ β + θ 1 j = 1 D φ l n y j D φ l n e 1 + β 1 j = 1 D φ l n 1 y j θ + β j = 1 D φ ( 1 y j θ ) + R i j = 1 D φ l n 1 e e 1 + e 1 e 1 y j θ β + γ φ T τ ; θ , β ,
where
γ φ T 1 ; θ , β = R c 1 + 1 * l n 1 e e 1 1 e 1 e 1 T 1 θ β , γ 2 T τ ; θ , β = 1 , and γ 3 T 2 ; θ , β = R c 2 + 1 * l n 1 e e 1 1 e 1 e 1 T 2 θ β .
By partially differentiating Equation (7) with reference to θ ^ and β ^ , the subsequent two findings are produced. After being equal to zero, likelihood equations must be simultaneously solved in order to create the MLEs.
l φ θ = D φ θ + j = 1 D φ l n y j β 1 j = 1 D φ y j θ l n y j 1 y j θ β j = 1 D φ y j θ l n y j R i j = 1 D φ e 1 y j θ β y j θ β 1 y j θ β 1 l n y j e 1 e e 1 + e 1 e 1 y j θ β + γ φ T τ ; θ , β θ ,  
and
l φ β = D φ β + 2 j = 1 D φ l n 1 y j θ + R i j = 1 D φ e 1 y j θ β 1 y j θ β l n 1 y j θ e 1 e e 1 + e 1 e 1 y j θ β + γ φ T τ ; θ , β β ,  
where φ = 1 ,   3 and τ = 1 ,   2 , respectively, we have
γ φ T τ ; θ , β θ = R c τ + 1 * e 1 T τ θ β β T τ θ l n T τ 1 e e 1 + e 1 e 1 T τ θ β , γ φ T τ ; θ , β β = R c τ + 1 * e 1 T τ θ β β l n 1 T τ θ e 1 e e 1 + e 1 e 1 T τ θ β .
According to Equations (8) and (9), it is necessary to simultaneously satisfy a system of two nonlinear equations in order to derive the MLEs of θ   and β   in the KMKu model. As a result, for θ   and β , there is not, and cannot be computed, an analytical closed-form solution. Thus, it may be estimated for each specific GPHCS-T2 data set using numerical techniques like the Newton-Raphson iterative method. When the estimates of θ   and β are derived by replacing them with θ ^   and β ^ , the MLEs R ^ ( t ) and h ^ ( t ) , respectively, may be easily computed.

3. Bayes Estimator

The HPD intervals for the Bayes estimators of θ   , β , R(t), and h(t) are developed using the SE loss function. To do this, it is assumed that the KMKu parameters θ   and β , respectively, have independent gamma priors of the forms ω ( ν 1 ,   ν 2 ) and ω ( ν 3 ,   ν 4 ) .
The normal distribution can be a standard choice for data if the domain of that distribution is from −∞ to ∞, and the beta distribution can be a standard choice for data if the domain of that distribution is from 0 to 1. Similarly, the gamma distribution can be a standard choice for non-negative continuous data if the domain of the gamma distribution is from 0 to ∞. This is one of the most important reasons, but there are other reasons as follows:
  • We believe the main motivation for the gamma prior is usually to constrain the random variables to positive values.
  • The gamma distribution is considered one of the most important and well-known statistical distributions because it is compatible with many engineering, mathematical, statistical, and medical applications.
  • The gamma distribution is one of the most famous distributions that is used in mathematical solutions (integrations), especially when the data are from 0 to ∞.
  • In previous studies, the gamma distribution was the most popular prior distribution and was associated with the best statistical results.
Gamma priors should be considered for a variety of reasons, including the fact that they are (1) adjustable, (2) offer diverse shapes based on parameter values, and (3) fairly basic and brief and might not generate a solution to a challenging estimation problem. Then, the combined previous density of θ and β is determined; for more details on this topic, see [15,16].
π θ , β θ ν 1 1 β ν 3 1 e θ ν 2 + β ν 4
If it is anticipated that for i = 1 ,   2 ,   3 , 4 , ν i > 0 are known. The joint posterior pdf of θ   and β , Equations (6) and (10), when combined, results.
π φ θ , β y _ θ D φ + ν 1 1 β D φ + ν 3 1 e θ ν 2 + β ν 4 j = 1 D φ θ β y j θ 1 e 1 1 y j θ β 1 e 1 y j θ β 1 e e 1 + e 1 e 1 y θ β R i ψ φ T τ ; θ , β
The Bayes estimate, η ~ θ , β , of θ and β respectively, under SE loss, η θ , β   is what is meant by the posterior expectation of Equation (11), which is given.
η ~ θ , β = 0 0 η θ , β π φ θ , β y _ d θ d β .
It is clear from Equation (11), that it is impossible to explicitly express the marginal pdfs of θ   and β . In order to accomplish this, we recommend creating samples from Equation (11) utilizing Bayes MCMC methods to calculate the joint Bayes estimates and supplying their HPD intervals. The complete conditional pdfs of θ   and β are provided for the MCMC sampler from Equation (11) to be performed as intended.
π φ θ θ β ,   y _ θ D φ + ν 1 1 e θ ν 2 j = 1 D φ θ β y j θ 1 e 1 1 y j θ β 1 e 1 y j θ β 1 e e 1 + e 1 e 1 y θ β R i ψ φ T τ ; θ , β ,  
and
π φ β β θ ,   y _ β D φ + ν 3 1 e β ν 4 j = 1 D φ θ β y j θ 1 e 1 1 y j θ β 1 e 1 y j θ β 1 e e 1 + e 1 e 1 y θ β R i ψ φ T τ ; θ , β .
The Metropolis-Hastings (M-H) approach is considered to be the best solution to this problem because no analytical method exists to reduce the posterior pdfs of θ and β in Equations (12) and (13), respectively, to any known distribution (for further information, see [17,18]. The sampling method of the M-H algorithm is implemented according to:
First, establish the starting points, θ 0 = θ ^ and β 0 = β ^ .
Set S = 1 after that.
Thirdly, from N ( μ ^ 1 , σ ^ 1 ) and N ( μ ^ 2 , σ ^ 2 ) , respectively, create θ * and β * .
The fourth step: Obtaining ϱ θ = m i n 1 , π φ θ θ * β s 1 ; y _ π φ θ θ s 1 β s 1 ; y _ and ϱ β = m i n 1 , π φ β β * θ s ; y _ π φ β β s 1 θ s ; y _ .
Fifth, use the uniform U ( 0,1 ) distribution to generate the samples u 1 and u 2 .
Sixth: Set θ S = θ * and β S = β * ,   respectively, if u 1 and u 2 are both smaller than ϱ θ and ϱ β , respectively. Set θ S = θ S 1   and β S = β S 1 , correspondingly, if not.
Seventh: Establish that S equals S + 1.
Eighth: Repeating steps three through seven a number of times B will give you the values for θ S   and β S   for S = 1 ,   2 , ,   B .
Ninth: To calculate the RF in Equation (4) and hrf in Equation (5), use θ S   and β S for S = 1 ,   2 , ,   B , respectively, for a given mission period t > 0 .
R S t = 1 e e 1 1 e 1 e 1 y θ S β S , y > 0 ,
and
h S t = θ S β S y θ S 1 e 1 1 y θ S β S 1 1 e 1 y θ S θ S , y > 0 .
The convergence of the MCMC sampler must be ensured, and starting, θ 0 and β 0   values must be eliminated. The first simulated variants, let us say B 0 , are removed as burn-ins. Therefore, using the remaining B B 0 samples of θ , β ,   R ( t ) , or h t , (let us suppose η ), the Bayesian estimates are computed. On the basis of the SE loss function, the Bayes MCMC estimates of η are shown.
η ~ = 1 B B 0 S = B 0 + 1 Β η S

4. Interval Estimators

The HPD interval estimators in this section are based on acquired MCMC-simulated variations, as opposed to the approximative confidence estimators of θ , β ,   R ( t ) , or h t that are based on observed Fisher information.

4.1. Asymptotic Intervals

To compute the ACIs for θ and β , the Fisher information matrix must first be inverted to produce the asymptotic variance-covariance (AVC) matrix. According to certain regularity criteria, ( θ ^ , β ^ ) is nearly normal with a mean ( θ , β ) and variance I   1   θ , β . In agreement with [19], we estimate I   1   θ , β   by I   1   θ ^ , β ^ , replacing θ ^ and β ^ for θ and β .
I   1   θ ^ , β ^ = a 11 a 12 a 21 a 22 1
where
a 11 = 2 l φ θ 2 = D ρ θ 2 β 1 j = 1 D φ y j θ l n y j 1 + l n y j 1 y j θ 1 y j θ 2 β j = 1 D φ y j θ l n y j 2 + R i j = 1 D φ β y j θ l n y j 2 e 1 y j θ β 1 y j θ β 1 1 + β y j θ 1 y j θ β + y j θ β 1 1 y j θ e 1 e e 1 + e 1 e 1 y j θ β 2 + 2 γ φ T τ ; θ , β θ 2
a 12 = a 21 = 2 l φ θ β = j = 1 D φ y j θ l n y j 1 y j θ j = 1 D φ y j θ l n y j R i j = 1 D φ y j θ l n y j e 1 y j θ β 1 y j θ β 1 1 + β 1 y j θ β l n 1 y j θ + β l n 1 y j θ 1 y j θ e 1 e e 1 + e 1 e 1 y j θ β 2 + 2 γ φ T τ ; θ , β θ β ,
a 22 = 2 l φ β 2 = D φ β 2 + R i j = 1 D φ l n 1 y j θ 2 e 1 y j θ β 1 y j θ β + 1 y j θ 2 β e 1 e e 1 + e 1 e 1 y j θ β 2 + 2 γ φ T τ ; θ , β β 2 ,
where
2 γ φ T τ ; θ , β θ 2 = R c τ + 1 * β l n T τ T τ θ e 1 T τ θ β 1 e e 1 + e 1 e 1 T τ θ β 1 T τ θ 1 T τ θ β 1 + T τ 2 θ e 1 β 1 T τ θ β 1 1 e e 1 + e 1 e 1 T τ θ β 2 ,
2 γ φ T τ ; θ , β θ β = R c 1 + 1 * e 1 e 1 T 1 θ β 1 T τ θ β l n 1 T τ θ 1 e e 1 1 e 1 e 1 T τ θ β ,
and
2 γ φ T τ ; θ , β β 2 = R c τ + 1 * l n 1 T τ θ e 1 T τ θ β e 1 e e 1 + e 1 e 1 T τ θ β 2 1 e e 1 + e 1 e 1 T τ θ β 1 + β 1 T τ θ β l n 1 T τ θ β e 1 1 T τ θ β l n 1 T τ θ e 1 T τ θ β .  
The two-sided 100 1 γ % ACIs are therefore given by
θ ^ ± Z γ 2 σ ^ 1 2 and β ^ ± Z γ 2 σ ^ 2 2 , for θ   and β , respectively, where Z γ 2 stands for the top γ 2 percentage points of the standard normal distribution, σ ^ 1 2 and σ ^ 2 2 are the primary diagonal elements of Equation (14). Furthermore, we employ the delta method to first establish the estimated variance of R ^ t and h ^ t (see [20]) before developing the ACIs of R t and h t as
σ ^ R ^ t 2 = ϵ R ^ T I 1 ε ^ ϵ R ^   and   σ ^ h ^ t 2 = ϵ h ^ T I 1 ε ^ ϵ h ^ ,
where ϵ R ^ T = R t θ R t β θ ^ ,   β ^ , and ϵ h ^ T = h t θ h t β θ ^ ,   β ^
Then, R(t) and h(t) both have two-sided 100 1 γ % ACIs that are supplied by R ^ t ± Z γ 2 σ ^ R ^ t 2 and h ^ t ± Z γ 2 σ ^ h ^ t 2 , respectively.
Adding bootstrapping techniques to improve estimators or create confidence intervals for θ , β , R t , or h t is easy.

4.2. HPD Intervals

The method put forward by [21] is used to create 100 1 γ % HPD interval estimates of θ ,   β , R t , or h t . First, we assign numerical values to the MCMC samples of ε j for j = B 0 + 1 ,   B 0 + 2 , ,   B   a s   ε B 0 + 1 ,   ε B 0 + 2 , ,   ε B correspondingly. The discovery is that the 100 1 γ % two-sided HPD interval of ε is supplied by ε j * ,   ε j * + 1 ϵ B B 0 , where j * = B 0 + 1 ,   B 0 + 2 , ,   B is selected so that ε j * + 1 ϵ B B 0 ε j * = m i n 1 j γ B B 0 ε j + 1 γ B B 0 ε j   ε j * ,   ε j * + 1 γ B B 0 .

5. Optimal PCS-T2 Designs

The experimenter may want to pick the “best” censoring scheme out of a collection of all accessible censoring schemes in order to provide the most details about the unknown parameters under investigation, especially in the context of dependability. First, [1] examined the problem of deciding which censoring strategy is most appropriate under various circumstances. However, a number of optimality criteria, R = R 1 ,   R 2 , ,   R s , where i = 1 s R i have been proposed, and several assessments of the top censoring strategies have been made. The precise values of n   (total test units), s (effective sample), and T i , i = 1 ,   2 (ideal test thresholds) are picked in advance according to the accessibility of the units, the accessibility of the experimental settings, and cost factors (see [22]). A number of articles in the literature have addressed the topic of contrasting two (or more) different censoring techniques. For examples, see [23,24]. To help us choose the best censoring strategy, O i , Table 2 offers a variety of widely used measures.
It is advised that the observed Fisher information, I 2 × 2 . values for O 1 , be maximized. For criterion O 2 and O 3 , we also wish to reduce the determinant and trace of I 2 × 2 . 1 . The best censoring strategy for multi-parameter distributions may be selected using scale-invariant criteria. While dealing with unknown multi-parameter distributions makes it more challenging to compare the two Fisher information matrices, dealing with single-parameter distributions allows for the use of scale-invariant criteria to compare a variety of criteria O 4 . The logarithmic MLE of the p t h quantile, l o g t ^ p , tends to have a variance that is minimized by the p-dependent criterion O 4 . As a result, the logarithm of the KMKu distribution for time t ^ p may be calculated using
l o g t ^ p = 1 ln e 1 p e 1 e 1 β 1 θ ,   0 < p < 1 ,
By using the delta technique to solve for Equation (4), the estimate of the variance for the l o g t ^ p   of the KMKu distribution is given as
V a r l o g t ^ p = l o g t ^ p T I 2 × 2 1 θ ^ ,   β ^ l o g t ^ p ,
where
l o g t ^ p T = θ l o g t ^ p , β l o g t ^ p θ = θ ^ , β = β ^
P R 1 = K 1 = n s K 1 r K 1 ( 1 r ) n s K 1 .
while i = 2 , 3 , , s 1 .   The maximum value of the O 1 criterion and the lowest value of O i ,   i = 2 , 3 , 4 , correspond to the best censoring. On the other hand, the greatest value of the O 1 criterion and the lowest value of the O i ,   i = 2 , 3 , 4 criterion correspond to the best censoring.

6. Simulation

Using different combinations of T i ; i = 1 , 2 (threshold points), n (sample size), s (size of censored sample), and R (censored removal), Monte-Carlo (MC) simulations were carried out to assess the true performance of the acquired point and interval estimators of θ , β ,   R T , and h T . To establish this goal, for KMKu(1.4, 1.5), KMKu(1.4, 0.5), and KMKu(0.4, 0.5), we replicated the GPHCS-T2 mechanism 1000 times. Taking T 1 , T 2 = 0.6 ,   0.85 , two different choices of n and s were used as (n = 30, 50, 100), and the choices of s were used as (s = 20, 25) at n = 30, (s = 35, 45) at n = 50, and (s = 70, 90) at n = 100. At T 1 = 0.6, the true values of R T 1 and h T 1 were 0.4278 and 1.4899, respectively. At T 2 = 0.85, the true values of R T 2 and h T 2 were 0.2526 and 3.3106, respectively.
Additionally, by utilizing the binomial elimination distribution and taking into account different censoring schemes for each combination of s and n , the following is conducted: according to the following probability mass function, the number of units removed at each failure time is expected to follow a binomial distribution.
P R i = K i R i 1 = K i 1 , , R 1 = K 1 = n s j = 1 i 1 K j K i r K i ( 1 r ) n s j = 1 i K j .
Additionally, assume that for any i , R i is independent of X i . In light of this, the likelihood function can be written as follows:
L x i , β ,   θ ,   r = L 1 ( x i , β ,   θ R = K ) P R = K ,
where
P R = K = P R 1 = K 1 , R 2 = K 2 , , R s 1 = K s 1 = P R s 1 = K s 1 R s 2 = K s 2 , , R 1 = K 1 ×           P R s 2 = K s 2 R s 3 = K s 3 , , R 1 = K 1 P R 2 = K 2 R 1 = K 1 ) P ( R 1 = K 1 .
That is,
P R = K = n s ! ( n s i = 1 s 1 K i ) ! i = 1 s 1 K i r i = 1 s 1 K i ( 1 r ) s 1 n s i = 1 s 1 s i K i ,
where the GPHCS-T2-based KMKu distribution’s parameters do not affect the binomial parameter r (Independent). We chose the binomial parameter r with varied values of 0.3 and 0.8.
The MLEs and 95% ACI estimates of θ , β ,   R t , and h t were assessed after 1000 GPHCS-T2 samples had been gathered using R 4.2.2 programming software and the “maxLik” library. We simulated 12,000 MCMC samples and omitted the first 2000 iterations as burn-in to obtain the Bayes point estimates along with their HPD interval estimates of the same unknown parameters using the “coda” library in the R 4.2.2 programming language. The estimates and their variances were equated with the Fisher information matrix of θ and β to produce the ML estimator, which it denoted as elective hyper-parameters, and this was contributed by [25]. This process allowed for the extraction of the hyper-parameters of the informative priors.
Some observations from Table 3, Table 4 and Table 5 include
  • The key general finding is that the suggested values for θ , β ,   R t , and h t performed well.
  • All estimations of θ , β , R(t), and h(t) functioned satisfactorily as n(or s) grew.
  • In most cases, the MSE, Bias, and WCI of all unknown parameters fell while their CPs grew as (T1, T2) increased.
  • Due to the gamma information, the Bayes estimates of θ , β ,   R t , and h t behaved more predictably than the other estimates. Regarding credible HPD intervals, the same statement might be made.
  • When the parameter of binomial r was increased, the proposed estimates of θ , β ,   R t , and h t performed better in most cases.

7. Application

The data set, which has been examined by [11], had 30 assessments of the tensile strength of polyester fibers. The following details are included in the data set: “0.023, 0.032, 0.054, 0.069, 0.081, 0.094, 0.105, 0.127, 0.148, 0.169, 0.188, 0.216, 0.255, 0.277, 0.311, 0.361, 0.376, 0.395, 0.432, 0.463, 0.481, 0.519, 0.529, 0.567, 0.642, 0.674, 0.752, 0.823, 0.887, 0.926”. For data on the strength of polyester fibers, where the Kolmogorov-Smirnov distance is 0.0569 with a p-value of 0.9999, [11] explores the MLE of this model using several measures of goodness-of-fit. The Kolmogorov-Smirnov test findings showed that the KMKu distribution fits the data on polyester fiber strength.
Two GPHCS-T2 samples with s = 20 and 25 were produced from the tensile strength of polyester fibers data in order to explain the proposed estimation methodology. The binomial removal has been used to obtain the GPHCS-T2 samples with different parameters of p = 0.2, 0.5, and 0.8. Table 6 lists the computed R(t) and h(t) at t = 0.6 and 0.85 by maximum likelihood estimates (MLE) and Bayesian estimation, respectively, along with their standard error (SE). By repeating the MCMC sampler 12,000 times and disregarding the first 2000 times as burn-in, the Bayes estimates (with their SE) were evaluated using incorrect gamma priors and are also provided in Table 4 because there was no prior knowledge about the unknown KMKu parameters θ , and β from the given data set. In order to estimate unknown hyperparameters for the computational logic, elective hyperparameters were employed. In terms of the minimum standard error and interval width values, it is evident from Table 6 that the MCMC estimates of θ , β , R t , and h t performed better than the others.
Figure 2, Figure 3 and Figure 4 were created to examine the maximum values of the estimators by profile likelihood as well as the existence and uniqueness of the log-likelihood function by contour plot with regard to different d and q options based on GPHCS-T2 samples with s = 20 and distinct p = 0.2, 0.5, and 0.8, respectively. Figure 5 clearly shows that the MCMC technique converged favorably and that the recommended size of the burn-in sample was adequate to completely nullify the impact of the recommended beginning values. Figure 5 demonstrates that the estimated estimates of θ and β were roughly symmetrical for each sample when s = 20.
Figure 6, Figure 7 and Figure 8 were created to examine the maximum values of the estimators by profile likelihood as well as the existence and uniqueness of the log-likelihood function by contour plot with regard to different d and q options based on GPHCS-T2 samples with s = 25 and distinct p = 0.2, 0.5, and 0.8, respectively. Figure 9 clearly shows that the MCMC technique converged favorably and that the recommended size of the burn-in sample was adequate to completely nullify the impact of the recommended beginning values. Figure 9 demonstrates that the estimated estimates of θ and β were roughly symmetrical for each sample when s = 25.

8. Conclusions and Discussion

This paper examines the reliability analysis of the unknown parameters, reliability, and hazard rate functions for the generalized type-II progressive hybrid censoring-based KMKu model. The “maxLik” package of the R programming language was used to compute the frequentist estimates with their asymptotic confidence intervals for the unknown parameters and any function of them. Since the likelihood function was produced in complex form, the posterior density function was obtained in nonlinear form. Consequently, the Bayesian estimates and the related HPD intervals were created using the Metropolis-Hastings technique and accounting for the squared error loss function. Numerous simulation experiments were run utilizing various total test unit choices, observed failure data, threshold times, and progressive censoring schemes in order to compare the behavior of the collected estimates. The outcomes demonstrated that the Bayes–MCMC strategy performed substantially better than the frequentist approach. Under generalized type-II progressive hybrid censoring, it was suggested to estimate the KMKu distribution’s parameters, reliability, and hazard functions using the Bayesian MCMC paradigm. We believe that the technique and results described here will be helpful to reliability practitioners and that they will be used to inform future censoring tactics. The 30 assessments of the tensile strength of polyester fibers are used to demonstrate how the recommended strategies may be applied in real-world circumstances. The most important results can be summarized in the following points:
  • The key general finding is that the suggested values for θ , β ,   R t , and h t performed well.
  • All estimations of θ ,   β , R(t), and h(t) functioned satisfactorily as n (or s) grew.
  • In most cases, the MSE, Bias, and WCI of all unknown parameters fell while their CPs grew as (T1, T2) increased.
  • Due to the gamma information, the Bayes estimates of θ , β , R t , and h t behaved more predictably than the other estimates. Regarding credible HPD intervals, the same statement might be made.
  • In most cases, the proposed estimates of θ , β , R t , and h t performed better when the parameter of binomial r was increased.
  • The MLE has a unique solution and a maximum value of log-likelihood.

Author Contributions

Conceptualization, R.A., E.M.A. and H.R.; methodology, R.A., E.M.A. and H.R.; writing—original draft preparation, R.A., E.M.A. and H.R.; writing—review and editing, R.A., E.M.A. and H.R. All the authors have the same contribution and agree to publish the manuscript in Symmetry Journal. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R50), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R50), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Balakrishnan, N.; Aggarwala, R. Progressive Censoring: Theory, Methods, and Applications; Springer Science & Business Media: Berlin, Germany, 2000. [Google Scholar]
  2. Balakrishnan, N. Progressive censoring methodology: An appraisal. TEST 2007, 16, 211–259. [Google Scholar] [CrossRef]
  3. Balakrishnan, N.; Cramer, E. The Art of Progressive Censoring; Springer: New York, NY, USA, 2014. [Google Scholar]
  4. Kundu, D.; Joarder, A. Analysis of Type-II progressively hybrid censored data. Comput. Stat. Data Anal. 2006, 50, 2509–2528. [Google Scholar] [CrossRef]
  5. Childs, A.; Chandrasekar, B.; Balakrishnan, N. Exact likelihood inference for an exponential parameter under progressive hybrid censoring schemes. In Statistical Models and Methods for Biomedical and Technical Systems; Vonta, F., Nikulin, M., Limnios, N., Huber-Carol, C., Eds.; Birkhäuser: Boston, MA, USA, 2008; pp. 319–330. [Google Scholar]
  6. Lee, K.; Sun, H.; Cho, Y. Exact likelihood inference of the exponential parameter under generalized Type II progressive hybrid censoring. J. Korean Stat. Soc. 2016, 45, 123–136. [Google Scholar] [CrossRef]
  7. Ateya, S.; Mohammed, H. Prediction under Burr-XII distribution based on generalized Type-II progressive hybrid censoring scheme. J. Egypt. Math. Soc. 2018, 26, 491–508. [Google Scholar]
  8. Seo, J.I. Objective Bayesian analysis for the Weibull distribution with partial information under the generalized Type-II progressive hybrid censoring scheme. Commun. Stat.-Simul. Comput. 2020, 51, 5157–5173. [Google Scholar] [CrossRef]
  9. Cho, S.; Lee, K. Exact likelihood inference for a competing risks model with generalized Type-II progressive hybrid censored exponential data. Symmetry 2021, 13, 887. [Google Scholar] [CrossRef]
  10. Nagy, M.; Bakr, M.E.; Alrasheedi, A.F. Analysis with applications of the generalized Type-II progressive hybrid censoring sample from Burr Type-XII model. Math. Probl. Eng. 2022, 2022, 1241303. [Google Scholar] [CrossRef]
  11. Alotaibi, N.; Elbatal, I.; Shrahili, M.; Al-Moisheer, A.S.; Elgarhy, M.; Almetwally, E.M. Statistical Inference for the Kavya–Manoharan Kumaraswamy Model under Ranked Set Sampling with Applications. Symmetry 2023, 15, 587. [Google Scholar] [CrossRef]
  12. Henningsen, A.; Toomet, O. maxlik: A package for maximum likelihood estimation in R. Comput Stat. 2011, 26, 443–458. [Google Scholar] [CrossRef]
  13. Plummer, M.; Best, N.; Cowles, K.; Vines, K. CODA: Convergence diagnosis and output analysis for MCMC. RNews 2006, 6, 7–11. [Google Scholar]
  14. Luo, C.; Shen, L.; Xu, A. Modelling and estimation of system reliability under dynamic operating environments and lifetime ordering constraints. Reliab. Eng. Syst. Saf. 2022, 218, 108136. [Google Scholar] [CrossRef]
  15. Shirong, Z.; Ancha, X.; Yincai, T.; Lijuan, S. Fast Bayesian Inference of Reparameterized Gamma Process with Random Effects. IEEE Trans. Reliab. 2023, 1–14. [Google Scholar] [CrossRef]
  16. Zhuang, L.; Xu, A.; Wang, X.-L. A prognostic driven predictive maintenance framework based on Bayesian deep learning. Reliab. Eng. Syst. Saf. 2023, 234, 109181. [Google Scholar] [CrossRef]
  17. Gelman, A.; Carlin, J.B.; Stern, H.S.; Rubin, D.B. Bayesian Data Analysis, 2nd ed.; Chapman and Hall/CRC: Boca Raton, FL, USA, 2004. [Google Scholar]
  18. Lynch, S.M. Introduction to Applied Bayesian Statistics and Estimation for Social Scientists; Springer: New York, NY, USA, 2007. [Google Scholar]
  19. Lawless, J.F. Statistical Models and Methods for Lifetime Data, 2nd ed.; John Wiley and Sons: Hoboken, NJ, USA, 1982. [Google Scholar]
  20. Greene, W.H. Econometric Analysis, 4th ed.; Prentice-Hall: New York, NY, USA, 2000. [Google Scholar]
  21. Chen, M.H.; Shao, Q.M. Monte Carlo estimation of Bayesian credible and HPD intervals. J. Comput. Graph. Stat. 1999, 8, 69–92. [Google Scholar]
  22. Ng, H.K.T.; Chan, P.S.; Balakrishnan, N. Optimal progressive censoring plans for the Weibull distribution. Technometrics 2004, 46, 470–481. [Google Scholar] [CrossRef]
  23. Pradhan, B.; Kundu, D. Inference and optimal censoring schemes for progressively censored Birnbaum–Saunders distribution. J. Stat. Plan. Inference 2013, 143, 1098–1108. [Google Scholar] [CrossRef]
  24. Alotaibi, R.; Mutairi, A.; Almetwally, E.M.; Park, C.; Rezk, H. Optimal Design for a Bivariate Step-Stress Accelerated Life Test with Alpha Power Exponential Distribution Based on Type-I Progressive Censored Samples. Symmetry 2022, 14, 830. [Google Scholar] [CrossRef]
  25. Dey, S.; Singh, S.; Tripathi, Y.M.; Asgharzadeh, A. Estimation and prediction for a progressively censored generalized inverted exponential distribution. Stat. Methodol. 2016, 32, 185–202. [Google Scholar]
Figure 1. Generalized type-II progressive hybrid cases.
Figure 1. Generalized type-II progressive hybrid cases.
Symmetry 15 01671 g001
Figure 2. Profile likelihood and contour plot s = 20 p = 0.2.
Figure 2. Profile likelihood and contour plot s = 20 p = 0.2.
Symmetry 15 01671 g002
Figure 3. Profile likelihood and contour plot s = 20 p = 0.5.
Figure 3. Profile likelihood and contour plot s = 20 p = 0.5.
Symmetry 15 01671 g003
Figure 4. Profile likelihood and contour plot s = 20 p = 0.8.
Figure 4. Profile likelihood and contour plot s = 20 p = 0.8.
Symmetry 15 01671 g004
Figure 5. Bayesian checks when s = 20.
Figure 5. Bayesian checks when s = 20.
Symmetry 15 01671 g005
Figure 6. Profile likelihood and contour plot s = 25; p = 0.2.
Figure 6. Profile likelihood and contour plot s = 25; p = 0.2.
Symmetry 15 01671 g006
Figure 7. Profile likelihood and contour plot s = 25; p = 0.5.
Figure 7. Profile likelihood and contour plot s = 25; p = 0.5.
Symmetry 15 01671 g007
Figure 8. Profile likelihood and contour plot s = 25; p = 0.8.
Figure 8. Profile likelihood and contour plot s = 25; p = 0.8.
Symmetry 15 01671 g008
Figure 9. Bayesian checks when s = 25.
Figure 9. Bayesian checks when s = 25.
Symmetry 15 01671 g009
Table 1. The notations of the GPHCS-T2.
Table 1. The notations of the GPHCS-T2.
φ C φ D φ ψ φ T τ ; β R c τ + 1 *
1 j = 1 c 1 i = j s R i + 1 c 1 1 F T 1 R c 1 + 1 * n c 1 i = 1 s 1 R i
2 j = 1 s i = j s R i + 1 s10
3 j = 1 c 2 i = j s R i + 1 c 2 1 F T 2 R c 2 + 1 * n c 2 i = 1 c 2 R i
Table 2. Illustrations of numerous helpful censoring methods and best practices.
Table 2. Illustrations of numerous helpful censoring methods and best practices.
CriterionMethod
O 1 Maximize trace I 2 × 2 .
O 2 Minimize trace I 2 × 2 . 1
O 3 Minimize det I 2 × 2 . 1
O 4 Minimize V a r l o g t ^ p ,   0 < p < 1
Table 3. Bias, MSE, WCI, and CP for parameters and reliability measures: β   =   1.4 ,   θ   =   1.5 .
Table 3. Bias, MSE, WCI, and CP for parameters and reliability measures: β   =   1.4 ,   θ   =   1.5 .
MLEBayesian
nrs BiasMSEWACICPOptimalityBiasMSEWCCI
300.320 β 0.36470.960963.568795.8% O 1 0.8621320.04950.036580.9815
θ 0.30830.307411.807495.3% O 2 0.0504310.06520.01780.6268
R(0.6)−1.07300.010213.568795.8% O 3 24.39956−1.09290.001590.9815
H(0.6)2.60112.155031.807495.3% O 4 0.3972592.27940.158970.6268
R(0.85)−1.31900.003113.568795.8% −1.32690.000260.9815
H(0.85)9.890528.989231.807495.3% 8.25791.306690.6268
25 β 0.31330.685183.004996.0% O 1 0.6032380.03480.011820.6606
θ 0.19310.205281.607595.7% O 2 0.0286220.01960.001630.4069
R(0.6)−1.09620.006443.004996.0% O 3 29.16632−1.10160.000500.6606
H(0.6)2.61771.643011.607595.7% O 4 0.6372782.27920.054060.4069
R(0.85)−1.32780.001903.004996.0% −1.32950.000090.6606
H(0.85)9.668320.976841.607595.7% 8.18520.423060.4069
0.820 β 0.44781.112703.745895.6% O 1 0.8975040.06860.043681.0098
θ 0.28810.299771.826095.3% O 2 0.0545450.05410.008530.6131
R(0.6)−1.09250.008273.745895.6% O 3 23.77178−1.09940.001471.0098
H(0.6)2.77042.385741.826095.3% O 4 0.4109672.32690.178860.6131
R(0.85)−1.32850.002173.745895.6% −1.32950.000251.0098
H(0.85)10.376433.316161.826095.3% 8.37461.535290.6131
25 β 0.30050.659952.960196.2% O 1 0.5884280.03120.008970.6494
θ 0.18920.201031.594295.3% O 2 0.028220.01990.001620.4063
R(0.6)−1.09540.006342.960196.2% O 3 28.83162−1.10090.000410.6494
H(0.6)2.59611.543401.594295.3% O 4 0.2903272.27130.041120.4063
R(0.85)−1.32730.001912.960196.2% −1.32930.000080.6494
H(0.85)9.593920.021471.594295.3% 8.16340.318890.4063
500.335 β 0.23210.338552.092695.2% O 1 0.35240.02190.006820.2936
θ 0.22060.163851.331094.8% O 2 0.01090.03890.004040.1817
R(0.6)0.01500.004910.268494.9% O 3 40.43910.00440.000560.0984
H(0.6)0.26600.867753.501395.4% O 4 0.40100.02300.034280.7479
R(0.85)0.00190.001510.152295.8% 0.00050.000100.0418
H(0.85)1.223310.3147211.646495.4% 0.11730.241801.8061
45 β 0.10090.179891.615795.3% O 1 0.23080.01270.001850.1390
θ 0.06880.072921.024195.5% O 2 0.00510.00800.000430.0697
R(0.6)0.00380.003060.216395.2% O 3 52.6779−0.00090.000130.0468
H(0.6)0.12370.559472.893195.4% O 4 0.38500.02280.009440.3526
R(0.85)0.00230.000990.122895.9% −0.00080.000020.0202
H(0.85)0.54015.799629.204495.3% 0.07290.066400.8489
0.835 β 0.22920.331202.070495.7% O 1 0.34130.03020.007620.2866
θ 0.17100.133531.266694.8% O 2 0.01030.02890.002770.1549
R(0.6)0.00430.004130.251694.5% O 3 40.87150.00030.000490.0905
H(0.6)0.29840.882223.492995.3% O 4 0.55100.04720.036680.6994
R(0.85)−0.00130.001320.142495.8% −0.00100.000090.0382
H(0.85)1.232810.2712711.602395.7% 0.16970.269681.7342
45 β 0.10240.170881.652795.8% O 1 0.23130.01320.001830.1361
θ 0.06140.070911.076795.4% O 2 0.00520.00770.000470.0735
R(0.6)0.00200.003010.219695.0% O 3 52.8830−0.00110.000130.0450
H(0.6)0.12880.556912.915395.4% O 4 0.39530.02420.009340.3352
R(0.85)0.00180.001000.124095.6% −0.00090.000020.0195
H(0.85)0.54935.701399.373595.9% 0.07620.065410.8280
1000.370 β 0.14280.138841.349895.3% O 1 0.14690.01290.002570.1755
θ 0.12690.065020.867495.3% O 2 0.00200.01890.001210.1107
R(0.6)0.00450.002180.182394.6% O 3 80.62940.00160.000210.0575
H(0.6)0.19010.404642.380795.0% O 4 0.61150.01620.013320.4378
R(0.85)−0.00170.000650.099996.1% −0.00010.000040.0249
H(0.85)0.77294.382127.630095.3% 0.07020.091841.0646
90 β 0.06680.079461.074095.7% O 1 0.10650.00620.000540.0852
θ 0.04970.037110.729995.5% O 2 0.00110.00450.000200.0495
R(0.6)0.00130.001620.157795.0% O 3 102.5854−0.00030.000050.0275
H(0.6)0.08980.256121.953394.7% O 4 0.88140.01070.002990.2093
R(0.85)0.00000.000500.088095.0% −0.00040.000010.0120
H(0.85)0.36252.560876.113195.6% 0.03510.019670.5192
0.870 β 0.13090.133181.336195.4% O 1 0.14420.01460.002610.1649
θ 0.10460.060100.869694.9% O 2 0.00200.01570.001070.1047
R(0.6)0.00200.002200.183894.7% O 3 80.59250.00050.000210.0568
H(0.6)0.18030.402882.386895.0% O 4 0.60900.02200.013610.4150
R(0.85)−0.00180.000650.100195.9% −0.00050.000040.0241
H(0.85)0.71254.253017.590195.2% 0.08160.093671.0046
90 β 0.05050.068591.007895.4% O 1 0.10390.00570.000470.0761
θ 0.03180.034040.712794.5% O 2 0.00110.00360.000190.0488
R(0.6)0.00040.001550.154495.2% O 3 103.7077−0.00040.000050.0271
H(0.6)0.06820.230491.863895.7% O 4 0.70990.01030.002660.1933
R(0.85)0.00030.000470.085395.4% −0.00040.000010.0113
H(0.85)0.27432.236595.765995.4% 0.03270.017190.4665
Table 4. Bias, MSE, WCI and CP for parameters and Reliability measures: β   =   1.4 ,   θ   =   0.5 .
Table 4. Bias, MSE, WCI and CP for parameters and Reliability measures: β   =   1.4 ,   θ   =   0.5 .
MLEBayesian
nrs BiasMSEWACICPOptimalityBiasMSEWCCI
300.320 β 0.28830.664022.989295.1% O 1 0.55290.03730.01720.9061
θ 0.09230.030170.577196.5% O 2 0.00450.01910.00100.2005
R(0.6)−1.31420.002732.989295.1% O 3 156.4684−1.32230.00020.9061
H(0.6)4.41983.850750.577196.5% O 4 3.66853.85340.11990.2005
R(0.85)−1.37800.000512.989295.1% −1.38340.00000.9061
H(0.85)11.246228.572080.577196.5% 9.61020.77610.2005
25 β 0.23300.415962.358795.6% O 1 0.36940.02420.00440.6032
θ 0.05650.019120.495095.3% O 2 0.00240.00580.00010.1311
R(0.6)−1.32220.001662.358795.6% O 3 187.5492−1.32440.00010.6032
H(0.6)4.31402.509740.495095.3% O 4 4.25333.82490.03140.1311
R(0.85)−1.38120.000292.358795.6% −1.38390.00000.6032
H(0.85)10.898018.073720.495095.3% 9.52590.20030.1311
0.820 β 0.31120.594762.767495.4% O 1 0.51830.04460.01760.9133
θ 0.08620.029090.577395.3% O 2 0.00420.01680.00090.1972
R(0.6)−1.32070.002182.767495.4% O 3 156.8018−1.32360.00020.9133
H(0.6)4.49023.511180.577395.3% O 4 1.98393.87410.12420.1972
R(0.85)−1.38080.000362.767495.4% −1.38380.00000.9133
H(0.85)11.403525.721570.577395.3% 9.66010.79880.1972
25 β 0.27600.469802.460694.4% O 1 0.38560.03040.00580.6125
θ 0.05600.019510.501995.5% O 2 0.00250.00540.00010.1304
R(0.6)−1.32540.001732.460694.4% O 3 188.0939−1.32510.00010.6125
H(0.6)4.42692.854140.501995.5% O 4 10.41863.84160.04130.1304
R(0.85)−1.38220.000272.460694.4% −1.38420.00000.6125
H(0.85)11.188920.490500.501995.5% 9.56780.26370.1304
500.335 β 0.19280.265871.875595.3% O 1 0.23000.02090.006250.2453
θ 0.06360.014640.403895.4% O 2 0.00100.01110.000340.0531
R(0.6)0.00220.001380.145496.0% O 3 266.35100.00030.000090.0379
H(0.6)0.45331.644294.704495.0% O 4 4.78860.05040.044670.6665
R(0.85)0.00170.000220.058395.2% −0.00010.000010.0118
H(0.85)1.263111.5929712.400895.4% 0.13650.283601.6561
45 β 0.10250.153571.483495.2% O 1 0.15280.01110.001360.1177
θ 0.02520.008560.349295.8% O 2 0.00050.00270.000050.0238
R(0.6)0.00100.000900.117795.9% O 3 332.8660−0.00070.000020.0179
H(0.6)0.24380.975873.754595.7% O 4 9.70830.02840.009810.3228
R(0.85)0.00150.000150.047395.2% −0.00030.000000.0057
H(0.85)0.67696.726019.818995.4% 0.07370.061770.7925
0.835 β 0.22160.300421.966195.3% O 1 0.23600.02860.007330.2601
θ 0.05910.014480.411195.5% O 2 0.00100.00970.000310.0514
R(0.6)−0.00140.001280.140495.7% O 3 264.3238−0.00090.000090.0380
H(0.6)0.52961.843314.902995.3% O 4 6.82330.07150.052060.7134
R(0.85)0.00050.000200.055594.9% −0.00040.000010.0117
H(0.85)1.457113.0984612.993195.3% 0.18860.332391.7468
45 β 0.12510.175001.565694.3% O 1 0.15850.01470.001890.1420
θ 0.02440.007820.333495.3% O 2 0.00050.00230.000040.0228
R(0.6)−0.00080.000970.122496.3% O 3 331.8112−0.00110.000030.0200
H(0.6)0.30531.123583.981194.4% O 4 3.44150.03800.013550.3838
R(0.85)0.00100.000150.048595.8% −0.00040.000000.0064
H(0.85)0.83137.7114210.391694.4% 0.09770.085630.9557
1000.370 β 0.10200.103121.194196.0% O 1 0.09300.00710.001760.1407
θ 0.03900.006420.274595.2% O 2 0.00020.00630.000130.0334
R(0.6)0.00120.000650.099794.9% O 3 526.49830.00060.000030.0220
H(0.6)0.24020.651293.021795.9% O 4 3.01620.01580.012760.3863
R(0.85)0.00070.000100.038695.3% 0.00010.000020.0070
H(0.85)0.66674.498847.897096.1% 0.04510.079910.9471
90 β 0.06610.068700.994895.5% O 1 0.06880.00550.000460.0784
θ 0.01640.003690.229295.4% O 2 0.00010.00140.000020.0160
R(0.6)−0.00050.000480.085895.7% O 3 654.6192−0.00030.000010.0117
H(0.6)0.16050.445072.539795.2% O 4 6.73610.01390.003380.2125
R(0.85)0.00030.000070.033295.6% −0.00010.000010.0038
H(0.85)0.43703.013816.589495.3% 0.03610.021090.5289
0.870 β 0.10360.104801.202895.8% O 1 0.09290.01150.002040.1556
θ 0.03200.006090.279195.3% O 2 0.00020.00500.000110.0339
R(0.6)−0.00020.000640.099495.7% O 3 525.8864−0.00020.000040.0232
H(0.6)0.24840.669143.056795.8% O 4 4.79020.02830.014720.4259
R(0.85)0.00030.000090.038196.1% −0.00020.000020.0073
H(0.85)0.68124.589277.965895.9% 0.07540.092501.0493
90 β 0.06520.072661.025895.7% O 1 0.06890.00600.000520.0794
θ 0.01610.003730.231095.1% O 2 0.00010.00130.000020.0161
R(0.6)−0.00010.000520.089494.8% O 3 650.3684−0.00040.000010.0124
H(0.6)0.15830.473382.626095.6% O 4 6.40730.01530.003790.2188
R(0.85)0.00050.000080.034594.8% −0.00020.000010.0040
H(0.85)0.43173.194736.802595.6% 0.03960.023620.5400
Table 5. Bias, MSE, WCI and CP for parameters and Reliability measures: β   =   0.4 ,   θ   =   0.5 .
Table 5. Bias, MSE, WCI and CP for parameters and Reliability measures: β   =   0.4 ,   θ   =   0.5 .
MLEBayesian
nrs BiasMSEWACICPOptimalityBiasMSEWCCI
300.320 β 0.04770.035720.717396.0% O 1 0.08220.00690.00090.2488
θ 0.16930.102281.064295.4% O 2 0.00110.04540.00660.3352
R(0.6)0.06530.011340.717396.0% O 3 0.08220.03920.00140.2488
H(0.6)1.07360.262081.064295.4% O 4 0.00111.00390.00890.3352
R(0.85)−0.12420.009770.717396.0% −0.14180.00080.2488
H(0.85)3.10861.579951.064295.4% 2.86250.04090.3352
25 β 0.04200.027390.627995.8% O 1 0.05910.00510.00030.1640
θ 0.09380.055900.851194.9% O 2 0.00060.01390.00090.2075
R(0.6)0.04130.008540.627995.8% O 3 170.12430.02910.00030.1640
H(0.6)1.07970.214440.851194.9% O 4 5.01561.00380.00290.2075
R(0.85)−0.13950.007290.627995.8% −0.14780.00020.1640
H(0.85)3.07471.225710.851194.9% 2.84700.01340.2075
0.820 β 0.07470.041260.740896.2% O 1 0.08180.01160.00120.2625
θ 0.16330.094371.020496.1% O 2 0.00110.03780.00510.3251
R(0.6)0.04430.010030.740896.2% O 3 131.11600.03300.00120.2625
H(0.6)1.15780.295061.020496.1% O 4 3.22921.01970.01100.3251
R(0.85)−0.14360.008640.740896.2% −0.14670.00070.2625
H(0.85)3.29181.797181.020496.1% 2.89330.05270.3251
25 β 0.05890.027820.612095.7% O 1 0.05920.00650.00030.1698
θ 0.09910.049710.783495.8% O 2 0.00060.01250.00060.2058
R(0.6)0.02960.007430.612095.7% O 3 159.31680.02750.00030.1698
H(0.6)1.13300.215050.783495.8% O 4 2.22421.00840.00340.2058
R(0.85)−0.15170.006180.612095.7% −0.14910.00020.1698
H(0.85)3.19241.234570.783495.8% 2.85640.01570.2058
500.335 β 0.03690.018390.511895.6% O 1 0.03870.00440.000320.0627
θ 0.10550.043980.710895.9% O 2 0.00020.02290.001650.1006
R(0.6)0.01830.005810.290395.5% O 3 230.42770.00520.000440.0790
H(0.6)0.07950.146271.467295.2% O 4 2.93440.01020.003260.2111
R(0.85)0.00860.005230.281795.6% 0.00210.000260.0627
H(0.85)0.24010.841283.471995.8% 0.03280.014980.4164
45 β 0.02920.012870.430095.8% O 1 0.02840.00310.000090.0311
θ 0.05230.022650.553495.8% O 2 0.00010.00650.000190.0406
R(0.6)0.00300.004140.252094.7% O 3 287.71030.00000.000100.0390
H(0.6)0.07020.108131.259995.7% O 4 3.09680.00880.000920.1034
R(0.85)−0.00030.003740.239794.6% −0.00080.000060.0317
H(0.85)0.18850.594572.932495.7% 0.02200.004230.2099
0.835 β 0.04510.018270.499795.6% O 1 0.03750.00650.000390.0608
θ 0.09160.034380.632395.2% O 2 0.00020.01810.001040.0866
R(0.6)0.00710.005240.282695.7% O 3 223.97690.00180.000380.0752
H(0.6)0.10930.146841.440595.6% O 4 1.67210.01770.003780.2016
R(0.85)−0.00040.004840.272795.5% −0.00050.000240.0604
H(0.85)0.29650.832953.385395.6% 0.04680.017800.4059
45 β 0.02950.013410.439295.7% O 1 0.02800.00340.000100.0323
θ 0.04480.020860.538696.2% O 2 0.00010.00540.000160.0383
R(0.6)0.00030.004150.252593.7% O 3 288.9783−0.00060.000090.0395
H(0.6)0.07230.111321.277496.0% O 4 3.51880.00980.000970.1069
R(0.85)−0.00180.003830.242694.8% −0.00120.000060.0328
H(0.85)0.18960.614932.984295.8% 0.02360.004510.2164
1000.370 β 0.02670.008550.347395.7% O 1 0.01710.00240.000120.0376
θ 0.06900.017830.448495.0% O 2 0.00000.01170.000430.0568
R(0.6)0.00900.002840.206195.8% O 3 442.45350.00260.000160.0489
H(0.6)0.06450.072901.028395.1% O 4 2.99470.00570.001250.1240
R(0.85)0.00200.002630.200895.7% 0.00100.000100.0387
H(0.85)0.17970.401802.384195.5% 0.01770.005680.2530
90 β 0.01290.005360.282794.7% O 1 0.01280.00130.000030.0181
θ 0.02560.009070.359795.0% O 2 0.00000.00290.000060.0233
R(0.6)0.00240.002050.177195.8% O 3 572.98850.00000.000030.0222
H(0.6)0.03100.048030.850894.7% O 4 2.39730.00380.000310.0595
R(0.85)0.00070.001870.169694.9% −0.00030.000020.0179
H(0.85)0.08380.254371.950694.9% 0.00960.001370.1225
0.870 β 0.02740.007950.332895.3% O 1 0.01680.00340.000130.0392
θ 0.05650.015000.426195.1% O 2 0.00000.00980.000340.0558
R(0.6)0.00350.002800.207094.8% O 3 441.28200.00110.000150.0484
H(0.6)0.06930.068690.991394.9% O 4 3.29560.00920.001330.1304
R(0.85)−0.00160.002540.197695.5% −0.00020.000100.0390
H(0.85)0.18340.372932.284594.9% 0.02430.006100.2614
90 β 0.01090.005150.278295.1% O 1 0.01270.00130.000030.0183
θ 0.02120.009040.363595.5% O 2 0.00000.00260.000050.0248
R(0.6)0.00220.002120.180295.2% O 3 577.12180.00000.000030.0236
H(0.6)0.02530.046440.839394.9% O 4 2.56520.00370.000280.0612
R(0.85)0.00110.001910.171294.9% −0.00040.000020.0190
H(0.85)0.06930.245081.922494.9% 0.00920.001250.1231
Table 6. MLE and Bayesian estimation.
Table 6. MLE and Bayesian estimation.
MLEBayesian
EstimatesSER(0.6)R(0.85)EstimatesSER(0.6)R(0.85)
sp H(0.6)H(0.85)H(0.6)H(0.85)
200.2 β 0.88840.52640.32530.11890.98840.32960.30470.1017
θ 1.00330.25142.74806.48841.05560.23652.98117.1010
0.5 β 1.53760.57600.21110.04361.35770.32430.25050.0608
θ 1.22930.26924.191910.42351.23990.17613.77599.3192
0.8 β 1.54040.57510.20910.04311.69210.38770.18000.0324
θ 1.22310.26774.202210.44391.20900.16814.552111.3874
250.2 β 1.49280.49250.18580.03921.55200.30810.17950.0360
θ 1.07760.22034.181910.21391.09660.14644.307410.5745
0.5 β 1.52310.50010.18840.03891.55720.31570.18110.0362
θ 1.11430.22944.229910.38591.10850.15074.312110.6013
0.8 β 1.52210.49960.18830.03891.45110.30530.20440.0451
θ 1.11290.22934.228510.38041.12520.15234.05719.9345
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alotaibi, R.; Almetwally, E.M.; Rezk, H. Reliability Analysis of Kavya Manoharan Kumaraswamy Distribution under Generalized Progressive Hybrid Data. Symmetry 2023, 15, 1671. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15091671

AMA Style

Alotaibi R, Almetwally EM, Rezk H. Reliability Analysis of Kavya Manoharan Kumaraswamy Distribution under Generalized Progressive Hybrid Data. Symmetry. 2023; 15(9):1671. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15091671

Chicago/Turabian Style

Alotaibi, Refah, Ehab M. Almetwally, and Hoda Rezk. 2023. "Reliability Analysis of Kavya Manoharan Kumaraswamy Distribution under Generalized Progressive Hybrid Data" Symmetry 15, no. 9: 1671. https://0-doi-org.brum.beds.ac.uk/10.3390/sym15091671

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop