Estimation of parameters for trend-renewal processes

Methods of estimating unknown parameters of a trend function for trend-renewal processes are investigated in the case when the renewal distribution function is unknown. If the renewal distribution is unknown, then the likelihood function of the trend-renewal process is unknown and consequently the maximum likelihood method cannot be used. In such a situation we propose three other methods of estimating the trend parameters. The methods proposed can also be used to predict future occurrence times. The performance of the estimators based on these methods is illustrated numerically for some trend-renewal processes for which the statistical inference is analytically intractable.     

Estimation of Weibull distribution parameters for irregular interval group failure data with unknown failure times

SUMMARY This paper presents three methods for estimating Weibull distribution parameters for the case of irregular interval group failure data with unknown failure times. The methods are based on the concepts of the piecewise linear distribution function (PLDF), an average interval failure rate (AIFR) and sequential updating of the distribution function (SUDF), and use an analytical approach similar to that of Ackoff and Sasieni for regular interval group data. Results from a large number of simulated case problems generated with specified values of Weibull distribution parameters have been presented, which clearly indicate that the SUDF method produces near-perfect parameter estimates for all types of failure pattern. The performances of the PLDF and AIFR methods have been evaluated by goodness-of-fit testing and statistical confidence limits on the shape parameter. It has been found that, while the PLDF method produces acceptable parameter estimates, the AIFR method may fail for low and high shape parameter values that represent the cases of random and wear-out types of failure. A real-life application of the proposed methods is also presented, by analyzing failures of hydrogen make-up compressor valves in a petroleum refinery.     

Numerical Comparison Between Different Empirical Prediction Intervals Under the Fay-Herriot Model

Recently, an empirical best linear unbiased predictor is widely used as a practical approach to small area inference. It is also of interest to construct empirical prediction intervals. However, we do not know which method should be used from among the several existing prediction intervals. In this article, we first obtain an empirical prediction interval by using the residual maximum likelihood method for estimating unknown model variance parameters. Then we compare the later with other intervals with the residual maximum likelihood method. Additionally, some different parametric bootstrap methods for constructing empirical prediction intervals are also compared in a simulation study.     

Interval estimation for the generalized inverted exponential distribution under progressive first failure censoring

In this paper, we consider the inferential procedures for the generalized inverted exponential distribution under progressive first failure censoring. The exact confidence interval for the scale parameter is derived. The generalized confidence intervals (GCIs) for the shape parameter and some commonly used reliability metrics such as the quantile and the reliability function are explored. Then the proposed procedure is extended to the prediction interval for the future measurement. The GCIs for the reliability of the stress-strength model are discussed under both equal scale and unequal scale scenarios. Extensive simulations are used to demonstrate the performance of the proposed GCIs and prediction interval. Finally, an example is used to illustrate the proposed methods.     

Inversion of the Renewal Density with Dead-time

Stationary renewal point processes are defined by the probability distribution of the distances between successive points (lifetimes) that are independent and identically distributed random variables. For some applications it is also interesting to define the properties of a renewal process by using the renewal density. There are well-known expressions of this density in terms of the probability density of the lifetimes. It is more difficult to solve the inverse problem consisting in the determination of the density of the lifetimes in terms of the renewal density. Theoretical expressions between their Laplace transforms are available but the inversion of these transforms is often very difficult to obtain in closed form. We show that this is possible for renewal processes presenting a dead-time property characterized by the fact that the renewal density is zero in an interval including the origin. We present the principle of a recursive method allowing the solution of this problem and we apply this method to the case of some processes with input dead-time. Computer simulations on Poisson and Erlang (2) processes show quite good agreement between theoretical calculations and experimental measurements on simulated data.     

Robust Bayesian estimation of a bathtub-shaped distribution under progressive Type-II censoring

The bathtub-shaped failure rate function has been used for modeling the life spans of a number of electronic and mechanical products, as well as for modeling the life spans of humans, especially when some of the data are censored. This article addresses robust methods for the estimation of unknown parameters in a two-parameter distribution with a bathtub-shaped failure rate function based on progressive Type-II censored samples. Here, a class of flexible priors is considered by using the hierarchical structure of a conjugate prior distribution, and corresponding posterior distributions are obtained in a closed-form. Then, based on the square error loss function, Bayes estimators of unknown parameters are derived, which depend on hyperparameters as parameters of the conjugate prior. In order to eliminate the hyperparameters, hierarchical Bayesian estimation methods are proposed, and these proposed estimators are compared to one another based on the mean squared error, through Monte Carlo simulations for various progressively Type-II censoring schemes. Finally, a real dataset is presented for the purpose of illustration.     

Estimation and prediction for Chen distribution with bathtub shape under progressive censoring

We consider estimation of the unknown parameters of Chen distribution [Chen Z. A new two-parameter lifetime distribution with bathtub shape or increasing failure rate function. Statist Probab Lett. 2000;49:155–161] with bathtub shape using progressive-censored samples. We obtain maximum likelihood estimates by making use of an expectation–maximization algorithm. Different Bayes estimates are derived under squared error and balanced squared error loss functions. It is observed that the associated posterior distribution appears in an intractable form. So we have used an approximation method to compute these estimates. A Metropolis–Hasting algorithm is also proposed and some more approximate Bayes estimates are obtained. Asymptotic confidence interval is constructed using observed Fisher information matrix. Bootstrap intervals are proposed as well. Sample generated from MH algorithm are further used in the construction of HPD intervals. Finally, we have obtained prediction intervals and estimates for future observations in one- and two-sample situations. A numerical study is conducted to compare the performance of proposed methods using simulations. Finally, we analyse real data sets for illustration purposes.     

Bayesian prediction of waiting times in stochastic models

The authors show how saddlepoint techniques lead to highly accurate approximations for Bayesian predictive densities and cumulative distribution functions in stochastic model settings where the prior is tractable, but not necessarily the likelihood or the predictand distribution. They consider more specifically models involving predictions associated with waiting times for semi‐Markov processes whose distributions are indexed by an unknown parameter θ. Bayesian prediction for such processes when they are not stationary is also addressed and the inverse‐Gaussian based saddlepoint approximation of Wood, Booth & Butler (1993) is shown to accurately deal with the nonstationarity whereas the normal‐based Lugannani & Rice (1980) approximation cannot, Their methods are illustrated by predicting various waiting times associated with M/M/q and M/G/1 queues. They also discuss modifications to the matrix renewal theory needed for computing the moment generating functions that are used in the saddlepoint methods.     

Analysis of progressively type-I interval censored competing risks data for a class of an exponential distribution

ABSTRACT

In this paper, we consider some problems of point estimation and point prediction when the competing risks data from a class of exponential distribution are progressive type-I interval censored. The maximum likelihood estimation and mid-point approximation method are proposed for the estimations of parameters. Also several point predictors of censored units such as the maximum likelihood predictor, the best unbiased predictor and the conditional median predictor are obtained. The methods discussed here are applied when the lifetime distributions of the latent failure times are independent and Weibull-distributed. Finally a simulation study is given by using Monte-Carlo simulations to compare the performances of the different methods and one data analysis has been presented for illustrative purposes.     

Two-dimensional Renewal Function Approximation

Two-dimensional renewal functions, which are naturally extensions of one-dimensional renewal functions, have wide applicability in areas where two random variables are needed to characterize the underlying process. These functions satisfy the renewal equation, which is not amenable for analytical solutions. This paper proposes a simple approximation for the computation of the two- dimensional renewal function based only on the first two moments and the correlation coefficient of the variables. The approximation yields exact values of renewal function for bivariate exponential distribution function. Illustrations are presented to compare our approximation with that of Iskandar (1991) who provided a computational procedure which requires the use of the bivariate distribution function of the two variables. A two-dimensional warranty model is used to illustrate the approximation.     

On spatiotemporal prediction for on-line monitoring data

Spatiotemporal prediction is of interest in many areas of applied statistics, especially in environmental monitoring with on-line data information. At first, this article reviews the approaches for spatiotemporal modeling in the context of stochastic processes and then introduces the new class of spatiotemporal dynamic linear models. Further, the methods for linear spatial data analysis, universal kriging and trend surface prediction, are related to the method of spatial linear Bayesian analysis. The Kalman filter is the preferred method for temporal linear Bayesian inferences. By combining the Kalman filter recursions with the trend surface predictor and universal kriging predictor, the prior and posterior spatiotemporal predictors for the observational process are derived, which form the main result of this article. The problem of spatiotemporal linear prediction in the case of unknown first and second order moments is treated as well.     

Bayes inference for the modulated power law process

The Modulated Power Law process has been recently proposed as a suitable model for describing the failure pattern of repairable systems when both renewal-type behaviour and time trend are present. Unfortunately, the maximum likelihood method provides neither accurate confidence intervals on the model parameters for small or moderate sample sizes nor predictive intervals on future observations.

This paper proposes a Bayes approach, based on both non-informative and vague prior, as an alternative to the classical method. Point and interval estimation of the parameters, as well as point and interval prediction of future failure times, are given. Monte Carlo simulation studies show that the Bayes estimation and prediction possess good statistical properties in a frequentist context and, thus, are a valid alternative to the maximum likelihood approach.

Numerical examples illustrate the estimation and prediction procedures.     

Estimation and prediction for a Burr type-III distribution with progressive censoring

We consider estimation of unknown parameters and reliability characteristics of a Burr type-III distribution under progressive censoring. Predictive estimates for censored observations and the associated prediction intervals are also obtained. We derive maximum-likelihood estimators of unknown quantities using the EM algorithm and then also obtain the observed Fisher information matrix. We provide various Bayes estimators for unknown parameters under the squared error loss function. Highest posterior density and asymptotic intervals are also constructed. We evaluate performance of proposed methods using simulations. Finally, an illustrative example is presented in support of the methods discussed.     

Shortest confidence and prediction intervals for the log-normal

We provide the shortest prediction interval for X, and the shortest confidence interval for the median of X, when X has the log-normal distribution for both the case σ², the variance of log X, known and unknown. Tables are given to assist the practitioner in constructing these intervals. A real-life example is provided to illustrate the results.     

PREDICTION FOR EXPONENTIAL LIFETIMES BASED ON STEP-STRESS TESTING

ABSTRACT

This paper presents methods for constructing prediction limits for a step-stress model in accelerated life testing. An exponential life distribution with a mean that is a log-linear function of stress, and a cumulative exposure model are assumed. Two prediction problems are discussed. One concerns the prediction of the life at a design stress, and the other concerns the prediction of a future life during the step-stress testing. Both predictions require the knowledge of some model parameters. When estimates for the model parameters are available, a calibration method based on simulations is proposed for correcting the prediction intervals (regions) obtained by treating the parameter estimates as the true parameter values. Finally, a numerical example is given to illustrate the prediction procedure.     

Incorporating diagnostic accuracy into the estimation of discrete survival function

Empirical distribution function (EDF) is a commonly used estimator of population cumulative distribution function. Survival function is estimated as the complement of EDF. However, clinical diagnosis of an event is often subjected to misclassification, by which the outcome is given with some uncertainty. In the presence of such errors, the true distribution of the time to first event is unknown. We develop a method to estimate the true survival distribution by incorporating negative predictive values and positive predictive values of the prediction process into a product-limit style construction. This will allow us to quantify the bias of the EDF estimates due to the presence of misclassified events in the observed data. We present an unbiased estimator of the true survival rates and its variance. Asymptotic properties of the proposed estimators are provided and these properties are examined through simulations. We evaluate our methods using data from the VIRAHEP-C study.     

Comparison of renewal function estimators for censored data

The computation of the renewal function when the distribution function is completely known has received much attention in the literature. However, in many cases the form of the distribution function is unknown and has to be estimated nonparametrically. A nonparametric estimator for the renewal function for complete data was suggested by Frees (1986). In many cases, however, censoring of the lifetime might occur. We shall present parametric and nonparametric estimators of the renewal function based on censored data. In a simulation study we compare the nonparametric estimators with parametric estimators for the Weibull and lognormal distribution. The study suggests that the nonparametric estimator is a viable alternative to the parametric estimators when the lifetime distribution is unknown. Also, the nonparametric estimator is computationally simpler than the parametric estimator.     

Prediction intervals for future Weibull residual data

In reliability theory, risk analysis, renewal processes and actuarial studies, the residual lifetimes data play an important essential role in studying the conditional tail of the lifetime data. In this paper, based on some observed ordered residual Weibull data, we introduce different prediction methods for obtaining prediction intervals (PIs) of future residual lifetimes including likelihood, Wald, moments, parametric bootstrap, and highest conditional methods. Monte Carlo simulations are performed to compare the performances of the so obtained PIs and one data analysis is performed for illustration purposes.     

Statistical inference based on generalized Lindley record values

This paper addresses the problems of frequentist and Bayesian estimation for the unknown parameters of generalized Lindley distribution based on lower record values. We first derive the exact explicit expressions for the single and product moments of lower record values, and then use these results to compute the means, variances and covariance between two lower record values. We next obtain the maximum likelihood estimators and associated asymptotic confidence intervals. Furthermore, we obtain Bayes estimators under the assumption of gamma priors on both the shape and the scale parameters of the generalized Lindley distribution, and associated the highest posterior density interval estimates. The Bayesian estimation is studied with respect to both symmetric (squared error) and asymmetric (linear-exponential (LINEX)) loss functions. Finally, we compute Bayesian predictive estimates and predictive interval estimates for the future record values. To illustrate the findings, one real data set is analyzed, and Monte Carlo simulations are performed to compare the performances of the proposed methods of estimation and prediction.     

Estimation of the generalized exponential renewal function

When the shape parameter is a non-integer of the generalized exponential (GE) distribution, the analytical renewal function (RF) usually is not tractable. To overcome this, the approximation method has been used in this paper. In the proposed model, the n-fold convolution of the GE cumulative distribution function (CDF) is approximated by n-fold convolutions of gamma and normal CDFs. We obtain the GE RF by a series approximation model. The method is very simple in the computation. Numerical examples have shown that the approximate models are accurate and robust. When the parameters are unknown, we present the asymptotic confidence interval of the RF. The validity of the asymptotic confidence interval is checked via numerical experiments.