Semi-supervised approach to event time annotation using longitudinal electronic health records

Large clinical datasets derived from insurance claims and electronic health record (EHR) systems are valuable sources for precision medicine research. These datasets can be used to develop models for personalized prediction of risk or treatment response. Efficiently deriving prediction models using real world data, however, faces practical and methodological challenges. Precise information on important clinical outcomes such as time to cancer progression are not readily available in these databases. The true clinical event times typically cannot be approximated well based on simple extracts of billing or procedure codes. Whereas, annotating event times manually is time and resource prohibitive. In this paper, we propose a two-step semi-supervised multi-modal automated time annotation (MATA) method leveraging multi-dimensional longitudinal EHR encounter records. In step I, we employ a functional principal component analysis approach to estimate the underlying intensity functions based on observed point processes from the unlabeled patients. In step II, we fit a penalized proportional odds model to the event time outcomes with features derived in step I in the labeled data where the non-parametric baseline function is approximated using B-splines. Under regularity conditions, the resulting estimator of the feature effect vector is shown as root-n consistent. We demonstrate the superiority of our approach relative to existing approaches through simulations and a real data example on annotating lung cancer recurrence in an EHR cohort of lung cancer patients from Veteran Health Administration.

Predicting outcomes of phase III oncology trials with Bayesian mediation modeling of tumor response

Pivotal cancer trials often fail to yield evidence in support of new therapies thought to offer promising alternatives to standards-of-care. Conducting randomized controlled trials in oncology tends to be considerably more expensive than studies of other diseases with comparable sample size. Moreover, phase III trial design often takes place with a paucity of survival data for experimental therapies. Experts have explained the failures on the basis of design flaws which produce studies with unrealistic expectations. This article presents a framework for predicting outcomes of phase III oncology trials using Bayesian mediation models. Predictions, which arise from interim analyses, derive from multivariate modeling of the relationships among treatment, tumor response, and their conjoint effects on survival. Acting as a safeguard against inaccurate pre-trial design assumptions, the methodology may better facilitate rapid closure of negative studies. Additionally the models can be used to inform re-estimations of sample size for under-powered trials that demonstrate survival benefit via tumor response mediation. The methods are applied to predict the outcomes of two colorectal cancer studies. Simulation is used to evaluate and compare models in the absence versus presence of reliable surrogate markers of survival.

Semiparametric regression analysis of partly interval-censored failure time data with application to an AIDS clinical trial

Failure time data subject to various types of censoring commonly arise in epidemiological and biomedical studies. Motivated by an AIDS clinical trial, we consider regression analysis of failure time data that include exact and left-, interval-, and/or right-censored observations, which are often referred to as partly interval-censored failure time data. We study the effects of potentially time-dependent covariates on partly interval-censored failure time via a class of semiparametric transformation models that includes the widely used proportional hazards model and the proportional odds model as special cases. We propose an EM algorithm for the nonparametric maximum likelihood estimation and show that it unifies some existing approaches developed for traditional right-censored data or purely interval-censored data. In particular, the proposed method reduces to the partial likelihood approach in the case of right-censored data under the proportional hazards model. We establish that the resulting estimator is consistent and asymptotically normal. In addition, we investigate the proposed method via simulation studies and apply it to the motivating AIDS clinical trial.

Regression analysis of arbitrarily censored survival data under the proportional odds model

Arbitrarily censored data are referred to as the survival data that contain a mixture of exactly observed, left-censored, interval-censored, and right-censored observations. Existing research work on regression analysis on arbitrarily censored data is relatively sparse and mainly focused on the proportional hazards model and the accelerated failure time model. This article studies the proportional odds (PO) model and proposes a novel estimation approach through an expectation-maximization (EM) algorithm for analyzing such data. The proposed EM algorithm has many appealing properties such as being robust to initial values, easy to implement, converging fast, and providing the variance estimate of the regression parameter estimate in closed form. An informal diagnosis plot is developed for checking the PO model assumption. Our method has shown excellent performance in estimating the regression parameters as well as the baseline survival function in a simulation study. A real-life dataset about metastatic colorectal cancer is analyzed for illustration. An R package regPO has been created for practitioners to implement our method.

Semiparametric regression analysis of interval-censored data with informative dropout

Interval-censored data arise when the event time of interest can only be ascertained through periodic examinations. In medical studies, subjects may not complete the examination schedule for reasons related to the event of interest. In this article, we develop a semiparametric approach to adjust for such informative dropout in regression analysis of interval-censored data. Specifically, we propose a broad class of joint models, under which the event time of interest follows a transformation model with a random effect and the dropout time follows a different transformation model but with the same random effect. We consider nonparametric maximum likelihood estimation and develop an EM algorithm that involves simple and stable calculations. We prove that the resulting estimators of the regression parameters are consistent, asymptotically normal, and asymptotically efficient with a covariance matrix that can be consistently estimated through profile likelihood. In addition, we show how to consistently estimate the survival function when dropout represents voluntary withdrawal and the cumulative incidence function when dropout is an unavoidable terminal event. Furthermore, we assess the performance of the proposed numerical and inferential procedures through extensive simulation studies. Finally, we provide an application to data on the incidence of diabetes from a major epidemiological cohort study.

Case-cohort studies with interval-censored failure time data

The case-cohort design has been widely used as a means of cost reduction in assembling or measuring expensive covariates in large cohort studies. The existing literature on the case-cohort design is mainly focused on right-censored data. In practice, however, the failure time is often subject to interval-censoring; it is known only to fall within some random time interval. In this paper, we consider the case-cohort study design for interval-censored failure time and develop a sieve semiparametric likelihood approach for analyzing data from this design under the proportional hazards model. We construct the likelihood function using inverse probability weighting and build the sieves with Bernstein polynomials. The consistency and asymptotic normality of the resulting regression parameter estimator are established and a weighted bootstrap procedure is considered for variance estimation. Simulations show that the proposed method works well for practical situations, and an application to real data is provided.

An Expectation Maximization algorithm for fitting the generalized odds-rate model to interval censored data

The generalized odds-rate model is a class of semiparametric regression models, which includes the proportional hazards and proportional odds models as special cases. There are few works on estimation of the generalized odds-rate model with interval censored data because of the challenges in maximizing the complex likelihood function. In this paper, we propose a gamma-Poisson data augmentation approach to develop an Expectation Maximization algorithm, which can be used to fit the generalized odds-rate model to interval censored data. The proposed Expectation Maximization algorithm is easy to implement and is computationally efficient. The performance of the proposed method is evaluated by comprehensive simulation studies and illustrated through applications to datasets from breast cancer and hemophilia studies. In order to make the proposed method easy to use in practice, an R package 'ICGOR' was developed. Copyright © 2016 John Wiley & Sons, Ltd.

Maximum likelihood estimation for semiparametric transformation models with interval-censored data

Interval censoring arises frequently in clinical, epidemiological, financial and sociological studies, where the event or failure of interest is known only to occur within an interval induced by periodic monitoring. We formulate the effects of potentially time-dependent covariates on the interval-censored failure time through a broad class of semiparametric transformation models that encompasses proportional hazards and proportional odds models. We consider nonparametric maximum likelihood estimation for this class of models with an arbitrary number of monitoring times for each subject. We devise an EM-type algorithm that converges stably, even in the presence of time-dependent covariates, and show that the estimators for the regression parameters are consistent, asymptotically normal, and asymptotically efficient with an easily estimated covariance matrix. Finally, we demonstrate the performance of our procedures through simulation studies and application to an HIV/AIDS study conducted in Thailand.

Fitting Cox Models with Doubly Censored Data Using Spline-Based Sieve Marginal Likelihood

In some applications, the failure time of interest is the time from an originating event to a failure event, while both event times are interval censored. We propose fitting Cox proportional hazards models to this type of data using a spline-based sieve maximum marginal likelihood, where the time to the originating event is integrated out in the empirical likelihood function of the failure time of interest. This greatly reduces the complexity of the objective function compared with the fully semiparametric likelihood. The dependence of the time of interest on time to the originating event is induced by including the latter as a covariate in the proportional hazards model for the failure time of interest. The use of splines results in a higher rate of convergence of the estimator of the baseline hazard function compared with the usual nonparametric estimator. The computation of the estimator is facilitated by a multiple imputation approach. Asymptotic theory is established and a simulation study is conducted to assess its finite sample performance. It is also applied to analyzing a real data set on AIDS incubation time.

Semiparametric odds rate model for modeling short-term and long-term effects with application to a breast cancer genetic study

The proportional odds model is commonly used in the analysis of failure time data. The assumption of constant odds ratios over time in the proportional odds model, however, can be violated in some applications. Motivated by a genetic study with breast cancer patients, we propose a novel semiparametric odds rate model for the analysis of right-censored survival data. The proposed model incorporates the short-term and long-term covariate effects on the failure time data and includes the proportional odds model as a nested model. We develop efficient likelihood-based inference procedures and establish the large sample properties of the proposed nonparametric maximum likelihood estimators. Simulation studies demonstrate that the proposed methods perform well in practical settings. An application to the motivating example is provided.

Semiparametric bayes' proportional odds models for current status data with underreporting

Current status data are a type of interval-censored event time data in which all the individuals are either left or right censored. For example, our motivation is drawn from a cross-sectional study, which measured whether or not fibroid onset had occurred by the age of an ultrasound exam for each woman. We propose a semiparametric Bayesian proportional odds model in which the baseline event time distribution is estimated nonparametrically by using adaptive monotone splines in a logistic regression model and the potential risk factors are included in the parametric part of the mean structure. The proposed approach has the advantage of being straightforward to implement using a simple and efficient Gibbs sampler, whereas alternative semiparametric Bayes' event time models encounter problems for current status data. The model is generalized to allow systematic underreporting in a subset of the data, and the methods are applied to an epidemiologic study of uterine fibroids.

Interval censoring

Interval-censored failure time data occur in many medical investigations as well as other studies such as demographical and sociological studies. They include the usual right-censored failure time data as a special case but provide much more complex structure and less relevant information than the right-censored data. This article reviews some basic concepts, issues and the corresponding statistical approaches related to the analysis of interval-censored data as well as recent advances. In particular, we discuss estimation of a survival function, comparison of several treatments and regression analysis as well as competing risks analysis and truncation in the presence of interval censoring. A well-known example of interval-censored data is described and analysed to illustrate some of the statistical procedures discussed.

A semiparametric probit model for case 2 interval-censored failure time data

Interval-censored data occur naturally in many fields and the main feature is that the failure time of interest is not observed exactly, but is known to fall within some interval. In this paper, we propose a semiparametric probit model for analyzing case 2 interval-censored data as an alternative to the existing semiparametric models in the literature. Specifically, we propose to approximate the unknown nonparametric nondecreasing function in the probit model with a linear combination of monotone splines, leading to only a finite number of parameters to estimate. Both the maximum likelihood and the Bayesian estimation methods are proposed. For each method, regression parameters and the baseline survival function are estimated jointly. The proposed methods make no assumptions about the observation process and can be applicable to any interval-censored data with easy implementation. The methods are evaluated by simulation studies and are illustrated by two real-life interval-censored data applications.

Linear transformation models for interval-censored data

In statistical analysis, when the value of a random variable is only known to be between two bounds, we say that this random variable is interval censored. This complicated censoring pattern is a common problem in research fields such as clinical trials or actuarial studies and raises challenges for statistical analysis. In this paper, we focus on regression analysis of case 2 interval-censored data. We first briefly review existing regression methods and an estimation approach under the class of linear transformation models developed by Zhang et al. We then propose a method for survival probability prediction via generalized estimating equations. We also consider a graphical model checking technique and a model selection tool. Some theoretical properties are established and the performance of our procedures is evaluated and illustrated by numerical studies including a real-life data analysis.

Semi-parametric accelerated failure time regression analysis with application to interval-censored HIV/AIDS data

This paper demonstrates a way to investigate a potentially non-linear relationship between an interval-censored response variable and a continuously distributed explanatory variable. A potentially non-linear effect of a continuous explanatory variable on the response is incorporated into an accelerated failure time model, forming a partial linear model. A sieve maximum likelihood estimator (MLE) is suggested to simultaneously estimate all the parameters. The sieve MLE is shown to be asymptotically efficient and normally distributed. Simulation studies show that the proposed estimators for the scale and regression parameters are robust and efficient, and the estimator for the non-linear function is able to capture the shape of a variety of smooth non-linear functions. The model is applied to observational HIV data, where the response variable is the time to suppression of HIV viral load after initiation of antiretroviral therapy, and baseline viral load is investigated as a potentially non-linear effect.

A baseline-free procedure for transformation models under interval censorship

An important property of Cox regression model is that the estimation of regression parameters using the partial likelihood procedure does not depend on its baseline survival function. We call such a procedure baseline-free. Using marginal likelihood, we show that an baseline-free procedure can be derived for a class of general transformation models under interval censoring framework. The baseline-free procedure results a simplified and stable computation algorithm for some complicated and important semiparametric models, such as frailty models and heteroscedastic hazard/rank regression models, where the estimation procedures so far available involve estimation of the infinite dimensional baseline function. A detailed computational algorithm using Markov Chain Monte Carlo stochastic approximation is presented. The proposed procedure is demonstrated through extensive simulation studies, showing the validity of asymptotic consistency and normality. We also illustrate the procedure with a real data set from a study of breast cancer. A heuristic argument showing that the score function is a mean zero martingale is provided.

Flexible parametric proportional-hazards and proportional-odds models for censored survival data, with application to prognostic modelling and estimation of treatment effects

Modelling of censored survival data is almost always done by Cox proportional-hazards regression. However, use of parametric models for such data may have some advantages. For example, non-proportional hazards, a potential difficulty with Cox models, may sometimes be handled in a simple way, and visualization of the hazard function is much easier. Extensions of the Weibull and log-logistic models are proposed in which natural cubic splines are used to smooth the baseline log cumulative hazard and log cumulative odds of failure functions. Further extensions to allow non-proportional effects of some or all of the covariates are introduced. A hypothesis test of the appropriateness of the scale chosen for covariate effects (such as of treatment) is proposed. The new models are applied to two data sets in cancer. The results throw interesting light on the behaviour of both the hazard function and the hazard ratio over time. The tools described here may be a step towards providing greater insight into the natural history of the disease and into possible underlying causes of clinical events. We illustrate these aspects by using the two examples in cancer.

Modeling multivariate survival data by a semiparametric random effects proportional odds model

In this article, the focus is on the analysis of multivariate survival time data with various types of dependence structures. Examples of multivariate survival data include clustered data and repeated measurements from the same subject, such as the interrecurrence times of cancer tumors. A random effect semiparametric proportional odds model is proposed as an alternative to the proportional hazards model. The distribution of the random effects is assumed to be multivariate normal and the random effect is assumed to act additively to the baseline log-odds function. This class of models, which includes the usual shared random effects model, the additive variance components model, and the dynamic random effects model as special cases, is highly flexible and is capable of modeling a wide range of multivariate survival data. A unified estimation procedure is proposed to estimate the regression and dependence parameters simultaneously by means of a marginal-likelihood approach. Unlike the fully parametric case, the regression parameter estimate is not sensitive to the choice of correlation structure of the random effects. The marginal likelihood is approximated by the Monte Carlo method. Simulation studies are carried out to investigate the performance of the proposed method. The proposed method is applied to two well-known data sets, including clustered data and recurrent event times data.

A semiparametric estimate of treatment effects with censored data

A semiparametric estimate of an average regression effect with right-censored failure time data has recently been proposed under the Cox-type model where the regression effect beta(t) is allowed to vary with time. In this article, we derive a simple algebraic relationship between this average regression effect and a measurement of group differences in k-sample transformation models when the random error belongs to the G(rho) family of Harrington and Fleming (1982, Biometrika 69, 553-566), the latter being equivalent to the conditional regression effect in a gamma frailty model. The models considered here are suitable for the attenuating hazard ratios that often arise in practice. The results reveal an interesting connection among the above three classes of models as alternatives to the proportional hazards assumption and add to our understanding of the behavior of the partial likelihood estimate under nonproportional hazards. The algebraic relationship provides a simple estimator under the transformation model. We develop a variance estimator based on the empirical influence function that is much easier to compute than the previously suggested resampling methods. When there is truncation in the right tail of the failure times, we propose a method of bias correction to improve the coverage properties of the confidence intervals. The estimate, its estimated variance, and the bias correction term can all be calculated with minor modifications to standard software for proportional hazards regression.