New control charts for monitoring the Weibull percentiles under complete data and Type‐II censoring

In this paper, we propose 3 new control charts for monitoring the lower Weibull percentiles under complete data and Type‐II censoring. In transforming the Weibull distribution to the smallest extreme value distribution, Pascaul et al (2017) presented an exponentially weighted moving average (EWMA) control chart, hereafter referred to as EWMA‐SEV‐Q, based on a pivotal quantity conditioned on ancillary statistics. We extended their concept to construct a cumulative sum (CUSUM) control chart denoted by CUSUM‐SEV‐Q. We provide more insights of the statistical properties of the monitoring statistic. Additionally, in transforming a Weibull distribution to a standard normal distribution, we propose EWMA and CUSUM control charts, denoted as EWMA‐YP and CUSUM‐YP, respectively, based on a pivotal quantity for monitoring the Weibull percentiles with complete data. With complete data, the EWMA‐YP and CUSUM‐YP control charts perform better than the EWMA‐SEV‐Q and CUSUM‐SEV‐Q control charts in terms of average run length. In Type‐II censoring, the EWMA‐SEV‐Q chart is slightly better than the CUSUM‐SEV‐Q chart in terms of average run length. Two numerical examples are used to illustrate the applications of the proposed control charts.     

Some Percentile Estimators for Weibull Parameters

A percentile estimator for the shape parameter of the Weibull distribution, based on the 17th and 97th sample percentiles, is proposed which is asymptotically about 66% efficient when compared with the MLE (maximum likelihood estimator). A two-observation percentile estimator, based on the 40th and 82nd sample percentiles, for the scale parameter when the shape parameter is unknown is asymptotically about 82y0 efficient when compared with the MLE. The 24th and 93rd sample percentiles yield asymptotically about 41ye jointly efficient percentile estimators for both the scale and shape parameters in a class of two-observation percentile estimators when compared with their MLEs. Some other simple percentile estimators for these parameters are also briefly discussed. Finally, asymptotic properties of these estimators are investigated and their application in statistical inference problems is mentioned.     

Robust Estimation for Weibull Distribution in Partially Accelerated Life Tests with Early Failures

Maximum likelihood estimation (MLE) is a frequently used method for estimating distribution parameters in constant stress partially accelerated life tests (CS‐PALTs). However, using the MLE to estimate the parameters for a Weibull distribution may be problematic in CS‐PALTs. First, the equation for the shape parameter estimator derived from the log‐likelihood function is difficult to solve for the occurrence of nonlinear equations. Second, the sample size is typically not large in life tests. The MLE, a typical large‐sample inference method, may be unsuitable. Test items unsuitable for stress conditions may become early failures, which have extremely short lifetimes. The early failures may cause parameter estimate bias. For addressing early failures in the Weibull distribution in CS‐PALTs, we propose an M‐estimation method based on a Weibull Probability Plot (WPP) framework, which leads a closed‐form expression for the shape parameter estimator. We conducted a simulation study to compare the M‐estimation method with the MLE method. The results show that, with early‐failure samples, the M‐estimation method performs better than the MLE does. Copyright © 2015 John Wiley & Sons, Ltd.     

More on the mis‐specification of the shape parameter with Weibull‐to‐exponential transformation

When lifetimes follow Weibull distribution with known shape parameter, a simple power transformation could be used to transform the data to the case of exponential distribution, which is much easier to analyze. Usually, the shape parameter cannot be known exactly and it is important to investigate the effect of mis‐specification of this parameter. In a recent article, it was suggested that the Weibull‐to‐exponential transformation approach should not be used as the confidence interval for the scale parameter has very poor statistical property. However, it would be of interest to study the use of Weibull‐to‐exponential transformation when the mean time to failure or reliability is to be estimated, which is a more common question. In this paper, the effect of mis‐specification of Weibull shape parameters on these quantities is investigated. For reliability‐related quantities such as mean time to failure, percentile lifetime and mission reliability, the Weibull‐to‐exponential transformation approach is generally acceptable. For the cases when the data are highly censored or when small tail probability is concerned, further studies are needed, but these are known to be difficult statistical problems for which there are no standard solutions. Copyright © 2000 John Wiley & Sons, Ltd.     

Reliability analysis of wind turbine subassemblies based on the 3-P Weibull model via an ergodic artificial bee colony algorithm

The Weibull distribution is the most widely used model for the reliability evaluation of wind turbine subassemblies. Considering the important role of the location parameter in the three-parameter (3-P) Weibull model and its rare application in wind turbines, this study conducted a reliability analysis of wind turbine subassemblies based on field data that obeyed the 3-P Weibull distribution model via maximum likelihood estimation (MLE). An improved ergodic artificial bee colony algorithm (ErgoABC) was proposed by introducing the chaos search theory, global best solution, and Lévy flights strategy into the classical artificial bee colony (ABC) algorithm to determine the maximum likelihood estimates of the Weibull distribution parameters. This was validated against simulation calculations and proved to be efficient for high-dimensional function optimization and parameter estimation of the 3-P Weibull distribution. Finally, reliability analyses of the wind turbine subassemblies based on different types of field failure data were conducted using ErgoABC. The results show that the 3-P Weibull model can reasonably evaluate the lifetime distribution of critical wind turbine subassemblies, such as generator slip rings and main shafts, on which the location parameter has a significant effect.     

Comparison of Methods to Estimate the Time‐to‐failure Distribution in Degradation Tests

Degradation tests are alternative approaches to lifetime tests and accelerated lifetime tests in reliability studies. Based on a degradation process of a product quality characteristic over time, degradation tests provide enough information to estimate the time‐to‐failure distribution. Some estimation methods, such as analytical, the numerical or the approximated, can be used to obtain the time‐to‐failure distribution. They are chosen according to the complexity of the degradation model used in the data analysis. An example of the application and analysis of degradation tests is presented in this paper to characterize the durability of a product and compare the various estimation methods of the time‐to‐failure distribution. The example refers to a degradation process related to an automobile's tyre, and was carried out to estimate its average distance covered and some percentiles of interest. Copyright © 2004 John Wiley & Sons, Ltd.     

A Bootstrap Control Chart for Weibull Percentiles

The problem of detecting a shift of a percentile of a Weibull population in a process monitoring situation is considered. The parametric bootstrap method is used to establish lower and upper control limits for monitoring percentiles when process measurements have a Weibull distribution. Small percentiles are of importance when observing tensile strength and it is desirable to detect their downward shift. The performance of the proposed bootstrap percentile charts is considered based on computer simulations, and some comparisons are made with an existing Weibull percentile chart. The new bootstrap chart indicates a shift in the process percentile substantially quicker than the previously existing chart, while maintaining comparable average run lengths when the process is in control. An illustrative example concerning the tensile strength of carbon fibers is presented. Copyright © 2005 John Wiley & Sons, Ltd.     

Analysis of field failure data when modeled Weibull slope increases with time

Weibull time‐to‐fail distributions cannot be correctly estimated from field data when manufacturing populations from different vintages have different failure modes. To investigate the pitfalls of ongoing Weibull parameter estimation, two cases, based upon real events, were analyzed. First, a time‐to‐fail distribution was generated assuming the same Weibull shape parameter representing an increasing failure rate for each monthly batch or vintage of production. The shape parameter was estimated from simulated field data at regular periods as the population accumulated service time. Estimates of the shape parameter were not constant, but gradually decreased (as had occurred in a real system) with added service time. In the second case, field reliability performance was modeled to match the actual historical data for one product from a disk drive manufacturer. The actual data was proprietary and was not directly available for analysis. A production schedule was modeled with a mix of two failure characteristics. The population reaching the field in the first 12 months had a low, constant failure rate. For the second and third years of production, higher volumes were introduced that had the higher, increasing failure rates of the first case. Assessment of the mixed population at each month of calendar time resulted in an increasing Weibull shape parameter estimate at each assessment. When the two populations were separated and estimated properly, a better fit with more accurate estimates of Weibull shape parameters resulted. Copyright © 2010 John Wiley & Sons, Ltd.     

Reliability analysis using exponentiated Weibull distribution and inverse power law

Today in reliability analysis, the most used distribution to describe the behavior of devices is the Weibull distribution. Nonetheless, the Weibull distribution does not provide an excellent fit to lifetime datasets that exhibit bathtub shaped or upside‐down bathtub shaped (unimodal) failure rates, which are often encountered in the performance of products such as electronic devices (ED). In this paper, a reliability model based on the exponentiated Weibull distribution and the inverse power law model is proposed, this new model provides a better approach to model the performance and fit of the lifetimes of electronic devices. A case study based on the lifetime of a surface‐mounted electrolytic capacitor is presented in this paper. Besides, it was found that the estimation of the proposed model differs from the Weibull classical model and that affects the mean time to failure (MTTF) of the capacitor under analysis.     

The omega probability distribution and its applications in reliability theory

A new three‐parameter probability distribution called the omega probability distribution is introduced, and its connection with the Weibull distribution is discussed. We show that the asymptotic omega distribution is just the Weibull distribution and point out that the mathematical properties of the novel distribution allow us to model bathtub‐shaped hazard functions in two ways. On the one hand, we demonstrate that the curve of the omega hazard function with special parameter settings is bathtub shaped and so it can be utilized to describe a complete bathtub‐shaped hazard curve. On the other hand, the omega probability distribution can be applied in the same way as the Weibull probability distribution to model each phase of a bathtub‐shaped hazard function. Here, we also propose two approaches for practical statistical estimation of distribution parameters. From a practical perspective, there are two notable properties of the novel distribution, namely, its simplicity and flexibility. Also, both the cumulative distribution function and the hazard function are composed of power functions, which on the basis of the results from analyses of real failure data, can be applied quite effectively in modeling bathtub‐shaped hazard curves.     

Reliability analysis for electronic devices using beta‐Weibull distribution

Today, in reliability analysis, the most used distribution to describe the behavior of electronic products under voltage profiles is the Weibull distribution. Nevertheless, the Weibull distribution does not provide a good fit to lifetime datasets that exhibit bathtub‐shaped or upside‐down bathtub–shaped (unimodal) failure rates, which are often encountered in the reliability analysis of electronic devices. In this paper, a reliability model based on the beta‐Weibull distribution and the inverse power law is proposed. This new model provides a better approach to model the performance and fit of the lifetimes of electronic devices. To estimate the parameters of the proposed model, a Bayesian analysis is used. A case study based on the lifetime of a surface mounted electrolytic capacitor is presented, the results showed that the estimation of the proposed model differs from the inverse power law–Weibull and that it affects directly the mean time to failure, the failure rate, the behavior, and the performance of the capacitor under analysis.     

Parameter estimation for bivariate Weibull distribution using generalized moment method for reliability evaluation

Bivariate Weibull distribution can address the life of a system exhibiting 2‐dimensional characteristics in risk and reliability engineering. The applicability of bivariate Weibull distribution has been hindered by its difficulty with parameter estimation, as the number of parameters in bivariate Weibull distribution is more than those in univariate Weibull distribution. Considering a particular structure of a bivariate Weibull distribution model, this paper proposes a generalized moment method (GMM) for parameter estimation. This GMM method is simple, and it has proved to be efficient. The GMM can guarantee the existence and the uniqueness of the solution. A confidence interval for each estimator is derived from the moments of the bivariate distribution. The paper presents a simulation case and 2 real cases to demonstrate the proposed methods.     

On estimating percentiles of the Weibull distribution by the linear regression method

Probability estimators developed previously by the authors have been used to obtain unbiased estimates of the Weibull parameters by the linear regression method. Using these unbiased estimators, percentiles of the Weibull distribution have been estimated. Since these percentiles are determined from the estimated parameters, they also have distributions and subsequently are determined for five sample sizes. Analysis has shown that the distributions of these estimated percentiles are neither normal, lognormal, three-parameter Weibull nor three-parameter log-Weibull. A new methodology to estimate the percentile with a specified level of confidence has been introduced. The step-by-step use of the methodology is demonstrated by examples in this paper.     

Point and Interval Estimation Procedures for the Two-Parameter Weibull and Extreme-Value Distributions

Point estimators of parameters of the first asymptotic distribution of smallest (extreme) values, or, the extreme-value distribution, are surveyed and compared. Those investigated are maximum-likelihood and moment estimators, inefficient estimators based on only a few ordered observations, and various linear estimation methods. A combination of Monte Carlo approximations and exact small-sample and asymptotic results has been used to compare the expected loss (with loss equal to squared error) of these various point estimators. Since the logarithms of variates having the two-parameter Weibull distribution are variates from the extreme-value distribution, the investigation is applicable to the estimation of Weibull parameters. Interval estimation procedures are also discussed.     

A sequential constant‐stress accelerated life testing scheme and its Bayesian inference

In the analysis of accelerated life testing (ALT) data, some stress‐life model is typically used to relate results obtained at stressed conditions to those at use condition. For example, the Arrhenius model has been widely used for accelerated testing involving high temperature. Motivated by the fact that some prior knowledge of particular model parameters is usually available, this paper proposes a sequential constant‐stress ALT scheme and its Bayesian inference. Under this scheme, test at the highest stress is firstly conducted to quickly generate failures. Then, using the proposed Bayesian inference method, information obtained at the highest stress is used to construct prior distributions for data analysis at lower stress levels. In this paper, two frameworks of the Bayesian inference method are presented, namely, the all‐at‐one prior distribution construction and the full sequential prior distribution construction. Assuming Weibull failure times, we (1) derive the closed‐form expression for estimating the smallest extreme value location parameter at each stress level, (2) compare the performance of the proposed Bayesian inference with that of MLE by simulations, and (3) assess the risk of including empirical engineering knowledge into ALT data analysis under the proposed framework. Step‐by‐step illustrations of both frameworks are presented using a real‐life ALT data set. Copyright © 2008 John Wiley & Sons, Ltd.     

Monitoring the Weibull Shape Parameter with Type II Censored Data

In this article, we introduce a method for monitoring the Weibull shape parameter β with type II (failure) censored data. The control limits depend on the sample size, the number of censored observations, the target average run length, and the stable value of β. The method assumes that the scale parameter α is constant during each sampling period, which is true under rational subgrouping. The proposed method utilizes the relationship between Weibull and smallest extreme value distribution. We propose an unbiased estimator of σ = 1/β as the monitoring statistic. We derive the control limits for one‐sided and two‐sided charts for several stable process average run lengths. We discuss two schemes, namely, the control‐limits‐only scheme and the control‐limits‐with‐warning‐lines scheme. The stable process average run length performance of the proposed charts is studied and compared with those of other charts for monitoring β under similar assumptions. Copyright © 2014 John Wiley & Sons, Ltd.     

EWMA chart based on Bayes‐conditional pivotal quantity for Weibull percentiles under complete data and type‐II censoring

Bayes‐conditional control chart has been used for monitoring the Weibull percentiles with complete data and type‐II censoring. Firstly, the Weibull data are transformed to the smallest extreme value (SEV) distribution. Secondly, the posterior median of quantiles is used as a monitoring statistic. Finally, a pivotal quantity based on the monitoring statistic with its conditional distribution function is derived for obtaining the control limits. This control chart is denoted as Shewhart‐SEV‐ . In this study, we extend this work based on an exponential weighted moving average model named exponential weighted moving average‐SEV‐ for monitoring the Weibull percentiles. We provide the statistical properties of the monitoring statistic. The average run length and the standard deviation of run lengths, computed by the integral equation approach, are used as performance measures. The results indicate that the proposed chart performs better than the Shewhart‐SEV‐ . The breaking strength of carbon fibers is used to illustrate the application of the proposed control chart.     

Robust Regression using Probability Plots for Estimating the Weibull Shape Parameter

The Weibull shape parameter is important in reliability estimation as it characterizes the ageing property of the system. Hence, this parameter has to be estimated accurately. This paper presents a study of the efficiency of using robust regression methods over the ordinary least‐squares regression method based on a Weibull probability plot. The emphasis is on the estimation of the shape parameter of the two‐parameter Weibull distribution. Both the case of small data sets with outliers and the case of data sets with multiple‐censoring are considered. Maximum‐likelihood estimation is also compared with linear regression methods. Simulation results show that robust regression is an effective method in reducing bias and it performs well in most cases. Copyright © 2006 John Wiley & Sons, Ltd.     

A two-parameter parsimonious bathtub model for the analysis of system failure times with illustrations to reliability data

In this study, a two-parameter, upper-bounded probability distribution called the tau distribution is introduced and its applications in reliability engineering are presented. Each of the parameters of the tau distribution has a clear semantic meaning. Namely, one of them determines the upper bound of the distribution, while the value of the other parameter influences the shape of the cumulative distribution function. A remarkable property of this new probability distribution is that its probability density function, survival function, hazard rate function (HRF), and quantile function can all be expressed in terms of its cumulative distribution function. The HRF of the proposed probability distribution can exhibit an increasing trend and various bathtub shapes with or without a low and long-flat phase (useful time phase), which makes this new distribution suitable for modeling a wide range of real-world problems. The constraint maximum likelihood estimation, percentile estimation, approximate Bayesian computation, and approximate quantile estimation computation are proposed to calculate the unknown parameters of the model. The suitability of the estimation methods is verified with the aid of simulation and real-world data results. The modeling capability of the tau distribution was compared with that of some well-known two- and three-parameter probability distributions using two data sets known from the literature of reliability engineering: time between failures data of a machining center, and time to failure of data acquisition system cards. Based on empirical results, the new distribution may be viewed as a viable competitor to the Weibull, Gamma, Chen, and modified Weibull distributions.     

Monitoring type‐2 censored reliability data in multistage processes

In this paper, we propose control charts to monitor the Weibull scale parameter of type‐2 censored reliability data in multistage processes. A cumulative sum control chart and 2 exponentially weighted moving average control charts based on conditional expected values are devised to detect decreases in the mean level of reliability‐related quality characteristic. The proposed control schemes are based on standard smallest extreme value distributions derived from Weibull processes to effectively account for the cascade property, which is the main characteristic of multistage processes. Subsequently, simulation study is conducted to evaluate the performance of the control charts using average run length criterion. Extra quadratic loss, performance comparison index, and relative average run length are also used to compare the detect ability of our proposed monitoring procedures. Moreover, sensitivity analysis is done to study the impact of failure number in the sample size and to investigate the robustness of the proposed monitoring procedures against the shift in the previous stage. Finally, a real case study in a glass bottle–making company is investigated to illustrate the performance of the competing control charts. The results reveal the superiority of the cumulative sum control chart.