Accelerated Degradation Models for Failure Based on Geometric Brownian Motion and Gamma Processes

Based on a generalized cumulative damage approach with a stochastic process describing degradation, new accelerated life test models are presented in which both observed failures and degradation measures can be considered for parametric inference of system lifetime. Incorporating an accelerated test variable, we provide several new accelerated degradation models for failure based on the geometric Brownian motion or gamma process. It is shown that in most cases, our models for failure can be approximated closely by accelerated test versions of Birnbaum–Saunders and inverse Gaussian distributions. Estimation of model parameters and a model selection procedure are discussed, and two illustrative examples using real data for carbon-film resistors and fatigue crack size are presented.     

Variables acceptance reliability sampling plan for items subject to inverse Gaussian degradation process

Until now, in the literature, a variety of acceptance reliability sampling plans have been developed based on different life test plans. In most of the reliability sampling plans, the decision procedures to accept or reject the corresponding lot are developed based on the lifetimes of the items observed on tests, or the number of failures observed during a pre-specified testing time. However, frequently, the items are subject to degradation phenomena and, in these cases, the observed degradation level of the item can be used as a decision statistic. In this paper, we develop a variables acceptance sampling plan based on the information on the degradation process of the items, assuming that the degradation process follows the inverse Gaussian process. It is shown that the developed sampling plan improves the reliability performance of the items conditional on the acceptance in the test and that the lifetimes of items after the reliability sampling test are stochastically larger than those before the test. A study comparing the proposed degradation-based sampling plan with the conventional sampling plan which is based on a life test is also performed.KEYWORDS: Variables sampling plan, degradation test, inverse Gaussian process, mixture distribution, stochastic ordering     

Optimal step-stress accelerated degradation test plans for inverse Gaussian process based on proportional degradation rate model

The purpose of this paper is to address the optimal design of the step-stress accelerated degradation test (SSADT) issue when the degradation process of a product follows the inverse Gaussian (IG) process. For this design problem, an important task is to construct a link model to connect the degradation magnitudes at different stress levels. In this paper, a proportional degradation rate model is proposed to link the degradation paths of the SSADT with stress levels, in which the average degradation rate is proportional to an exponential function of the stress level. Two optimization problems about the asymptotic variances of the lifetime characteristics' estimators are investigated. The optimal settings including sample size, measurement frequency and the number of measurements for each stress level are determined by minimizing the two objective functions within a given budget constraint. As an example, the sliding metal wear data are used to illustrate the proposed model.     

Modelling Accelerated Degradation Data Using Wiener Diffusion With A Time Scale Transformation

Engineering degradation tests allow industry to assess the potential life span of long-life products that do not fail readily under accelerated conditions in life tests. A general statistical model is presented here for performance degradation of an item of equipment. The degradation process in the model is taken to be a Wiener diffusion process with a time scale transformation. The model incorporates Arrhenius extrapolation for high stress testing. The lifetime of an item is defined as the time until performance deteriorates to a specified failure threshold. The model can be used to predict the lifetime of an item or the extent of degradation of an item at a specified future time. Inference methods for the model parameters, based on accelerated degradation test data, are presented. The model and inference methods are illustrated with a case application involving self-regulating heating cables. The paper also discusses a number of practical issues encountered in applications. This revised version was published online in July 2006 with corrections to the Cover Date.     

Normal Inverse Gaussian Distributions and Stochastic Volatility Modelling

The normal inverse Gaussian distribution is defined as a variance-mean mixture of a normal distribution with the inverse Gaussian as the mixing distribution. The distribution determines an homogeneous Lévy process, and this process is representable through subordination of Brownian motion by the inverse Gaussian process. The canonical, Lévy type, decomposition of the process is determined. As a preparation for developments in the latter part of the paper the connection of the normal inverse Gaussian distribution to the classes of generalized hyperbolic and inverse Gaussian distributions is briefly reviewed. Then a discussion is begun of the potential of the normal inverse Gaussian distribution and Lévy process for modelling and analysing statistical data, with particular reference to extensive sets of observations from turbulence and from finance. These areas of application imply a need for extending the inverse Gaussian Lévy process so as to accommodate certain, frequently observed, temporal dependence structures. Some extensions, of the stochastic volatility type, are constructed via an observation-driven approach to state space modelling. At the end of the paper generalizations to multivariate settings are indicated.     

Inverse Gaussian process models for bivariate degradation analysis: A Bayesian perspective

This article conducts a Bayesian analysis for bivariate degradation models based on the inverse Gaussian (IG) process. Assume that a product has two quality characteristics (QCs) and each of the QCs is governed by an IG process. The dependence of the QCs is described by a copula function. A bivariate simple IG process model and three bivariate IG process models with random effects are investigated by using Bayesian method. In addition, a simulation example is given to illustrate the effectiveness of the proposed methods. Finally, an example about heavy machine tools is presented to validate the proposed models.     

Longitudinal conditional models with intermittent missingness: SAS code and applications

This work provides a set of macros performed with SAS (Statistical Analysis System) for Windows, which can be used to fit conditional models under intermittent missingness in longitudinal data. A formalized transition model, including random effects for individuals and measurement error, is presented. Model fitting is based on the missing completely at random or missing at random assumptions, and the separability condition. The problem translates to maximization of the marginal observed data density only, which for Gaussian data is again Gaussian, meaning that the likelihood can be expressed in terms of the mean and covariance matrix of the observed data vector. A simulation study is presented and misspecification issues are considered. A practical application is also given, where conditional models are fitted to the data from a clinical trial that assessed the effect of a Cuban medicine on a disease of the respiratory system.     

Geostatistical Modelling Using Non‐Gaussian Matérn Fields

This work provides a class of non‐Gaussian spatial Matérn fields which are useful for analysing geostatistical data. The models are constructed as solutions to stochastic partial differential equations driven by generalized hyperbolic noise and are incorporated in a standard geostatistical setting with irregularly spaced observations, measurement errors and covariates. A maximum likelihood estimation technique based on the Monte Carlo expectation‐maximization algorithm is presented, and a Monte Carlo method for spatial prediction is derived. Finally, an application to precipitation data is presented, and the performance of the non‐Gaussian models is compared with standard Gaussian and transformed Gaussian models through cross‐validation.     

A Non‐Parametric Estimator of the Spectral Density of a Continuous‐Time Gaussian Process Observed at Random Times

Abstract. In numerous applications data are observed at random times and an estimated graph of the spectral density may be relevant for characterizing and explaining phenomena. By using a wavelet analysis, one derives a non‐parametric estimator of the spectral density of a Gaussian process with stationary increments (or a stationary Gaussian process) from the observation of one path at random discrete times. For every positive frequency, this estimator is proved to satisfy a central limit theorem with a convergence rate depending on the roughness of the process and the moment of random durations between successive observations. In the case of stationary Gaussian processes, one can compare this estimator with estimators based on the empirical periodogram. Both estimators reach the same optimal rate of convergence, but the estimator based on wavelet analysis converges for a different class of random times. Simulation examples and an application to biological data are also provided.     

Failure Inference and Optimization for Step Stress Model Based on Bivariate Wiener Model

In this article, we consider the situation under a life test, in which the failure time of the test units are not related deterministically to an observable stochastic time varying covariate. In such a case, the joint distribution of failure time and a marker value would be useful for modeling the step stress life test. The problem of accelerating such an experiment is considered as the main aim of this article. We present a step stress accelerated model based on a bivariate Wiener process with one component as the latent (unobservable) degradation process, which determines the failure times and the other as a marker process, the degradation values of which are recorded at times of failure. Parametric inference based on the proposed model is discussed and the optimization procedure for obtaining the optimal time for changing the stress level is presented. The optimization criterion is to minimize the approximate variance of the maximum likelihood estimator of a percentile of the products’ lifetime distribution.     

Functional Autoregression for Sparsely Sampled Data

We develop a hierarchical Gaussian process model for forecasting and inference of functional time series data. Unlike existing methods, our approach is especially suited for sparsely or irregularly sampled curves and for curves sampled with nonnegligible measurement error. The latent process is dynamically modeled as a functional autoregression (FAR) with Gaussian process innovations. We propose a fully nonparametric dynamic functional factor model for the dynamic innovation process, with broader applicability and improved computational efficiency over standard Gaussian process models. We prove finite-sample forecasting and interpolation optimality properties of the proposed model, which remain valid with the Gaussian assumption relaxed. An efficient Gibbs sampling algorithm is developed for estimation, inference, and forecasting, with extensions for FAR(p) models with model averaging over the lag p. Extensive simulations demonstrate substantial improvements in forecasting performance and recovery of the autoregressive surface over competing methods, especially under sparse designs. We apply the proposed methods to forecast nominal and real yield curves using daily U.S. data. Real yields are observed more sparsely than nominal yields, yet the proposed methods are highly competitive in both settings. Supplementary materials, including R code and the yield curve data, are available online.     

Estimation for binary models generated by Gaussian autoregressive processes

Regression-type and partial likelihood models are presented for binary data obtained by clipping a Gaussian autoregressive process. Five methods for estimating parameters of the model are proposed and compared via a simulation study. A real data analysis is also presented.     

A Bayesian prediction using the skew Gaussian distribution

A model based on the skew Gaussian distribution is presented to handle skewed spatial data. It extends the results of popular Gaussian process models. Markov chain Monte Carlo techniques are used to generate samples from the posterior distributions of the parameters. Finally, this model is applied in the spatial prediction of weekly rainfall. Cross-validation shows that the predictive performance of our model compares favorably with several kriging variants.     

Analysis of Bivariate Survival Data using Shared Inverse Gaussian Frailty Model

In this article, we consider shared frailty model with inverse Gaussian distribution as frailty distribution and log-logistic distribution (LLD) as baseline distribution for bivariate survival times. We fit this model to three real-life bivariate survival data sets. The problem of analyzing and estimating parameters of shared inverse Gaussian frailty is the interest of this article and then compare the results with shared gamma frailty model under the same baseline for considered three data sets. Data are analyzed using Bayesian approach to the analysis of clustered survival data in which there is a dependence of failure time observations within the same group. The variance component estimation provides the estimated dispersion of the random effects. We carried out a test for frailty (or heterogeneity) using Bayes factor. Model comparison is made using information criteria and Bayes factor. We observed that the shared inverse Gaussian frailty model with LLD as baseline is the better fit for all three bivariate data sets.     

Estimating degradation by a wiener diffusion process subject to measurement error

Most materials and components degrade physically before they fail. Engineering degradation tests are designed to measure these degradation processes. Measurements in the tests reflect the inherent randomness of degradation itself as well as measurement errors created by imperfect instruments, procedures and environments. This paper describes a statistical model for measured degradation data that takes both sources of variation into account. The degradation process in the model is taken to be a Wiener diffusion process. The measurement errors are assumed to be independent normal random outcomes that are independent of the degradation process. The paper describes inference procedures for the model and discusses some practical issues that must be considered in dealing with the statistical problem. A case study is presented.     

Optimal Design for Step-Stress Accelerated Degradation Test with Multiple Performance Characteristics Based on Gamma Processes

Step-stress accelerated degradation test (SSADT) plays an important role in assessing the lifetime distribution of highly reliable products under normal operating conditions when there are not enough test units available for testing purposes. Recently, the optimal SSADT plans are presented based on an underlying assumption that there is only one performance characteristic. However, many highly reliable products usually have complex structure, with their reliability being evaluated by two or more performance characteristics. At the same time, the degradation of these performance characteristics would be always positive and strictly increasing. In such a case, the gamma process is usually considered as a degradation process due to its independent and nonnegative increments properties. Therefore, it is of great interest to design an efficient SSADT plan for the products with multiple performance characteristics based on gamma processes. In this work, we first introduce reliability model of the degradation products with two performance characteristics based on gamma processes, and then present the corresponding SSADT model. Next, under the constraint of total experimental cost, the optimal settings such as sample size, measurement times, and measurement frequency are obtained by minimizing the asymptotic variance of the estimated 100 qth percentile of the product’s lifetime distribution. Finally, a numerical example is given to illustrate the proposed procedure.     

Inverse gaussian accelerated test models based on cumulative damage

The failure of a system under environmental stress often can be described by an accelerated test model which incorporates the environmental variable L. Here, the failure of such a system at environmental level L is modeled as the first passage of accumulated damage to a critical threshold value. Assuming a discrete additive damage model leads to a Birnbaum–Saunders-type distribution for the failure time which can be closely approximated by an inverse Gaussian-type model. However, if a continuous damage model based on a Gaussian process is assumed, a more general family of inverse Gaussian accelerated test models is obtained. Three sets of failure data are discussed to illustrate the usefulness of this general family.     

Testing for linearity in generalized linear models using kernel smoothing

A test statistic proposed by Li (1999) for testing the adequacy of heteroscedastic nonlinear regression models using nonparametric kernel smoothers is applied to testing for linearity in generalized linear models. Simulation results for models with centered gamma and inverse Gaussian errors are presented to illustrate the performance of the resulting test compared with log-likelihood ratio tests for specific parametric alternatives. The test is applied to a data set of coronary heart disease status (Hosmer and Lemeshow, (1990).     

Modeling binary familial data using Gaussian copula

Modeling binary familial data has been a challenging task due to the dependence among family members and the constraints imposed on the joint probability distribution of the binary responses. This paper investigates some useful familial dependence structures and proposes analyzing binary familial data using Gaussian copula model. Advantages of this approach are discussed as well as some computational details. An numerical example is also presented with an aim to show the capability of Gaussian copula model in more sophisticated data analysis.     

On a new goodness‐of‐fit process for families of copulas

A goodness‐of‐fit procedure is proposed for parametric families of copulas. The new test statistics are functionals of an empirical process based on the theoretical and sample versions of Spearman's dependence function. Conditions under which this empirical process converges weakly are seen to hold for many families including the Gaussian, Frank, and generalized Farlie–Gumbel–Morgenstern systems of distributions, as well as the models with singular components described by Durante [Durante ( 2007 ) Comptes Rendus Mathématique. Académie des Sciences. Paris, 344, 195–198]. Thanks to a parametric bootstrap method that allows to compute valid P‐values, it is shown empirically that tests based on Cramér–von Mises distances keep their size under the null hypothesis. Simulations attesting the power of the newly proposed tests, comparisons with competing procedures and complete analyses of real hydrological and financial data sets are presented. The Canadian Journal of Statistics 37: 80‐101; 2009 © 2009 Statistical Society of Canada