A simple linear regression method for quantitative trait loci linkage analysis with censored observations

Standard quantitative trait loci (QTL) mapping techniques commonly assume that the trait is both fully observed and normally distributed. When considering survival or age-at-onset traits these assumptions are often incorrect. Methods have been developed to map QTL for survival traits; however, they are both computationally intensive and not available in standard genome analysis software packages. We propose a grouped linear regression method for the analysis of continuous survival data. Using simulation we compare this method to both the Cox and Weibull proportional hazards models and a standard linear regression method that ignores censoring. The grouped linear regression method is of equivalent power to both the Cox and Weibull proportional hazards methods and is significantly better than the standard linear regression method when censored observations are present. The method is also robust to the proportion of censored individuals and the underlying distribution of the trait. On the basis of linear regression methodology, the grouped linear regression model is computationally simple and fast and can be implemented readily in freely available statistical software.     

Bayesian analysis of survival data with missing censoring indicators

In some large clinical studies, it may be impractical to perform the physical examination to every subject at his/her last monitoring time in order to diagnose the occurrence of the event of interest. This gives rise to survival data with missing censoring indicators where the probability of missing may depend on time of last monitoring and some covariates. We present a fully Bayesian semi‐parametric method for such survival data to estimate regression parameters of the proportional hazards model of Cox. Theoretical investigation and simulation studies show that our method performs better than competing methods. We apply the proposed method to analyze the survival data with missing censoring indicators from the Orofacial Pain: Prospective Evaluation and Risk Assessment study.     

Cox regression model with randomly censored covariates

This paper deals with a Cox proportional hazards regression model, where some covariates of interest are randomly right‐censored. While methods for censored outcomes have become ubiquitous in the literature, methods for censored covariates have thus far received little attention and, for the most part, dealt with the issue of limit‐of‐detection. For randomly censored covariates, an often‐used method is the inefficient complete‐case analysis (CCA) which consists in deleting censored observations in the data analysis. When censoring is not completely independent, the CCA leads to biased and spurious results. Methods for missing covariate data, including type I and type II covariate censoring as well as limit‐of‐detection do not readily apply due to the fundamentally different nature of randomly censored covariates. We develop a novel method for censored covariates using a conditional mean imputation based on either Kaplan–Meier estimates or a Cox proportional hazards model to estimate the effects of these covariates on a time‐to‐event outcome. We evaluate the performance of the proposed method through simulation studies and show that it provides good bias reduction and statistical efficiency. Finally, we illustrate the method using data from the Framingham Heart Study to assess the relationship between offspring and parental age of onset of cardiovascular events.     

Semiparametric regression analysis of interval-censored data

We propose a semiparametric approach to the proportional hazards regression analysis of interval-censored data. An EM algorithm based on an approximate likelihood leads to an M-step that involves maximizing a standard Cox partial likelihood to estimate regression coefficients and then using the Breslow estimator for the unknown baseline hazards. The E-step takes a particularly simple form because all incomplete data appear as linear terms in the complete-data log likelihood. The algorithm of Turnbull (1976, Journal of the Royal Statistical Society, Series B 38, 290-295) is used to determine times at which the hazard can take positive mass. We found multiple imputation to yield an easily computed variance estimate that appears to be more reliable than asymptotic methods with small to moderately sized data sets. In the right-censored survival setting, the approach reduces to the standard Cox proportional hazards analysis, while the algorithm reduces to the one suggested by Clayton and Cuzick (1985, Applied Statistics 34, 148-156). The method is illustrated on data from the breast cancer cosmetics trial, previously analyzed by Finkelstein (1986, Biometrics 42, 845-854) and several subsequent authors.     

A statistical measure of tissue heterogeneity with application to 3D PET sarcoma data

In vivo measurement of local tissue characteristics by modern bioimaging techniques such as positron emission tomography (PET) provides the opportunity to analyze quantitatively the role that tissue heterogeneity may play in understanding biological function. This paper develops a statistical measure of the heterogeneity of a tissue characteristic that is based on the deviation of the distribution of the tissue characteristic from a unimodal elliptically contoured spatial pattern. An efficient algorithm is developed for computation of the measure based on volumetric region of interest data. The technique is illustrated by application to data from PET imaging studies of fluorodeoxyglucose utilization in human sarcomas. A set of 74 sarcoma patients (with five-year follow-up survival information) were evaluated for heterogeneity as well as a number of other potential prognostic indicators of survival. A Cox proportional hazards analysis of these data shows that the degree of heterogeneity of the sarcoma is the major risk factor associated with patient death. Some theory is developed to analyze the asymptotic statistical behavior of the heterogeneity estimator. In the context of data arising from Poisson deconvolution (PET being the prime example), the heterogeneity estimator, which is a non-linear functional of the PET image data, is consistent and converges at a rate that is parametric in the injected dose.     

Correlated versus uncorrelated frailty Cox models: A comparison of different estimation procedures

In many studies in medicine, including clinical trials and epidemiological investigations, data are clustered into groups such as health centers or herds in veterinary medicine. Such data are usually analyzed by hierarchical regression models to account for possible variation between groups. When such variation is large, it is of potential interest to explore whether additionally the effect of a within‐group predictor varies between groups. In survival analysis, this may be investigated by including two frailty terms at group level in a Cox proportional hazards model. Several estimation methods have been proposed to estimate this type of frailty Cox models. We review four of these methods, apply them to real data from veterinary medicine, and compare them using a simulation study.     

Approximate case influence for the proportional hazards regression model with censored data

A method is presented for approximating the influence of individual cases upon regression coefficient estimates obtained from the Cox proportional hazards model. Observations can thus be identified which may greatly influence statistical inferences regarding the effects of prognostic factors upon survival time. An example from a cancer clinical trial is given.     

Regression analysis of grouped survival data with application to breast cancer data

Use of the proportional hazards regression model (Cox 1972) substantially liberalized the analysis of censored survival data with covariates. Available procedures for estimation of the relative risk parameter, however, do not adequately handle grouped survival data, or large data sets with many tied failure times. The grouped data version of the proportional hazards model is proposed here for such estimation. Asymptotic likelihood results are given, both for the estimation of the regression coefficient and the survivor function. Some special results are given for testing the hypothesis of a zero regression coefficient which leads, for example, to a generalization of the log-rank test for the comparison of several survival curves. Application to breast cancer data, from the National Cancer Institute-sponsored End Results Group, indicates that previously noted race differences in breast cancer survival times are explained to a large extent by differences in disease extent and other demographic characteristics at diagnosis.     

Application of the Principal Components Method and the Proportional Hazards Regression Model to Analysis of Survival Data

One of factor analysis techniques, viz. the principal components method, and the proportional hazards regression model (Cox, 1972) are applied in this work to study the significance of various factors characterizing the patient, the disease, and the method of treatment in the survival. The application of these methods to analysis of survival data for cervical cancer patients has shown, in particular, the tumor growth rate to be the crucial factor in distribution of the patients survival time and to be even more important than the therapy characteristics.     

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.     

Cause‐Specific Cumulative Incidence Estimation and the Fine and Gray Model Under Both Left Truncation and Right Censoring

Summary The standard estimator for the cause‐specific cumulative incidence function in a competing risks setting with left truncated and/or right censored data can be written in two alternative forms. One is a weighted empirical cumulative distribution function and the other a product‐limit estimator. This equivalence suggests an alternative view of the analysis of time‐to‐event data with left truncation and right censoring: individuals who are still at risk or experienced an earlier competing event receive weights from the censoring and truncation mechanisms. As a consequence, inference on the cumulative scale can be performed using weighted versions of standard procedures. This holds for estimation of the cause‐specific cumulative incidence function as well as for estimation of the regression parameters in the Fine and Gray proportional subdistribution hazards model. We show that, with the appropriate filtration, a martingale property holds that allows deriving asymptotic results for the proportional subdistribution hazards model in the same way as for the standard Cox proportional hazards model. Estimation of the cause‐specific cumulative incidence function and regression on the subdistribution hazard can be performed using standard software for survival analysis if the software allows for inclusion of time‐dependent weights. We show the implementation in the R statistical package. The proportional subdistribution hazards model is used to investigate the effect of calendar period as a deterministic external time varying covariate, which can be seen as a special case of left truncation, on AIDS related and non‐AIDS related cumulative mortality.     

A prognostic five long–noncoding RNA signature for patients with rectal cancer

This study aimed to identify prognostic long noncoding RNAs (lncRNAs) signature for predicting the prognosis of patients with rectal cancer. LncRNA-sequencing data and clinicopathological data of patients with rectal cancer were retrieved from The Cancer Genome Atlas database. Univariate and multivariate Cox proportional hazards regression analysis, the least absolute shrinkage, and selection operator analysis and the Kaplan-Meier curve method were employed to identify prognostic lncRNAs and construct multi-lncRNA signature. Finally, five lncRNAs (AC079789.1, AC106900.2, AL121987.1, AP004609.1, and LINC02163) were identified to construct a five-lncRNA signature. According to the five-lncRNA signature, patients with rectal cancer were divided into a high-risk group and low-risk group. Patients with rectal cancer had significantly poorer overall survival in the high-risk group than in the low-risk group. We used a time-dependent receiver operating characteristic curve to assess the power of the five-lncRNA signature by calculating the area under the curve (AUC). The AUCs for predicting 3-year survival and 5-year survival were 0.742 and 0.935, respectively, which indicated a good performance of the five-lncRNA signature. The five-lncRNA signature was independently associated with the prognosis of patients with rectal cancer through using univariate and multivariate Cox regression analysis. The biological function of the five lncRNAs was enriched in some cancer-related biological processes and pathways by performing functional enrichment analysis of their correlated protein-coding genes. In conclusion, we developed a five-lncRNA signature as a potential indicator for rectal cancer.     

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.     

Applying a data duplication technique in linear regression analysis of waiting time to pregnancy

This analysis demonstrates the application of a data duplication technique in linear regression with censored observations of the waiting time to third pregnancy ending in two outcome types, using data from Malaysia. The linear model not only confirmed the results obtained by the Cox proportional hazards model, but also identified two additional significant factors. The method provides a useful alternative when Cox proportionality assumption of the hazards is violated.     

Bayesian cure rate frailty models with application to a root canal therapy study

Due to natural or artificial clustering, multivariate survival data often arise in biomedical studies, for example, a dental study involving multiple teeth from each subject. A certain proportion of subjects in the population who are not expected to experience the event of interest are considered to be "cured" or insusceptible. To model correlated or clustered failure time data incorporating a surviving fraction, we propose two forms of cure rate frailty models. One model naturally introduces frailty based on biological considerations while the other is motivated from the Cox proportional hazards frailty model. We formulate the likelihood functions based on piecewise constant hazards and derive the full conditional distributions for Gibbs sampling in the Bayesian paradigm. As opposed to the Cox frailty model, the proposed methods demonstrate great potential in modeling multivariate survival data with a cure fraction. We illustrate the cure rate frailty models with a root canal therapy data set.     

Variance Calculations for Direct Adjusted Survival Curves,with Applications to Testing for No Treatment Effect

Several investigators have recently constructed survival curves adjusted for imbalances in prognostic factors by a method which we call direct adjustment. We present methods for calculating variances of these direct adjusted survival curves and their differences. Estimates of the adjusted curves, their variances, and the variances of their differences are compared for non-parametric (Kaplan-Meier), semi-parametric (Cox) and parametric (Weibull) models applied to censored exponential data. Semi-parametric proportional hazards models were nearly fully efficient for estimating differences in adjusted curves, but parametric estimates of individual adjusted curves may be substantially more precise. Standardized differences between direct adjusted survival curves may be used to test the null hypothesis of no treatment effect. This procedure may prove especially useful when the proportional hazards assumption is questionable.     

Unlocking the potential of survival data for model organisms through a new database and online analysis platform: SurvCurv

Lifespan measurements, also called survival records, are a key phenotype in research on aging. If external hazards are excluded, aging alone determines the mortality in a population of model organisms. Understanding the biology of aging is highly desirable because of the benefits for the wide range of aging‐related diseases. However, it is also extremely challenging because of the underlying complexity. Here, we describe SurvCurv, a new database and online resource focused on model organisms collating survival data for storage and analysis. All data in SurvCurv are manually curated and annotated. The database, available at www.ebi.ac.uk/thornton-srv/databases/SurvCurv/ , offers various functions including plotting, Cox proportional hazards analysis, mathematical mortality models and statistical tests. It facilitates reanalysis and allows users to analyse their own data and compare it with the largest repository of model‐organism data from published experiments, thus unlocking the potential of survival data and demographics in model organisms.     

The Mizon-Richard encompassing test for the Cox and Aalen additive hazards models.

The Cox hazards model (Cox, 1972, Journal of the Royal Statistical Society, Series B34, 187-220) for survival data is routinely used in many applied fields, sometimes, however, with too little emphasis on the fit of the model. A useful alternative to the Cox model is the Aalen additive hazards model (Aalen, 1980, in Lecture Notes in Statistics-2, 1-25) that can easily accommodate time changing covariate effects. It is of interest to decide which of the two models that are most appropriate to apply in a given application. This is a nontrivial problem as these two classes of models are nonnested except only for special cases. In this article we explore the Mizon-Richard encompassing test for this particular problem. It turns out that it corresponds to fitting of the Aalen model to the martingale residuals obtained from the Cox regression analysis. We also consider a variant of this method, which relates to the proportional excess model (Martinussen and Scheike, 2002, Biometrika 89, 283-298). Large sample properties of the suggested methods under the two rival models are derived. The finite-sample properties of the proposed procedures are assessed through a simulation study. The methods are further applied to the well-known primary biliary cirrhosis data set.     

Joint Analysis of Survival Time and Longitudinal Categorical Outcomes

In biomedical or public health research, it is common for both survival time and longitudinal categorical outcomes to be collected for a subject, along with the subject’s characteristics or risk factors. Investigators are often interested in finding important variables for predicting both survival time and longitudinal outcomes which could be correlated within the same subject. Existing approaches for such joint analyses deal with continuous longitudinal outcomes. New statistical methods need to be developed for categorical longitudinal outcomes. We propose to simultaneously model the survival time with a stratified Cox proportional hazards model and the longitudinal categorical outcomes with a generalized linear mixed model. Random effects are introduced to account for the dependence between survival time and longitudinal outcomes due to unobserved factors. The Expectation–Maximization (EM) algorithm is used to derive the point estimates for the model parameters, and the observed information matrix is adopted to estimate their asymptotic variances. Asymptotic properties for our proposed maximum likelihood estimators are established using the theory of empirical processes. The method is demonstrated to perform well in finite samples via simulation studies. We illustrate our approach with data from the Carolina Head and Neck Cancer Study (CHANCE) and compare the results based on our simultaneous analysis and the separately conducted analyses using the generalized linear mixed model and the Cox proportional hazards model. Our proposed method identifies more predictors than by separate analyses.