Statistical inference for dependent competing risks data under adaptive Type-II progressive hybrid censoring

In this article, we consider statistical inference based on dependent competing risks data from Marshall-Olkin bivariate Weibull distribution. The maximum likelihood estimates of the unknown model parameters have been computed by using the Newton-Raphson method under adaptive Type II progressive hybrid censoring with partially observed failure causes. The existence and uniqueness of maximum likelihood estimates are derived. Approximate confidence intervals have been constructed via the observed Fisher information matrix using the asymptotic normality property of the maximum likelihood estimates. Bayes estimates and highest posterior density credible intervals have been calculated under gamma-Dirichlet prior distribution by using the Markov chain Monte Carlo technique. Convergence of Markov chain Monte Carlo samples is tested. In addition, a Monte Carlo simulation is carried out to compare the effectiveness of the proposed methods. Further, three different optimality criteria have been taken into account to obtain the most effective censoring plans. Finally, a real-life data set has been analyzed to illustrate the operability and applicability of the proposed methods

Optimal Test Plan of Step Stress Partially Accelerated Life Testing for Alpha Power Inverse Weibull Distribution under Adaptive Progressive Hybrid Censored Data and Different Loss Functions

Accelerated life tests are used to explore the lifetime of extremely reliable items by subjecting them to elevated stress levels from stressors to cause early failures, such as temperature, voltage, pressure, and so on. The alpha power inverse Weibull (APIW) distribution is of great significance and practical applications due to its appealing characteristics, such as its flexibilities in the probability density function and the hazard rate function. We analyze the step stress partially accelerated life testing model with samples from the APIW distribution under adaptive type II progressively hybrid censoring. We first obtain the maximum likelihood estimates and two types of approximate confidence intervals of the distributional parameters and then derive Bayes estimates of the unknownparameters under different loss functions. Furthermore, we analyze three probable optimum test techniques for identifying the best censoring under different optimality criteria methods. We conduct simulation studies to assess the finite sample performance of the proposed methodology. Finally, we provide a real data example to further demonstrate the proposed technique

On the maximum likelihood estimation for progressively censored lifetimes from constant-stress and step-stress accelerated tests

In order to gather the information about the lifetime distribution of a product, a standard life testing method at normal operating conditions is not practical when the product has an extremely long lifespan. Accelerated life testing solves this difficult issue by subjecting the test units at higher stress levels than normal for quicker and more failure data. The lifetime at the design stress is then estimated through extrapolation using an appropriate regression model. Estimation of the regression parameters based on exponentially distributed lifetimes from accelerated life tests has been considered by a number of authors using numerical methods but without systematic or analytical validation. In this article, we propose an alternative approach based on a simple and easy-to-apply graphical method, which also establishes the existence and uniqueness of the maximum likelihood estimates for constant-stress and step-stress accelerated life tests under progressive censorings

Optimal Experimental Planning of Reliability Experiments Based on Coherent Systems

In industrial engineering and manufacturing, assessing the reliability of a product or system is an important topic. Life-testing and reliability experiments are commonly used reliability assessment methods to gain sound knowledge about product or system lifetime distributions. Usually, a sample of items of interest is subjected to stresses and environmental conditions that characterize the normal operating conditions. During the life-test, successive times to failure are recorded and lifetime data are collected. Life-testing is useful in many industrial environments, including the automobile, materials, telecommunications, and electronics industries. There are different kinds of life-testing experiments that can be applied for different purposes. For instance, accelerated life tests (ALTs) and censored life tests are commonly used to acquire information in reliability and life-testing experiments with the presence of time and resource limitations. Statistical inference based on the data obtained from a life test and effectively planning a life-testing experiment subject to some constraints are two important problems statisticians are interested in. The experimental design problem for a life test has long been studied; however, the experimental planning considering putting the experimental units into systems for a life-test has not been studied. In this thesis, we study the optimal experimental planning problem in multiple stress levels life-testing experiments and progressively Type-II censored life-testing experiments when the test units can be put into coherent systems for the experiment. Based on the notion of system signature, a tool in structure reliability to represent the structure of a coherent system, under different experimental settings, models and assumptions, we derive the maximum likelihood estimators of the model parameters and the expected Fisher information matrix. Then, we use the expected Fisher information matrix to obtain the asymptotic variance-covariance matrix of the maximum likelihood estimators when $n$-component coherent systems are used in the life-testing experiment. Based on different optimality criteria, such as $D$-optimality, $A$-optimality and $V$-optimality, we obtain the optimal experimental plans under different settings. Numerical and Monte Carlo simulation studies are used to demonstrate the advantages and disadvantages of using systems in life-testing experiments

Order-statistics-based inferences for censored lifetime data and financial risk analysis

This thesis was submitted for the degree of Doctor of Philosophy and awarded by Brunel University.This thesis focuses on applying order-statistics-based inferences on lifetime analysis and financial risk measurement. The first problem is raised from fitting the Weibull distribution to progressively censored and accelerated life-test data. A new orderstatistics- based inference is proposed for both parameter and con dence interval estimation. The second problem can be summarised as adopting the inference used in the first problem for fitting the generalised Pareto distribution, especially when sample size is small. With some modifications, the proposed inference is compared with classical methods and several relatively new methods emerged from recent literature. The third problem studies a distribution free approach for forecasting financial volatility, which is essentially the standard deviation of financial returns. Classical models of this approach use the interval between two symmetric extreme quantiles of the return distribution as a proxy of volatility. Two new models are proposed, which use intervals of expected shortfalls and expectiles, instead of interval of quantiles. Different models are compared with empirical stock indices data. Finally, attentions are drawn towards the heteroskedasticity quantile regression. The proposed joint modelling approach, which makes use of the parametric link between the quantile regression and the asymmetric Laplace distribution, can provide estimations of the regression quantile and of the log linear heteroskedastic scale simultaneously. Furthermore, the use of the expectation of the check function as a measure of quantile deviation is discussed

Stress-strength reliability estimation for the inverted exponentiated Rayleigh distribution under unified progressive hybrid censoring with application

In this paper, we studied the estimation of a stress-strength reliability model (R = P(X>Y)) based on inverted exponentiated Rayleigh distribution under the unified progressive hybrid censoring scheme (unified PHCS). The maximum likelihood estimates of the unknown parameters were obtained using the stochastic expectation-maximization algorithm (stochastic EMA). The asymptotic confidence intervals were also created. Under squared error and Linex and generalized entropy loss functions, the Gibbs sampler together with Metropolis-Hastings algorithm was provided to compute the Bayes estimates and the credible intervals. Extensive simulations were performed to see the effectiveness of the proposed estimation methods. Also, parallel to the development of reliability studies, it is necessary to study its application in different sciences such as engineering. Therefore, droplet splashing data under two nozzle pressures were proposed to exemplify the theoretical outcomes

Inference and optimal design for the k-level step-stress accelerated life test based on progressive Type-I interval censored power Rayleigh data

In this paper, a new generalization of the one parameter Rayleigh distribution called the Power Rayleigh (PRD) was employed to model the life of the tested units in the step-stress accelerated life test. Under progressive Type-I interval censored data, the cumulative exposure distribution was considered to formulate the life model, assuming the scale parameter of PRD has the inverse power function at each stress level. Point estimates of the model parameters were obtained via the maximum likelihood estimation method, while interval estimates were obtained using the asymptotic normality of the derived estimators and the bootstrap resampling method. An extensive simulation study of $k = 4$ levels of stress in different combinations of the life test under different progressive censoring schemes was conducted to investigate the performance of the obtained point and interval estimates. Simulation results indicated that point estimates of the model parameters are closest to their initial true values and have relatively small mean squared errors. Accordingly, the interval estimates have small lengths and their coverage probabilities are almost convergent to the 95% significance level. Based on the Fisher information matrix, the D-optimality and the A-optimality criteria are implemented to determine the optimal design of the life test by obtaining the optimum inspection times and optimum stress levels that improve the estimation procedures and give more efficient estimates of the model parameters. Finally, the developed inferential procedures were also applied to a real dataset

Regression methods for survival and multistate models.

A common research interest in medical, biological, and engineering research is determining whether certain independent variables are correlated with the survival or failure times. Standard statistical techniques cannot usually be applied for failure-time data due to the lack of complete data or in other word, due to censoring. From a statistical perspective, the study of time to event data is even more challenging when further complexities such as high dimensionality or multivariablity is added to the model. In this dissertation, we consider the predicating patient survival from proteomic profile of patient serum using matrix-assisted laser desorption/ionization time-of-flight (MALDI-TOF) data of non-small cell lung cancer patients. Due to much larger dimension of features in a mass spectrum compared to the study sample size, traditional linear regression modeling of survival times with high number of proteomic features is not feasible. Hence, we consider latent factor and regularized/penalized methods for fitting such models in order to predict patient survival from the mass spectrometry features. Extensive numerical studies involving both simulated as well as real mass spectrometry data are used to compare four popular regression methods, namely, partial least squares (PLS), sparse partial least square (SPLS), least absolute shrinkage and selection operator (LASSO) and elastic net regularization, on processed spectra. Right censoring is handled through a residual based multiple imputation. Overall, more complex methods such as the elastic net and SPLS result in better performances provided the operational parameters are chosen carefully via cross validation. For survival time prediction, we recommend using the elastic net based on a selected set of features. As a type of multivariate survival data, multistate models have a wide range of applications. Most of the existing regression approaches to analyze such data are based on parametric and semi-parametric procedures in which one should rely on specific model structures. In this dissertation, we construct non-parametric regression estimators of a number of temporal functions in a multistate system based on a univariate continuous baseline covariate. These estimators include state occupation probabilities, state entry, exit and waiting (sojourn) times distribution functions of a general progressive (e.g. acyclic) multistate model. The data are subject to right censoring and the censoring mechanism is explainable by observable covariates that could be time dependent. The resulting estimators are valid even if the multistate process is non-Markov. The performance of the estimators is studied using a detailed simulation. We illustrate our estimators using a data set on bone marrow transplant patients. Finally, some extension of the proposed methods to more general case with multivariate covariates are presented along with plans for future developments