Bayesian Analysis of Masked Series System Lifetime Data

The problem of analyzing series system lifetime data with masked or partial information on cause of failure is recent, compared to that of the standard competing risks model. A generic Gibbs sampling scheme is developed in this article towards a Bayesian analysis for a general parametric competing risks model with masked cause of failure data. The masking probabilities are not subjected to the symmetry assumption and independent Dirichlet priors are used to marginalize these nuisance parameters. The developed methodology is illustrated for the case where the components of a series system have independent log-Normal life distributions by employing independent Normal-Gamma priors for these component lifetime parameters. The Gibbs sampling scheme developed for the required analysis can also be used to provide a Bayesian analysis of data arising from the conventional competing risks model of independent log-Normals, which interestingly has so far remained by and large neglected in the literature. The developed methodology is deployed to analyze a masked lifetime data of PS/2 computer systems.     

Maximum likelihood analysis of masked series system lifetime data

Component lifetime parameters of a series system are estimated from system lifetimes and masked cause of failure observations. The time and cause of system failures are assumed to follow a competing risks model. The masking probabilities of the minimum random subsets are not subjected to the symmetry assumption. Sufficient regularity conditions are provided, justifying the maximum likelihood analysis. Maximum likelihood estimates of both the lifetime parameters and masking probabilities are generically computed via an EM algorithm. An appropriate set of asymptotically pivotal quantities are also derived. Such maximum likelihood based estimates are further refined by bootstrap. The developed techniques are illustrated by numerical examples of independent Weibull component lifetimes with distinct scale and shape parameters.     

Consistency of the Generalized MLE with Interval-Censored and Masked Competing Risks Data

We consider nonparametric estimation based on interval-censored competing risks data with masked failure cause. The generalized maximum likelihood estimator of the joint survival function of the failure time and the failure cause is studied under mixed case interval censorship and random partition masking. Strong consistency in the L ₁(μ)-topology is established for some finite measure μ which is derived from the joint censoring and masking distribution. Under additional regularity assumptions we also establish the strong consistencies in the topologies of weak convergence, point-wise convergence, and uniform convergence.     

The random partition masking model for interval-censored and masked competing risks data

We study the nonparametric maximum likelihood estimate (NPMLE) of the cdf or sub-distribution functions of the failure time for the failure causes in a series system. The study is motivated by a cancer research data (from the Memorial Sloan-Kettering Cancer Center) with interval-censored time and masked failure cause. The NPMLE based on this data set suggests that the existing masking models are not appropriate. We propose a new model called the random partition masking model, which does not rely on the commonly used symmetry assumption (namely, given the failure cause, the probability of observing the masked failure causes is independent of the failure time; see Flehinger et al. Inference about defects in the presence of masking, Technometrics 38 (1996), pp. 247–255). The RPM model is easier to implement in simulation studies than the existing models. We discuss the algorithms for computing the NPMLE and study its asymptotic properties. Our simulation and data analysis indicate that the NPMLE is feasible for a moderate sample size.     

ANALYSIS OF MASKED DATA IN A DEPENDENT COMPETING RISKS MODEL UNDER UNKNOWN ENVIRONMENT

ABSTRACT

System failure data is often analyzed to estimate component reliabilities. Due to cost and time constraints, the exact component causing the failure of the system cannot be identified in some cases. This phenomenon is called masking. Further, it is sometimes necessary for us to take account of the influence of the operating environment. Here we consider a series system, operating under unknown environment, of two components whose failure times follow the Marshall-Olkin bivariate exponential distribution. We present a maximum likelihood approach for obtaining estimators from the masked data for this system. From a simulation study, we found that the relative errors of the estimates are almost well behaved even for small or moderate expected number of systems whose cause of failure is identified.     

Bayesian analysis of competing risks with partially masked cause of failure

Summary. Bayesian analysis of system failure data from engineering applications under a competing risks framework is considered when the cause of failure may not have been exactly identified but has only been narrowed down to a subset of all potential risks. In statistical literature, such data are termed masked failure data. In addition to masking, failure times could be right censored owing to the removal of prototypes at a prespecified time or could be interval censored in the case of periodically acquired readings. In this setting, a general Bayesian formulation is investigated that includes most commonly used parametric lifetime distributions and that is sufficiently flexible to handle complex forms of censoring. The methodology is illustrated in two engineering applications with a special focus on model comparison issues.     

Control chart for monitoring the Weibull shape parameter under two competing risks

ABSTRACT

In this paper, we propose a control chart to monitor the Weibull shape parameter where the observations are censored due to competing risks. We assume that the failure occurs due to two competing risks that are independent and follow Weibull distribution with different shape and scale parameters. The control charts are proposed to monitor one or both of the shape parameters of competing risk distributions and established based on the conditional expected values. The proposed control chart for both shape parameters is used in certain situations and allows to monitor both shape parameters in only one chart. The control limits depend on the sample size, number of failures due to each risk and the desired stable average run length (ARL). We also consider the estimation problem of the target parameters when the Phase I sample is incomplete. We assumed that some of the products that fail during the life testing have a cause of failure that is only known to belong to a certain subset of all possible failures. This case is known as masking. In the presence of masking, the expectation-maximization (EM) algorithm is proposed to estimate the parameters. For both cases, with and without masking, the behaviour of ARLs of charts is studied through the numerical methods. The influence of masking on the performance of proposed charts is also studied through a simulation study. An example illustrates the applicability of the proposed charts.     

Parametric Modeling for Survival with Competing Risks and Masked Failure Causes

We consider a life testing situation in which systems are subject to failure from independent competing risks. Following a failure, immediate (stage-1) procedures are used in an attempt to reach a definitive diagnosis. If these procedures fail to result in a diagnosis, this phenomenon is called masking. Stage-2 procedures, such as failure analysis or autopsy, provide definitive diagnosis for a sample of the masked cases. We show how stage-1 and stage-2 information can be combined to provide statistical inference about (a) survival functions of the individual risks, (b) the proportions of failures associated with individual risks and (c) probability, for a specified masked case, that each of the masked competing risks is responsible for the failure. Our development is based on parametric distributional assumptions and the special case for which the failure times for the competing risks have a Weibull distribution is discussed in detail.     

Bayesian Analysis of Masked Data in Step-stress Accelerated Life Testing

This article considers a k level step-stress accelerated life testing (ALT) on series system products, where independent Weibull-distributed lifetimes are assumed for the components. Due to cost considerations or environmental restrictions, causes of system failures are masked and type-I censored observations might occur in the collected data. Bayesian approach combined with auxiliary variables is developed for estimating the parameters of the model. Further, the reliability and hazard rate functions of the system and components are estimated at a specified time at use stress level. The proposed method is illustrated through a numerical example based on two priors and various masking probabilities.     

Inverse probability weighted estimators for single-index models with missing covariates

Abstract

In this article, we consider the inverse probability weighted estimators for a single-index model with missing covariates when the selection probabilities are known or unknown. It is shown that the estimator for the index parameter by using estimated selection probabilities has a smaller asymptotic variance than that with true selection probabilities, thus is more efficient. Therefore, the important Horvitz-Thompson property is verified for the index parameter in single index model. However, this difference disappears for the estimators of the link function. Some numerical examples and a real data application are also conducted to illustrate the performances of the estimators.     

Dependent masking and system life data analysis: Bayesian inference for two-component systems

Data from field operations of a system is often used to estimate the reliability of components. Under ideal circumstances, this system field data contains the time to failure along with information on the exact component responsible for the system failure. However, in many cases, the exact component causing the failure of the system cannot be identified, and is considered to be masked. Previously developed models for estimation of component reliability from masked system life data have been based upon the assumption that masking occurs independently of the true cause of system failure. In this paper we develop a Bayesian methodology for estimating component reliabilities from masked system life data when the probability of masking is dependent upon the true cause of system failure. The Bayesian approach is illustrated for the case of a two-component system of exponentially distributed components.     

Robust Bayesian analysis for parallel system with masked data under inverse Weibull lifetime distribution

Abstract

This paper considers the statistical analysis of masked data in a parallel system with inverse Weibull distributed components under type II censoring. Based on Gamma conjugate prior, the Bayesian estimation as well as the hierarchical Bayesian estimation for the parameters and the reliability function of system are obtained by using the Bayesian theory and the hierarchical Bayesian method. Finally, Monte Carlo simulations are provided to compare the performances of the estimates under different masking probabilities and effective sample sizes.     

Assessing component reliability using lifetime data from systems

ABSTRACT

In this paper, we introduce a competing risks model for the lifetimes of components that differs from the classical competing risks models by the fact that it is not directly observable which component has failed. We propose two statistical methods for estimating the reliability of components from failure data on a system. Our methods are applied to simulated failure data, in order to illustrate the performance of the methods.     

A Note on the Lower Bounds on Maximum Number of Clear Two-Factor Interactions for (2 m−p III) and (2 m−p IV) Designs

ABSTRACT

In the reliability analysis of mechanical repairable equipment subjected to reliability deterioration with operating time, two forms of the non-homogeneous Poisson processes, namely the Power-Law (PL) and the Log-Linear (LL) model, have found general acceptance in the literature. Inferential procedures, conditioned on the assumption of the PL or LL model, underestimate the overall uncertainty about a quantity of interest because the PL and LL models can provide different estimates of the quantity of interest, even when both of them adequately fit the observed data. In this paper, a composite estimation procedure, which uses the PL and LL models as competing models, is proposed in the framework of Bayesian statistics, thus allowing the uncertainty involved in model selection to be considered. A model-free approach is then proposed for incorporating technical information on the failure mechanism into the inferential procedure. Such an approach, which is based on two model-free quantities defined irrespectively of the functional form of the failure model, prevents that the prior information on the failure mechanism can improperly introduce prior probabilities on the adequacy of each model to fit the observed data. Finally, numerical applications are provided to illustrate the proposed procedures.     

Shared frailty models with baseline generalized Pareto distribution

Abstract

In this article, we have considered three different shared frailty models under the assumption of generalized Pareto Distribution as baseline distribution. Frailty models have been used in the survival analysis to account for the unobserved heterogeneity in an individual risks to disease and death. These three frailty models are with gamma frailty, inverse Gaussian frailty and positive stable frailty. Then we introduce the Bayesian estimation procedure using Markov chain Monte Carlo (MCMC) technique to estimate the parameters. We applied these three models to a kidney infection data and find the best fitted model for kidney infection data. We present a simulation study to compare true value of the parameters with the estimated values. Model comparison is made using Bayesian model selection criterion and a well-fitted model is suggested for the kidney infection data.     

Constant-partially accelerated life tests for Marshall–Olkin exponential series system with dependent masked data

This paper considers the constant-partially accelerated life tests for series system products, where dependent M-O bivariate exponential distribution is assumed for the components.

Based on progressive type-II censored and masked data, the maximum likelihood estimates for the parameters and acceleration factors are obtained by using the decomposition approach. In addition, this method can also be applied to the Bayes estimates, which are too complex to obtain as usual way. Finally, a Monte Carlo simulation study is carried out to verify the accuracy of the methods under different masking probabilities and censoring schemes.     

Perturbation analysis of Markov modulated fluid models

ABSTRACT

We consider perturbations of positive recurrent Markov modulated fluid models. In addition to the infinitesimal generator of the phases, we also perturb the rate matrix, and analyze the effect of those perturbations on the matrix of first return probabilities to the initial level. Our main contribution is the construction of a substitute for the matrix of first return probabilities, which enables us to analyze the effect of the perturbation under consideration.     

Randomized response model versus mixture model in the crossover design for clinical trials

ABSTRACT

The clinical trials are usually designed with the implicit assumption that data analysis will occur only after the trial is completed. It is a challenging problem if the sponsor wishes to evaluate the drug efficacy in the middle of the study without breaking the randomization codes. In this article, the randomized response model and mixture model are introduced to analyze the data, masking the randomization codes of the crossover design. Given the probability of treatment sequence, the test of mixture model provides higher power than the test of randomized response model, which is inadequate in the example. The paired t-test has higher powers than both models if the randomization codes are broken. The sponsor may stop the trial early to claim the effectiveness of the study drug if the mixture model concludes a positive result.     

Inference for a series system with dependent masked data under progressive interval censoring

This paper considers the statistical analysis of masked data in a series system with Burr-XII distributed components. Based on progressively Type-I interval censored sample, the maximum likelihood estimators for the parameters are obtained by using the expectation maximization algorithm, and the associated approximate confidence intervals are also derived. In addition, Gibbs sampling procedure using important sampling is applied for obtaining the Bayesian estimates of the parameters, and Monte Carlo method is employed to construct the credible intervals. Finally, a simulation study is proposed to illustrate the efficiency of the methods under different removal schemes and masking probabilities.     

Fitting time series models for longitudinal survey data under informative sampling

The purpose of this paper is to account for informative sampling in fitting time series models, and in particular an autoregressive model of order one, for longitudinal survey data. The idea behind the proposed approach is to extract the model holding for the sample data as a function of the model in the population and the first-order inclusion probabilities, and then fit the sample model using maximum-likelihood, pseudo-maximum-likelihood and estimating equations methods. A new test for sampling ignorability is proposed based on the Kullback–Leibler information measure. Also, we investigate the issue of the sensitivity of the sample model to incorrect specification of the conditional expectations of the sample inclusion probabilities. The simulation study carried out shows that the sample-likelihood-based method produces better estimators than the pseudo-maximum-likelihood method, and that sensitivity to departures from the assumed model is low. Also, we find that both the conventional t-statistic and the Kullback–Leibler information statistic for testing of sampling ignorability perform well under both informative and noninformative sampling designs.