Approximate Identifiability, Moments and Censored Data

It is well known that the joint distribution of a pair of random variables ( X,Y ) is not identifiable on the basis of the joint distribution of the function (min ( X,Y ), 1[ X < Y ]). This paper introduces the concept of approximate identifiability and studies its relevance to the function (min ( X,Y ), Y ). It shows that the distribution of ( X,Y ) is approximately identifiable on the basis of the distribution of (min ( X,Y ), Y ). The identification is explicitly executed by a method of moments. The method is applied to the analysis of censored distributions arising in the theory of clinical trials and is compared to the standard method of Kaplan and Meier.     

A Test for Independence of Competing Risks with Discrete Failure Times

The identifiability problem in Competing Risks is well known. In particular, it implies that independent action or otherwise of the risks cannot be inferred from data alone. However, Crowder (1996) showed that, in the case of purely discrete failure times, an inference can be made. An algebraic criterion was derived which bears essentially on the independence in question. The condition was presented in a theoretical setting but it was pointed out that the quantities involved can be estimated from data and that, therefore, there is the potential to develop practical tests for the hypothesis of independence. It is the purpose of this paper to construct such a test. This revised version was published online in July 2006 with corrections to the Cover Date.     

Non-parametric Quantile Regression with Censored Data

Abstract.  Censored regression models have received a great deal of attention in both the theoretical and applied statistics literature. Here, we consider a model in which the response variable is censored but not the covariates. We propose a new estimator of the conditional quantiles based on the local linear method, and give an algorithm for its numerical implementation. We study its asymptotic properties and evaluate its performance on simulated data sets.     

Nonparametric Quantile Estimation with Correlated Failure Time Data

In biomedical studies, correlated failure time data arise often. Although point and confidence interval estimation for quantiles with independent censored failure time data have been extensively studied, estimation for quantiles with correlated failure time data has not been developed. In this article, we propose a nonparametric estimation method for quantiles with correlated failure time data. We derive the asymptotic properties of the quantile estimator and propose confidence interval estimators based on the bootstrap and kernel smoothing methods. Simulation studies are carried out to investigate the finite sample properties of the proposed estimators. Finally, we illustrate the proposed method with a data set from a study of patients with otitis media.     

Nonparametric Test for Doubly Interval-Censored Failure Time Data

This paper considers comparison of discrete failure time distributions when the survival time of interest measures elapsed time between two related events and observations on the occurrences of both events could be interval-censored. This kind of data is often referred to as doubly interval-censored failure time data. If the occurrence of the first event defining the survival time can be exactly observed, the data are usually referred to as interval-censored data. For the comparison problem based on interval-censored failure time data, Sun (1996) proposed a nonparametric test procedure. In this paper we generalize the procedure given in Sun (1996) to doubly interval-censored data case and the generalized test is evaluated by simulations.     

Marginal Analysis For Clustered Failure Time Data

Clustered failure time data are commonly encountered in biomedical research where the study subjects from the same cluster (e.g., family) share the common genetic and/or environmental factors such that the failure times within the same cluster are correlated. Two approaches that are commonly used to account for the intra-cluster association are frailty models and marginal models. In this paper, we study the marginal proportional hazards model, where the structure of dependence between individuals within a cluster is unspecified. An estimation procedure is developed based on a pseudo-likelihood approach, and a risk set sampling method is proposed for the formulation of the pseudo-likelihood. The asymptotic properties of the proposed estimators are studied, and the related issues regarding the statistical efficiencies are discussed. The performances of the proposed estimator are demonstrated by the simulation studies. A data example from a child vitamin A supplementation trial in Nepal (Nepal Nutrition Intervention Project-Sarlahi, or NNIPS) is used to illustrate this methodology.     

On assessing independence of Competing Risks when failure times are discrete

The traditional approach to modelling for Competing Risks, via a multivariate distribution of latent failure times, is very natural for many applications but suffers from a well-documented problem of identifiability. However, the demonstrations of this problem in the literature apply to essentially continuous latent failure times where any atoms of probability in their distributions are not too intrusive. It is shown in this paper that for discrete failure times the classic results on the identifiability problem concerning the existence of equivalent independent risks are incomplete.     

Statistical Analysis of Linear Degradation and Failure Time Data with Multiple Failure Modes

The paper considers linear degradation and failure time models with multiple failure modes. Dependence of traumatic failure intensities on the degradation level are included into the models. Estimators of traumatic event cumulative intensities, and of various reliability characteristics are proposed. Prediction of residual reliability characteristics given a degradation value at a given moment is discussed. Non-parametric, semiparametric and parametric estimation methods are given. Theorems on simultaneous asymptotic distribution of random functions characterising degradation and intensities of traumatic events are proposed. Asymptotic properties of unconditional and residual reliability characteristics estimators are given. Real tire wear and failure time data are analysed.     

A Discrete Time Parametric Model for the Analysis of Failure Time Data

This paper introduces a parametric discrete failure time model which allows a variety of smooth hazard function shapes, including shapes which are not readily available with continuous failure time models. The model is easy to fit, and statistical inference is simple. Further, it is readily extended to allow for differences between subjects while retaining the ease of fit and simplicity of statistical inference. The performance of the discrete time analysis is demonstrated by application to several data sets.     

Weighted Estimator for the Linear Transformation Models with Multivariate Failure Time Data

In this article, a simple and efficient weighted method is proposed to improve the estimation efficiency for the linear transformation models with multivariate failure time data. Asymptotic properties of the estimators with a closed-form variance-covariance matrix are established. In addition, a goodness-of-fit test is developed to evaluate the adequacy of the model. The performance of proposed method and the comparison on the efficiency between the proposed method and the working independence method (Lu, 2005) are conducted in finite-sample situation by simulation studies. Finally a real data set from the Busselton Population Health Surveys is illustrated to validate the proposed methodology. The related proofs of the theorems are given in the Appendix.     

Characterizations of Competing Risks in Terms of Independent-Risks Proxy Models

In competing risks a failure time T and a cause C , one of p possible, are observed. A traditional representation is via a vector ( T ₁, ..., T_p ) of latent failure times such that T = min( T ₁, ..., T_p ); C is defined by T = T_C in the basic situation of failure from a single cause. There are several results in the literature to the effect that a joint distribution for ( T ₁, ..., T_p ), in which the T_j are independent, can always be constructed to yield any given bivariate distribution for ( C , T ). For this reason the prevailing wisdom is that independence cannot be assessed from competing risks data, not even with arbitrarily large sample sizes (e.g. Prentice et al. , 1978). A result was given by Crowder (1996) which shows that, under certain circumstances, independence can be assessed. The various results will be drawn together and a complete characterization can now be given in terms of independent-risks proxy models.     

Composite Estimating Equation Method for the Accelerated Failure Time Model with Length‐biased Sampling Data

Length‐biased sampling data are often encountered in the studies of economics, industrial reliability, epidemiology, genetics and cancer screening. The complication of this type of data is due to the fact that the observed lifetimes suffer from left truncation and right censoring, where the left truncation variable has a uniform distribution. In the Cox proportional hazards model, Huang & Qin (Journal of the American Statistical Association, 107, 2012, p. 107) proposed a composite partial likelihood method which not only has the simplicity of the popular partial likelihood estimator, but also can be easily performed by the standard statistical software. The accelerated failure time model has become a useful alternative to the Cox proportional hazards model. In this paper, by using the composite partial likelihood technique, we study this model with length‐biased sampling data. The proposed method has a very simple form and is robust when the assumption that the censoring time is independent of the covariate is violated. To ease the difficulty of calculations when solving the non‐smooth estimating equation, we use a kernel smoothed estimation method (Heller; Journal of the American Statistical Association, 102, 2007, p. 552). Large sample results and a re‐sampling method for the variance estimation are discussed. Some simulation studies are conducted to compare the performance of the proposed method with other existing methods. A real data set is used for illustration.     

Estimation of Parameters of Mixed Failure Time Distribution Based on an Extended Modified Sampling Scheme

Abstract

The problem of estimation of parameters of a mixture of degenerate (at zero) and exponential distribution is considered by Dixit and Prasad [Dixit, V. U. (Nee: Jayade, V. D.), Prasad, M. S. (1990 Dixit, V. U. and Prasad, M. S. 1990. Estimation of parameters of mixed failure time distribution. Commun.in Statist.-Theory Meth, 19(12): 4667–4678. (Nee: Jayade, V. D.) [Google Scholar]). Estimation of parameters of mixed failure time distribution. Commun.in Statist.-Theory Meth., 19(12):4667–4678]. The sampling scheme proposed in it is extended to k positive observations in Dixit [Dixit, V. U. (1993 Dixit, V. U. 1993. “Statistical Inference for AR (1) Process with Mixed Errors”. Kolhapur, , India: Shivaji University. Unpublished Ph.D. thesis [Google Scholar]). Statistical Inference for AR (1) Process with Mixed Errors. Unpublished Ph.D. thesis, Shivaji University Kolhapur, India] and moment estimator, MLE and UMVUE based on it are obtained and their finite sample and asymptotic properties are studied. These results are presented in this paper. It is interesting to mention that the sampling scheme proposed by Shinde and Shanubhogue [Shinde, R. L., Shanubhogue, A. (2000 Shinde, R. L. and Shanubhogue, A. 2000. Estimation of parameters and the mean life of a mixed failure time distribution. Commun. Statist.-Theory Meth, 29(11): 2621–2642. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). Estimation of parameters and the mean life of a mixed failure time distribution. Commun. Statist.-Theory Meth. 29(11):2621–2642] is a particular case of the sampling scheme proposed in Dixit [Dixit, V. U. (1993 Dixit, V. U. 1993. “Statistical Inference for AR (1) Process with Mixed Errors”. Kolhapur, , India: Shivaji University. Unpublished Ph.D. thesis [Google Scholar]). Statistical Inference for AR (1) Process with Mixed Errors. Unpublished Ph.D. thesis, Shivaji University Kolhapur, India] for n = k.     

Estimated Generalized Estimating Equation for Correlated Failure Time Data with Auxiliary Covariates

In this article, we study a marginal hazard model with common baseline hazard for correlated failure time data. We assume that the true covariate is measured precisely in a subset of the whole study cohort, whereas an auxiliary information for the true covariate is available for the whole cohort. We first estimate the relative risk function empirically. Then we obtain the estimator for the regression parameter by replacing the relative risk function with its estimator in a generalized estimating equation (GEE) proposed by Cai (1992 Cai , J. ( 1992 ). Generalized estimation equations for censored multivariate failure time data. Ph.D. thesis, University of Washington, Seattle, Washington . [Google Scholar]). A key feature of this method is that it is nonparametric with respect to the association between the missing covariate and the observed auxiliary covariate. The proposed estimator is shown to be consistent and asymptotically normal. Furthermore, we present a corrected Breslow-type estimator for the cumulative hazard function. Simulation studies are conducted to evaluate the proposed method.     

Model Selection For A Family Of Discrete Failure Time Distributions

In recent years, a large number of new discrete distributions have appeared in the literature. However, flexible discrete models which, at the same time, allow for easy statistical inference, are still an exception. This paper makes a detailed analysis of a family of discrete failure time distributions which meets both requirements. It examines the maximum likelihood estimation of the unknown parameters and presents a goodness-of-fit test for this model. The test is used for the selection of an appropriate model for datasets of frequencies of the duration of atmospheric circulation patterns.     

Estimation of the marginal survival time in the presence of dependent competing risks using inverse probability of censoring weighted (IPCW) methods

In medical studies, there is interest in inferring the marginal distribution of a survival time subject to competing risks. The Kyushu Lipid Intervention Study (KLIS) was a clinical study for hypercholesterolemia, where pravastatin treatment was compared with conventional treatment. The primary endpoint was time to events of coronary heart disease (CHD). In this study, however, some subjects died from causes other than CHD or were censored due to loss to follow-up. Because the treatments were targeted to reduce CHD events, the investigators were interested in the effect of the treatment on CHD events in the absence of causes of death or events other than CHD. In this paper, we present a method for estimating treatment group-specific marginal survival curves of time-to-event data in the presence of dependent competing risks. The proposed method is a straightforward extension of the Inverse Probability of Censoring Weighted (IPCW) method to settings with more than one reason for censoring. The results of our analysis showed that the IPCW marginal incidence for CHD was almost the same as the lower bound for which subjects with competing events were assumed to be censored at the end of all follow-up. This result provided reassurance that the results in KLIS were robust to competing risks.     

Rank Regression Analysis of Multivariate Failure Time Data Based on Marginal Linear Models

Abstract.  Multivariate failure time data arises when each study subject can potentially ex-perience several types of failures or recurrences of a certain phenomenon, or when failure times are sampled in clusters. We formulate the marginal distributions of such multivariate data with semiparametric accelerated failure time models (i.e. linear regression models for log-transformed failure times with arbitrary error distributions) while leaving the dependence structures for related failure times completely unspecified. We develop rank-based monotone estimating functions for the regression parameters of these marginal models based on right-censored observations. The estimating equations can be easily solved via linear programming. The resultant estimators are consistent and asymptotically normal. The limiting covariance matrices can be readily estimated by a novel resampling approach, which does not involve non-parametric density estimation or evaluation of numerical derivatives. The proposed estimators represent consistent roots to the potentially non-monotone estimating equations based on weighted log-rank statistics. Simulation studies show that the new inference procedures perform well in small samples. Illustrations with real medical data are provided.     

The Current Duration Approach to Estimating Time to Pregnancy

Abstract. Time‐to‐pregnancy (TTP) is the duration from the time a couple starts trying to become pregnant until they succeed. It is considered one of the most direct methods to measure natural fecundity in humans. Statistical tools for designing and analysing time to pregnancy studies belong to survival analysis, but several features require special attention. Prospective designs are difficult to carry out and retrospective (pregnancy‐based) designs, being widely used in this area, do not allow efficiently including couples remaining childless. A third possible design starts from a cross‐sectional sample of couples currently trying to become pregnant, using current duration (backward recurrence time) as basis for the estimation of TTP. Regression analysis is then most conveniently carried out in the accelerated failure time model. This paper surveys some practical and technical‐statistical issues in implementing this approach in a large telephone‐based survey, the Epidemiological Observatory of Fecundity in France (Obseff).     

Semiparametric Likelihood Estimation in the Clayton–Oakes Failure Time Model

Multivariate failure time data arise when the sample consists of clusters and each cluster contains several possibly dependent failure times. The Clayton–Oakes model (Clayton, 1978; Oakes, 1982) for multivariate failure times characterizes the intracluster dependence parametrically but allows arbitrary specification of the marginal distributions. In this paper, we discuss estimation in the Clayton–Oakes model when the marginal distributions are modeled to follow the Cox (1972) proportional hazards regression model. Parameter estimation is based on an approximate generalized maximum likelihood estimator. We illustrate the model's application with example datasets.     

Statistical Analysis of General Degradation Path Model and Failure Time Data with Multiple Failure Modes

ABSTRACT

This article considers degradation and failure time models with multiple failure modes which used to study the problem of longevity and aging in survival analysis and reliability. Degradation process is modeled using general nonparametric, nonlinear path models. Semi-parametric models for the intensities of the traumatic failures are used supposing that these intensities depend on degradation level. Semi-parametric estimators of various reliability characteristics are proposed and asymptotic properties of the estimators are obtained. The theoretical results are illustrated using simulated data.