Bayesian accelerated life testing under competing log-location-scale family of causes of failure

This article provides Bayesian analyses of data arising from multi-stress accelerated life testing of series systems. The component log-lifetimes are assumed to independently belong to some log-concave location-scale family of distributions. The location parameters are assumed to depend on the stress variables through a linear stress translation function. Bayesian analyses and associated predictive inference of reliability characteristics at usage stresses are performed using Gibbs sampling from the joint posterior. The developed methodology is numerically illustrated by analyzing a real data set through Bayesian model averaging of the two popular cases of Weibull and log-normal, with the later getting a special focus in this article as a slightly easier example of the log-location-scale family. A detailed simulation study is also carried out to compare the performance of various Bayesian point estimators for the log-normal case.     

Bayesian optimal life‐testing plan under the balanced two sample type‐II progressive censoring scheme

Joint progressive censoring schemes are quite useful to conduct comparative life‐testing experiment of different competing products. Recently, Mondal and Kundu (“A New Two Sample Type‐II Progressive Censoring Scheme,” Commun Stat‐Theory Methods; 2018) introduced a joint progressive censoring scheme on two samples known as the balanced joint progressive censoring (BJPC) scheme. Optimal planning of such progressive censoring scheme is an important issue to the experimenter. This article considers optimal life‐testing plan under the BJPC scheme using the Bayesian precision and D‐optimality criteria, assuming that the lifetimes follow Weibull distribution. In order to obtain the optimal BJPC life‐testing plans, one needs to carry out an exhaustive search within the set of all admissible plans under the BJPC scheme. However, for large sample size, determination of the optimal life‐testing plan is difficult by exhaustive search technique. A metaheuristic algorithm based on the variable neighborhood search method is employed for computation of the optimal life‐testing plan. Optimal plans are provided under different scenarios. The optimal plans depend upon the values of the hyperparameters of the prior distribution. The effect of different prior information on optimal scheme is studied.     

Data analysis for accelerated life tests via Weibull-gamma frailty regression models

In this article, we study data analysis methods for accelerated life test (ALT) with blocking. Unlike the previous assumption of normal distribution for random block effects, we advocate the use of Weibull regression model with gamma random effects for making statistical inference of ALT data. To estimate the unknown parameters in the proposed model, maximum likelihood estimation and Bayesian estimation methods are provided. We illustrate the proposed methods using real data examples and simulation examples. Numerical results suggest that distribution of random effects has minimal impact on the estimation of fixed effects in the Weibull regression models. Furthermore, to demonstrate the advantage of our proposed model, we also provide methods to compare ALT plans and thus identify the optimal ALT plans.     

一般寿命分布和定时截尾的Bayes变量抽样方案

林（1994）研究了指数分布和定时截尾的变量抽样方案．本文将讨论一般寿命分布和定时截尾的一次抽样方案．在多项式损失函数的假设下,我们讨论了Weibull分布、双参数指数分布和-分布三种情形,并着重讨论Weibull分布的情形．本文还提出了一个可用于近似地确定最优抽样方案的有报算法,并且进行了灵敏度分析,还同林较早的模型（1990,1994）做了比较．     

A load share model for non-identical components of a k-out-of-m system

We consider a k-out-of-m load sharing system when the lifetimes of the components are not necessarily identically distributed random variables. For such systems, a model for the load sharing phenomenon through the exponentiated conditional survival functions of ordered failure times is proposed. This model is more general than the load sharing model with identically distributed component lifetimes and leads to a different family of distributions for ordered random variables. A general expression for the reliability of the system is given. The computations of the reliability for a two component parallel load sharing system corresponding to the exponential and Weibull distributions are discussed. For illustrative purpose, we discuss the inference procedures for a two component parallel load sharing system corresponding to the exponential distributions. A simulation study is carried out to assess the proposed estimation and testing procedures. The applicability of the proposed load sharing model is shown through two data sets.     

Two sample Bayesian prediction intervals for order statistics based on the inverse exponential-type distributions using right censored sample

In this paper, two sample Bayesian prediction intervals for order statistics (OS) are obtained. This prediction is based on a certain class of the inverse exponential-type distributions using a right censored sample. A general class of prior density functions is used and the predictive cumulative function is obtained in the two samples case. The class of the inverse exponential-type distributions includes several important distributions such the inverse Weibull distribution, the inverse Burr distribution, the loglogistic distribution, the inverse Pareto distribution and the inverse paralogistic distribution. Special cases of the inverse Weibull model such as the inverse exponential model and the inverse Rayleigh model are considered.     

Bayesian planning of step‐stress accelerated degradation tests under various optimality criteria

Step‐stress accelerated degradation testing (SSADT) has become a common approach to predicting lifetime for highly reliable products that are unlikely to fail in a reasonable time under use conditions or even elevated stress conditions. In literature, the planning of SSADT has been widely investigated for stochastic degradation processes, such as Wiener processes and gamma processes. In this paper, we model the optimal SSADT planning problem from a Bayesian perspective and optimize test plans by determining both stress levels and the allocation of inspections. Large‐sample approximation is used to derive the asymptotic Bayesian utility functions under 3 planning criteria. A revisited LED lamp example is presented to illustrate our method. The comparison with optimal plans from previous studies demonstrates the necessity of considering the stress levels and inspection allocations simultaneously.     

A stochastic expectation-maximization algorithm for the analysis of system lifetime data with known signature

Statistical estimation of the model parameters of component lifetime distribution based on system lifetime data with known system structure is discussed here. We propose the use of stochastic expectation-maximization (SEM) algorithm for obtaining the maximum likelihood estimates of model parameters based on complete and censored system lifetimes. Different ways of implementing the SEM algorithm are also studied. We have shown that the proposed methods are feasible and are easy to implement for various families of component lifetime distributions. The proposed methodologies are then illustrated with two popular lifetime models—the Weibull and Birnbaum-Saunders distributions. Monte Carlo simulation is then used to compare the performance of the proposed methods with the corresponding estimation by direct maximization. Finally, two illustrative examples are presented along with some concluding remarks.     

A Comparison of Two Sequential Metropolis-Hastings Algorithms with Standard Simulation Techniques in Bayesian Inference in Reliability Models Involving a Generalized Gamma Distribution

In this paper we consider the generalized gamma distribution as introduced in Gåsemyr and Natvig (1998). This distribution enters naturally in Bayesian inference in exponential survival models with left censoring. In the paper mentioned above it is shown that the weighted sum of products of generalized gamma distributions is a conjugate prior for the parameters of component lifetimes, having autopsy data in a Marshall-Olkin shock model. A corresponding result is shown in Gåsemyr and Natvig (1999) for independent, exponentially distributed component lifetimes in a model with partial monitoring of components with applications to preventive system maintenance. A discussion in the present paper strongly indicates that expressing the posterior distribution in terms of the generalized gamma distribution is computationally efficient compared to using the ordinary gamma distribution in such models. Furthermore, we present two types of sequential Metropolis-Hastings algorithms that may be used in Bayesian inference in situations where exact methods are intractable. Finally these types of algorithms are compared with standard simulation techniques and analytical results in arriving at the posterior distribution of the parameters of component lifetimes in special cases of the mentioned models. It seems that one of these types of algorithms may be very favorable when prior assessments are updated by several data sets and when there are significant discrepancies between the prior assessments and the data.     

Quadrature Methods for Bayesian Optimal Design of Experiments With Nonnormal Prior Distributions

Many optimal experimental designs depend on one or more unknown model parameters. In such cases, it is common to use Bayesian optimal design procedures to seek designs that perform well over an entire prior distribution of the unknown model parameter(s). Generally, Bayesian optimal design procedures are viewed as computationally intensive. This is because they require numerical integration techniques to approximate the Bayesian optimality criterion at hand. The most common numerical integration technique involves pseudo Monte Carlo draws from the prior distribution(s). For a good approximation of the Bayesian optimality criterion, a large number of pseudo Monte Carlo draws is required. This results in long computation times. As an alternative to the pseudo Monte Carlo approach, we propose using computationally efficient Gaussian quadrature techniques. Since, for normal prior distributions, suitable quadrature techniques have already been used in the context of optimal experimental design, we focus on quadrature techniques for nonnormal prior distributions. Such prior distributions are appropriate for variance components, correlation coefficients, and any other parameters that are strictly positive or have upper and lower bounds. In this article, we demonstrate the added value of the quadrature techniques we advocate by means of the Bayesian D-optimality criterion in the context of split-plot experiments, but we want to stress that the techniques can be applied to other optimality criteria and other types of experimental designs as well. Supplementary materials for this article are available online.     

Bayesian single and double variable sampling plans for the Weibull distribution with censoring

The sampling inspection problem is one of the main research topics in quality control. In this paper, we employ Bayesian decision theory to study single and double variable sampling plans, for the Weibull distribution, with Type II censoring. A general loss function which includes the sampling cost, the time-consuming cost, the salvage value, and the after-sales cost is proposed to determine the Bayes risk and the corresponding optimal sampling plan. Explicit expressions for the Bayes risks for both single and double sampling plans are derived, respectively. Numerical examples are given to illustrate the effectiveness of the proposed method. Comparisons between single and double sampling plans are made, and sensitivity analysis is performed.     

Weibull分布恒定应力加速寿命试验的Bayes估计

对定时和定数截尾样本情形CE模型下Weibull分布场合恒定应力加速寿命试验进行了Bayes统计分析,利用Laplace方法给出了该模型的近似Bayes估计.最后通过模拟实例表明该Bayes估计是有效的.     

Bayesian planning of optimal step-stress accelerated life test for log-location-scale distributions

This paper introduces some Bayesian optimal design methods for step-stress accelerated life test planning with one accelerating variable, when the acceleration model is linear in the accelerated variable or its function, based on censored data from a log-location-scale distributions. In order to find the optimal plan, we propose different Monte Carlo simulation algorithms for different Bayesian optimal criteria. We present an example using the lognormal life distribution with Type-I censoring to illustrate the different Bayesian methods and to examine the effects of the prior distribution and sample size. By comparing the different Bayesian methods we suggest that when the data have large(small) sample size B₁(τ) (B₂(τ)) method is adopted. Finally, the Bayesian optimal plans are compared with the plan obtained by maximum likelihood method.     

CE模型下Weibull分布序加试验的Bayes分析

序进应力加速寿命试验是一种最为有效而经济的寿命试验方法,随着其理论的日趋成熟,在实践中开始得到应用和推广．本文给出了逆幂律模型下Weilbull分布定时和定数场合序进应力加速寿命试验的一种Bayes统计分析,并利用Gibbs抽样方法解决了分布的形状参数取为连续先验时各参数的Bayes估计．这种先验意义更明确,实例表明这是一种非常有效的方法．     

Statistical inference under adaptive progressive censoring scheme

In this paper, a general exponential form of the underlying distribution and a general conjugate prior are used to discuss the maximum likelihood and Bayesian estimation based on an adaptive progressive censored sample. A general procedure for deriving the point and interval Bayesian prediction of the future progressive censored from the same sample as well as that from an unobserved future sample is also developed. The Weibull, Pareto, and Burr Type-XII distributions are then used as illustrative examples. Finally, two numerical examples are presented for illustrating all the inferential procedures developed here.     

Statistical inference for coherent systems with Weibull distributed component lifetimes under complete and incomplete information

Point estimators for the parameters of the component lifetime distribution in coherent systems are evolved assuming to be independently and identically Weibull distributed component lifetimes. We study both complete and incomplete information under continuous monitoring of the essential component lifetimes. First, we prove that the maximum likelihood estimator (MLE) under complete information based on progressively Type‐II censored system lifetimes uniquely exists and we present two approaches to compute the estimates. Furthermore, we consider an ad hoc estimator, a max‐probability plan estimator and the MLE for the parameters under incomplete information. In order to compute the MLEs, we consider a direct maximization of the likelihood and an EM‐algorithm–type approach, respectively. In all cases, we illustrate the results by simulations of the five‐component bridge system and the 10‐component parallel system, respectively.     

Design of progressively censored group sampling plans for Weibull distributions: An optimization problem

Optimization algorithms provides efficient solutions to many statistical problems. Essentially, the design of sampling plans for lot acceptance purposes is an optimization problem with several constraints, usually related to the quality levels required by the producer and the consumer. An optimal acceptance sampling plan is developed in this paper for the Weibull distribution with unknown scale parameter. The proposed plan combines grouping of items, sudden death testing in each group and progressive group removals, and its decision criterion is based on the uniformly most powerful life test. A mixed integer programming problem is first solved for determining the minimum number of failures required and the corresponding acceptance constant. The optimal number of groups is then obtained by minimizing a balanced estimation of the expected test cost. Excellent approximately optimal solutions are also provided in closed-forms. The sampling plan is considerably flexible and allows to save experimental time and cost. In general, our methodology achieves solutions that are quite robust to small variations in the Weibull shape parameter. A numerical example about a manufacturing process of gyroscopes is included for illustration.     

Classical versus Bayesian risks in acceptance sampling: a sensitivity analysis

Assuming a beta prior distribution on the fraction defective, $p$ , failure-censored sampling plans for Weibull lifetime models using classical (or average) and Bayesian (or posterior) producer’s and consumer’s risks are designed to determine the acceptability of lots of a given product. The average risk criterion provides a certain assurance that good (bad) lots will be accepted (rejected), whereas the posterior risk criterion provides a determined confidence that an accepted (rejected) lot is indeed good (bad). The performance of classical and Bayesian risks are analyzed in developing sampling plans when the lifetime variable follows the Weibull distribution. Several figures and tables illustrate the sensitivity of the risks and optimal sample sizes for selected censoring levels and specifications according to the available prior information on $p$ . The analysis clarifies the distinction among the different risks for a given sampling plan, and the effect of the prior knowledge on the required sample size. The study shows that, under uncertainty in the prior variance of $p$ , the designs using Bayesian risks are more appropriate.     

Some results on information properties of coherent systems

This paper considers information properties of coherent systems when component lifetimes are independent and identically distributed. Some results on the entropy of coherent systems in terms of ordering properties of component distributions are proposed. Moreover, various sufficient conditions are given under which the entropy order among systems as well as the corresponding dual systems hold. Specifically, it is proved that under some conditions, the entropy order among component lifetimes is preserved under coherent system formations. The findings are based on system signatures as a useful measure from comparison purposes. Furthermore, some results on the system's entropy are derived when lifetimes of components are dependent and identically distributed. Several illustrative examples are also given.     

Bayesian analysis of fixed-interval preventive-maintenance models

This article develops methods for making accurate decisionswhen scheduling preventive maintenance in systems where inter-eventtimes can be modelled by a delayed renewal process or delayedalternating renewal process. A practical application, relatingto the reliability and maintenance of a relatively low-levelcomponent (valve) in a continuous-process industry over a periodof six years is presented to demonstrate and compare the differentapproaches. Our analyses indicate a cost-effective recommendationfor maintenance practice in this context. Our main thrust relates to the use of Bayesian methodology inorder to obtain rational, admissible decisions. Particular advancesover previous research allow for informative prior distributions,better approximations which lead to improved accuracy, non-negligibledowntimes, and general lifetime distributions. General analyticsolutions are sought for the simpler models, in order to achieveaccuracy and insight. The resulting integrals can only be solvedto give an infinite series and one approximation to the requiredsolution is obtained by truncating this series. Two other approximationsare developed, based on expansions of the prior predictive andlog-posterior distributions. A simulation approach is also developed to include prior informationand hence provide alternative approximations to these optimaldecisions. With exponential lifetime distributions, the relevantposterior lifetime distributions are non-central Pareto. Thissimulation is simple to program, compared to the approximations,but requires more computing time. It is accurate and extendseasily to situations involving greater complexity. We considertwo such extensions, the inclusion of downtimes and Weibulllifetime distributions.