Small sample study on estimators of life distributions and of mean survival times using randomly censored data

Whereas large-sample properties of the estimators of survival distributions using censored data have been studied by many authors, exact results for small samples have been difficult to obtain. In this paper we obtain the exact expression for the ath moment (a > 0) of the Bayes estimator of survival distribution using the censored data under proportional hazard model. Using the exact expression we compute the exact mean, variance and MSE of the Bayes estimator. Also two estimators ofthe mean survival time based on the Kaplan-Meier estimator and the Bayes estimator are compared for small samples under proportional hazards.     

Nonparametric bayes estimation of the survival function from a record of failures and follow-ups

In some observational studies, we have random censoring model. However, the data available may be partially observable censored data consisting of the observed failure times and only those nonfailure times which are subject to follow-up. Suzuki (1985) discussed the problem of nonparametric estimation of the survival function from such partially observable censored data. In this article, we derive a nonparametric Bayes estimator of the survival function for such data of failures and follow-ups under a Dirichlet process prior and squared error loss. The limiting properties such as the mean square consistency, weak convergence and strong consistency of the Bayes estimator are studied. Finally, the procedures developed are illustrated by means of an example.     

A semi-parametric estimation of a survival function from incomplete and doubly censored data

The purpose of this paper is to present a semi-parametric estimation of a survival function when analyzing incomplete and doubly censored data. Under the assumption that the chance of censoring is not related to the individual's survivorship, we propose a consistent estimation of survival. The derived estimator treats the uncensored observations nonparametrically and uses parametric models for both right and left censored data. Some asymptotic properties and simulation studies are also presented in order to analyze the behavior of the proposed estimator.     

Aalen's linear model for doubly censored data

Double censoring often occurs in registry studies when left censoring is present in addition to right censoring. In this work, we examine estimation of Aalen's nonparametric regression coefficients based on doubly censored data. We propose two estimation techniques. The first type of estimators, including ordinary least squared (OLS) estimator and weighted least squared (WLS) estimators, are obtained using martingale arguments. The second type of estimator, the maximum likelihood estimator (MLE), is obtained via expectation-maximization (EM) algorithms that treat the survival times of left censored observations as missing. Asymptotic properties, including the uniform consistency and weak convergence, are established for the MLE. Simulation results demonstrate that the MLE is more efficient than the OLS and WLS estimators.     

The estimation of treatment effect for the censored bivariate data under the univariate censoring

This paper establishes a nonparametric estimator for the treatment effect on censored bivariate data under unvariate censoring. This proposed estimator is based on the one from Lin and Ying(1993)'s nonparametric bivariate survival function estimator, which is itself a generalized version of Park and Park(1995)' quantile estimator. A Bahadur type representation of quantile functions were obtained from the marginal survival distribution estimator of Lin and Ying' model. The asymptotic property of this estimator is shown below and the simulation studies are also given     

A Nonparametric Estimator of Survival Functions for Arbitrarily Truncated and Censored Data

It is well-known that the nonparametric maximum likelihood estimator (NPMLE) may severely under-estimate the survival function with left truncated data. Based on the Nelson estimator (for right censored data) and self-consistency we suggest a nonparametric estimator of the survival function, the iterative Nelson estimator (INE), for arbitrarily truncated and censored data, where only few nonparametric estimators are available. By simulation we show that the INE does well in overcoming the under-estimation of the survival function from the NPMLE for left-truncated and interval-censored data. An interesting application of the INE is as a diagnostic tool for other estimators, such as the monotone MLE or parametric MLEs. The methodology is illustrated by application to two real world problems: the Channing House and the Massachusetts Health Care Panel Study data sets.     

Parametric Estimation of Menarcheal Age Distribution Based on Recall Data

Menarche, the onset of menstruation, is an important maturational event of female childhood. Most of the studies of age at menarche make use of dichotomous (status quo) data. More information can be harnessed from recall data, but such data are often censored in a informative way. We show that the usual maximum likelihood estimator based on interval censored data, which ignores the informative nature of censoring, can be biased and inconsistent. We propose a parametric estimator of the menarcheal age distribution on the basis of a realistic model of the recall phenomenon. We identify the additional information contained in the recall data and demonstrate theoretically as well as through simulations the advantage of the maximum likelihood estimator based on recall data over that based on status quo data.     

Robust and efficient parameter estimation based on censored data with stochastic covariates

Analysis of random censored life-time data along with some related stochastic covariables is of great importance in many applied sciences. The parametric estimation technique commonly used under this set-up is based on the efficient but non-robust likelihood approach. In this paper, we propose a robust parametric estimator for censored data with stochastic covariates based on the minimum density power divergence approach. The resulting estimator also has competitive efficiency with respect to the maximum likelihood estimator under pure data. The strong robustness property of the proposed estimator with respect to the presence of outliers is examined and illustrated through an appropriate real data example and simulation studies. Further, the theoretical asymptotic properties of the proposed estimator are also derived in terms of a general class of M-estimators based on the estimating equation.     

Estimation of the association for bivariate doubly censored data

We study the problem of estimating the association between two related survival variables when they follow a copula model and the bivariate doubly censored data is available. A two-stage estimation procedure is proposed and the asymptotic properties of the proposed estimator are established. Simulation studies are conducted to investigate the finite sample properties of the proposed estimate.     

Estimation of the maximal moment exponent with censored data

Heavy-tailed distributions have been used to model phenomena in which extreme events occur with high probability. In these type of occurrences, it is likely that extreme events are not observable after a certain threshold. Appropriate estimators are needed to deal with this type of censored data. We show that the well-known Hill-Hall estimator is unable to deal with censored data and yields highly biased estimates. We propose and study an unbiased modified maximum likelihood estimator, as well as a truncated tail regression estimator. We assess the expected value and the variance of these estimators in the cases of stable- and Pareto-distributed data.     

Estimation Of A Survival Function With Doubly Censored Data And Dirichlet Process Prior Knowledge On The Observable Variable

This article presents a nonparametric Bayesian procedure for estimating a survival curve in a proportional hazard model when some of the data are censored on the left and some are censored on the right. The method works under the assumption that there is a Dirichlet process prior knowledge on the observable variable. Strong consistency of the estimator is proved and an example is given. To finish some simulation is presented to analyze the estimator.     

A copula model for bivariate hybrid censored survival data with application to the MACS study

A copula model for bivariate survival data with hybrid censoring is proposed to study the association between survival time of individuals infected with HIV and persistence time of infection with an additional virus. Survival with HIV is right censored and the persistence time of the additional virus is subject to interval censoring case 1. A pseudo-likelihood method is developed to study the association between the two event times under such hybrid censoring. Asymptotic consistency and normality of the pseudo-likelihood estimator are established based on empirical process theory. Simulation studies indicate good performance of the estimator with moderate sample size. The method is applied to a motivating HIV study which investigates the effect of GB virus type C (GBV-C) co-infection on survival time of HIV infected individuals.     

Inverse censoring weighted median regression

We implement semiparametric random censorship model aided inference for censored median regression models. This is based on the idea that, when the censoring is specified by a common distribution, a semiparametric survival function estimator acts as an improved weight in the so-called inverse censoring weighted estimating function. We show that the proposed method will always produce estimates of the model parameters that are as good as or better than an existing estimator based on the traditional Kaplan–Meier weights. We also provide an illustration of the method through an analysis of a lung cancer data set.     

Simple non-parametric estimators for unemployment duration analysis

Summary.  We consider an extension of conventional univariate Kaplan–Meier-type estimators for the hazard rate and the survivor function to multivariate censored data with a censored random regressor. It is an Akritas-type estimator which adapts the non-parametric conditional hazard rate estimator of Beran to more typical data situations in applied analysis. We show with simulations that the estimator has nice finite sample properties and our implementation appears to be fast. As an application we estimate non-parametric conditional quantile functions with German administrative unemployment duration data.     

Inference for proportional hazard model with propensity score

Since the publication of the seminal paper by Cox (1972), proportional hazard model has become very popular in regression analysis for right censored data. In observational studies, treatment assignment may depend on observed covariates. If these confounding variables are not accounted for properly, the inference based on the Cox proportional hazard model may perform poorly. As shown in Rosenbaum and Rubin (1983), under the strongly ignorable treatment assignment assumption, conditioning on the propensity score yields valid causal effect estimates. Therefore we incorporate the propensity score into the Cox model for causal inference with survival data. We derive the asymptotic property of the maximum partial likelihood estimator when the model is correctly specified. Simulation results show that our method performs quite well for observational data. The approach is applied to a real dataset on the time of readmission of trauma patients. We also derive the asymptotic property of the maximum partial likelihood estimator with a robust variance estimator, when the model is incorrectly specified.     

Survivor Function Estimators Under Group Sequential Monitoring Based on the Logrank Statistic

In this paper we investigate a group sequential analysis of censored survival data with staggered entry, in which the trial is monitored using the logrank test while comparisons of treatment and control Kaplan-Meier curves at various time points are performed at the end of the trial. We concentrate on two-sample tests under local alternatives. We describe the relationship of the asymptotic bias of Kaplan-Meier curves between the two groups. We show that even if the asymptotic bias of the Kaplan-Meier curve is negligible relative to the true survival, this is not the case for the difference between the curves of the two arms of the trial. A corrected estimator for the difference between the survival curves is presented and by simulations we show that the corrected estimator reduced the bias dramatically and has a smaller variance. The methods of estimation are applied to the Beta-Blocker Heart Attack Trial (1982), a well-known group sequential trial.     

The Gini concentration test for survival data

We apply the well known Gini index to the measurement of concentration in survival times within groups of patients, and as a way to compare the distribution of survival times across groups of patients in clinical studies. In particular, we propose an estimator of a restricted version of the index from right censored data. We derive the asymptotic distribution of the resulting Gini statistic, and construct an estimator for its asymptotic variance. We use these results to propose a novel test for differences in the heterogeneity of survival distributions, which may suggest the presence of a differential treatment effect for some groups of patients. We focus in particular on traditional and generalized cure rate models, i.e., mixture models with a distribution of the lifetimes of the cured patients that is either degenerate at infinity or has a density. Results from a simulation study suggest that the Gini index is useful in some situations, and that it should be considered together with existing tests (in particular, the Log-rank, Wilcoxon, and Gray–Tsiatis tests). Use of the test is illustrated on the classic data arising from the Eastern Cooperative Oncology Group melanoma clinical trial E1690.     

Estimation of the bivariate distribution function for censored gap times

In many medical studies, patients may experience several events during follow-up. The times between consecutive events (gap times) are often of interest and lead to problems that have received much attention recently. In this work, we consider the estimation of the bivariate distribution function for censored gap times. Some related problems such as the estimation of the marginal distribution of the second gap time and the conditional distribution are also discussed. In this article, we introduce a nonparametric estimator of the bivariate distribution function based on Bayes’ theorem and Kaplan–Meier survival function and explore the behavior of the four estimators through simulations. Real data illustration is included.     

Kernel Survival Function Estimation Based on Doubly Censored Data

This work concerns the estimation of a smooth survival function based on doubly censored data. We establish strong consistency and asymptotic normality for a kernel estimator. Moreover, we also obtain an asymptotic expression for the mean integrated squared error, which yields an optimum bandwidth in terms of readily estimable quantities.     

Double robust estimation in longitudinal marginal structural models

In this article we consider estimation of causal parameters in a marginal structural model for the discrete intensity of the treatment specific counting process (e.g. hazard of a treatment specific survival time) based on longitudinal observational data on treatment, covariates and survival. We define three estimators: the inverse probability of treatment weighted (IPTW) estimator, the maximum likelihood estimator (MLE), and a double robust (DR) estimator. The DR estimator is obtained by following a general methodology for constructing double robust estimating functions in censored data models as described in van der Laan and Robins (Unified Methods for Censored Longitudinal Data and Causality, 2002). The double-robust estimator is consistent and asymptotically linear when either the treatment mechanism or the partial likelihood of the observed data is consistently estimated. We illustrate the superiority of the DR estimator relative to the IPTW and ML estimators in a simulation study. The proposed methodology is also applied to estimate the causal effect of exercise on physical functioning in a longitudinal study of seniors in Sonoma County.