Nothing Special   »   [go: up one dir, main page]
Next Article in Journal
Physics-Informed Neural Networks for the Heat Equation with Source Term under Various Boundary Conditions
Next Article in Special Issue
Implementation Aspects in Regularized Structural Equation Models
Previous Article in Journal
A Framework for Determining Collision Likelihood Using Continuous Friction Values in a Connected Vehicle Environment
Previous Article in Special Issue
A Mathematical Study on a Fractional-Order SEIR Mpox Model: Analysis and Vaccination Influence
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Algorithms
Volume 16
Issue 9
10.3390/a16090427
algorithms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Confidence Regions for Weibull Parameters and Quantiles under Progressive Censoring

by
Arturo J. Fernández
Departamento de Matemáticas, Estadística e Investigación Operativa and Instituto de Matemáticas y Aplicaciones (IMAULL), Universidad de La Laguna (ULL), 38071 Santa Cruz de Tenerife, Spain
Algorithms 2023, 16(9), 427; https://doi.org/10.3390/a16090427
Submission received: 8 August 2023 / Revised: 21 August 2023 / Accepted: 4 September 2023 / Published: 6 September 2023
(This article belongs to the Special Issue Mathematical Models and Their Applications IV)
Browse Figures
Figure 1
<p>Posterior pdf of <math display="inline"><semantics> <mrow> <mo>(</mo> <mi>θ</mi> <mo>,</mo> <mi>α</mi> <mo>)</mo> </mrow> </semantics></math> in <a href="#sec5-algorithms-16-00427" class="html-sec">Section 5</a>.</p> "> Figure 2
<p>Posterior pdf of <math display="inline"><semantics> <mi>θ</mi> </semantics></math> in <a href="#sec5-algorithms-16-00427" class="html-sec">Section 5</a>.</p> "> Figure 3
<p>Posterior pdf of <math display="inline"><semantics> <mi>α</mi> </semantics></math> in <a href="#sec5-algorithms-16-00427" class="html-sec">Section 5</a>.</p> "> Figure 4
<p>Minimum-area 0.95-confidence region (solid) for <math display="inline"><semantics> <mrow> <mo>(</mo> <mi>θ</mi> <mo>,</mo> <mi>α</mi> <mo>)</mo> </mrow> </semantics></math> and the corresponding region (dashed) proposed by Wu (2002) in the example considered in <a href="#sec5-algorithms-16-00427" class="html-sec">Section 5</a>.</p> "> Figure 5
<p>Smallest-size 0.90-confidence/credibility region for <math display="inline"><semantics> <mrow> <mo>(</mo> <msub> <mi>q</mi> <mrow> <mn>0.05</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>q</mi> <mrow> <mn>0.50</mn> </mrow> </msub> <mo>)</mo> </mrow> </semantics></math> in the example considered in <a href="#sec5-algorithms-16-00427" class="html-sec">Section 5</a>.</p> ">
Versions Notes

Abstract

:
Confidence regions for the Weibull parameters with minimum areas among all those based on the Conditionality Principle are constructed using an equivalent diffuse Bayesian approach. The process is valid for scenarios involving standard failure and progressive censorship, and complete data. Optimal conditional confidence sets for two Weibull quantiles are also derived. Simulation-based algorithms are provided for computing the smallest-area regions with fixed confidence levels. Importantly, the proposed confidence sets satisfy the Sufficiency, Likelihood and Conditionality Principles in contrast to the unconditional regions based on maximum likelihood estimators and other insufficient statistics. The suggested perspective can be applied to parametric estimation and hypothesis testing, as well as to the determination of minimum-size confidence sets for other invariantly estimable functions of the Weibull parameters. A dataset concerning failure times of an insulating fluid is studied for illustrative and comparative purposes.

1. Introduction

Life-testing and reliability experiments are typically terminated before all sample items fail. Under standard failure censoring, the test concludes when a specified number of units have failed; see, e.g., Bhattacharyya [1], LaRiccia [2], Schneider and Weissfeld [3], Fernández [4] and Jaheen and Okasha [5] and references therein. Progressive censorship is an extension of failure censoring in which a certain number of live units can be excluded from the study (i.e., censored) at each failure time. This pattern of censorship offers substantial versatility to the researcher, and also permits the collection of degradation or deterioration data with the objective of analyzing the aging mechanism. The progressive censoring scheme has been widely analyzed in recent decades. Papers by Kemaloglu and Gebizlioglu [6], Wang et al. [7], Lee et al. [8], Almongy et al. [9], Chen and Gui [10] and Abu-Moussa et al. [11] are just a sample. Comprehensive analyses of the state of the art on progressive censorship are provided in the works of Balakrishnan and Aggarwala [12] and Balakrishnan and Cramer [13].
The Weibull distribution with scale parameter θ and shape parameter α is a flexible log-location–scale model for the analysis of time-to-event data that is valuable in many disciplines, including economics, biometry, management, engineering and the actuarial, social and environmental sciences. The Weibull W ( θ , α ) distribution plays a relevant role in many survival and reliability analyses, and has been successfully applied to describe the reliability of both components and equipment in industrial engineering, as well as human and animal disease mortality. Various studies and applications of the Weibull model can be found in Thoman et al. [14], Meeker and Escobar [15], Nordman and Meeker [16], Chen et al. [17], Tsai et al. [18], Fernández [19], Roy [20], Algarni [21], Boult et al. [22], Li et al. [23] and Yu et al. [24]. This model reduces to the exponential distribution when α = 1 ; see, e.g., Fernández et al. [25], Lee et al. [26], Fernández [27], Yousef et al. [28] and Tanackov et al. [29].
The development of confidence sets for Weibull parameters and quantiles from test samples is of great interest and practical relevance in many experimental analyses. In particular, these regions are useful for model selection, parametric estimation and hypothesis testing. In practice, joint confidence sets for the Weibull scale and shape parameters, θ and α , are often based on the unconditional distribution of pivotal quantities related to the maximum likelihood estimator (MLE) of ( θ , α ) , denoted as ( θ ^ , α ^ ) . However, ( θ ^ , α ^ ) does not contain all the sample information. In the Weibull case, the whole sample is minimally sufficient. Confidence regions for ( θ , α ) can also be derived by using the conditional distribution of ( θ ^ , α ^ ) , given the observed values of the ancillary statistics. From a strictly logical point of view, the conditional approach seems more appropriate because it obeys the Sufficient Principle. Uniformly most accurate confidence sets do not exist in the Weibull case. Region size minimization is an alternative optimality criterion that is commonly used to select the best confidence set; see, for example, Casella and Berger [30], Lehmann and Romano [31] and Fernández [32]. At first glance, as smaller confidence sets contain fewer points, they are less likely to cover false values. Interval estimation from progressively censored data has been considered by many authors, including Viveros and Balakrishnan [33], Wu [34], Lawless [35] and Fernández [36].
This paper deals with the construction of minimum-area confidence regions for Weibull parameters and quantiles based on progressively censored data when the principle of conditioning on ancillary statistics is adopted. Our approach is based on the general results of Hora and Buehler [37] for location–scale parameter problems, which are also valid for log-location–scale models, such as the Weibull distribution, because the logarithmic transformation is strictly monotonic. According to Hora and Buehler [37], conditional frequentist confidence sets and Bayesian credibility sets for invariantly estimable functions are numerically equivalent when the analyst assumes independent Lebesgue measure priors on the location parameter and the logarithm of the scale parameter. In our setting, the Weibull parameters and quantiles correspond to invariantly estimable functions, and the proposed optimal confidence regions would include the points with the highest posterior density (HPD) assuming independent flat priors for log ( θ ) and log ( α ) . The above diffuse Bayesian approach satisfies the Conditionality, Likelihood and Sufficiency Principles, and often allows us to substantially reduce the areas of confidence regions.
The remainder of this paper is structured as follows. Given a progressively censored sample from the Weibull model, the next section presents the likelihood function, as well as a diffuse improper prior density for the Weibull parameters and the corresponding posterior density function. Minimum-area confidence regions for the Weibull parameters based on the Conditionality Principle, which coincide with the Bayesian HPD credibility sets in the diffuse case, are derived in Section 3, whereas Section 4 is concerned with the determination of the smallest-size joint confidence set for two arbitrary Weibull quantiles. Algorithms to find the optimal confidence regions via simulation, as well as applications to hypotheses testing, are also suggested. A numerical example regarding failure times for an insulating fluid between two electrodes is considered in Section 5 for illustrative and comparative purposes. Finally, Section 6 offers some concluding remarks.

2. Weibull Models and Progressive Censoring

Suppose that the lifetime X of a certain device follows the Weibull distribution with scale parameter θ > 0 and shape parameter α > 0 , which is denoted as X ∼ W ( θ , α ) . The probability density function (pdf) and cumulative distribution function (cdf) of X are then given by
f ( x ∣ θ , α ) = α x α − 1 / θ α exp x / θ α and F ( x ∣ θ , α ) = 1 − exp − x θ α , x > 0 ,
respectively. Moreover, the reliability or survival function of X is defined as S ( x ∣ θ , α ) = exp − x / θ α and the failure rate is given by h ( x ∣ θ , α ) = α / θ α x α − 1 for x > 0 .
The k-th moment of X ∼ W ( θ , α ) is obtained to be E [ X k ∣ θ , α ] = θ k Γ ( 1 + k / α ) ,   k = 1 , 2 , … , where Γ ( · ) is the well-known gamma function. The parameter θ determines the scaling of the density, whereas the parameter α controls its shape. In many practical cases, the survival of populations with increasing ( α > 1 ) , decreasing ( α < 1 ) , or constant ( α = 1 ) hazard risks can be modeled by Weibull distributions.
Assume now that n randomly selected units from a W ( θ , α ) population with unknown parameters θ and α are put on life test under a progressive censoring scheme r = r 1 , … , r s , and also that x = x 1 , … , x s is the the observed realization of the random sample of failure times X = X 1 , … , X s . That is, n 0 = n units are simultaneously placed on test at time zero in the life-testing experiment; for i = 1 , … , s − 1 ,   r i randomly selected living units are retired from the study at the i-th observed failure time, x i ; so, prior to the i + 1 -th failure, there are n i = n − i − ∑ j = 1 i r j units on inspection; finally, at the time of the s-th observed failure, x s , the test is concluded, i.e., the remaining r s units are removed from the analysis. The constants s and r i ,   i = 1 , … , s , are prefixed integers which must satisfy the assumptions: 1 ≤ s ≤ n ,   0 ≤ r i ≤ n i − 1 − 1 for i = 1 , … , s − 1 , and r s = n s − 1 − 1 .
The likelihood function for ( θ , α ) given x , r , is then defined by
L ( θ , α ∣ x , r ) = ∏ i = 1 s n i − 1 f ( x i ∣ θ , α ) 1 − F ( x i ∣ θ , α ) r i , θ , α > 0 .
In accordance with (1) and (2), the likelihood becomes
L ( θ , α ∣ x , r ) = N w α − 1 α s θ − s α exp v α / θ α , θ , α > 0 ,
where w = W ( x ) and v α = V α ( x , r ) are the observed values of the random quantities
W ≡ W ( X ) = ∏ i = 1 s X i and V α ≡ V α ( X , r ) = ∑ i = 1 s r i + 1 X i α ,
respectively, and N = ∏ i = 1 s n i − 1 . Given the censoring scheme r , it is clear that the observed sample X is minimal sufficient for ( θ , α ) .
Hereafter, it will be assumed that s > 2 and x 1 < x s . Note that the probability that X 1 ≥ X s is zero when s ≥ 2 . In such a case, the unique MLE of ( θ , α ) , denoted by ( θ ^ , α ^ ) , can be derived by solving the equations ∂ log L ( θ , α ∣ x , r ) / ∂ θ = 0 and ∂ log L ( θ , α ∣ x , r ) / ∂ α = 0 via iterative procedures. In our situation, the MLE of ( θ , α ) is an insufficient statistic. Hence, the confidence regions for ( θ , α ) based on ( θ ^ , α ^ ) do not satisfy the Sufficiency Principle.
As mentioned earlier, the Weibull distribution is a log-location-scale parameter model. Specifically, the random variable Y = − log ( X ) has a Gumbel G μ , σ distribution with location and scale parameters μ = − log ( θ ) and σ = 1 / α , respectively. From Hora and Buehler [37], it follows for any level ε ∈ ( 0 , 1 ) that the conditional ε -confidence sets for Gumbel parameters and quantiles are numerically equivalent to the corresponding ε -credibility sets obtained with the improper prior density π 0 ( μ , σ ) ∝   1 / σ ; see also Lawless ([35], p. 565, Property 2). Since the logarithmic transformation is strictly monotonic, the above results are also valid for Weibull parameters and quantiles when the prior pdf is defined by π ( θ , α ) ∝ 1 / ( θ α ) for θ , α > 0 . Assuming this diffuse prior model, the posterior pdf of ( θ , α ) given x , r can be expressed as
π * ( θ , α ∣ x , r ) = w α − 1 θ − s α − 1 α s − 1 K exp v α / θ α , θ , α > 0 ,
where
K ≡ K ( x , r ) = ∫ 0 ∞ ( ∏ i = 1 s x i ) α − 1 α s − 2 ∑ i = 1 s r i + 1 x i α s d α = ∫ 0 ∞ w α − 1 α s − 2 v α s d α .
The posterior pdf of ( θ , α ) given x , r is unimodal, which implies that the HPD estimate (or posterior mode) of ( θ , α ) , denoted by ( θ ˜ , α ˜ ) , is unique. The posterior pdf of α given x , r is given by
π * ( α ∣ x , r ) = w α − 1 α s − 2 K v α s , α > 0 ,
whereas the posterior pdf of t = 1 / θ α conditional on α given x , r is defined as
π * ( t ∣ α ; x , r ) = v α s t s − 1 exp − v α t ( s − 1 ) ! , t > 0 .
Therefore, the posterior distribution of 2 v α / θ α given α and ( x , r ) is chi-square with 2 s degrees of freedom, i.e., 2 v α / θ α ∣ α , ( x , r ) ∼ χ 2 s 2 . The posterior cdf of α conditional to x , r is then given by
G ( z ) = ∫ 0 z w α − 1 α s − 2 K v α s d α , z > 0 .
As graphical illustrations, Figure 1 shows the posterior pdf for ( θ , α ) conditional to x , r associated with the case to be analyzed in Section 5. The corresponding posterior pdfs for θ and α are displayed in Figure 2 and Figure 3.

3. Optimal Regions for the Weibull Parameters

Assume that ε ∈ ( 0 , 1 ) denotes the confidence or credibility level and also that c ε represents the ( 1 − ε ) -quantile of the random posterior density of ( θ , α ) given x , r . The Bayesian HPD ε -credibility region for ( θ , α ) is then defined as
R ε ≡ R ε ( x , r ) = ( θ , α ) : π * ( θ , α ∣ x , r ) ≥ c ε ,
where
Pr ( θ , α ) ∈ R ε ∣ x , r = ∬ R ε π * ( θ , α ∣ x , r ) d θ d α = ε .
The region R ε is also the conditional (frequentist) ε -confidence region for ( θ , α ) with minimum area. The Bayesian credibility degree ε coincides with the frequency-based confidence level of the random region R ε ( X , r ) . Therefore, the smallest-size 100 ε % confidence set for ( θ , α ) based on the Conditionality Principle given ( x , r ) may be defined as
R ε ≡ R ε ( x , r ) = ( θ , α ) : w α − 1 θ − s α − 1 α s − 1 exp − v α / θ α ≥ K c ε .
Since the joint posterior pdf of θ and α derived in (3) is unimodal, it is clear that R ε is a simply connected region. Hence, the set R ε is bounded by a single curve C ε , which does not intersect itself, i.e., the region limited by the contour C ε results in the required ε -confidence set.
An approximate value of c ε can be obtained through simulation. A simple algorithm for determining a random sample of size m from the posterior distribution of ( θ , α ) conditional to x , r can be sketched as follows: Given a large integer number m , for i = 1 , … , m , simulate an observation u i from the uniform distribution U ( 0 , 1 ) and another value d i from the χ 2 s 2 distribution, and then determine α i = G − 1 ( u i ) and θ i = ( 2 v α i / d i ) 1 / α i . In such a case, ( θ 1 , α 1 ) , … , ( θ m , α m ) constitute a random sample of size m from the posterior distribution of ( θ , α ) . Since c ε is the ( 1 − ε ) -quantile of the random variable π * ( θ , α ∣ x , r ) , an approximation of c ε is given by the ( 1 − ε ) -quantile of the simulated sample π * ( θ i , α i ∣ x , r ) i = 1 m .
For interested readers, Thomopoulos [38] focuses on the fundamentals of Monte Carlo methods using basic computer simulation techniques.
The smallest confidence regions presented in this paper can also be applied in hypotheses testing. For instance, if x is the observed value of X and r is the progressive censoring scheme, the p-value associated to the test of the null hypothesis H 0 : ( θ , α ) = ( θ 0 , α 0 ) versus the alternative hypothesis H 1 : ( θ , α ) ≠ ( θ 0 , α 0 ) based on the smallest confidence sets for ( θ , α ) would be defined by p = 1 − ε 0 , where
ε 0 = min ε : θ 0 , α 0 ∈ R ε .
It is easy to show that
R ε 0 = θ , α : π * ( θ , α ∣ x , r ) ≥ π * ( θ 0 , α 0 ∣ x , r ) ,
i.e., the constant c ε 0 equals π * ( θ 0 , α 0 ∣ x , r ) . An approximation of ε 0 is given by the proportion of the simulated sample data π * ( θ i , α i ∣ x , r ) i = 1 m that are at least π * ( θ 0 , α 0 ∣ x , r ) . That is, the p-value is approximately given by the proportion of simulated sample data π * ( θ i , α i ∣ x , r ) i = 1 m that are less than π * ( θ 0 , α 0 ∣ x , r ) . More formally,
p ≈ 1 m ∑ i = 1 m I π * ( θ i , α i ∣ x , r ) < π * ( θ 0 , α 0 ∣ x , r ) ,
where I · denotes the indicator function.

4. Optimal Regions for Two Weibull Quantiles

Given 0 < u < 1 , the Weibull u-quantile is defined as q u = θ − log ( 1 − u ) 1 / α . Our goal in this section is to construct the smallest-size confidence region for the pair of Weibull quantiles ( q a , q b ) , where 0 < a < b < 1 .
The posterior pdf of ( q a , q b ) given x , r can be expressed as
π a , b * ( q a , q b ∣ x , r ) = π * ( θ , α ∣ x , r ) ∂ ( θ , α ) ∂ ( q a , q b ) , 0 < q a < q b ,
where ∂ ( θ , α ) / ∂ ( q a , q b ) denotes the Jacobian (determinant) for the change of variables from ( θ , α ) to ( q a , q b ) .
The Jacobian is defined by
∂ ( θ , α ) ∂ ( q a , q b ) = ∂ θ ∂ q a ∂ θ ∂ q b ∂ α ∂ q a ∂ α ∂ q b = ∂ θ ∂ q a ∂ α ∂ q b − ∂ θ ∂ q b ∂ α ∂ q a ,
where q a = θ − log ( 1 − a ) 1 / α and q b = θ − log ( 1 − b ) 1 / α . After some calculations, it is derived that
∂ ( θ , α ) ∂ ( q a , q b ) = ( k a − k b ) q b k b / ( k a − k b ) q a k a / ( k a − k b ) log ( q a / q b ) 2 ,
where k j = log { − log ( 1 − j ) } for j = a , b , because
∂ θ ∂ q a = − k b k a − k b q b q a k a / ( k a − k b ) , ∂ θ ∂ q b = k a k a − k b q b q a k b / ( k a − k b )
and
∂ α ∂ q a = k b − k a q a log ( q a / q b ) 2 , ∂ α ∂ q b = k a − k b q b log ( q a / q b ) 2 .
As a consequence of the above results, the posterior pdf of ( q a , q b ) given x , r is defined as
π a , b * ( q a , q b ∣ x , r ) = π * q b k a / ( k a − k b ) q a k b / ( k a − k b ) , k a − k b log ( q a / q b ) ∣ x , r ( k b − k a ) q b k b / ( k a − k b ) q a k a / ( k a − k b ) log ( q a / q b ) 2
for 0 < q a < q b . Obviously, ( q a ; i , q b ; i ) i = 1 m , where
q a ; i = θ i − log ( 1 − a ) 1 / α i and q b ; i = θ i − log ( 1 − b ) 1 / α i
for i = 1 , … , m , constitute a random sample of size m from the posterior distribution of ( q a , q b ) conditional to x , r .
In this case, the minimum-area ε -confidence region for ( q a , q b ) , denoted by R a , b ; ε ≡ R a , b ; ε ( x , r ) , is defined as
R a , b ; ε = ( q a , q b ) : π a , b * ( q a , q b ∣ x , r ) ≥ d ε ,
where
Pr ( q a , q b ) ∈ R a , b ; ε ∣ x , r = ∬ R a , b ; ε π a , b * ( q a , q b ∣ x , r ) d q a d q b = ε .
As discussed, the smallest-size confidence regions are relevant to practitioners because they are less likely to contain spurious parameter values.
An approximation of the constant d ε is given by the ( 1 − ε ) -quantile of the simulated sample
π a , b * ( q a ; i , q b ; i ∣ x , r ) i = 1 m
because d ε is precisely the ( 1 − ε ) -quantile of the random variable π a , b * ( q a , q b ∣ x , r ) .
In our situation, the p-value associated to the test H 0 : ( q a , q b ) = ( q a 0 , q b 0 ) against H 1 : ( q a , q b ) ≠ ( q a 0 , q b 0 ) based on the smallest confidence sets for ( q a , q b ) would be defined by p = 1 − ε 0 , where
ε 0 = min ε : ( q a 0 , q b 0 ) ∈ R a , b ; ε .
Hence,
R a , b ; ε 0 = q a , q b : π a , b * ( q a , q b ∣ x , r ) ≥ d ε 0 ,
where d ε 0 = π a , b * ( q a 0 , q b 0 ∣ x , r ) . Furthermore, if an analyst uses the above random sample, it is clear that the p-value for testing H 0 versus H 1 is approximately given by
p ≈ 1 m ∑ i = 1 m I π a , b * ( q a ; i , q b ; i ∣ x , r ) < π a , b * ( q a 0 , q b 0 ∣ x , r ) .
Note that testing H 0 : q a = q a 0 and q b = q b 0 is equivalent to checking H 0 ′ : S ( q a 0 ) = 1 − a and S ( q b 0 ) = 1 − b , which implies that the corresponding reliabilities of the device in study at times q a 0 and q b 0 are 1 − a and 1 − b . For example, H 0 : q 0.01 = 1 and q 0.05 = 2 is identical to the null hypothesis H 0 ′ : S ( 1 ) = 0.99 and S ( 2 ) = 0.95 .

5. Illustrative Applications

A progressively censored sample studied by Balakrishnan and Cramer ([13], p. 9) is considered in this section to illustrate the results developed above. This sample is based on the data reported by Nelson ([39], p. 105) concerning times to breakdown (in minutes) of an insulating fluid between two electrodes subject to a voltage of 34 kV. According to engineering considerations, for a fixed voltage level, the time to breakdown, X , follows a Weibull distribution.
In our case, the progressive censoring scheme is r = ( 0 , 0 , 3 , 0 , 3 , 0 , 0 , 5 ) and the sample of observed failure times is given by:
x = ( 0.19 , 0.78 , 0.96 , 1.31 , 2.78 , 4.85 , 6.50 , 7.35 ) .
Therefore, s = 8 ,   n = 19 and w = W ( x ) = ∏ i = 1 s x i = 120.054 . It can be shown that the maximum likelihood estimates of θ and α are θ ^ = 9.22542 and α ^ = 0.974323 , respectively. Moreover, the value of K ≡ K ( x , r ) is obtained to be K = exp ( − 34.4158 ) .
The posterior pdfs for ( θ , α ) ,   θ and α given x , r are plotted in Figure 1, Figure 2 and Figure 3, respectively. The HPD estimate (or posterior mode) of ( θ , α ) is given by ( θ ˜ , α ˜ ) = ( 8.28923 , 0.930714 ) .
Balakrishnan and Cramer ([13], p. 387) derived that the joint 95% confidence region for θ and α suggested by Wu [34], which is denoted by W 0.95 ≡ W 0.95 ( x , r ) , is defined by
W 0.95 = ( θ , α ) : 0.2807 ≤ α ≤ 1.9648 , 2 v α 31.2070 1 / α ≤ θ ≤ 2 v α 6.0684 1 / α .
Assuming that the confidence level is ε = 0.95 , the optimal (minimum area) 95% confidence region for ( θ , α ) ,   R 0.95 ≡ R 0.95 ( x , r ) , proposed in this paper is given by
R 0.95 = ( θ , α ) : π * ( θ , α ∣ x , r ) ≥ c 0.95 ,
where the 0.05-quantile of the random posterior density of ( θ , α ) given x , r is c 0.95 = 0.00268513 . This set is also the Bayesian HPD 0.95 -credibility region for ( θ , α ) in the noninformative case. For illustrative and comparative purposes, the 95% confidence regions W 0.95 and R 0.95 are depicted in Figure 4.
The area of the optimal region is A r e a [ R 0.95 ] = 33.1901 , whereas A r e a [ W 0.95 ] = 148.137 . Note also that, if α is small, the values of θ such that ( θ , α ) ∈ W 0.95 could be very large. For instance, the point ( 0.2807 , 2014.25 ) is contained in W 0.95 , which is clearly unrealistic because θ ^ = 9.22542 is too small compared to 2014.25 . In general, our approach can greatly reduce the areas of the confidence regions for the Weibull parameters. In the above situation, A r e a [ W 0.95 ] / A r e a [ R 0.95 ] = 4.46412 .
Consider that a reliability engineer aims to check whether ( θ , α ) = ( 10 , 1 ) is reasonable or not. Since the p-value for testing H 0 : ( θ , α ) = ( 10 , 1 ) versus H 1 : ( θ , α ) ≠ ( 10 , 1 ) is calculated to be p = 0.840618 , the values θ = 10 and α = 1 are quite admissible. In contrast, θ = 10 and α = 2 are not reasonable because the p-value for testing H 0 : ( θ , α ) = ( 10 , 2 ) against H 1 : ( θ , α ) ≠ ( 10 , 2 ) is only p = 0.0134487 .
Suppose now that an analyst seeks to determine the smallest-size 90% confidence region for the pair of Weibull quantiles ( q a , q b ) , where a = 0.05 and b = 0.50 . In this case, the minimum-area 90% confidence region for ( q 0.05 , q 0.50 ) ,   R a , b ; 0.90 ≡ R a , b ; 0.90 ( x , r ) , is defined as
R a , b ; 0.90 = ( q a , q b ) : π a , b * ( q a , q b ∣ x , r ) ≥ d 0.90 ,
where d 0.90 = 0.115945 . This set, which is also the diffuse Bayesian HPD 0.90 -credibility region for ( q 0.05 , q 0.50 ) , is depicted in Figure 5.
Evidently, the proposed confidence region can be used to perform hypothesis tests about ( q 0.05 , q 0.50 ) . In particular, the null hypothesis H 0 : ( q 0.05 , q 0.50 ) = ( 0.5 , 7 ) cannot be rejected when ε = 0.90 is the level of confidence. Specifically, the p-value is obtained to be 0.159557 . In contrast, H 0 : ( q 0.05 , q 0.50 ) = ( 0.5 , 8 ) is not acceptable because the p-value is now only 0.0416269 .

6. Concluding Remarks

Optimal joint confidence regions for the scale and shape parameters of the Weibull distribution and two Weibull quantiles are presented in this paper when available data are progressively censored. The proposed confidence sets have minimum area among all those which are based on the Conditionality Principle, and they numerically coincide with the Bayesian highest posterior density credibility sets in the noninformative case.
Smallest-area confidence regions are found by using simulation methods and numerical integration. The suggested approach is valid for both standard failure and progressive censoring, as well as for uncensored samples, and is also applicable to hypothesis testing.
Our methodology obeys the Conditionality, Sufficiency and Likelihood Principles. In contrast, the unconditional methods based on the MLEs and other insufficient statistics violate these principles. In our view, reducing available sample information to insufficient statistics is not appropriate. Moreover, in terms of area, the optimal confidence regions offer appreciable gains over the existing confidence sets. Furthermore, the reduction in area is overwhelming in some cases. In addition, our perspective allows us to construct minimum-size confidence sets for other invariantly estimable functions of the Weibull parameters.

Funding

This research was partially supported by MCIN/AEI/10.13039/501100011033 under grant PID2019-110442GB-I00.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Bhattacharyya, G.K. On asymptotics of maximum likelihood and related estimators based on Type II censored data. J. Am. Stat. Assoc. 1985, 80, 398–404. [Google Scholar] [CrossRef]
  2. LaRiccia, V.N. Asymptotically chi-squared distributed tests of normality for Type II censored samples. J. Am. Stat. Assoc. 1986, 81, 1026–1031. [Google Scholar] [CrossRef]
  3. Schneider, H.; Weissfeld, L. Inference based on Type II censored samples. Biometrics 1986, 42, 531–536. [Google Scholar] [CrossRef] [PubMed]
  4. Fernández, A.J. Highest posterior density estimation from multiply censored Pareto data. Stat. Pap. 2008, 49, 333–341. [Google Scholar] [CrossRef]
  5. Jaheen, Z.F.; Okasha, H.M. E-Bayesian estimation for the Burr type XII model based on type-2 censoring. Appl. Math. Model. 2011, 35, 4730–4737. [Google Scholar] [CrossRef]
  6. Kemaloglu, S.A.; Gebizlioglu, O.L. Risk analysis under progressive type II censoring with binomial claim numbers. J. Comput. Appl. Math. 2009, 233, 61–72. [Google Scholar] [CrossRef]
  7. Wang, B.X.; Yu, K.; Jones, M. Inference under progressively Type II right censored sampling for certain lifetime distributions. Technometrics 2010, 52, 453–460. [Google Scholar] [CrossRef]
  8. Lee, W.-C.; Wu, J.-W.; Hong, M.-L.; Lin, L.-S.; Chan, R.-L. Assessing the lifetime performance index of Rayleigh products based on the Bayesian estimation under progressive type II right censored samples. J. Comput. Appl. Math. 2011, 235, 1676–1688. [Google Scholar] [CrossRef]
  9. Almongy, H.M.; Alshenawy, F.Y.; Almetwally, E.M.; Abdo, D.A. Applying Transformer Insulation Using Weibull Extended Distribution Based on Progressive Censoring Scheme. Axioms 2021, 10, 100. [Google Scholar] [CrossRef]
  10. Chen, S.; Gui, W. Estimation of Unknown Parameters of Truncated Normal Distribution under Adaptive Progressive Type-II Censoring Scheme. Mathematics 2021, 9, 49. [Google Scholar] [CrossRef]
  11. Abu-Moussa, M.H.; Alsadat, N.; Sharawy, A. On Estimation of Reliability Functions for the Extended Rayleigh Distribution under Progressive First-Failure Censoring Model. Axioms 2023, 12, 680. [Google Scholar] [CrossRef]
  12. Balakrishnan, N.; Aggarwala, R. Progressive Censoring: Theory, Methods, and Applications; Springer Science Business Media: Berlin/Heidelberg, Germany; Birkhäuser, Boston: Cambridge, MA, USA, 2000. [Google Scholar]
  13. Balakrishnan, N.; Cramer, E. The Art of Progressive Censoring: Applications to Reliability and Quality; Birkhäuser, Springer: New York, NY, USA, 2014. [Google Scholar]
  14. Thoman, D.R.; Bain, L.J.; Antle, C.E. Maximum likelihood estimation, exact confidence intervals for reliability and tolerance limits in the Weibull distribution. Technometrics 1970, 12, 363–373. [Google Scholar] [CrossRef]
  15. Meeker, W.Q.; Escobar, L.A. Statistical Methods for Reliability Data; Wiley: New York, NY, USA, 1998. [Google Scholar]
  16. Nordman, D.J.; Meeker, W.Q. Weibull prediction intervals for a future number of failures. Technometrics 2002, 44, 15–23. [Google Scholar] [CrossRef]
  17. Chen, J.W.; Li, K.H.; Lam, Y. Bayesian single and double variable sampling plans for the Weibull distribution with censoring. Eur. J. Oper. Res. 2007, 177, 1062–1073. [Google Scholar] [CrossRef]
  18. Tsai, T.-R.; Lu, Y.-T.; Wu, S.-J. Reliability sampling plans for Weibull distribution with limited capacity of test facility. Comput. Ind. Eng. 2008, 55, 721–728. [Google Scholar] [CrossRef]
  19. Fernández, A.J. Optimum attributes component test plans for k-out-of-n:F Weibull systems using prior information. Eur. J. Oper. Res. 2015, 240, 688–696. [Google Scholar] [CrossRef]
  20. Roy, S. Bayesian accelerated life test plans for series systems with Weibull component lifetimes. Appl. Math. Model. 2018, 62, 383–403. [Google Scholar] [CrossRef]
  21. Algarni, A. Group Acceptance Sampling Plan Based on New Compounded Three-Parameter Weibull Model. Axioms 2022, 11, 438. [Google Scholar] [CrossRef]
  22. Boult, T.E.; Windesheim, N.M.; Zhou, S.; Pereyda, C.; Holder, L.B. Weibull-Open-World (WOW) Multi-Type Novelty Detection in CartPole3D. Algorithms 2022, 15, 381. [Google Scholar] [CrossRef]
  23. Li, C.; Günther, M.; Dhamija, A.R.; Cruz, S.; Jafarzadeh, M.; Ahmad, T.; Boult, T.E. Agglomerative Clustering with Threshold Optimization via Extreme Value Theory. Algorithms 2022, 15, 170. [Google Scholar] [CrossRef]
  24. Yu, Q.; Carlson, O.; Sagitov, S. Optimal Maintenance Schedule for a Wind Power Turbine with Aging Components. Algorithms 2023, 16, 334. [Google Scholar] [CrossRef]
  25. Fernández, A.J.; Bravo, J.I.; De Fuentes, Í. Computing maximum likelihood estimates from Type II doubly censored exponential data. Stat. Methods Appl. 2002, 11, 187–200. [Google Scholar] [CrossRef]
  26. Lee, W.-C.; Wu, J.-W.; Hong, C.-W. Assessing the lifetime performance index of products with the exponential distribution under progressively type II right censored samples. J. Comput. Appl. Math. 2009, 231, 648–656. [Google Scholar] [CrossRef]
  27. Fernández, A.J. On calculating generalized confidence intervals for the two-parameter exponential reliability function. Statistics 2007, 41, 129–135. [Google Scholar] [CrossRef]
  28. Yousef, M.M.; Hassan, A.S.; Alshanbari, H.M.; El-Bagoury, A.-A.H.; Almetwally, E.M. Bayesian and Non-Bayesian Analysis of Exponentiated Exponential Stress-–Strength Model Based on Generalized Progressive Hybrid Censoring Process. Axioms 2022, 11, 455. [Google Scholar] [CrossRef]
  29. Tanackov, I.; Prentkovskis, O.; Jevtić, Ž.; Stojić, G.; Ercegovac, P. A New Method for Markovian Adaptation of the Non-Markovian Queueing System Using the Hidden Markov Model. Algorithms 2019, 12, 133. [Google Scholar] [CrossRef]
  30. Casella, G.; Berger, R.L. Statistical Inference; Duxbury/Thomson Learning: Pacific Grove, CA, USA, 2002. [Google Scholar]
  31. Lehmann, E.L.; Romano, J.P. Testing Statistical Hypotheses; Springer: New York, NY, USA, 2005. [Google Scholar]
  32. Fernández, A.J. Smallest Pareto confidence regions and applications. Comput. Stat. Data Anal. 2013, 62, 11–25. [Google Scholar] [CrossRef]
  33. Viveros, R.; Balakrishnan, N. Interval estimation of parameters of life from progressively censored data. Technometrics 1994, 36, 84–91. [Google Scholar] [CrossRef]
  34. Wu, S.J. Estimations of the parameters of the Weibull distribution with progressively censored data. J. Jpn. Stat. Soc. 2002, 32, 155–163. [Google Scholar] [CrossRef]
  35. Lawless, J.F. Statistical Models and Methods for Lifetime Data; Wiley: New York, NY, USA, 2003. [Google Scholar]
  36. Fernández, A.J. Computing optimal confidence sets for Pareto models under progressive censoring. J. Comput. Appl. Math. 2014, 258, 168–180. [Google Scholar] [CrossRef]
  37. Hora, R.B.; Buehler, R.J. Fiducial theory and invariant estimation. Ann. Math. Stat. 1966, 37, 643–656. [Google Scholar] [CrossRef]
  38. Thomopoulos, N.T. Essentials of Monte Carlo Simulation; Springer: New York, NY, USA, 2013. [Google Scholar]
  39. Nelson, W. Applied Life Data Analysis; Wiley: New York, NY, USA, 1982. [Google Scholar]
Algorithms 16 00427 g001
Figure 1. Posterior pdf of ( θ , α ) in Section 5.
Figure 1. Posterior pdf of ( θ , α ) in Section 5.
Algorithms 16 00427 g001
Algorithms 16 00427 g002
Figure 2. Posterior pdf of θ in Section 5.
Figure 2. Posterior pdf of θ in Section 5.
Algorithms 16 00427 g002
Algorithms 16 00427 g003
Figure 3. Posterior pdf of α in Section 5.
Figure 3. Posterior pdf of α in Section 5.
Algorithms 16 00427 g003
Algorithms 16 00427 g004
Figure 4. Minimum-area 0.95-confidence region (solid) for ( θ , α ) and the corresponding region (dashed) proposed by Wu (2002) in the example considered in Section 5.
Figure 4. Minimum-area 0.95-confidence region (solid) for ( θ , α ) and the corresponding region (dashed) proposed by Wu (2002) in the example considered in Section 5.
Algorithms 16 00427 g004
Algorithms 16 00427 g005
Figure 5. Smallest-size 0.90-confidence/credibility region for ( q 0.05 , q 0.50 ) in the example considered in Section 5.
Figure 5. Smallest-size 0.90-confidence/credibility region for ( q 0.05 , q 0.50 ) in the example considered in Section 5.
Algorithms 16 00427 g005
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

Fernández, A.J. Optimal Confidence Regions for Weibull Parameters and Quantiles under Progressive Censoring. Algorithms 2023, 16, 427. https://doi.org/10.3390/a16090427

AMA Style

Fernández AJ. Optimal Confidence Regions for Weibull Parameters and Quantiles under Progressive Censoring. Algorithms. 2023; 16(9):427. https://doi.org/10.3390/a16090427

Chicago/Turabian Style

Fernández, Arturo J. 2023. "Optimal Confidence Regions for Weibull Parameters and Quantiles under Progressive Censoring" Algorithms 16, no. 9: 427. https://doi.org/10.3390/a16090427

APA Style

Fernández, A. J. (2023). Optimal Confidence Regions for Weibull Parameters and Quantiles under Progressive Censoring. Algorithms, 16(9), 427. https://doi.org/10.3390/a16090427

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