A modified chi-square test for Bertholon model with censored data

In this work, we propose the construction of a chi-squared goodness-of-fit test in censored data case, for Bertholon model which can analyse various competing risks of failure or death. This test is based on a modification of the Nikulin-Rao-Robson (NRR) statistic proposed by Bagdonavicius and Nikulin (2011a Bagdonavicius, V., Nikulin, M. (2011a). Chi-squared tests for general composite hypotheses from censored samples. Comptes Rendus Mathématiques: Series I 349(3–4):219–223. [Google Scholar], 2011b Bagdonavicius, V., Nikulin, M. (2011b). Chi-squared goodness-of-fit test for right censored data. International Journal of Applied Mathematics and Statistics 24:30–50. [Google Scholar]) for censored data. We applied this test to numerical examples from simulated samples and real data.     

Estimation of the Parameters of the Birnbaum–Saunders Distribution

Some alternative estimators to the maximum likelihood estimators of the two parameters of the Birnbaum–Saunders distribution are proposed. Most have high efficiencies as measured by root mean square error and are robust to departure from the model as well as to outliers. In addition, the proposed estimators are easy to compute. Both complete and right-censored data are discussed. Simulation studies are provided to compare the performance of the estimators.     

On the Additive Weibull Distribution

We derive explicit expressions for the moments, incomplete moments, quantile function and generating function of the additive Weibull model pioneered by Xie and Lai (1995 Xie, M., Lai, C.D. (1995). Reliability analysis using an additive Weibull model with bathtub-shaped failure rate function. Reliab. Eng. Syst. Safety 52:87–93.[Crossref], [Web of Science ®] , [Google Scholar]), which is a quite flexible distribution for fitting lifetime data with bathtub-shaped failure rate function. In addition, we estimate the model parameters by maximum likelihood and determine the observed information matrix. The flexibility of the additive Weibull distribution is illustrated by means of one application to real data.     

On Maximum Likelihood Estimation for the Two Parameter Burr XII Distribution

This article considers three related aspects of maximum likelihood estimation of parameters in the two-parameter Burr XII distribution. Specifically, we first provide further clarification to some limiting results in Wingo (1993 Wingo , D. R. ( 1993 ). Maximum likelihood estimation of Burr XII distribution parameters under Type II censoring . Microelectron. Reliab. 33 : 1251 – 1257 .[Crossref], [Web of Science ®] , [Google Scholar]). We then focus on details in a proof of the uniqueness of the maximum likelihood estimators. Finally, we consider using the likelihood approach for data which does not satisfy Wingo's criterion, and show that this results in fitting either a Pareto distribution or an intuitively sensible degenerate distribution to the data. The discussion here is completely general, and not restricted to data obtained under Type II censoring.     

Using BBPSO Algorithm to Estimate the Weibull Parameters with Censored Data

This article proposes the maximum likelihood estimates based on bare bones particle swarm optimization (BBPSO) algorithm for estimating the parameters of Weibull distribution with censored data, which is widely used in lifetime data analysis. This approach can produce more accuracy of the parameter estimation for the Weibull distribution. Additionally, the confidence intervals for the estimators are obtained. The simulation results show that the BB PSO algorithm outperforms the Newton–Raphson method in most cases in terms of bias, root mean square of errors, and coverage rate. Two examples are used to demonstrate the performance of the proposed approach. The results show that the maximum likelihood estimates via BBPSO algorithm perform well for estimating the Weibull parameters with censored data.     

Partially Accelerated Life Tests for the Weibull Distribution Under Multiply Censored Data

This article aims to estimate the parameters of the Weibull distribution in step-stress partially accelerated life tests under multiply censored data. The step partially acceleration life test is that all test units are first run simultaneously under normal conditions for a pre-specified time, and the surviving units are then run under accelerated conditions until a predetermined censoring time. The maximum likelihood estimates are used to obtaining the parameters of the Weibull distribution and the acceleration factor under multiply censored data. Additionally, the confidence intervals for the estimators are obtained. Simulation results show that the maximum likelihood estimates perform well in most cases in terms of the mean bias, errors in the root mean square and the coverage rate. An example is used to illustrate the performance of the proposed approach.     

Parameter estimation of the censored δ-shock model on uniform interval

The censored δ-shock model is a special kind of shock model and it has very important research values in the reliability theory. In this paper, we discuss the parameter estimation of the censored δ-shock model when the inter-arrival times between two successive shocks follows uniform distribution on [a, b]. With the maximum likelihood estimation, we obtain the parameter estimator and expectation of estimator and lifetime of the model. By numerical simulation, we get the empirical result on the estimator of δ and the relationship between δ and other parameters.     

Improved likelihood inference for the shape parameter in Weibull regression

We obtain adjustments to the profile likelihood function in Weibull regression models with and without censoring. Specifically, we consider two different modified profile likelihoods: (i) the one proposed by Cox and Reid [Cox, D.R. and Reid, N., 1987, Parameter orthogonality and approximate conditional inference. Journal of the Royal Statistical Society B, 49, 1–39.], and (ii) an approximation to the one proposed by Barndorff–Nielsen [Barndorff–Nielsen, O.E., 1983, On a formula for the distribution of the maximum likelihood estimator. Biometrika, 70, 343–365.], the approximation having been obtained using the results by Fraser and Reid [Fraser, D.A.S. and Reid, N., 1995, Ancillaries and third-order significance. Utilitas Mathematica, 47, 33–53.] and by Fraser et al. [Fraser, D.A.S., Reid, N. and Wu, J., 1999, A simple formula for tail probabilities for frequentist and Bayesian inference. Biometrika, 86, 655–661.]. We focus on point estimation and likelihood ratio tests on the shape parameter in the class of Weibull regression models. We derive some distributional properties of the different maximum likelihood estimators and likelihood ratio tests. The numerical evidence presented in the paper favors the approximation to Barndorff–Nielsen's adjustment.     

Optimal experimental plan for multi-level stress testing with Weibull regression under progressive Type-II extremal censoring

In the design of constant-stress life-testing experiments, the optimal allocation in a multi-level stress test with Type-I or Type-II censoring based on the Weibull regression model has been studied in the literature. Conventional Type-I and Type-II censoring schemes restrict our ability to observe extreme failures in the experiment and these extreme failures are important in the estimation of upper quantiles and understanding of the tail behaviors of the lifetime distribution. For this reason, we propose the use of progressive extremal censoring at each stress level, whereas the conventional Type-II censoring is a special case. The proposed experimental scheme allows some extreme failures to be observed. The maximum likelihood estimators of the model parameters, the Fisher information, and asymptotic variance–covariance matrices of the maximum likelihood estimates are derived. We consider the optimal experimental planning problem by looking at four different optimality criteria. To avoid the computational burden in searching for the optimal allocation, a simple search procedure is suggested. Optimal allocation of units for two- and four-stress-level situations is determined numerically. The asymptotic Fisher information matrix and the asymptotic optimal allocation problem are also studied and the results are compared with optimal allocations with specified sample sizes. Finally, conclusions and some practical recommendations are provided.     

On a Selection Weibull Distribution

In this article, a selection Weibull distribution is investigated. First, some properties and representations of the model with some plots of the density and hazard rate functions are illustrated. Second, some simple relations of this model with some distributions discussed. In addition, maximum likelihood estimators obtained with numerical methods, and compared by three sub-models with a data set that shows the performance of our model. Finally, a simulation study presented for all parameters.     

Accelerated Life Test Data Analysis for Repairable Systems

In this article, the imperfect maintenance model and proportional intensity model are used to analyze failure data of repairable systems in accelerated life testing. In the design and development phase of products, we should collect and analyze failure data quickly with small proto-type products. Thus, we test the products under accelerated conditions and if the products fail, then we repair and use those continuously in the life testing. Acceleration and repair models are needed to analyze the failure data. An age reduction model (Brown et al.'s Brown et al. 1983 Brown , J. F. , Mahoney , J. F. , Sivazlian , B. D. ( 1983 ). Hysteresis repair in discounted replacement problems . IIE Trans. 15 : 156 – 165 . [CSA] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar], model) and relationship between scale parameter and stress level are assumed. The stress acts multiplicatively on the baseline cumulative intensity. The maximum likelihood method is used, the log-likelihood function is derived, and a maximizing procedure is proposed. In simulation studies, we investigate the accuracy and trends of the maximum likelihood estimator.     

Generalized Inverse Gamma Distribution and its Application in Reliability

In this article, we introduce a new reliability model of inverse gamma distribution referred to as the generalized inverse gamma distribution (GIG). A generalization of inverse gamma distribution is defined based on the exact form of generalized gamma function of Kobayashi (1991). This function is useful in many problems of diffraction theory and corrosion problems in new machines. The new distribution has a number of lifetime special sub-models. For this model, some of its statistical properties are studied. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is derived. We also demonstrate the usefulness of this distribution on a real data set.     

On an Extension of the Exponentiated Weibull Distribution

In this article, the exponentiated Weibull distribution is extended by the Marshall-Olkin family. Our new four-parameter family has a hazard rate function with various desired shapes depending on the choice of its parameters and, thus, it is very flexible in data modeling. It also contains two mixed distributions with applications to series and parallel systems in reliability and also contains several previously known lifetime distributions. We shall study some basic distributional properties of the new distribution. Some closed forms are derived for its moment generating function and moments as well as moments of its order statistics. The model parameters are estimated by the maximum likelihood method. The stress–strength parameter and its estimation are also investigated. Finally, an application of the new model is illustrated using two real datasets.     

On Parameter Inference for Step-Stress Accelerated Life Test with Geometric Distribution

In this article, we present the parameter inference in step-stress accelerated life tests under the tampered failure rate model with geometric distribution. We deal with Type-II censoring scheme involved in experimental data, and provide the maximum likelihood estimate and confidence interval of the parameters of interest. With the help of the Monte-Carlo simulation technique, a comparison of precision of the confidence limits is demonstrated for our method, the Bootstrap method, and the large-sample based procedure. The application of two industrial real datasets shows the proposed method efficiency and feasibility.     

The Marshall-Olkin logistic-exponential distribution

We introduce a three-parameter extension of the exponential distribution which contains as sub-models the exponential, logistic-exponential and Marshall-Olkin exponential distributions. The new model is very flexible and its associated density function can be decreasing or unimodal. Further, it can produce all of the four major shapes of the hazard rate, that is, increasing, decreasing, bathtub and upside-down bathtub. Given that closed-form expressions are available for the survival and hazard rate functions, the new distribution is quite tractable. It can be used to analyze various types of observations including censored data. Computable representations of the quantile function, ordinary and incomplete moments, generating function and probability density function of order statistics are obtained. The maximum likelihood method is utilized to estimate the model parameters. A simulation study is carried out to assess the performance of the maximum likelihood estimators. Two actual data sets are used to illustrate the applicability of the proposed model.     

Maximum likelihood estimation for the proportional odds model with mixed interval-censored failure time data

This article discusses regression analysis of mixed interval-censored failure time data. Such data frequently occur across a variety of settings, including clinical trials, epidemiologic investigations, and many other biomedical studies with a follow-up component. For example, mixed failure times are commonly found in the two largest studies of long-term survivorship after childhood cancer, the datasets that motivated this work. However, most existing methods for failure time data consider only right-censored or only interval-censored failure times, not the more general case where times may be mixed. Additionally, among regression models developed for mixed interval-censored failure times, the proportional hazards formulation is generally assumed. It is well-known that the proportional hazards model may be inappropriate in certain situations, and alternatives are needed to analyze mixed failure time data in such cases. To fill this need, we develop a maximum likelihood estimation procedure for the proportional odds regression model with mixed interval-censored data. We show that the resulting estimators are consistent and asymptotically Gaussian. An extensive simulation study is performed to assess the finite-sample properties of the method, and this investigation indicates that the proposed method works well for many practical situations. We then apply our approach to examine the impact of age at cranial radiation therapy on risk of growth hormone deficiency in long-term survivors of childhood cancer.     

Empirical likelihood method for the multivariate accelerated failure time models

In applications, multivariate failure time data appears when each study subject may potentially experience several types of failures or recurrences of a certain phenomenon, or failure times may be clustered. Three types of marginal accelerated failure time models dealing with multiple events data, recurrent events data and clustered events data are considered. We propose a unified empirical likelihood inferential procedure for the three types of models based on rank estimation method. The resulting log-empirical likelihood ratios are shown to possess chi-squared limiting distributions. The properties can be applied to do tests and construct confidence regions without the need to solve the rank estimating equations nor to estimate the limiting variance-covariance matrices. The related computation is easy to implement. The proposed method is illustrated by extensive simulation studies and a real example.     

On Jointly Estimating Parameters and Missing Data by Maximizing the Complete-Data Likelihood

One approach to handling incomplete data occasionally encountered in the literature is to treat the missing data as parameters and to maximize the complete-data likelihood over the missing data and parameters. This article points out that although this approach can be useful in particular problems, it is not a generally reliable approach to the analysis of incomplete data. In particular, it does not share the optimal properties of maximum likelihood estimation, except under the trivial asymptotics in which the proportion of missing data goes to zero as the sample size increases.     

Uniqueness of the Maximum Likelihood Estimate of the Weibull Distribution Tampered Failure Rate Model

Abstract

Bhattacharyya and Soejoeti (Bhattacharyya, G. K., Soejoeti, Z. A. (1989 Bhattacharyya, G. K. and Soejoeti, Z. A. 1989. Tampered failure rate model for step-stress accelerated life test. Commun. Statist.—Theory Meth., 18(5): 1627–1643. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Tampered failure rate model for step-stress accelerated life test. Commun. Statist.—Theory Meth. 18(5):1627–1643.) pro- posed the TFR model for step-stress accelerated life tests. Under the TFR model, this article proves that the maximum likelihood estimate of the shape parameters is unique for the Weibull distribution in a multiple step-stress accelerated life test, and investigates the accuracy of the maximum likelihood estimate using the Monte-Carlo simulation.     

A Finite Mixture Three-Parameter Weibull Model for the Analysis of Wind Speed Data

This article presents a mixture three-parameter Weibull distribution to model wind speed data. The parameters are estimated by using maximum likelihood (ML) method in which the maximization problem is regarded as a nonlinear programming with only inequality constraints and is solved numerically by the interior-point method. By applying this model to four lattice-point wind speed sequences extracted from National Centers for Environmental Prediction (NCEP) reanalysis data, it is observed that the mixture three-parameter Weibull distribution model proposed in this paper provides a better fit than the existing Weibull models for the analysis of wind speed data under study.