Regression analysis of multivariate interval-censored failure time data with application to tumorigenicity experiments

This paper discusses multivariate interval-censored failure time data that occur when there exist several correlated survival times of interest and only interval-censored data are available for each survival time. Such data occur in many fields. One is tumorigenicity experiments, which usually concern different types of tumors, tumors occurring in different locations of animals, or together. For regression analysis of such data, we develop a marginal inference approach using the additive hazards model and apply it to a set of bivariate interval-censored data arising from a tumorigenicity experiment. Simulation studies are conducted for the evaluation of the presented approach and suggest that the approach performs well for practical situations.     

Maximum likelihood estimation for interval-censored data using a Weibull-based accelerated failure time model.

The accelerated failure time regression model is most commonly used with right-censored survival data. This report studies the use of a Weibull-based accelerated failure time regression model when left- and interval-censored data are also observed. Two alternative methods of analysis are considered. First, the maximum likelihood estimates (MLEs) for the observed censoring pattern are computed. These are compared with estimates where midpoints are substituted for left- and interval-censored data (midpoint estimator, or MDE). Simulation studies indicate that for relatively large samples there are many instances when the MLE is superior to the MDE. For samples where the hazard rate is flat or nearly so, or where the percentage of interval-censored data is small, the MDE is adequate. An example using Framingham Heart Study data is discussed.     

Bayesian nonparametric analysis of restricted mean survival time

The restricted mean survival time (RMST) evaluates the expectation of survival time truncated by a prespecified time point, because the mean survival time in the presence of censoring is typically not estimable. The frequentist inference procedure for RMST has been widely advocated for comparison of two survival curves, while research from the Bayesian perspective is rather limited. For the RMST of both right- and interval-censored data, we propose Bayesian nonparametric estimation and inference procedures. By assigning a mixture of Dirichlet processes (MDP) prior to the distribution function, we can estimate the posterior distribution of RMST. We also explore another Bayesian nonparametric approach using the Dirichlet process mixture model and make comparisons with the frequentist nonparametric method. Simulation studies demonstrate that the Bayesian nonparametric RMST under diffuse MDP priors leads to robust estimation and under informative priors it can incorporate prior knowledge into the nonparametric estimator. Analysis of real trial examples demonstrates the flexibility and interpretability of the Bayesian nonparametric RMST for both right- and interval-censored data.     

Fast approximations of pseudo-observations in the context of right censoring and interval censoring

In the context of right-censored and interval-censored data, we develop asymptotic formulas to compute pseudo-observations for the survival function and the restricted mean survival time (RMST). These formulas are based on the original estimators and do not involve computation of the jackknife estimators. For right-censored data, Von Mises expansions of the Kaplan–Meier estimator are used to derive the pseudo-observations. For interval-censored data, a general class of parametric models for the survival function is studied. An asymptotic representation of the pseudo-observations is derived involving the Hessian matrix and the score vector. Theoretical results that justify the use of pseudo-observations in regression are also derived. The formula is illustrated on the piecewise-constant-hazard model for the RMST. The proposed approximations are extremely accurate, even for small sample sizes, as illustrated by Monte Carlo simulations and real data. We also study the gain in terms of computation time, as compared to the original jackknife method, which can be substantial for a large dataset.     

A goodness-of-fit test for the marginal Cox model for correlated interval-censored failure time data

The marginal Cox model approach is perhaps the most commonly used method in the analysis of correlated failure time data (Cai, 1999; Cai and Prentice, 1995; Lin, 1994; Wei, Lin and Weissfeld, 1989). It assumes that the marginal distributions for the correlated failure times can be described by the Cox model and leaves the dependence structure completely unspecified. This paper discusses the assessment of the marginal Cox model for correlated interval-censored data and a goodness-of-fit test is presented for the problem. The method is applied to a set of correlated interval-censored data arising from an AIDS clinical trial.     

Residuals for proportional hazards models with interval-censored survival data

We develop diagnostic tools for use with proportional hazards models for interval-censored survival data. We propose counterparts to the Cox-Snell, Lagakos (or martingale), deviance, and Schoenfeld residuals. Many of the properties of these residuals carry over to the interval-censored case. In particular, the interval-censored versions of the Lagakos and Schoenfeld residuals may be derived as components of suitable score statistics. The Lagakos residuals may be used to check regression relationships, while the Schoenfeld residuals can help to detect nonproportional hazards in semiparametric models. The methods apply to parametric models and to the semiparametric model with discrete observation times.     

Sieve Estimation for the Cox Model with Clustered Interval-Censored Failure Time Data

Clustered interval-censored failure time data occur when the failure times of interest are clustered into small groups and known only to lie in certain intervals. A number of methods have been proposed for regression analysis of clustered failure time data, but most of them apply only to clustered right-censored data. In this paper, a sieve estimation procedure is proposed for fitting a Cox frailty model to clustered interval-censored failure time data. In particular, a two-step algorithm for parameter estimation is developed and the asymptotic properties of the resulting sieve maximum likelihood estimators are established. The finite sample properties of the proposed estimators are investigated through a simulation study and the method is illustrated by the data arising from a lymphatic filariasis study.     

Parametric survival models for interval-censored data with time-dependent covariates

We present a parametric family of regression models for interval-censored event-time (survival) data that accomodates both fixed (e.g. baseline) and time-dependent covariates. The model employs a three-parameter family of survival distributions that includes the Weibull, negative binomial, and log-logistic distributions as special cases, and can be applied to data with left, right, interval, or non-censored event times. Standard methods, such as Newton-Raphson, can be employed to estimate the model and the resulting estimates have an asymptotically normal distribution about the true values with a covariance matrix that is consistently estimated by the information function. The deviance function is described to assess model fit and a robust sandwich estimate of the covariance may also be employed to provide asymptotically robust inferences when the model assumptions do not apply. Spline functions may also be employed to allow for non-linear covariates. The model is applied to data from a long-term study of type 1 diabetes to describe the effects of longitudinal measures of glycemia (HbA1c) over time (the time-dependent covariate) on the risk of progression of diabetic retinopathy (eye disease), an interval-censored event-time outcome.     

A proportional hazards model for multivariate interval-censored failure time data

This paper focuses on the methodology developed for analyzing a multivariate interval-censored data set from an AIDS observational study. A purpose of the study was to determine the natural history of the opportunistic infection cytomeglovirus (CMV) in an HIV-infected individual. For this observational study, laboratory tests were performed at scheduled clinic visits to test for the presence of the CMV virus in the blood and in the urine (called CMV shedding in the blood and urine). The study investigators were interested in determining whether the stage of HIV disease at study entry was predictive of an increased risk for CMV shedding in either the blood or the urine. If all patients had made each clinic visit, the data would be multivariate grouped failure time data and published methods could be used. However, many patients missed several visits, and when they returned, their lab tests indicated a change in their blood and/or urine CMV shedding status, resulting in interval-censored failure time data. This paper outlines a method for applying the proportional hazards model to the analysis of multivariate interval-censored failure time data from a study of CMV in HIV-infected patients.     

A semiparametric model for regression analysis of interval-censored failure time data

Left-, right-, and interval-censored response time data arise in a variety of settings, including the analyses of data from laboratory animal carcinogenicity experiments, clinical trials, and longitudinal studies. For such incomplete data, the usual regression techniques such as the Cox (1972, Journal of the Royal Statistical Society, Series B 34, 187-220) proportional hazards model are inapplicable. In this paper, we present a method for regression analysis which accommodates interval-censored data. We present applications of this methodology to data sets from a study of breast cancer patients who were followed for cosmetic response to therapy, a small animal tumorigenicity study, and a clinical trial.     

Regression analysis of doubly censored failure time data with frailty

In doubly censored failure time data, the survival time of interest is defined as the elapsed time between an initial event and a subsequent event, and the occurrences of both events cannot be observed exactly. Instead, only right- or interval-censored observations on the occurrence times are available. For the analysis of such data, a number of methods have been proposed under the assumption that the survival time of interest is independent of the occurrence time of the initial event. This article investigates a different situation where the independence may not be true with the focus on regression analysis of doubly censored data. Cox frailty models are applied to describe the effects of covariates and an EM algorithm is developed for estimation. Simulation studies are performed to investigate finite sample properties of the proposed method and an illustrative example from an acquired immune deficiency syndrome (AIDS) cohort study is provided.     

Survival analysis in clinical trials: past developments and future directions

The field of survival analysis emerged in the 20th century and experienced tremendous growth during the latter half of the century. The developments in this field that have had the most profound impact on clinical trials are the Kaplan-Meier (1958, Journal of the American Statistical Association 53, 457-481) method for estimating the survival function, the log-rank statistic (Mantel, 1966, Cancer Chemotherapy Report 50, 163-170) for comparing two survival distributions, and the Cox (1972, Journal of the Royal Statistical Society, Series B 34, 187-220) proportional hazards model for quantifying the effects of covariates on the survival time. The counting-process martingale theory pioneered by Aalen (1975, Statistical inference for a family of counting processes, Ph.D. dissertation, University of California, Berkeley) provides a unified framework for studying the small- and large-sample properties of survival analysis statistics. Significant progress has been achieved and further developments are expected in many other areas, including the accelerated failure time model, multivariate failure time data, interval-censored data, dependent censoring, dynamic treatment regimes and causal inference, joint modeling of failure time and longitudinal data, and Baysian methods.     

Analysis of doubly-censored survival data, with application to AIDS

This paper proposes nonparametric and weakly structured parametric methods for analyzing survival data in which both the time origin and the failure event can be right- or interval-censored. Such data arise in clinical investigations of the human immunodeficiency virus (HIV) when the infection and clinical status of patients are observed only at several time points. The proposed methods generalize the self-consistency algorithm proposed by Turnbull (1976, Journal of the Royal Statistical Society, Series B 38, 290-295) for singly-censored univariate data, and are illustrated with the results from a study of hemophiliacs who were infected with HIV by contaminated blood factor.     

Inference for set-based effects in genetic association studies with interval-censored outcomes

The rapid acceleration of genetic data collection in biomedical settings has recently resulted in the rise of genetic compendiums filled with rich longitudinal disease data. One common feature of these data sets is their plethora of interval-censored outcomes. However, very few tools are available for the analysis of genetic data sets with interval-censored outcomes, and in particular, there is a lack of methodology available for set-based inference. Set-based inference is used to associate a gene, biological pathway, or other genetic construct with outcomes and is one of the most popular strategies in genetics research. This work develops three such tests for interval-censored settings beginning with a variance components test for interval-censored outcomes, the interval-censored sequence kernel association test (ICSKAT). We also provide the interval-censored version of the Burden test, and then we integrate ICSKAT and Burden to construct the interval censored sequence kernel association test—optimal (ICSKATO) combination. These tests unlock set-based analysis of interval-censored data sets with analogs of three highly popular set-based tools commonly applied to continuous and binary outcomes. Simulation studies illustrate the advantages of the developed methods over ad hoc alternatives, including protection of the type I error rate at very low levels and increased power. The proposed approaches are applied to the investigation that motivated this study, an examination of the genes associated with bone mineral density deficiency and fracture risk.     

Nonparametric estimation of the joint distribution of a survival time subject to interval censoring and a continuous mark variable

This article considers three nonparametric estimators of the joint distribution function for a survival time and a continuous mark variable when the survival time is interval censored and the mark variable may be missing for interval-censored observations. Finite and large sample properties are described for the nonparametric maximum likelihood estimator (NPMLE) as well as estimators based on midpoint imputation (MIDMLE) and coarsening the mark variable (CMLE). The estimators are compared using data from a simulation study and a recent phase III HIV vaccine efficacy trial where the survival time is the time from enrollment to infection and the mark variable is the genetic distance from the infecting HIV sequence to the HIV sequence in the vaccine. Theoretical and empirical evidence are presented indicating the NPMLE and MIDMLE are inconsistent. Conversely, the CMLE is shown to be consistent in general and thus is preferred.     

Exploring Migration Data Using Interval-Censored Time-to-Event Models

Abstract: Ecologists and wildlife biologists rely on periodic observation of radiocollared animals to study habitat use, survival, movement, and migration, resulting in response times (e.g., mortality and migration) known only to occur within an interval of time. We illustrate methods for analyzing interval-censored data using data on the timing of fall migration (from spring-summer-fall to winter ranges) for white-tailed deer (Odocoileus virginianus) in northern Minnesota, USA, during years 1991–1992 to 2005–2006. We compare both nonparametric and parametric methods for estimating the cumulative distribution function of migration times, and we suggest a parametric (cure rate) model that accounts for conditional (facultative) migrators as a potential alternative to traditional parametric models. Lastly, we illustrate methods for exploring the effect of environmental covariates on migration timing. Models with time-dependent covariates (snow depth, temp) were sensitive to the treatment of the data (as interval-censored or known event times), suggesting the need to account for interval-censoring when modeling the effect of these covariates.     

A proportional hazards model for interval-censored failure time data

This paper develops a method for fitting the proportional hazards regression model when the data contain left-, right-, or interval-censored observations. Results given for testing the hypothesis of a zero regression coefficient lead to a generalization of the log-rank test for comparison of several survival curves. The method is used to analyze data from an animal tumorigenicity study and also a clinical trial.     

A simple local sensitivity analysis tool for nonignorable coarsening: application to dependent censoring

Right- and interval-censored data are common special cases of coarsened data (Heitjan and Rubin, 1991, Annals of Statistics19, 2244-2253). As with missing data, standard statistical methods that ignore the random nature of the coarsening mechanism may lead to incorrect inferences. We extend a simple sensitivity analysis tool, the index of local sensitivity to nonignorability (Troxel, Ma, and Heitjan, 2004, Statistica Sinica14, 1221-1237), to the evaluation of nonignorability of the coarsening process in the general coarse-data model. By converting this index into a simple graphical display one can easily assess the sensitivity of key inferences to nonignorable coarsening. We illustrate the validity of the method with a simulated example, and apply it to right-censored data from an observational study of cardiac transplantation and to interval-censored data on time to detectable viral load from a clinical trial in HIV disease.     

A modification of Peto's nonparametric estimation of survival curves for interval-censored data

Peto (1973, Applied Statistics, 22, 86-91) gave a nonparametric generalized maximum-likelihood estimate of the survival function for interval-censored data. His method has a tendency to concentrate probability masses at the endpoints of the intervals, even for the ordinary grouped data, instead of spreading them through the intervals, as one might expect them to be in the underlying distribution. We describe a modification that overcomes this. The new estimate reduces to the standard binomial estimate when applied to grouped data. It also reduces to the Kaplan-Meier estimate when applied to survival data that consist of only exact or right-censored observations. Both estimates are maximum-likelihood estimates but are based on different interpretations of the endpoints of the intervals.     

Nonparametric screening and feature selection for ultrahigh-dimensional Case II interval-censored failure time data

For the analysis of ultrahigh-dimensional data, the first step is often to perform screening and feature selection to effectively reduce the dimensionality while retaining all the active or relevant variables with high probability. For this, many methods have been developed under various frameworks but most of them only apply to complete data. In this paper, we consider an incomplete data situation, case II interval-censored failure time data, for which there seems to be no screening procedure. Basing on the idea of cumulative residual, a model-free or nonparametric method is developed and shown to have the sure independent screening property. In particular, the approach is shown to tend to rank the active variables above the inactive ones in terms of their association with the failure time of interest. A simulation study is conducted to demonstrate the usefulness of the proposed method and, in particular, indicates that it works well with general survival models and is capable of capturing the nonlinear covariates with interactions. Also the approach is applied to a childhood cancer survivor study that motivated this investigation.