Semiparametric Joint Model Research Articles

Semiparametric multivariate joint model for skewed-longitudinal and survival data: A Bayesian approach.

Joint models and statistical inference for longitudinal and survival data have been an active area of statistical research and have mostly coupled a longitudinal biomarker-based mixed-effects model with normal distribution and an event time-based survival model. In practice, however, the following issues may standout: (i) Normality of model error in longitudinal models is a routine assumption, but it may be unrealistically violating data features of subject variations. (ii) Data collected are often featured by the mixed types of multiple longitudinal outcomes which are significantly correlated, ignoring their correlation may lead to biased estimation. Additionally, a parametric model specification may be inflexible to capture the complicated patterns of longitudinal data. (iii) Missing observations in the longitudinal data are often encountered; the missing measures are likely to be informative (nonignorable) and ignoring this phenomenon may result in inaccurate inference. Multilevel item response theory (MLIRT) models have been increasingly used to analyze the multiple longitudinal data of mixed types (ie, continuous and categorical) in clinical studies. In this article, we develop an MLIRT-based semiparametric joint model with skew-t distribution that consists of an extended MLIRT model for the mixed types of multiple longitudinal data and a Cox proportional hazards model, linked through random-effects. A Bayesian approach is employed for joint modeling. Simulation studies are conducted to assess performance of the proposed models and method. A real example from primary biliary cirrhosis clinical study is analyzed to estimate parameters in the joint model and also evaluate sensitivity of parameter estimates for various plausible nonignorable missing data mechanisms.

Bayesian semiparametric joint model of multivariate longitudinal and survival data with dependent censoring.

We consider a novel class of semiparametric joint models for multivariate longitudinal and survival data with dependent censoring. In these models, unknown-fashion cumulative baseline hazard functions are fitted by a novel class of penalized-splines (P-splines) with linear constraints. The dependence between the failure time of interest and censoring time is accommodated by a normal transformation model, where both nonparametric marginal survival function and censoring function are transformed to standard normal random variables with bivariate normal joint distribution. Based on a hybrid algorithm together with the Metropolis-Hastings algorithm within the Gibbs sampler, we propose a feasible Bayesian method to simultaneously estimate unknown parameters of interest, and to fit baseline survival and censoring functions. Intensive simulation studies are conducted to assess the performance of the proposed method. The use of the proposed method is also illustrated in the analysis of a data set from the International Breast Cancer Study Group.

Copula modeling for the estimation of measures of marker classification and predictiveness performance with survival outcomes.

The discriminative and predictive power of a continuous-valued marker for survival outcomes can be summarized using the receiver operating characteristic and predictiveness curves, respectively. In this paper, fully parametric and semi-parametric copula-based constructions of the joint model of the marker and the survival time are developed for characterizing, plotting, and analyzing both curves along with other underlying performance measures. The formulations require a copula function, a parametric specification for the margin of the marker, and either a parametric distribution or a non-parametric estimator for the margin of the time to event, to respectively characterize the fully parametric and semi-parametric joint models. Estimation is carried out using maximum likelihood and a two-stage procedure for the parametric and semi-parametric models, respectively. Resampling-based methods are used for computing standard errors and confidence bounds for the various parameters, curves, and associated measures. Graphical inspection of residuals from each conditional distribution is employed as a guide for choosing a copula from a set of candidates. The performance of the estimators of various classification and predictiveness measures is assessed in simulation studies, assuming different copula and censoring scenarios. The methods are illustrated with the analysis of two markers using the familiar primary biliary cirrhosis data set.

Causal inference with longitudinal data subject to irregular assessment times.

Data collected in the context of usual care present a rich source of longitudinal data for research, but often require analyses that simultaneously enable causal inferences from observational data while handling irregular and informative assessment times. An inverse-weighting approach to this was recently proposed, and handles the case where the assessment times are at random (ie, conditionally independent of the outcome process given the observed history). In this paper, we extend the inverse-weighting approach to handle a special case of assessment not at random, where assessment and outcome processes are conditionally independent given past observed covariates and random effects. We use multiple outputation to accomplish the same purpose as inverse-weighting, and apply it to the Liang semi-parametric joint model. Moreover, we develop an alternative joint model that does not require covariates for the outcome model to be known at times where there is no assessment of the outcome. We examine the performance of these methods through simulation and illustrate them through a study of the causal effect of wheezing on time spent playing outdoors among children aged 2-9 years and enrolled in the TargetKids! study.

Bayesian semiparametric joint modeling of a count outcome and inconveniently timed longitudinal predictors.

The Women's Health Initiative (WHI) Life and Longevity After Cancer (LILAC) study is an excellent resource for studying the quality of life following breast cancer treatment. At study entry, women were asked about new symptoms that appeared following their initial cancer treatment. In this article, we were interested in using regression modeling to estimate associations of clinical and lifestyle factors at cancer diagnosis (independent variables) with the number of new symptoms (dependent variable). Although clinical and lifestyle data were collected longitudinally, few measurements were obtained at diagnosis or at a consistent timepoint prior to diagnosis, which complicates the analysis. Furthermore, parametric count models, such as the Poisson and negative binomial, do not fit the symptom data well. Thus, motivated by the issues encountered in LILAC, we propose two Bayesian joint models for longitudinal data and a count outcome. Our two models differ according to the assumption on the outcome distribution: one uses a negative binomial (NB) distribution and the other a nonparametric rounded mixture of Gaussians (RMG). The mean of each count distribution is dependent on imputed values of continuous, binary, and ordinal variables at a time point of interest (e.g. diagnosis). To facilitate imputation, longitudinal variables are modeled jointly using a linear mixed model for a latent underlying normal random variable, and a Dirichlet process prior is assigned to the random subject-specific effects to relax distribution assumptions. In simulation studies, the RMG joint model exhibited superior power and predictive accuracy over the NB model when the data were not NB. The RMG joint model also outperformed an RMG model containing predictors imputed using the last value carried forward, which generated estimates that were biased toward the null. We used our models to examine the relationship between sleep health at diagnosis and the number of new symptoms following breast cancer treatment in LILAC.

A Semiparametric Bayesian Joint Modelling of Skewed Longitudinal and Competing Risks Failure Time Data: With Application to Chronic Kidney Disease

In clinical and epidemiological studies, when the time-to-event(s) and the longitudinal outcomes are associated, modelling them separately may give biased estimates. A joint modelling approach is required to obtain unbiased results and to evaluate their association. In the joint model, a subject may be exposed to more than one type of failure event (competing risks). Considering the competing event as an independent censoring of the time-to-event process may underestimate the true survival probability and give biased results. Within the joint model, longitudinal outcomes may have nonlinear (irregular) trajectories over time and exhibit skewness with heavy tails. Accordingly, fully parametric mixed-effect models may not be flexible enough to model this type of complex longitudinal data. In addition, assuming a Gaussian distribution for model errors may be too restrictive to adequately represent within-individual variations and may lack robustness against deviation from distributional assumptions. To simultaneously overcome these issues, in this paper, we presented semiparametric joint models for competing risks failure time and skewed-longitudinal data by using a smoothing spline approach and a multivariate skew-t distribution. We also considered different parameterization approaches in the formulation of joint models and used a Bayesian approach to make the statistical inference. We illustrated the proposed methods by analyzing real data on a chronic kidney disease. To evaluate the performance of the methods, we also carried out simulation studies. The results of both the application and simulation studies revealed that the joint modelling approach proposed in this study performed well when the semiparametric, random-effects parameterization, and skew-t distribution specifications were taken into account.

A semiparametric joint model for cluster size and subunit-specific interval-censored outcomes.

Clustered data frequently arise in biomedical studies, where observations, or subunits, measured within a cluster are associated. The cluster size is said to be informative, if the outcome variable is associated with the number of subunits in a cluster. In most existing work, the informative cluster size issue is handled by marginal approaches based on within-cluster resampling, or cluster-weighted generalized estimating equations. Although these approaches yield consistent estimation of the marginal models, they do not allow estimation of within-cluster associations and are generally inefficient. In this paper, we propose a semiparametric joint model for clustered interval-censored event time data with informative cluster size. We use a random effect to account for the association among event times of the same cluster as well as the association between event times and the cluster size. For estimation, we propose a sieve maximum likelihood approach and devise a computationally-efficient expectation-maximization algorithm for implementation. The estimators are shown to be strongly consistent, with the Euclidean components being asymptotically normal and achieving semiparametric efficiency. Extensive simulation studies are conducted to evaluate the finite-sample performance, efficiency and robustness of the proposed method. We also illustrate our method via application to a motivating periodontal diseasedataset.

Efficient Algorithms and Implementation of a Semiparametric Joint Model for Longitudinal and Competing Risk Data: With Applications to Massive Biobank Data.

Semiparametric joint models of longitudinal and competing risk data are computationally costly, and their current implementations do not scale well to massive biobank data. This paper identifies and addresses some key computational barriers in a semiparametric joint model for longitudinal and competing risk survival data. By developing and implementing customized linear scan algorithms, we reduce the computational complexities from O(n2) or O(n3) to O(n) in various steps including numerical integration, risk set calculation, and standard error estimation, where n is the number of subjects. Using both simulated and real-world biobank data, we demonstrate that these linear scan algorithms can speed up the existing methods by a factor of up to hundreds of thousands when n > 104, often reducing the runtime from days to minutes. We have developed an R package, FastJM, based on the proposed algorithms for joint modeling of longitudinal and competing risk time-to-event data and made it publicly available on the Comprehensive R Archive Network (CRAN).

Bayesian semiparametric joint modeling of longitudinal explanatory variables of mixed types and a binary outcome.

Many prospective biomedical studies collect longitudinal clinical and lifestyle data that are both continuous and discrete. In some studies, there is interest in the association between a binary outcome and the values of these longitudinal measurements at a specific time point. A common problem in these studies is inconsistency in timing of measurements and missing follow-ups which can lead to few measurements at the time of interest. Some methods have been developed to address this problem, but are only applicable to continuous measurements. To address this limitation, we propose a new class of joint models for a binary outcome and longitudinal explanatory variables of mixed types. The longitudinal model uses a latent normal random variable construction with regression splines to model time-dependent trends in mean with a Dirichlet Process prior assigned to random effects to relax distribution assumptions. We also standardize timing of the explanatory variables by relating the binary outcome to imputed longitudinal values at a set time point. The proposed model is evaluated through simulation studies and applied to data from a cancer survivor study of participants in the Women's Health Initiative.

Clustered longitudinal data subject to irregular observation.

Data collected longitudinally as part of usual health care is becoming increasingly available for research, and is often available across several centres. Because the frequency of follow-up is typically determined by the patient's health, the timing of measurements may be related to the outcome of interest. Failure to account for the informative nature of the observation process can result in biased inferences. While methods for accounting for the association between observation frequency and outcome are available, they do not currently account for clustering within centres. We formulate a semi-parametric joint model to include random effects for centres as well as subjects. We also show how inverse-intensity weighted GEEs can be adapted to account for clustering, comparing stratification, frailty models, and covariate adjustment to account for clustering in the observation process. The finite-sample performance of the proposed methods is evaluated through simulation and the methods illustrated using a study of the relationship between outdoor play and air quality in children aged 2-9 living in the Greater Toronto Area.

Semiparametric Joint Models to Assess the Effect of Timing of the Indication for Treatment on Survival with Application to Cancer Screening

Semiparametric Joint Models to Assess the Effect of Timing of the Indication for Treatment on Survival with Application to Cancer Screening

Covariate-Adjusted Regression for Distorted Longitudinal Data With Informative Observation Times

ABSTRACTIn many longitudinal studies, repeated response and predictors are not directly observed, but can be treated as distorted by unknown functions of a common confounding covariate. Moreover, longitudinal data involve an observation process which may be informative with a longitudinal response process in practice. To deal with such complex data, we propose a class of flexible semiparametric covariate-adjusted joint models. The new models not only allow for the longitudinal response to be correlated with observation times through latent variables and completely unspecified link functions, but they also characterize distorted longitudinal response and predictors by unknown multiplicative factors depending on time and a confounding covariate. For estimation of regression parameters in the proposed models, we develop a novel covariate-adjusted estimating equation approach which does not rely on forms of link functions and distributions of frailties. The asymptotic properties of resulting parameter estimators are established and examined by simulation studies. A longitudinal data example containing calcium absorption and intake measurements is provided for illustration. Supplementary materials for this article are available online.

A joint model for recurrent events and a semi-competing risk in the presence of multi-level clustering.

Clinical trial designs often include multiple levels of clustering in which patients are nested within clinical sites and recurrent outcomes are nested within patients who may also experience a semi-competing risk. Traditional survival methods that analyze these processes separately may lead to erroneous inferences as they ignore possible dependencies. To account for the association between recurrent events and a semi-competing risk in the presence of two levels of clustering, we developed a semi-parametric joint model. The Gaussian quadrature with a piecewise constant baseline hazard was used to estimate the unspecified baseline hazards and the likelihood. Simulations showed that the proposed joint model has good statistical properties (i.e. <5% bias and 95% coverage) compared to the shared frailty and joint frailty models when informative censoring and multiple levels of clustering were present. The proposed method was applied to data from an AIDS clinical trial to investigate the impact of antiretroviral treatment on recurrent AIDS-defining events in the presence of a semi-competing risk of death and multi-level clustering and showed a significant dependency between AIDS-defining events and death at the patient level but not at the clinic level.

Flexible semiparametric joint modeling: an application to estimate individual lung function decline and risk of pulmonary exacerbations in cystic fibrosis

BackgroundEpidemiologic surveillance of lung function is key to clinical care of individuals with cystic fibrosis, but lung function decline is nonlinear and often impacted by acute respiratory events known as pulmonary exacerbations. Statistical models are needed to simultaneously estimate lung function decline while providing risk estimates for the onset of pulmonary exacerbations, in order to identify relevant predictors of declining lung function and understand how these associations could be used to predict the onset of pulmonary exacerbations.MethodsUsing longitudinal lung function (FEV1) measurements and time-to-event data on pulmonary exacerbations from individuals in the United States Cystic Fibrosis Registry, we implemented a flexible semiparametric joint model consisting of a mixed-effects submodel with regression splines to fit repeated FEV1 measurements and a time-to-event submodel for possibly censored data on pulmonary exacerbations. We contrasted this approach with methods currently used in epidemiological studies and highlight clinical implications.ResultsThe semiparametric joint model had the best fit of all models examined based on deviance information criterion. Higher starting FEV1 implied more rapid lung function decline in both separate and joint models; however, individualized risk estimates for pulmonary exacerbation differed depending upon model type. Based on shared parameter estimates from the joint model, which accounts for the nonlinear FEV1 trajectory, patients with more positive rates of change were less likely to experience a pulmonary exacerbation (HR per one standard deviation increase in FEV1 rate of change = 0.566, 95% CI 0.516–0.619), and having higher absolute FEV1 also corresponded to lower risk of having a pulmonary exacerbation (HR per one standard deviation increase in FEV1 = 0.856, 95% CI 0.781–0.937). At the population level, both submodels indicated significant effects of birth cohort, socioeconomic status and respiratory infections on FEV1 decline, as well as significant effects of gender, socioeconomic status and birth cohort on pulmonary exacerbation risk.ConclusionsThrough a flexible joint-modeling approach, we provide a means to simultaneously estimate lung function trajectories and the risk of pulmonary exacerbations for individual patients; we demonstrate how this approach offers additional insights into the clinical course of cystic fibrosis that were not possible using conventional approaches.

A joint model of cancer incidence, metastasis, and mortality.

Many diseases, especially cancer, are not static, but rather can be summarized by a series of events or stages (e.g. diagnosis, remission, recurrence, metastasis, death). Most available methods to analyze multi-stage data ignore intermediate events and focus on the terminal event or consider (time to) multiple events as independent. Competing-risk or semi-competing-risk models are often deficient in describing the complex relationship between disease progression events which are driven by a shared progression stochastic process. A multi-stage model can only examine two stages at a time and thus fails to capture the effect of one stage on the time spent between other stages. Moreover, most models do not account for latent stages. We propose a semi-parametric joint model of diagnosis, latent metastasis, and cancer death and use nonparametric maximum likelihood to estimate covariate effects on the risks of intermediate events and death and the dependence between them. We illustrate the model with Monte Carlo simulations and analysis of real data on prostate cancer from the SEER database.

A semiparametric joint model for terminal trend of quality of life and survival in palliative care research.

Palliative medicine is an interdisciplinary specialty focusing on improving quality of life (QOL) for patients with serious illness and their families. Palliative care programs are available or under development at over 80% of large US hospitals (300+ beds). Palliative care clinical trials present unique analytic challenges relative to evaluating the palliative care treatment efficacy which is to improve patients' diminishing QOL as disease progresses towards end of life (EOL). A unique feature of palliative care clinical trials is that patients will experience decreasing QOL during the trial despite potentially beneficial treatment. Often longitudinal QOL and survival data are highly correlated which, in the face of censoring, makes it challenging to properly analyze and interpret terminal QOL trend. To address these issues, we propose a novel semiparametric statistical approach to jointly model the terminal trend of QOL and survival data. There are two sub-models in our approach: a semiparametric mixed effects model for longitudinal QOL and a Cox model for survival. We use regression splines method to estimate the nonparametric curves and AIC to select knots. We assess the model performance through simulation to establish a novel modeling approach that could be used in future palliative care research trials. Application of our approach in a recently completed palliative care clinical trial is also presented.

Simultaneous inference for semiparametric mixed-effects joint models with skew distribution and covariate measurement error for longitudinal competing risks data analysis

ABSTRACTSemiparametric mixed-effects joint models are flexible for modeling complex longitudinal–competing risks data. Skew distributions are commonly observed for this type of data. Covariates in the joint models are usually measured with substantial errors. We propose a Bayesian method for semiparametric mixed-effects joint models with covariate measurement errors and skew distribution. The methods are illustrated with AIDS clinical data. Simulation results are conducted to validate the proposed methods.

Bayesian semiparametric mixed-effects joint models for analysis of longitudinal-competing risks data with skew distribution

Bayesian semiparametric mixed-effects joint models for analysis of longitudinal-competing risks data with skew distribution

Bayesian variable selection and estimation in semiparametric joint models of multivariate longitudinal and survival data.

This paper presents a novel semiparametric joint model for multivariate longitudinal and survival data (SJMLS) by relaxing the normality assumption of the longitudinal outcomes, leaving the baseline hazard functions unspecified and allowing the history of the longitudinal response having an effect on the risk of dropout. Using Bayesian penalized splines to approximate the unspecified baseline hazard function and combining the Gibbs sampler and the Metropolis-Hastings algorithm, we propose a Bayesian Lasso (BLasso) method to simultaneously estimate unknown parameters and select important covariates in SJMLS. Simulation studies are conducted to investigate the finite sample performance of the proposed techniques. An example from the International Breast Cancer Study Group (IBCSG) is used to illustrate the proposed methodologies.

Joint partially linear model for longitudinal data with informative drop-outs.

In biomedical research, a steep rise or decline in longitudinal biomarkers may indicate latent disease progression, which may subsequently cause patients to drop out of the study. Ignoring the informative drop-out can cause bias in estimation of the longitudinal model. In such cases, a full parametric specification may be insufficient to capture the complicated pattern of the longitudinal biomarkers. For these types of longitudinal data with the issue of informative drop-outs, we develop a joint partially linear model, with an aim to find the trajectory of the longitudinal biomarker. Specifically, an arbitrary function of time along with linear fixed and random covariate effects is proposed in the model for the biomarker, while a flexible semiparametric transformation model is used to describe the drop-out mechanism. Advantages of this semiparametric joint modeling approach are the following: 1) it provides an easier interpretation, compared to standard nonparametric regression models, and 2) it is a natural way to control for common (observable and unobservable) prognostic factors that may affect both the longitudinal trajectory and the drop-out process. We describe a sieve maximum likelihood estimation procedure using the EM algorithm, where the Akaike information criterion (AIC) and Bayesian information criterion (BIC) are considered to select the number of knots. We show that the proposed estimators achieve desirable asymptotic properties through empirical process theory. The proposed methods are evaluated by simulation studies and applied to prostate cancer data.