Impact of dependent left truncation in semiparametric competing risks methods: A simulation study

In this study, we investigated the robustness of the methods that account for independent left truncation when applied to competing risks settings with dependent left truncation. We specifically focused on the methods for the proportional cause-specific hazards model and the Fine–Gray model. Simulation experiments showed that these methods are not in general robust against dependent left truncation. The magnitude of the bias was analogous to the strength of the association between left truncation and failure times, the effect of the covariate on the competing cause of failure, and the baseline hazard of left truncation time.     

Analyzing Length-biased Data with Semiparametric Transformation and Accelerated Failure Time Models

Right-censored time-to-event data are often observed from a cohort of prevalent cases that are subject to length-biased sampling. Informative right censoring of data from the prevalent cohort within the population often makes it difficult to model risk factors on the unbiased failure times for the general population, because the observed failure times are length biased. In this paper, we consider two classes of flexible semiparametric models: the transformation models and the accelerated failure time models, to assess covariate effects on the population failure times by modeling the length-biased times. We develop unbiased estimating equation approaches to obtain the consistent estimators of the regression coefficients. Large sample properties for the estimators are derived. The methods are confirmed through simulations and illustrated by application to data from a study of a prevalent cohort of dementia patients.     

Failure time studies with intermittent observation and losses to follow-up

In health research interest often lies in modeling a failure time process but in many cohort studies failure status is only determined at scheduled assessment times. While the assessment times may be fixed upon study entry, individuals may become lost to follow-up and miss visits subsequent to the time of loss to follow-up. We consider a three-state model to characterize a joint failure and loss to follow-up process, and use it to investigate the impact of dependent loss to follow-up on standard parametric, nonparametric, and semiparametric analysis. The effect of dependent loss to follow-up is mitigated by fitting the joint model. The performance of standard methods is studied using the asymptotic theory of misspecified models, and the finite sample performance is examined for the standard and joint analyses through simulation studies. An application to data from a youth smoking prevention study is presented for illustration.     

Large sample tests for testing symmetry and independence in some bivariate exponential models

Bivariate Exponential Distribution (BVED) were introduced by Freund (1961), Marshall and Olkin (1967) and Block and Basu (1974) as models for the distributions of (X,Y) the failure times of dependent components (C₁,C₂). We study the structure of these models and observe that Freund model leads to a regular exponential family with a four dimensional orthogonal parameter. Marshall-Olkin model involving three parameters leads to a conditional or piece wise exponential family and Block-Basu model which also depends on three parameters is a sub-model of the Freund model and is a curved exponential family. We obtain a large sample tests for symmetry as well as independence of (X,Y) in each of these models by using the Generalized Likelihood Ratio Tests (GLRT) or tests basesd on MLE of the parameters and root n consistent estimators of their variance-covariance matrices.     

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.     

Comparing waiting times in a multi-stage model: A log-rank approach

A proper log-rank test for comparing two waiting (i.e. sojourn, gap) times under right censored data has been absent in the survival literature. The classical log-rank test provides a biased comparison even under independent right censoring since the censoring induced on the time since state entry depends on the entry time unless the hazards are semi-Markov. We develop test statistics for comparing K waiting time distributions from a multi-stage model in which censoring and waiting times may be dependent upon the transition history in the multi-stage model. To account for such dependent censoring, the proposed test statistics utilize an inverse probability of censoring weighted (IPCW) approach previously employed to define estimators for the cumulative hazard and survival function for waiting times in multi-stage models. We develop the test statistics as analogues to K-sample log-rank statistics for failure time data, and weak convergence to a Gaussian limit is demonstrated. A simulation study demonstrates the appropriateness of the test statistics in designs that violate typical independence assumptions for multi-stage models, under which naive test statistics for failure time data perform poorly, and illustrates the superiority of the test under proportional hazards alternatives to a Mann–Whitney type test. We apply the test statistics to an existing data set of burn patients.     

Optimal preventive policies for imperfect maintenance models with replacement first and last

ABSTRACT

This paper proposes preventive replacement policies for an operating system which may continuously works for N jobs with random working times and is imperfectly maintained upon failure. As a failure occurs, the system suffers one of the two types of failures based on some random mechanism: type-I (repairable or minor) failure is rectified by a minimal repair, or type-II (non repairable or catastrophic) failure is removed by a corrective replacement. A notation of preventive replacement last model is considered in which the system is replaced before any type-II failure at an operating time T or at number N of working times, whichever occurs last. Comparisons between such a preventive replacement last and the conventional replacement first are discussed in detail. For each model, the optimal schedule of preventive replacement that minimizes the mean cost rate is presented theoretically and determined numerically. Because the framework and analysis are general, the proposed models extend several existing results.     

An estimating equation for parametric shared frailty models with marginal additive hazards

Summary.  Multivariate failure time data arise when data consist of clusters in which the failure times may be dependent. A popular approach to such data is the marginal proportional hazards model with estimation under the working independence assumption. In some contexts, however, it may be more reasonable to use the marginal additive hazards model. We derive asymptotic properties of the Lin and Ying estimators for the marginal additive hazards model for multivariate failure time data. Furthermore we suggest estimating equations for the regression parameters and association parameters in parametric shared frailty models with marginal additive hazards by using the Lin and Ying estimators. We give the large sample properties of the estimators arising from these estimating equations and investigate their small sample properties by Monte Carlo simulation. A real example is provided for illustration.     

Frailty models for survival data

A frailty model is a random effects model for time variables, where the random effect (the frailty) has a multiplicative effect on the hazard. It can be used for univariate (independent) failure times, i.e. to describe the influence of unobserved covariates in a proportional hazards model. More interesting, however, is to consider multivariate (dependent) failure times generated as conditionally independent times given the frailty. This approach can be used both for survival times for individuals, like twins or family members, and for repeated events for the same individual. The standard assumption is to use a gamma distribution for the frailty, but this is a restriction that implies that the dependence is most important for late events. More generally, the distribution can be stable, inverse Gaussian, or follow a power variance function exponential family. Theoretically, large differences are seen between the choices. In practice, using the largest model makes it possible to allow for more general dependence structures, without making the formulas too complicated.This paper is a revised version of a review, which together with ten papers by the author made up a thesis for a Doctor of Science degree at the University of Copenhagen.     

Failure Inference From a Marker Process Based on a Bivariate Wiener Model

Many models have been proposed that relate failure times and stochastic time-varying covariates. In some of these models, failure occurs when a particular observable marker crosses a threshold level. We are interested in the more difficult, and often more realistic, situation where failure is not related deterministically to an observable marker. In this case, joint models for marker evolution and failure tend to lead to complicated calculations for characteristics such as the marginal distribution of failure time or the joint distribution of failure time and marker value at failure. This paper presents a model based on a bivariate Wiener process in which one component represents the marker and the second, which is latent (unobservable), determines the failure time. In particular, failure occurs when the latent component crosses a threshold level. The model yields reasonably simple expressions for the characteristics mentioned above and is easy to fit to commonly occurring data that involve the marker value at the censoring time for surviving cases and the marker value and failure time for failing cases. Parametric and predictive inference are discussed, as well as model checking. An extension of the model permits the construction of a composite marker from several candidate markers that may be available. The methodology is demonstrated by a simulated example and a case application.     

Finite time ruin probability and structural density properties in the presence of dependence in insurance risk model

In the present paper, we consider the classical compound Poisson risk model with dependence between claim sizes and claim inter-arrival time. We attempt to analyze the approximation of finite time ruin probability. The finite time ruin probabilities are plotted for fixed threshold value associated to the claim inter-arrival time and also for fixed dependence parameter in Nelsen (2006) copula separately. Additionally, a general form for joint density of the interclaim times and claim sizes is considered. With respect to the classical Gerber-Shiu's (1998) function, first some structural density properties of dependent collective risk model is obtained. Then the ladder height probability density function of claim sizes is computed and the dependency structure investigated for Erlang interclaim time. As the application, some dependent models of the interclaim times and claim sizes are studied.     

Modeling systems with two dependent components under bivariate shock models

Series and parallel systems consisting of two dependent components are studied under bivariate shock models. The random variables N₁ and N₂ that represent respectively the number of shocks until failure of component 1 and component 2 are assumed to be dependent and phase-type. The times between successive shocks are assumed to follow a continuous phase-type distribution, and survival functions and mean time to failure values of series and parallel systems are obtained in matrix forms. An upper bound for the joint survival function of the components is also provided under the particular case when the times between shocks follow exponential distribution.     

Inference from Accelerated Degradation and Failure Data Based on Gaussian Process Models

An important problem in reliability and survival analysis is that of modeling degradation together with any observed failures in a life test. Here, based on a continuous cumulative damage approach with a Gaussian process describing degradation, a general accelerated test model is presented in which failure times and degradation measures can be combined for inference about system lifetime. Some specific models when the drift of the Gaussian process depends on the acceleration variable are discussed in detail. Illustrative examples using simulated data as well as degradation data observed in carbon-film resistors are presented.     

A joint modelling approach for clustered recurrent events and death events

In dental implant research studies, events such as implant complications including pain or infection may be observed recurrently before failure events, i.e. the death of implants. It is natural to assume that recurrent events and failure events are correlated to each other, since they happen on the same implant (subject) and complication times have strong effects on the implant survival time. On the other hand, each patient may have more than one implant. Therefore these recurrent events or failure events are clustered since implant complication times or failure times within the same patient (cluster) are likely to be correlated. The overall implant survival times and recurrent complication times are both interesting to us. In this paper, a joint modelling approach is proposed for modelling complication events and dental implant survival times simultaneously. The proposed method uses a frailty process to model the correlation within cluster and the correlation within subjects. We use Bayesian methods to obtain estimates of the parameters. Performance of the joint models are shown via simulation studies and data analysis.     

Buckley–James-type estimation of quantile regression with recurrent gap time data

In longitudinal studies, an individual may potentially undergo a series of repeated recurrence events. The gap times, which are referred to as the times between successive recurrent events, are typically the outcome variables of interest. Various regression models have been developed in order to evaluate covariate effects on gap times based on recurrence event data. The proportional hazards model, additive hazards model, and the accelerated failure time model are all notable examples. Quantile regression is a useful alternative to the aforementioned models for survival analysis since it can provide great flexibility to assess covariate effects on the entire distribution of the gap time. In order to analyze recurrence gap time data, we must overcome the problem of the last gap time subjected to induced dependent censoring, when numbers of recurrent events exceed one time. In this paper, we adopt the Buckley–James-type estimation method in order to construct a weighted estimation equation for regression coefficients under the quantile model, and develop an iterative procedure to obtain the estimates. We use extensive simulation studies to evaluate the finite-sample performance of the proposed estimator. Finally, analysis of bladder cancer data is presented as an illustration of our proposed methodology.     

Progressively Type-II censored competing risks data for exponential distributions based on sequential order statistics

We consider the progressively Type-II censored competing risks model based on sequential order statistics. It is assumed that the latent failure times are independent and the failure of each unit influences the lifetime distributions of the latent failure times of surviving units. We provide explicit expressions for the likelihood function of the available data under the conditional proportional hazard rate (CPHR) and the power trend conditional proportional hazard rate (PTCPHR) models. Under CPHR and PTCPHR models and assumption that the baseline distributions of the latent failure times are exponential, classical and Bayesian estimates of the unknown parameters are provided. Monte Carlo simulations are then performed for illustrative purposes. Finally, two datasets are analyzed.     

Optimal scheduling replacement policies for a system with multiple random works

From the economical viewpoint in reliability theory, this paper addresses a scheduling replacement problem for a single operating system which works at random times for multiple jobs. The system is subject to stochastic failure which results the imperfect maintenance activity based on some random failure mechanism: minimal repair due to type-I (repairable) failure, or corrective replacement due to type-II (non-repairable) failure. Three scheduling models for the system with multiple jobs are considered: a single work, N tandem works, and N parallel works. To control the deterioration process, the preventive replacement is planned to undergo at a scheduling time T or the job's completion time of for each model. The objective is to determine the optimal scheduling parameters (T^* or N^*) that minimizes the mean cost rate function in a finite time horizon for each model. A numerical example is provided to illustrate the proposed analytical model. Because the framework and analysis are general, the proposed models extend several existing results.     

Estimating fecundability from data on waiting times to first conception

"This article tests assumptions invoked in the demographic literature to estimate the population distribution of fecundability from data on waiting times to first conception. In continuous time, the key assumption is that waiting times are realizations from a mixture of exponentials distribution. In discrete time, the key assumption is that waiting times are realizations from a mixture of geometrics distribution. The [U.S.] Hutterite data analyzed by Sheps (1965) are consistent with this assumption. Various models, however, have one representation in mixture of exponentials form. A fundamental identification problem plagues the conventional estimation procedure. Our analysis calls into question the conventional practice of checking model specification by using goodness-of-fit tests. The practical importance of the identification problem in duration models is demonstrated."     

Central Limit Theorem for ISE of Kernel Density Estimators in Censored Dependent Model

In some long-term studies, a series of dependent and possibly censored failure times may be observed. Suppose that the failure times have a common continuous distribution function F. A popular stochastic measure of the distance between the density function f of the failure times and its kernel estimate f _n is the integrated square error(ISE). In this article, we derive a central limit theorem for the integrated square error of the kernel density estimators under a censored dependent model.     

Extended Hazard Regression Model for Reliability and Survival Analysis

We propose an extended hazard regression model which allows the spread parameter to be dependent on covariates. This allows a broad class of models which includes the most common hazard models, such as the proportional hazards model, the accelerated failure time model and a proportional hazards/accelerated failure time hybrid model with constant spread parameter. Simulations based on sub-classes of this model suggest that maximum likelihood performs well even when only small or moderate-size data sets are available and the censoring pattern is heavy. The methodology provides a broad framework for analysis of reliability and survival data. Two numerical examples illustrate the results.