Joint two‐part Tobit models for longitudinal and time‐to‐event data

In this article, we show how Tobit models can address problems of identifying characteristics of subjects having left‐censored outcomes in the context of developing a method for jointly analyzing time‐to‐event and longitudinal data. There are some methods for handling these types of data separately, but they may not be appropriate when time to event is dependent on the longitudinal outcome, and a substantial portion of values are reported to be below the limits of detection. An alternative approach is to develop a joint model for the time‐to‐event outcome and a two‐part longitudinal outcome, linking them through random effects. This proposed approach is implemented to assess the association between the risk of decline of CD4/CD8 ratio and rates of change in viral load, along with discriminating between patients who are potentially progressors to AIDS from patients who do not. We develop a fully Bayesian approach for fitting joint two‐part Tobit models and illustrate the proposed methods on simulated and real data from an AIDS clinical study.     

Joint modeling of longitudinal and survival data with missing and left‐censored time‐varying covariates

We propose a joint model for longitudinal and survival data with time‐varying covariates subject to detection limits and intermittent missingness at random. The model is motivated by data from the Multicenter AIDS Cohort Study (MACS), in which HIV+ subjects have viral load and CD4 cell count measured at repeated visits along with survival data. We model the longitudinal component using a normal linear mixed model, modeling the trajectory of CD4 cell count by regressing on viral load, and other covariates. The viral load data are subject to both left censoring because of detection limits (17%) and intermittent missingness (27%). The survival component of the joint model is a Cox model with time‐dependent covariates for death because of AIDS. The longitudinal and survival models are linked using the trajectory function of the linear mixed model. A Bayesian analysis is conducted on the MACS data using the proposed joint model. The proposed method is shown to improve the precision of estimates when compared with alternative methods. Copyright © 2014 John Wiley & Sons, Ltd.     

Bayesian joint ordinal and survival modeling for breast cancer risk assessment

We propose a joint model to analyze the structure and intensity of the association between longitudinal measurements of an ordinal marker and time to a relevant event. The longitudinal process is defined in terms of a proportional‐odds cumulative logit model. Time‐to‐event is modeled through a left‐truncated proportional‐hazards model, which incorporates information of the longitudinal marker as well as baseline covariates. Both longitudinal and survival processes are connected by means of a common vector of random effects. General inferences are discussed under the Bayesian approach and include the posterior distribution of the probabilities associated to each longitudinal category and the assessment of the impact of the baseline covariates and the longitudinal marker on the hazard function. The flexibility provided by the joint model makes possible to dynamically estimate individual event‐free probabilities and predict future longitudinal marker values. The model is applied to the assessment of breast cancer risk in women attending a population‐based screening program. The longitudinal ordinal marker is mammographic breast density measured with the Breast Imaging Reporting and Data System (BI‐RADS) scale in biennial screening exams. © 2016 The Authors. Statistics in Medicine Published by John Wiley & Sons Ltd.     

Bayesian joint modeling of longitudinal and spatial survival AIDS data

Joint analysis of longitudinal and survival data has received increasing attention in the recent years, especially for analyzing cancer and AIDS data. As both repeated measurements (longitudinal) and time‐to‐event (survival) outcomes are observed in an individual, a joint modeling is more appropriate because it takes into account the dependence between the two types of responses, which are often analyzed separately. We propose a Bayesian hierarchical model for jointly modeling longitudinal and survival data considering functional time and spatial frailty effects, respectively. That is, the proposed model deals with non‐linear longitudinal effects and spatial survival effects accounting for the unobserved heterogeneity among individuals living in the same region. This joint approach is applied to a cohort study of patients with HIV/AIDS in Brazil during the years 2002–2006. Our Bayesian joint model presents considerable improvements in the estimation of survival times of the Brazilian HIV/AIDS patients when compared with those obtained through a separate survival model and shows that the spatial risk of death is the same across the different Brazilian states. Copyright © 2016 John Wiley & Sons, Ltd.     

Bayesian quantile regression‐based nonlinear mixed‐effects joint models for time‐to‐event and longitudinal data with multiple features

This article explores Bayesian joint models for a quantile of longitudinal response, mismeasured covariate and event time outcome with an attempt to (i) characterize the entire conditional distribution of the response variable based on quantile regression that may be more robust to outliers and misspecification of error distribution; (ii) tailor accuracy from measurement error, evaluate non‐ignorable missing observations, and adjust departures from normality in covariate; and (iii) overcome shortages of confidence in specifying a time‐to‐event model. When statistical inference is carried out for a longitudinal data set with non‐central location, non‐linearity, non‐normality, measurement error, and missing values as well as event time with being interval censored, it is important to account for the simultaneous treatment of these data features in order to obtain more reliable and robust inferential results. Toward this end, we develop Bayesian joint modeling approach to simultaneously estimating all parameters in the three models: quantile regression‐based nonlinear mixed‐effects model for response using asymmetric Laplace distribution, linear mixed‐effects model with skew‐t distribution for mismeasured covariate in the presence of informative missingness and accelerated failure time model with unspecified nonparametric distribution for event time. We apply the proposed modeling approach to analyzing an AIDS clinical data set and conduct simulation studies to assess the performance of the proposed joint models and method. Copyright © 2016 John Wiley & Sons, Ltd.     

Investigation of 2‐stage meta‐analysis methods for joint longitudinal and time‐to‐event data through simulation and real data application

Background

Joint modelling of longitudinal and time‐to‐event data is often preferred over separate longitudinal or time‐to‐event analyses as it can account for study dropout, error in longitudinally measured covariates, and correlation between longitudinal and time‐to‐event outcomes. The joint modelling literature focuses mainly on the analysis of single studies with no methods currently available for the meta‐analysis of joint model estimates from multiple studies.

Methods

We propose a 2‐stage method for meta‐analysis of joint model estimates. These methods are applied to the INDANA dataset to combine joint model estimates of systolic blood pressure with time to death, time to myocardial infarction, and time to stroke. Results are compared to meta‐analyses of separate longitudinal or time‐to‐event models. A simulation study is conducted to contrast separate versus joint analyses over a range of scenarios.

Results

Using the real dataset, similar results were obtained by using the separate and joint analyses. However, the simulation study indicated a benefit of use of joint rather than separate methods in a meta‐analytic setting where association exists between the longitudinal and time‐to‐event outcomes.

Conclusions

Where evidence of association between longitudinal and time‐to‐event outcomes exists, results from joint models over standalone analyses should be pooled in 2‐stage meta‐analyses.     

A joint logistic regression and covariate‐adjusted continuous‐time Markov chain model

The use of longitudinal measurements to predict a categorical outcome is an increasingly common goal in research studies. Joint models are commonly used to describe two or more models simultaneously by considering the correlated nature of their outcomes and the random error present in the longitudinal measurements. However, there is limited research on joint models with longitudinal predictors and categorical cross‐sectional outcomes. Perhaps the most challenging task is how to model the longitudinal predictor process such that it represents the true biological mechanism that dictates the association with the categorical response. We propose a joint logistic regression and Markov chain model to describe a binary cross‐sectional response, where the unobserved transition rates of a two‐state continuous‐time Markov chain are included as covariates. We use the method of maximum likelihood to estimate the parameters of our model. In a simulation study, coverage probabilities of about 95%, standard deviations close to standard errors, and low biases for the parameter values show that our estimation method is adequate. We apply the proposed joint model to a dataset of patients with traumatic brain injury to describe and predict a 6‐month outcome based on physiological data collected post‐injury and admission characteristics. Our analysis indicates that the information provided by physiological changes over time may help improve prediction of long‐term functional status of these severely ill subjects. Copyright © 2017 John Wiley & Sons, Ltd.     

A Bayesian mixture of semiparametric mixed‐effects joint models for skewed‐longitudinal and time‐to‐event data

In longitudinal studies, it is of interest to investigate how repeatedly measured markers in time are associated with a time to an event of interest, and in the mean time, the repeated measurements are often observed with the features of a heterogeneous population, non‐normality, and covariate measured with error because of longitudinal nature. Statistical analysis may complicate dramatically when one analyzes longitudinal–survival data with these features together. Recently, a mixture of skewed distributions has received increasing attention in the treatment of heterogeneous data involving asymmetric behaviors across subclasses, but there are relatively few studies accommodating heterogeneity, non‐normality, and measurement error in covariate simultaneously arose in longitudinal–survival data setting. Under the umbrella of Bayesian inference, this article explores a finite mixture of semiparametric mixed‐effects joint models with skewed distributions for longitudinal measures with an attempt to mediate homogeneous characteristics, adjust departures from normality, and tailor accuracy from measurement error in covariate as well as overcome shortages of confidence in specifying a time‐to‐event model. The Bayesian mixture of joint modeling offers an appropriate avenue to estimate not only all parameters of mixture joint models but also probabilities of class membership. Simulation studies are conducted to assess the performance of the proposed method, and a real example is analyzed to demonstrate the methodology. The results are reported by comparing potential models with various scenarios. Copyright © 2015 John Wiley & Sons, Ltd.     

A Bayesian joint model for zero‐inflated integers and left‐truncated event times with a time‐varying association: Applications to senior health care

Population aging in most industrialized societies has led to a dramatic increase in emergency medical demand among the elderly. In the context of private health care, an optimal allocation of the medical resources for seniors is commonly done by forecasting their life spans. Accounting for each subject's particularities is therefore indispensable, so the available data must be processed at an individual level. We use a large and unique dataset of insured parties aged 65 and older to appropriately relate the emergency care usage with mortality risk. Longitudinal and time‐to‐event processes are jointly modeled, and their underlying relationship can therefore be assessed. Such an application, however, requires some special features to also be considered. First, longitudinal demand for emergency services exhibits a nonnegative integer response with an excess of zeros due to the very nature of the data. These subject‐specific responses are handled by a zero‐inflated version of the hierarchical negative binomial model. Second, event times must account for the left truncation derived from the fact that policyholders must reach the age of 65 before they may begin to be observed. Consequently, a delayed entry bias arises for those individuals entering the study after this age threshold. Third, and as the main challenge of our analysis, the association parameter between both processes is expected to be age‐dependent, with an unspecified association structure. This is well‐approximated through a flexible functional specification provided by penalized B‐splines. The parameter estimation of the joint model is derived under a Bayesian scheme.     

Change‐point models to estimate the limit of detection

In many biological and environmental studies, measured data is subject to a limit of detection. The limit of detection is generally defined as the lowest concentration of analyte that can be differentiated from a blank sample with some certainty. Data falling below the limit of detection is left censored, falling below a level that is easily quantified by a measuring device. A great deal of interest lies in estimating the limit of detection for a particular measurement device. In this paper, we propose a change‐point model to estimate the limit of detection by using data from an experiment with known analyte concentrations. Estimation of the limit of detection proceeds by a two‐stage maximum likelihood method. Extensions are considered that allow for censored measurements and data from multiple experiments. A simulation study is conducted demonstrating that in some settings the change‐point model provides less biased estimates of the limit of detection than conventional methods. The proposed method is then applied to data from an HIV pilot study. Copyright © 2013 John Wiley & Sons, Ltd.     

Segmental modeling of viral load changes for HIV longitudinal data with skewness and detection limits

Although it is a common practice to analyze complex HIV longitudinal data using nonlinear mixed‐effects or nonparametric mixed‐effects models in literature, the following issues may standout. (i) In clinical practice, the profile of each subject's viral response may follow a ‘broken‐stick’‐like trajectory, indicating multiple phases of decline and increase in response. Such multiple phases (change points) may be an important indicator to help quantify treatment effect and improve management of patient care. To estimate change points, nonlinear mixed‐effects or nonparametric mixed‐effects models become a challenge because of complicated structures of model formulations. (ii) The commonly assumed distribution for model random errors is normal, but this assumption may unrealistically obscure important features of subject variations. (iii) The response observations (viral load) may be subject to left censoring due to a limit of detection. Inferential procedures can be complicated dramatically when data with asymmetric (skewed) characteristics and left censoring are observed in conjunction with change points as unknown parameters into models. There is relatively little work concerning all these features simultaneously. This article proposes segmental mixed‐effects models with skew distributions for the response process (with left censoring) under a Bayesian framework. A real data example is used to illustrate the proposed methods. Copyright © 2012 John Wiley & Sons, Ltd.     

A joint model for interval‐censored functional decline trajectories under informative observation

Multi‐state models are useful for modelling disease progression where the state space of the process is used to represent the discrete disease status of subjects. Often, the disease process is only observed at clinical visits, and the schedule of these visits can depend on the disease status of patients. In such situations, the frequency and timing of observations may depend on transition times that are themselves unobserved in an interval‐censored setting. There is a potential for bias if we model a disease process with informative observation times as a non‐informative observation scheme with pre‐specified examination times. In this paper, we develop a joint model for the disease and observation processes to ensure valid inference because the follow‐up process may itself contain information about the disease process. The transitions for each subject are modelled using a Markov process, where bivariate subject‐specific random effects are used to link the disease and observation models. Inference is based on a Bayesian framework, and we apply our joint model to the analysis of a large study examining functional decline trajectories of palliative care patients. Copyright © 2015 John Wiley & Sons, Ltd.     

Functional linear models for zero‐inflated count data with application to modeling hospitalizations in patients on dialysis

We propose functional linear models for zero‐inflated count data with a focus on the functional hurdle and functional zero‐inflated Poisson (ZIP) models. Although the hurdle model assumes the counts come from a mixture of a degenerate distribution at zero and a zero‐truncated Poisson distribution, the ZIP model considers a mixture of a degenerate distribution at zero and a standard Poisson distribution. We extend the generalized functional linear model framework with a functional predictor and multiple cross‐sectional predictors to model counts generated by a mixture distribution. We propose an estimation procedure for functional hurdle and ZIP models, called penalized reconstruction, geared towards error‐prone and sparsely observed longitudinal functional predictors. The approach relies on dimension reduction and pooling of information across subjects involving basis expansions and penalized maximum likelihood techniques. The developed functional hurdle model is applied to modeling hospitalizations within the first 2 years from initiation of dialysis, with a high percentage of zeros, in the Comprehensive Dialysis Study participants. Hospitalization counts are modeled as a function of sparse longitudinal measurements of serum albumin concentrations, patient demographics, and comorbidities. Simulation studies are used to study finite sample properties of the proposed method and include comparisons with an adaptation of standard principal components regression. Copyright © 2014 John Wiley & Sons, Ltd.     

Joint analysis of stochastic processes with application to smoking patterns and insomnia

This article proposes a joint modeling framework for longitudinal insomnia measurements and a stochastic smoking cessation process in the presence of a latent permanent quitting state (i.e., ‘cure’). We use a generalized linear mixed‐effects model and a stochastic mixed‐effects model for the longitudinal measurements of insomnia symptom and for the smoking cessation process, respectively. We link these two models together via the latent random effects. We develop a Bayesian framework and Markov Chain Monte Carlo algorithm to obtain the parameter estimates. We formulate and compute the likelihood functions involving time‐dependent covariates. We explore the within‐subject correlation between insomnia and smoking processes. We apply the proposed methodology to simulation studies and the motivating dataset, that is, the Alpha‐Tocopherol, Beta‐Carotene Lung Cancer Prevention study, a large longitudinal cohort study of smokers from Finland. Copyright © 2013 John Wiley & Sons, Ltd.     

Retracted: Bayesian inference on mixed‐effects location scale models with skew‐t distribution and mismeasured covariates for longitudinal data

In AIDS studies, heterogeneous between and within subject variations are often observed on longitudinal endpoints. To accommodate heteroscedasticity in the longitudinal data, statistical methods have been developed to model the mean and variance jointly. Most of these methods assume (conditional) normal distributions for random errors, which is not realistic in practice. In this article, we propose a Bayesian mixed‐effects location scale model with skew‐t distribution and mismeasured covariates for heterogeneous longitudinal data with skewness. The proposed model captures the between‐subject and within‐subject (WS) heterogeneity by modeling the between‐subject and WS variations with covariates as well as a random effect at subject level in the WS variance. Further, the proposed model also takes into account the covariate measurement errors, and commonly assumed normal distributions for model errors are substituted by skew‐t distribution to account for skewness. Parameter estimation is carried out in a Bayesian framework. The proposed method is illustrated with a Multicenter AIDS Cohort Study. Simulation studies are performed to assess the performance of the proposed method. Copyright © 2017 John Wiley & Sons, Ltd.     

Dynamic predictions using flexible joint models of longitudinal and time‐to‐event data

Joint models for longitudinal and time‐to‐event data are particularly relevant to many clinical studies where longitudinal biomarkers could be highly associated with a time‐to‐event outcome. A cutting‐edge research direction in this area is dynamic predictions of patient prognosis (e.g., survival probabilities) given all available biomarker information, recently boosted by the stratified/personalized medicine initiative. As these dynamic predictions are individualized, flexible models are desirable in order to appropriately characterize each individual longitudinal trajectory. In this paper, we propose a new joint model using individual‐level penalized splines (P‐splines) to flexibly characterize the coevolution of the longitudinal and time‐to‐event processes. An important feature of our approach is that dynamic predictions of the survival probabilities are straightforward as the posterior distribution of the random P‐spline coefficients given the observed data is a multivariate skew‐normal distribution. The proposed methods are illustrated with data from the HIV Epidemiology Research Study. Our simulation results demonstrate that our model has better dynamic prediction performance than other existing approaches. © 2017 The Authors. Statistics in Medicine Published by John Wiley & Sons Ltd.     

Frailty models for pneumonia to death with a left‐censored covariate

Frailty models are multiplicative hazard models for studying association between survival time and important clinical covariates. When some values of a clinical covariate are unobserved but known to be below a threshold called the limit of detection (LOD), naive approaches ignoring this problem, such as replacing the undetected value by the LOD or half of the LOD, often produce biased parameter estimate with larger mean squared error of the estimate. To address the LOD problem in a frailty model, we propose a flexible smooth nonparametric density estimator along with Simpson's numerical integration technique. This is an extension of an existing method in the likelihood framework for the estimation and inference of the model parameters. The proposed new method shows the estimators are asymptotically unbiased and gives smaller mean squared error of the estimates. Compared with the existing method, the proposed new method does not require distributional assumptions for the underlying covariates. Simulation studies were conducted to evaluate the performance of the new method in realistic scenarios. We illustrate the use of the proposed method with a data set from Genetic and Inflammatory Markers of Sepsis study in which interlekuin‐10 was subject to LOD. Copyright © 2015 John Wiley & Sons, Ltd.     

Statistical models for longitudinal zero‐inflated count data with applications to the substance abuse field

This study fills in the current knowledge gaps in statistical analysis of longitudinal zero‐inflated count data by providing a comprehensive review and comparison of the hurdle and zero‐inflated Poisson models in terms of the conceptual framework, computational advantage, and performance under different real data situations. The design of simulations represents the special features of a well‐known longitudinal study of alcoholism so that the results can be generalizable to the substance abuse field. When the hurdle model is more natural under the conceptual framework of the data, the zero‐inflated Poisson model tends to produce inaccurate estimates. Model performance improves with larger sample sizes, lower proportions of missing data, and lower correlations between covariates. The simulation also shows that the computational strength of the hurdle model disappears when random effects are included. Copyright © 2012 John Wiley & Sons, Ltd.     

Semiparametric Bayesian inference on skew–normal joint modeling of multivariate longitudinal and survival data

We propose a semiparametric multivariate skew–normal joint model for multivariate longitudinal and multivariate survival data. One main feature of the posited model is that we relax the commonly used normality assumption for random effects and within‐subject error by using a centered Dirichlet process prior to specify the random effects distribution and using a multivariate skew–normal distribution to specify the within‐subject error distribution and model trajectory functions of longitudinal responses semiparametrically. A Bayesian approach is proposed to simultaneously obtain Bayesian estimates of unknown parameters, random effects and nonparametric functions by combining the Gibbs sampler and the Metropolis–Hastings algorithm. Particularly, a Bayesian local influence approach is developed to assess the effect of minor perturbations to within‐subject measurement error and random effects. Several simulation studies and an example are presented to illustrate the proposed methodologies. Copyright © 2014 John Wiley & Sons, Ltd.     

A joint frailty model to estimate the recurrence process and the disease‐specific mortality process without needing the cause of death

In chronic diseases, such as cancer, recurrent events (such as relapses) are commonly observed; these could be interrupted by death. With such data, a joint analysis of recurrence and mortality processes is usually conducted with a frailty parameter shared by both processes. We examined a joint modeling of these processes considering death under two aspects: ‘death due to the disease under study' and ‘death due to other causes', which enables estimating the disease‐specific mortality hazard. The excess hazard model was used to overcome the difficulties in determining the causes of deaths (unavailability or unreliability); this model allows estimating the disease‐specific mortality hazard without needing the cause of death but using the mortality hazards observed in the general population. We propose an approach to model jointly recurrence and disease‐specific mortality processes within a parametric framework. A correlation between the two processes is taken into account through a shared frailty parameter. This approach allows estimating unbiased covariate effects on the hazards of recurrence and disease‐specific mortality. The performance of the approach was evaluated by simulations with different scenarios. The method is illustrated by an analysis of a population‐based dataset on colon cancer with observations of colon cancer recurrences and deaths. The benefits of the new approach are highlighted by comparison with the ‘classical' joint model of recurrence and overall mortality. Moreover, we assessed the goodness of fit of the proposed model. Comparisons between the conditional hazard and the marginal hazard of the disease‐specific mortality are shown, and differences in interpretation are discussed. Copyright © 2014 John Wiley & Sons, Ltd.