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.     

Analysis of times to repeated events in two-arm randomized trials with noncompliance and dependent censoring

This article develops randomization-based methods for times to repeated events in two-arm randomized trials with noncompliance and dependent censoring. Structural accelerated failure time models are assumed to capture causal effects on repeated event times and dependent censoring time, but the dependence structure among repeated event times and dependent censoring time is unspecified. Artificial censoring techniques to accommodate nonrandom noncompliance and dependent censoring are proposed. Estimation of the acceleration parameters are based on rank-based estimating functions. A simulation study is conducted to evaluate the performance of the developed methods. An illustration of the methods using data from an acute myeloid leukemia trial is provided.     

Semiparametric regression analysis of longitudinal data with informative drop-outs

Informative drop-out arises in longitudinal studies when the subject's follow-up time depends on the unobserved values of the response variable. We specify a semiparametric linear regression model for the repeatedly measured response variable and an accelerated failure time model for the time to informative drop-out. The error terms from the two models are assumed to have a common, but completely arbitrary joint distribution. Using a rank-based estimator for the accelerated failure time model and an artificial censoring device, we construct an asymptotically unbiased estimating function for the linear regression model. The resultant estimator is shown to be consistent and asymptotically normal. A resampling scheme is developed to estimate the limiting covariance matrix. Extensive simulation studies demonstrate that the proposed methods are suitable for practical use. Illustrations with data taken from two AIDS clinical trials are provided.     

Quality-adjusted survival estimation with periodic observations

Quality-adjusted survival is a measure that integrates both longevity and quality-of-life information. The analysis of quality-adjusted survival in a clinical study with data collected at periodic intervals encounters difficulties due to incomplete information. Based on observed time points, the time axis is partitioned into a set of disjoint time intervals, and under a Markovian assumption on patient's health status, the expected quality-adjusted survival is estimated as the summed product of the quality of life and its mean sojourn time of each health state within partitioned intervals. It is shown that the estimator is asymptotically normal with a simple variance calculation. A simulation study is conducted to investigate the behavior of the estimator, and a stroke study illustrates the use of the estimator.     

Sensitivity analysis for informative censoring in parametric survival models

Most statistical methods for censored survival data assume there is no dependence between the lifetime and censoring mechanisms, an assumption which is often doubtful in practice. In this paper we study a parametric model which allows for dependence in terms of a parameter delta and a bias function B(t, theta). We propose a sensitivity analysis on the estimate of the parameter of interest for small values of delta. This parameter measures the dependence between the lifetime and the censoring mechanisms. Its size can be interpreted in terms of a correlation coefficient between the two mechanisms. A medical example suggests that even a small degree of dependence between the failure and censoring processes can have a noticeable effect on the analysis.     

Testing for dependence between failure time and visit compliance with interval-censored data

Interval-censored failure-time data arise when subjects miss prescheduled visits at which the failure is to be assessed. The resulting intervals in which the failure is known to have occurred are overlapping. Most approaches to the analysis of these data assume that the visit-compliance process is ignorable with respect to likelihood analysis of the failure-time distribution. While this assumption offers considerable simplification, it is not always plausible. Here we test for dependence between the failure- and visit-compliance processes, applicable to studies in which data collection continues after the occurrence of the failure. We do not make any of the assumptions made by previous authors about the joint distribution of the visit-compliance process, a covariate process, and the failure time. Instead, we consider conditional models of the true failure history given the current visit compliance at each visit time, allowing for correlation across visit times. Because failure status is not known at some visit times due to missed visits, only models of the observed failure history given current visit compliance are estimable. We describe how the parameters from these models can be used to test for a negative association and how bounds on unestimable parameters provided by the observed data are needed additionally to infer a positive association. We illustrate the method with data from an AIDS study and we investigate the power of the test through a simulation study.     

Survival analysis of mortality in pre-weaning kids of Sirohi goat

Pre-weaning animals exit a flock through death induced by various reasons, causing significant economic losses to the goat producers. In this study, we investigated the survival from birth to weaning of Sirohi goat kids within framework of the survival analysis. Kid records were accessed from 1997 to 2017, with the information on 4417 pre-weaning animals of farmed Sirohi goat native to the Rajasthan State of India. A multivariable Cox regression was fitted to the data after checking the assumptions of regression. The explanatory variables were sex, type of birth, season of birth, birthweight, doe weight at kidding and year of birth. Model selection eliminated doe weight from the model, and sex, type of birth, season of birth, birthweight and year of birth were retained in the model. With model calibration also, these five covariates were retained in the model. The mortality on the first day after birth was 0.3%, constituting 3.5% of all pre-weaning mortality. The mortality until the end of weaning period was 7.8%. Regression analysis revealed that the higher birthweight at kidding was associated with reduced hazard of death among the kids. Male kids had higher hazards of death compared with female kids. The single-born kids had lower risks of death compared with twin-born kids after accounting for heterogeneity. The winter season had a very high adverse effect on the survival of the kids. With each passing year, risks of death decreased. The results of this study indicate that better survival of kids can be achieved by controlling both environmental and animal-related factors.     

Nonparametric inference for median costs with censored data

Summary Increasingly, estimations of health care costs are used to evaluate competing treatments or to assess the expected expenditures associated with certain diseases. In health policy and economics, the primary focus of these estimations has been on the mean cost, because the total cost can be derived directly from the mean cost, and because information about total resources utilized is highly relevant for policymakers. Yet, the median cost also could be important, both as an intuitive measure of central tendency in cost distribution and as a subject of interest to payers and consumers. In many prospective studies, cost data collection is sometimes incomplete for some subjects due to right censoring, which typically is caused by loss to follow-up or by limited study duration. Censoring poses a unique challenge for cost data analysis because of so-called induced informative censoring, in that traditional methods suited for survival data generally are invalid in censored cost estimation. In this article, we propose methods for estimating the median cost and its confidence interval (CI) when data are subject to right censoring. We also consider the estimation of the ratio and difference of two median costs and their CIs. These methods can be extended to the estimation of other quantiles and other informatively censored data. We conduct simulation and real data analysis in order to examine the performance of the proposed methods.     

An Asymptotically Optimal Adaptive Selection Procedure in the Proportional Hazards Model with Conditionally Independent Censoring

Assume k independent populations are given which are distributed according to R, …,R (ϑ_i ∈ Θ ⊆ R ). Taking samples of size n the population with the smallest ϑ-value is to be selected. Using the framework of Le Cam's decision theory (Le Cam , 1986; Strasser , 1985) under mild regularity assumptions, an asymptotically optimal selection procedure is derived for the sequence of localized models. In the proportional hazards model with conditionally independent censoring, an asymptotically optimal adaptive selection procedure is constructed by substituting the unknown nuisance parameter by a kernel estimator.     

Insider Censoring: Distortion of Data with Nondetects

Environmental data often include low-level concentrations below reporting limits. These data may be reported as “< RL,” where RL is one of several types of reporting limits. Some values also may be reported as a single number, but flagged with a qualifier (J-values) to indicate a difference in precision as compared to values above the RL. A currently used method for reporting censored environmental data called “insider censoring” produces a strong upward bias, while also distorting the shape of the data distribution. This results in inaccurate estimates of summary statistics and regression coefficients, distorts evaluations of whether data follow a normal distribution, and introduces inaccuracies into risk assessments and models. Insider censoring occurs when data measured as below the detection limit (< DL) are reported as less than the higher quantitation limit (< QL), whereas values between the DL and QL are reported as individual numbers. Three unbiased alternatives to insider censoring are presented so that laboratories and their data users can recognize, and remedy, this problem.     

Additive hazard regression with auxiliary covariates

We consider the additive hazard model when some of the truecovariates are measured only on a randomly selected validationset whereas auxiliary covariates are observed for all studysubjects. An updated pseudoscore estimation approach is proposedfor the parameters of the additive hazard model. It allows oneto fit the model with auxiliary covariates, while leaving thebaseline hazard unspecified. Asymptotic properties of the proposedestimators are established, and consistent standard error estimatorsare developed. Simulations demonstrate that the asymptotic approximationsof the proposed estimates are adequate for practical use. Areal example is used to illustrate the performance of the proposedmethod.     

Use of summary measures to adjust for informative missingness in repeated measures data with random effects

We discuss how to apply the conditional informative missing model of Wu and Bailey (1989, Biometrics 45, 939-955) to the setting where the probability of missing a visit depends on the random effects of the primary response in a time-dependent fashion. This includes the case where the probability of missing a visit depends on the true value of the primary response. Summary measures for missingness that are weighted sums of the indicators of missed visits are derived for these situations. These summary measures are then incorporated as covariates in a random effects model for the primary response. This approach is illustrated by analyzing data collected from a trial of heroin addicts where missed visits are informative about drug test results. Simulations of realistic experiments indicate that these time-dependent summary measures also work well under a variety of informative censoring models. These summary measures can achieve large reductions in estimation bias and mean squared errors relative to those obtained by using other summary measures.