- Research article
- Open Access
- Open Peer Review
Joint models for longitudinal and time-to-event data: a review of reporting quality with a view to meta-analysis
- Maria Sudell1Email authorView ORCID ID profile,
- Ruwanthi Kolamunnage-Dona†1 and
- Catrin Tudur-Smith†1
© The Author(s). 2016
- Received: 9 September 2016
- Accepted: 23 November 2016
- Published: 5 December 2016
The Correction to this article has been published in BMC Medical Research Methodology 2018 18:33
Joint models for longitudinal and time-to-event data are commonly used to simultaneously analyse correlated data in single study cases. Synthesis of evidence from multiple studies using meta-analysis is a natural next step but its feasibility depends heavily on the standard of reporting of joint models in the medical literature. During this review we aim to assess the current standard of reporting of joint models applied in the literature, and to determine whether current reporting standards would allow or hinder future aggregate data meta-analyses of model results.
We undertook a literature review of non-methodological studies that involved joint modelling of longitudinal and time-to-event medical data. Study characteristics were extracted and an assessment of whether separate meta-analyses for longitudinal, time-to-event and association parameters were possible was made.
The 65 studies identified used a wide range of joint modelling methods in a selection of software. Identified studies concerned a variety of disease areas. The majority of studies reported adequate information to conduct a meta-analysis (67.7% for longitudinal parameter aggregate data meta-analysis, 69.2% for time-to-event parameter aggregate data meta-analysis, 76.9% for association parameter aggregate data meta-analysis). In some cases model structure was difficult to ascertain from the published reports.
Whilst extraction of sufficient information to permit meta-analyses was possible in a majority of cases, the standard of reporting of joint models should be maintained and improved. Recommendations for future practice include clear statement of model structure, of values of estimated parameters, of software used and of statistical methods applied.
- Joint model
- Reporting standards
Joint modelling of longitudinal and time-to-event data is an area of increasing research [1–3], which allows the simultaneous modelling of a longitudinal (repeatedly measured over time) outcome such as weekly biomarker measurements, and a time-to-event (survival) outcome such as time to death. The model consists of two sub-models; a longitudinal sub-model (such as a linear mixed effects model) and a time-to-event sub-model (such as a cox proportional hazards models) which are linked using an association structure that quantifies the relationship between the outcomes of interest.
Within a single study, joint models have the potential to reduce parameter estimate bias, account for dropout in longitudinal studies and enable the inclusion of longitudinal covariates measured with error in time-to-event models [1, 4]. These qualities often make joint models preferable to separate longitudinal or time-to-event analyses. Joint models have been applied in the literature to investigate links between biomarkers and certain disease events (e.g. in cancer studies), and to account for informative study dropout.
Glass 1976  defined meta-analysis (MA) as the statistical analysis or pooling of results from separate studies. Such analyses can increase power and precision compared to original studies, or answer questions additional to those originally posed . MA can be performed on the original Individual Participant Data (IPD), or on the study level results (published in the literature or obtained from authors) termed Aggregate Data (AD). Overviews of MA methodology can be found in Whitehead 2002  and in the Cochrane handbook .
Whilst the benefits of joint modelling methods for individual studies are well established, little attention has been given towards the potential value of pooling estimates across similar studies in an aggregate data meta-analysis (AD-MA) of joint models. However before an AD-MA using published data can be undertaken, relevant studies must be identified, and the necessary information must be extracted. We aim to investigate the reporting of joint longitudinal and time-to-event models applied to real medical data in the literature to establish whether current reporting practices would allow sufficient data to be extracted to undertake AD-MA.
Identification of papers
We performed our systematic review in accordance to the guidelines of the Preferred Reporting of Items for Systematic Reviews and Meta-Analyses (PRISMA) . We searched the Medline, Pubmed and Scopus datasets for studies using joint models for longitudinal and time-to-event data to analyse medical data (search strategies available in the Additional file 1).
Papers mentioning joint models for longitudinal (or repeated measures over time) data and time-to-event (or survival, event time or event history) data were identified. Duplicates were identified and removed. Abstracts and keywords were then examined, and irrelevant papers were removed. Examples of disregarded papers include papers modelling body joints, papers discussing joint models as a future extension or alternative to methods used, or papers using two stage approaches rather than simultaneous estimation of the longitudinal and time-to-event sub-models. Papers not relating to medical or biostatistical datasets were discarded (e.g. data from plant or animal subjects except from modelling of human diseases input into animal hosts). Additionally papers involving repeated measures over space rather than time were discarded (e.g. repeated measures across tumour sites). If study relevancy was unclear from the abstract, the full text was obtained and viewed after which the study was included or discarded.
Any retained papers were sorted into an applied and a methodological group. Some methodology papers presented results from application to example datasets. These were considered reanalyses of data or demonstrations of methods rather than primary analyses to influence future practice. Also, methodological papers might be expected to better report results as their authors are experts in the area. The aim of this review was to assess how well joint models are reported in the general medical literature, so we focussed on the applied group only.
A blank data extraction form is presented in Additional file 2. During the investigation we refer to references identified as applying joint models to relevant datasets as studies. Other publications (e.g. those cited by studies) are referred to as papers. Information recorded from identified studies included publication year, author, journal, joint model type, sharing structure between the longitudinal and time-to-event sub-models, types of sub-models, Bayesian or frequentist methods, and software used. Disease area was recorded (with respect to the type of longitudinal and time-to-event data, for example studies modelling biomarkers in heart disease patients after a transplant operation were classed as transplant data).
The sources of the methods used were recorded. Specifically if the study developed methods specific to their dataset, “own methods developed” was recorded. If the study referenced specific papers as the source of the methods they used the papers referenced were recorded.
Availability of information required for a MA was also recorded (including the number of participants, significance level, and presence longitudinal, time-to-event and association parameters along with their precision estimates). The significance level used was identified through direct statement in the text, or specified confidence interval sizes on tables, graphs, or footnotes relating to the joint models fitted. For a MA to be considered possible, the number of participants and model coefficients had to be reported, with either a standard error, or a confidence interval with accompanying significance level.
We assume for an AD-MA of joint longitudinal and time-to-event models that a separate MA would be conducted for each of the longitudinal, the time-to-event and the association parameters from the identified studies. Consequently when identified studies were assessed for sufficient information to conduct a MA, they were assessed separately for longitudinal, for time-to-event and for association information. Ideally all three separate groups of MA would be conducted, however if insufficient information was reported only a subset of these MA might be feasible.
The reason for joint models use (see Henderson et al [9) may influence what information is presented in the study report. If joint models were used to account for informative dropout, the study might not report the time-to-event parameter information (although if the time-to-event endpoint is clinically defined then time-to-event estimates should be reported). Similarly if interest was to include a longitudinal variable measured with error as a time-variable covariate in a time-to-event model, the study may not clearly report longitudinal sub-model parameters. To investigate whether the reason for joint model use affected the proportion of possible MA for each of the longitudinal, time-to-event and association parameter components, we investigated the proportions of possible MA for studies using joint models to account for dropout, or to account for error in a time-varying covariate.
The aim of this review was not to perform any MA, solely to assess if MA were undertaken, what proportion of the identified studies could contribute.
Characteristics of identified studies
Year of publication
Full text or abstract
Full articles were obtained for 63 studies (96.9%) [9–21, 23–54, 56–73], with abstracts available for 2 studies (3.1%) [22, 55]. Of the identified studies, some individuals were listed as authors on multiple studies, suggesting that the group of individuals applying joint modelling methods may be limited.
Characteristics of identified studies
Full text or abstract available
Cancer related data
Patient status after transplants
Disability in the elderly
Heart related data
Statistics in Medicine
Journal of the Royal Statistical Society. Series C: Applied Statistics
Quality of Life Research
Journal of the American Geriatrics Society
Journal of the American Statistical Association
Journals of Gerontology - Series B Psychological Sciences and Social Sciences
Statistical Methods in Medical Research
Other (only one study per journal)
Reason for joint modelling use*
To investigate the link between longitudinal and time-to-event outcomes
To account for dropout
To include longitudinally measured variable in time-to-event model
To increase efficiency
To reduce bias
Easier to interpret
To use of all available data
The studies identified were published in a range of journals, with 8 journals occurring more than once (Table 1), indicating that there may not currently a preferred journal to present joint modelling studies in.
Reason for use of joint model
The reasons given for using joint models are summarised in Table 1 (with some studies providing multiple reasons giving a total greater than 65). The two most common reasons for joint model use were to investigate the link between the outcomes, or to account for dropout. Only 4 studies stated that they used a joint model to include a time varying covariate in the time-to-event sub-model.
Source of methods used
Methods used in identified studies
Source of joint modelling methods used
Own methods developed
Guo-Carlin 2004 
Rizopoulos 2010 JM R package 
Henderson et al 2000 
Tsiatis and Davidian 2004 
Rizopoulos 2012 
Wulfsohn and Tsiatis 1997 
Diggle et al 2008 
Crowther et al 2013 
Proust-Lima et al 2009 
Rizopoulos 2011 
Fixed and Random Effects
Current Value of Fixed and Random Effects
Current Slope (first derivative) of Fixed and Random Effects
Current Value of Fixed and Random Effects and Current Slope (first derivative) of Fixed and Random Effects
Fixed and random effects without covariates
Random Effects only
Random Effects with covariates
Random Effects without covariates
Random Effects unclear with or without covariates
Specialist sharing structure
Of the 65 included studies, 45 (69.2%) took a frequentist approach, 17 (26.2%) took a Bayesian approach, 1 (1.5%) took both (in separate model fits) and in 2 (3.1%) studies it was unclear (these were the two abstracts). The larger proportion of frequentist approaches may be attributable to the larger number of papers and based on frequentist methods. Additionally the main joint modelling textbook  deals with frequentist methods.
There were 21 unique model types recorded for the longitudinal sub-model (plus 1 study with unclear type). Linear mixed effects models were most common (35 studies, 53.8%), followed by mixed effect models with splines (6 studies, 9.2%) or mixed models with unspecified structure (5 studies, 7.7%). Other methods used included different mixed models dependent on latent class, non-linear models with or without splines, and models with change points.
The methods used for the time-to-event sub-model varied widely (and were unclear in 4 studies). The Cox proportional hazards (PH) model was most common (8 studies, 12.3%). Other methods included models with parametric baselines, such as a Weibull PH model (5 studies, 7.7%), a PH model with piecewise constant baseline (4 studies, 6.2%), or a spline modelled baseline hazard (2 studies, 3.1%). Parametric models included the Weibull (5 studies, 7.7%) and the exponential (1 study, 1.5%), and 1 study (1.5%) examined both Weibull and exponential models.
Sharing structure between longitudinal and time-to-event sub-models
The structures used to link the longitudinal and time-to-event sub-models are listed in Table 2. Some studies fitted multiple joint models, with varying sharing structures, allowing a total of more than 65 recorded sharing structures.
Any sharing structure (also termed association structure) involving a function of both fixed and random effects is designated “Fixed and Random effects” whereas those involving random effects but no fixed effects are termed “Random Effects only”. Fixed and random effects sharing structures (33 studies, 50.8%) model effects of aspects of the overall longitudinal outcome value on the time-to-event outcome. The random effects only sharing structures (27 studies, 41.5%) model the effect of individual deviation from the population mean longitudinal outcome on the time-to-event outcome. A description of fixed and random effects sharing structures is given in Rizopoulos 2012 , whilst Henderson et al  discuss random effects only sharing structures. Additionally Rizopoulos and Ghosh  and Gould  discuss a range of association structures.
The fixed and random effect sharing structures can be subdivided further. Current value refers to models inserting the current longitudinal trajectory value into the time-to-event sub-models, and is used when the current overall value of the longitudinal trajectory affects the risk of an event. The current slope or first derivative of the population trajectory could also be inserted into the time-to-event sub-model, in conjunction with the current value or alone, and is used to model the effect of rate of change of the longitudinal variable on the risk of an event. Another fixed and random effects sharing structure identified in one study, inserted the fixed and random effects coefficients of the longitudinal trajectory into the time-to-event sub-model without their related covariates.
Random effects only sharing structures can also be grouped. We define a random effects only sharing structure to contain covariates if it is of format such as α(U 0 + U 1 t) or α 1 U 0 + α 2 U 1 + α 3(U 0 + U 1 t) where the α terms are association parameters, the U terms are random effects, and t represents a covariate such as time. Alternatively, if the structure is similar to α(U 0 + U 1), where the random effect U 1 had a covariate t in the longitudinal sub-model, we define the random effects only sharing structure to not contain covariates (see Henderson et al  for further examples).
The specialist sharing structure group (4 studies (6.2%)) contained less common sharing structures such as associating the time-to-event and longitudinal sub-models through a multivariate distribution. Another option was the latent class structure [3, 82], used in at least one joint model in 3 studies (4.6%). Finally 4 (6.2%) studies, including the 2 abstracts, had unclear sharing structures.
We should note that choice of association structure should be driven by the data itself, and so it is expected to see a range of sharing structures given the range of disease areas of the identified studies.
Software used in joint model fits in included studies
R (JM) 
R (JMBayes) 
R (joineR) 
R (frailtypack) 
R (unspecified package)
R (own code developed, unclear if available)
SAS (PROC NLMIXED)
SAS (own code available)
JM Macro 
JMFit Macro 
Stata (stjm) 
WinBUGS (own code available)
WinBUGS (no available code)
OpenBUGS (no available code)
Fortran (code available, not study specific)
Fortran (own code developed)
Fortran (study states code available)
NONMEM (unspecified) 
C++ (own code unclear if available) 
Mplus (unspecified) 
More than one software listed/potentially used
R (JM) or SAS (unspecified)
R or SAS (unspecified)
WinBUGS and R (Directed Acyclic Graph provided)
WinBUGS and R (own code available)
Software used was not always stated by identified studies, a potential issue for future MA when determining exact modelling structure used. The most mentioned software was R, although SAS and Stata were also common. Some of the software identified requires more coding from users (such as C++ WinBUGs), which might explain the preference for software with specific joint modelling packages.
The current preference in R is for the JM package  (implements frequentist joint models that insert the fixed and random effects of the longitudinal sub-model into the time-to-event sub-model), and PROC NLMIXED in SAS (allows fitting of a range of non-linear mixed models). The popularity of JM might be explained by the availability of a textbook with worked examples of joint model implementation using the package .
In four studies (6.2%) more than one software was stated, it was unclear which implemented the joint model fit.
Did studies report sufficient information to contribute to meta-analyses?
Summary of information available to contribute to meta-analysis
Coefficients reported (%)
Precision reported (%)
Standard Errors reported (%)
Confidence Intervals (CI) reported (%)
Significance level reported (%)
Sample size reported (%)
MA possible given reported information (%)
All identified studies (N = 65)
Studies using joint models to account for dropout (N = 22)
Studies using joint models to include time varying covariate in time-to-event sub-model (N = 4)
The sample size was reported in the majority (98.5%) of the identified studies, with median sample size of 514 (IQR 277–1054.5, range 46–3814).
The association parameters were more commonly reported (51 studies (78.5%)) than those of the longitudinal and time-to-event sub-models (45 (69.2%) and 46 (70.8%) respectively). This could be attributable to the high proportion of identified studies that stated that joint models were used to investigate the link between longitudinal and time-to-event outcomes.
The number of studies where a precision measure was available was comparable to the number of studies that reported coefficients, across the three MA categories. Whilst a MA would be possible for each category if the parameter, the standard error and the sample size were reported, if the only available precision estimate was the confidence interval then the significance level was also required. In the studies we identified, the significance level was unclear for 8 (12.3%) studies, was 0.05 for 53 (81.5%) studies, 0.01 for 3 (4.6%) studies, 0.1 for 1 (1.5%) study.
Overall, a MA would be possible for the association parameter in 50 (76.9%) studies, for the longitudinal parameters in 44 (67.7%) studies, and for time-to-event parameters in 45 (69.2%) studies.
Ideally, a study would provide enough information to perform MA in all three separate groups. Sufficient information to allow all three MA to be undertaken was available from 38 (58.5%) of the studies. Only two MA were possible in 6 studies (9.2%), only 1 in 13 studies (20.0%) and in 8 studies (12.3%) there was not sufficient information to complete any MA.
The reasons for joint model use may affect the information stated in the study report. We re-examined the proportions of studies in which the three groups of MA could take place, dependent on the reason for joint model use (Table 4).
For the 22 studies which stated accounting for dropout as a reason for joint model use, we saw a much higher percentage (81.8%) for which MA of longitudinal sub-model parameters was possible, compared to for all studies (67.7%). However percentages of MA possible for the time-to-event coefficients or association parameters were smaller. This could be explained by studies using joint models to account for dropout being mainly interested in the parameters from the longitudinal sub-model.
Only 4 studies stated inclusion of a time varying covariate in a time-to-event model as one of their reasons for using joint models. There was a slight indication that the longitudinal coefficients for studies using joint models to include time varying covariates in time-to-event models are worse reported than for all identified studies, possibly because the longitudinal component of the joint model is of interest as a covariate rather than an outcome in these cases, however more information is needed before this relationship can be fully investigated.
Joint models for longitudinal and time-to-event data are often stated as beneficial compared to separate longitudinal or time-to-event analyses, as they can reduce bias and increase efficiency in model estimation (see Ibrahim et al 2010  for example). Additionally Powney et al 2014  discuss a study where joint models showed a significant difference between treatment groups that was not identified by separate analyses . These benefits of joint models reinforce the suggestion that in certain circumstances MA of joint models may be more appropriate than MA of separate models.
Recommendations for future reporting of joint models
1. Model structure (longitudinal and time-to-event sub-models, and association structure) be clearly stated
2. Estimates of all model coefficients and their precisions be available in main or supplementary material
3. Software (and packages) used to fit models be stated, or code available on request from author
4. Statistical methods should be mentioned in text accessible to search engine to aid identification of papers in meta-analyses
With the increasing use of joint models in the literature, ensuring they are well reported is vital so that the analyses can be interpreted fully and that the published data can be used in future evidence synthesis. We have identified that for the scenario of AD-MA of the results of joint models published in the literature, it was possible to perform MA from a high proportion but not all studies. We would recommend for future practice that regardless of the reason for joint model use, full model covariates with precision estimates be reported either in the study report or supplementary materials (Table 5) not only to aid interpretation within the study itself, but to ensure that future MA would be possible.
Additionally it is important that model structure (both sub-models and the sharing structures that link them) are clearly reported (Table 5). However in some cases model formulae were not reported, a particular issue identified for association parameters. We have noted the multiple association structures available to researchers. Due to the different interpretation of each, it may not make sense to pool association parameters from radically different sharing structures. Without clear statement of the model structure it may be difficult to conduct a future MA. It would be beneficial if studies applying joint models included statement of the model structure as standard to give clarity to the methods used.
Also, it is important especially in joint modelling to state the software and packages or functions used, as the software used could indicate the model structure, as well as methods used to obtain the parameter estimations (Table 5).
Consistently, the proportion of studies where association parameter MA was possible was higher than the proportion where MA was possible for longitudinal or time-to-event parameters. The association information may have been more commonly reported than the longitudinal or time-to-event information because the association parameters in shared random effect models quantify the link between the sub-models. It can be expected therefore that studies aiming to quantify the link between the sub-models report the association information more prominently than other model parameters.
Overall, this investigation has highlighted the need to fully report the coefficients and precision estimates in studies applying joint models to datasets. Whilst this review identified a limited number of studies that fulfilled our criteria, the range of disease areas covered by the studies was large with applications using a wide range of sub-models and association structures being published in an assortment of journals by a range of authors. It may be some years before there are sufficient studies published in one area to conduct a meta-analysis of joint models. Nevertheless our review has demonstrated that the use of joint models is increasing every year and the availability of software to fit a range of flexible models is likely to facilitate their application even further.
In the future the current standard of joint modelling should be maintained and improved upon, following the recommendations given in this paper, in order that information from published studies can be used for the purpose of MA.
This work was supported by the Health eResearch Centre (HeRC) funded by the Medical Research Council Grant MR/K006665/1.
Availability of data and material
The datasets analysed during the review are available from the corresponding author on reasonable request. The data collection form is included in the Additional file 2.
MS, RKD and CTS all participated in discussion of the design of the investigation. MS extracted data from the identified studies, conducted analyses and was a major contributor in writing the manuscript. All authors read and approved the final manuscript.
MS and CTS were authors on one of the studies identified in this review. RKD was involved in the original grant for the development of the R joint modelling package joineR.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Tsiatis AA, Davidian M. Joint modeling of longitudinal and time-to-event data: An overview. Stat Sinica. 2004;14(3):809–34.Google Scholar
- Gould AL, et al. Responses to discussants of 'Joint modeling of survival and longitudinal non-survival data: Current methods and issues. report of the DIA Bayesian joint modeling working group'. Stat Med. 2015;34(14):2202–3.PubMedPubMed CentralView ArticleGoogle Scholar
- Proust-Lima C, et al. Joint latent class models for longitudinal and time-to-event data: A review. Stat Methods Med Res. 2014;23(1):74–90.PubMedView ArticleGoogle Scholar
- Henderson R, Diggle P, Dobson A. Joint modelling of longitudinal measurements and event time data. Biostatistics. 2000;1(4):465–80.PubMedView ArticleGoogle Scholar
- Glass GV. Primary, Secondary, and Meta-Analysis of Research. Educ Res. 1976;5(10):3–8.View ArticleGoogle Scholar
- Higgins J, Green S. and (editors). Cochrane Handbook for Systematic Reviewsof Interventions Version 5.1.0 (updated March 2011). The Cochrane Collaboration, 2011. Available from http://handbook.cochrane.org/.
- Whitehead A. Meta-analysis of controlled clinical trials/Anne Whitehead. Statistics in practice. Chichester: John Wiley & Sons; 2002.View ArticleGoogle Scholar
- Moher D, et al. Methods of Systematic Reviews and Meta-Analysis: Preferred Reporting Items for Systematic Reviews and Meta-Analyses: The PRISMA Statement. J Clin Epidemiol. 2009;62:1006–12.PubMedView ArticleGoogle Scholar
- Andrinopoulou ER, et al. An introduction to mixed models and joint modeling: Analysis of valve function over time. Ann Thorac Surg. 2012;93(6):1765–72.PubMedView ArticleGoogle Scholar
- Andrinopoulou ER, et al. Joint modeling of two longitudinal outcomes and competing risk data. Stat Med. 2014;33(18):3167–78.PubMedView ArticleGoogle Scholar
- Argyropoulos C, et al. Dialyzer reuse and outcomes of high flux dialysis. PLoS ONE. 2015;10(6):e0129575.PubMedPubMed CentralView ArticleGoogle Scholar
- Artaud F, et al. Decline in fast gait speed as a predictor of disability in older adults. J Am Geriatr Soc. 2015;63(6):1129–36.PubMedView ArticleGoogle Scholar
- Berzuini C, Allemani C. Effectiveness of potent antiretroviral therapy on progression of human immunodeficiency virus: Bayesian modelling and model checking via counterfactual replicates. J Royal Stat Soc Series C Appl Stat. 2004;53(4):633–50.View ArticleGoogle Scholar
- Brombin C, Di Serio C, Rancoita PM. Joint modeling of HIV data in multicenter observational studies: A comparison among different approaches. Stat Methods Med Res. 2016;25(6):2472–87. Epub 2014 Mar 26.Google Scholar
- Brombin C, Di Serio C, Rancoita PM. Joint modeling of HIV data in multicenter observational studies: A comparison among different approaches. Stat Methods Med Res. 2016;25(6):2472-487. Epub 2014Mar 26.Google Scholar
- Chang CCH, et al. Minimizing attrition bias: a longitudinal study of depressive symptoms in an elderly cohort. Int Psychogeriatr. 2009;21(5):869–78.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen Q, et al. Joint modeling of longitudinal and survival data with missing and left-censored time-varying covariates. Stat Med. 2014;33(26):4560–76.PubMedPubMed CentralView ArticleGoogle Scholar
- Daher Abdi Z, et al. Impact of longitudinal exposure to mycophenolic acid on acute rejection in renal-transplant recipients using a joint modeling approach. Pharmacol Res. 2013;72:52–60.PubMedView ArticleGoogle Scholar
- Dantan E, et al. Joint model with latent state for longitudinal and multistate data. Biostatistics. 2011;12(4):723–36.PubMedView ArticleGoogle Scholar
- Deslandes E, Chevret S. Joint modeling of multivariate longitudinal data and the dropout process in a competing risk setting: Application to ICU data. BMC Med Res Methodol. 2010;10:69.PubMedPubMed CentralView ArticleGoogle Scholar
- Du H, Hahn EA, Cella D. The impact of missing data on estimation of health-related quality of life outcomes: An analysis of a randomized longitudinal clinical trial. Health Services and Outcomes Research Methodology. 2011;11(3-4):134–44.View ArticleGoogle Scholar
- Duvignac J, et al. CD4+ T-lymphocytes natural decrease in HAART-naive HIV-infected adults in Abidjan. HIV Clinical Trials. 2008;9(1):26–35.PubMedView ArticleGoogle Scholar
- Ediebah DE, et al. Joint modeling of longitudinal health-related quality of life data and survival. Qual Life Res. 2015;24(4):795–804.PubMedView ArticleGoogle Scholar
- Elhaj M, et al. Can serum surfactant protein D or CC-chemokine ligand 18 predict outcome of interstitial lung disease in patients with early systemic sclerosis? J Rheumatol. 2013;40(7):1114–20.PubMedPubMed CentralView ArticleGoogle Scholar
- Fairclough DL, et al. Handling missing quality of life data in HIV clinical trials: what is practical? Qual Life Res. 2008;17(1):61–73.PubMedView ArticleGoogle Scholar
- Garre FG, et al. A joint latent class changepoint model to improve the prediction of time to graft failure. J Royal Stat Soc Series A Stat Soc. 2008;171(1):299–308.Google Scholar
- Ghisletta P. Application of a joint multivariate longitudinal-survival analysis to examine the terminal decline hypothesis in the Swiss Interdisciplinary Longitudinal Study on the Oldest Old. J Gerontol B Psychol Sci Soc Sci. 2008;63(3):185–92.View ArticleGoogle Scholar
- Ghisletta P, McArdle JJ, Lindenberger U. Longitudinal cognition-survival relations in old and very old age: 13-year data from the Berlin aging study. Eur Psychol. 2006;11(3):204–23.View ArticleGoogle Scholar
- Goyal N, Gomeni R. A latent variable approach in simultaneous modeling of longitudinal and dropout data in schizophrenia trials. Eur Neuropsychopharmacol. 2013;23(11):1570–6.PubMedView ArticleGoogle Scholar
- Graham PL, Ryan LM, Luszcz MA. Joint modelling of survival and cognitive decline in the Australian Longitudinal Study of Ageing. J Royal Stat Soc Series C Appl Stat. 2011;60(2):221–38.View ArticleGoogle Scholar
- Guler I, et al. Joint modelling for longitudinal and time-to-event data: Application to liver transplantation data, in 14th International Conference on Computational Science and Its Applications, ICCSA 2014. Guimaraes: Springer Verlag; 2014. p. 580–93.Google Scholar
- Husser O, et al. Tumor marker carbohydrate antigen 125 predicts adverse outcome after transcatheter aortic valve implantation. JACC Cardiovasc Interv. 2013;6(5):487–96.PubMedView ArticleGoogle Scholar
- Ibrahim JG, Chen MH, Sinha D. Bayesian methods for joint modeling of longitudinal and survival data with applications to cancer vaccine trials. Stat Sin. 2004;14(3):863–83.Google Scholar
- Jacoby A, et al. Quality-of-life outcomes of initiating treatment with standard and newer antiepileptic drugs in adults with new-onset epilepsy: findings from the SANAD trial. Epilepsia. 2015;56(3):460–72.PubMedView ArticleGoogle Scholar
- Khoundabi B, et al. Acute kidney injury in ICU patients following non-cardiac surgery at Masih Daneshvari hospital: joint modeling application. Tanaffos. 2015;14(1):49–54.PubMedPubMed CentralGoogle Scholar
- Kupferman JC, et al. BP control and left ventricular hypertrophy regression in children with CKD. J Am Soc Nephrol. 2014;25(1):167–74.PubMedView ArticleGoogle Scholar
- Kypriotakis G, Vidrine DJ, Francis LE, Rose JH. The longitudinal relationship between quality of life and survival in advanced stage cancer. Psychooncology. 2016;25(2):225–31. doi:10.1002/pon.3846. Epub 2015 May 8.
- Lawson A.B, et al. Bayesian latent structure modeling of walking behavior in a physical activity intervention. Stat Methods Med Res. 2014Google Scholar
- Levine SZ, et al. Joint modeling of dropout and outcome in three pivotal clinical trials of schizophrenia. Schizophr Res. 2015;164(1-3):122–6.PubMedView ArticleGoogle Scholar
- Liu X, et al. Does C-reactive protein predict the long-term progression of interstitial lung disease and survival in patients with early systemic sclerosis? Arthritis Care and Research. 2013;65(8):1375–80.PubMedView ArticleGoogle Scholar
- Malehi AS, et al. Assessing the autoantibody levels in relation to recurrence of pemphigus: Joint modeling of longitudinal measurements and recurrent event times. Iran Red Crescent Med J. 2014;16(2):e13812.Google Scholar
- Marioni RE, et al. Cognitive lifestyle jointly predicts longitudinal cognitive decline and mortality risk. Eur J Epidemiol. 2014;29(3):211–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Mayeda ER1, Haan MN, Yaffe K, Kanaya AM, Neuhaus J. Does type 2 diabetes increase rate of cognitive decline in older Mexican Americans? Alzheimer Dis Assoc Disord. 2015;29(3):206-12. doi:10.1097/WAD.0000000000000083.
- McLain AC, Sundaram R, Buck Louis GM. Joint analysis of longitudinal and survival data measured on nested timescales by using shared parameter models: An application to fecundity data. J Royal Stat Soc Series C Appl Stat. 2015;64(2):339–57.View ArticleGoogle Scholar
- Medeiros FA, et al. Evaluation of progressive neuroretinal rim loss as a surrogate end point for development of visual field loss in glaucoma. Ophthalmology. 2014;121(1):100–9.PubMedView ArticleGoogle Scholar
- Meira-Freitas D, et al. Predicting progression in glaucoma suspects with longitudinal estimates of retinal ganglion cell counts. Investig Ophthalmol Vis Sci. 2013;54(6):4174–83.View ArticleGoogle Scholar
- Meira-Freitas D, et al. Predicting progression of glaucoma from rates of frequency doubling technology perimetry change. Ophthalmology. 2014;121(2):498–507.PubMedView ArticleGoogle Scholar
- Miki A, et al. Rates of retinal nerve fiber layer thinning in glaucoma suspect eyes. Ophthalmology. 2014;121(7):1350–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Murphy TE, et al. A longitudinal, observational study with many repeated measures demonstrated improved precision of individual survival curves using bayesian joint modeling of disability and survival. Exp Aging Res. 2015;41(3):221–39.PubMedPubMed CentralView ArticleGoogle Scholar
- Njagi EN, et al. A joint survival-longitudinal modelling approach for the dynamic prediction of rehospitalization in telemonitored chronic heart failure patients. Stat Model. 2013;13(3):179–98.View ArticleGoogle Scholar
- Núñez J, et al. Red blood cell distribution width is longitudinally associated with mortality and anemia in heart failure patients. Circ J. 2014;78(2):410–8.PubMedView ArticleGoogle Scholar
- Pan J, et al. Joint longitudinal and survival-cure models in tumour xenograft experiments. Stat Med. 2014;33(18):3229–40.PubMedView ArticleGoogle Scholar
- Pauler DK, Finkelstein DM. Predicting time to prostate cancer recurrence based on joint models for non-linear longitudinal biomarkers and event time outcomes. Stat Med. 2002;21(24):3897–911.PubMedView ArticleGoogle Scholar
- Paulsen JS, et al. Prediction of manifest Huntington's disease with clinical and imaging measures: a prospective observational study. Lancet Neurol. 2014;13(12):1193–201.PubMedPubMed CentralView ArticleGoogle Scholar
- Pike F, Weissfeld LA, Chang CCH. Joint modeling of multivariate censored longitudinal and event time data with application to the Genetic Markers of Inflammation Study. J Appl Stat. 2014;41(10):2178–91.View ArticleGoogle Scholar
- Pilla Reddy V, et al. Modelling and simulation of the positive and negative syndrome scale (PANSS) time course and dropout hazard in placebo Arms of schizophrenia clinical trials. Clin Pharmacokinet. 2012;51(4):261–75.PubMedView ArticleGoogle Scholar
- Sargent-Cox KA, Anstey KJ, Luszcz MA. Longitudinal change of self-perceptions of aging and mortality. J Gerontol B Psychol Sci Soc Sci. 2014;69(2):168–73.PubMedView ArticleGoogle Scholar
- Seid A, et al. Joint modeling of longitudinal CD4 cell counts and time-to-default from HAART treatment: A comparison of separate and joint models. Electron J Appl Stat Anal. 2014;7(2):292–314.Google Scholar
- Shen Y, et al. Joint modeling tumor burden and time to event data in oncology trials. Pharm Stat. 2014;13(5):286–93.PubMedView ArticleGoogle Scholar
- Skupien J, et al. Synergism between circulating tumor necrosis factor receptor 2 and HbA(1c) in determining renal decline during 5-18 years of follow-up in patients with type 1 diabetes and proteinuria. Diabetes Care. 2014;37(9):2601–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Su CT, et al. Changes in anthropometry and mortality in maintenance hemodialysis patients in the HEMO study. Am J Kidney Dis. 2013;62(6):1141–50.PubMedView ArticleGoogle Scholar
- Tang X, et al. Increasing Chimerism after Allogeneic Stem Cell Transplantation Is Associated with Longer Survival Time. Biol Blood Marrow Transplant. 2014;20(8):1139–44.PubMedPubMed CentralView ArticleGoogle Scholar
- Terrin N, Rodday AM, Parsons SK. Joint models for predicting transplant-related mortality from quality of life data. Qual Life Res. 2015;24(1):31–9.PubMedView ArticleGoogle Scholar
- Thabut G, et al. Survival benefit of lung transplant for cystic fibrosis since lung allocation score implementation. Am J Respir Crit Care Med. 2013;187(12):1335–40.PubMedPubMed CentralView ArticleGoogle Scholar
- Thiebaut R, et al. Joint modelling of bivariate longitudinal data with informative dropout and left-censoring, with application to the evolution of CD4+ cell count and HIV RNA viral load in response to treatment of HIV infection. Stat Med. 2005;24(1):65–82.PubMedView ArticleGoogle Scholar
- van den Hout A, Muniz-Terrera G. Joint models for discrete longitudinal outcomes in aging research. J. R. Stat. Soc. C. 2016;65:167–186. doi:10.1111/rssc.12114 Issue online: 17 December 2015.
- Van Der Linde D, et al. Natural history of discrete subaortic stenosis in adults: A multicentre study. Eur Heart J. 2013;34(21):1548–56.PubMedView ArticleGoogle Scholar
- Wang Y, Taylor JMG. Jointly Modeling Longitudinal and Event Time Data with Application to Acquired Immunodeficiency Syndrome. J Am Stat Assoc. 2001;96(455):895–905.View ArticleGoogle Scholar
- Watanabe H, et al. Factors affecting longitudinal functional decline and survival in amyotrophic lateral sclerosis patients. Amyotroph Lateral Scler Frontotemporal Degeneration. 2015;16(3-4):230–6.View ArticleGoogle Scholar
- Wu L, Liu W, Hu XJ. Joint inference on HIV viral dynamics and immune suppression in presence of measurement errors. Biometrics. 2010;66(2):327–35.PubMedView ArticleGoogle Scholar
- Xue QL, et al. Heterogeneity in rate of decline in grip, hip, and knee strength and the risk of all-cause mortality: The women's health and aging study II. J Am Geriatr Soc. 2010;58(11):2076–84.PubMedPubMed CentralView ArticleGoogle Scholar
- Yu M, Taylor JMG, Sandler HM. Individual prediction in prostate cancer studies using a joint longitudinal survival-cure model. J Am Stat Assoc. 2008;103(481):178–87.View ArticleGoogle Scholar
- Zhang JP, et al. Joint modeling of longitudinal changes in depressive symptoms and mortality in a sample of community-dwelling elderly people. Psychosom Med. 2009;71(7):704–14.PubMedPubMed CentralView ArticleGoogle Scholar
- Wulfsohn MS, Tsiatis AA. A joint model for survival and longitudinal data measured with error. Biometrics. 1997;53(1):330–9.PubMedView ArticleGoogle Scholar
- Guo X, Carlin BP. Separate and Joint Modeling of Longitudinal and Event Time Data Using Standard Computer Packages. Am Stat. 2004;58(1):16–24.View ArticleGoogle Scholar
- Rizopoulos D. JM: An R Package for the Joint Modelling of Longitudinal and Time-to-Event Data. J Stat Softw. 2010;35(9):1–33.View ArticleGoogle Scholar
- Rizopoulos D. Joint models for longitudinal and time-to-event data, with applications in R. Boca Raton: Chapman and Hall/CRC; 2012.View ArticleGoogle Scholar
- Crowther MJ, Abrams KR, Lambert PC. Flexible parametric joint modelling of longitudinal and survival data. Stat Med. 2012;31(30):4456–71.PubMedView ArticleGoogle Scholar
- Philipson, P, et al. joineR: Joint modelling of repeated measurements and time-to-event data. Comprehensive R Archive Network, 2012.Google Scholar
- Crowther, Michael J., (2013), STJM: Stata module to fit shared parameter joint models of longitudinal and survival data, http://EconPapers.repec.org/RePEc:boc:bocode:s457502.
- Crowther MJ, Abrams KR, Lambert PC. Joint modeling of longitudinal and survival data. Stata J. 2013;13(1):165–84.Google Scholar
- Proust-Lima C, et al. Joint modelling of multivariate longitudinal outcomes and a time-to-event: A nonlinear latent class approach. Comput Stat Data Anal. 2009;53(4):1142–54.View ArticleGoogle Scholar
- Rizopoulos D, Ghosh P. A Bayesian semiparametric multivariate joint model for multiple longitudinal outcomes and a time-to-event. Stat Med. 2011;30(12):1366–80. 15.PubMedView ArticleGoogle Scholar
- Lawrence Gould A, et al. Joint modeling of survival and longitudinal non-survival data: Current methods and issues. Report of the DIA Bayesian joint modeling working group. Stat Med. 2015;34(14):2181–95.PubMedView ArticleGoogle Scholar
- Ibrahim JG, Chu H, Chen LM. Basic concepts and methods for joint models of longitudinal and survival data. J Clin Oncol. 2010;28(16):2796–801.PubMedPubMed CentralView ArticleGoogle Scholar
- Powney M, et al. A review of the handling of missing longitudinal outcome data in clinical trials. Trials. 2014;15.Google Scholar
- Powell C, et al. MAGNEsium Trial In Children (MAGNETIC): a randomised, placebo-controlled trial and economic evaluation of nebulised magnesium sulphate in acute severe asthma in children. Health Technol Assess. 2013;17(32):1–216. 216p.PubMedPubMed CentralGoogle Scholar
- Wulfsohn MS, Tsiatis AA. A joint model for survival and longitudinal data measured with error. Int Biometric Soc. 1997;53(1):330–9.View ArticleGoogle Scholar
- Diggle PJ, Sousa I, Chetwynd AG. Joint modelling of repeated measurements and time-to-event outcomes: The fourth Armitage lecture. Stat Med. 2008;27(16):2981–98.PubMedView ArticleGoogle Scholar
- Rizopoulos D. Dynamic Predictions and Prospective Accuracy in Joint Models for Longitudinal and Time-to-Event Data. Biometrics. 2011;67(3):819–29.PubMedView ArticleGoogle Scholar
- R Core Team and R: A language and environment for statistical computing. In R Foundation for Statistical Computing, Vienna, Austria. 2015. https://www.R-project.org/.
- Rizopoulos D. The R Package JMbayes for Fitting Joint Models for Longitudinal and Time-to-Event Data using MCMC. 2014Google Scholar
- Rondeau V, Mazroui Y, Gonzalez JR. Frailtypack: An r package for the analysis of correlated survival data with frailty models using penalized likelihood estimation or parametrical estimation. J Stat Softw. 2012;47.Google Scholar
- SAS software. Copyright, SAS Institute Inc. SAS and all other SAS Institute Inc. product or service names are registered trademarks or trademarks of SAS Institute Inc., Cary, NC, USA.Google Scholar
- Garcia-Hernandez, A. and Rizopoulos, D. (2016). %JM: A SAS macro to fit jointly generalized mixed models for longitudinal data and time-to-event responses. Journal of Statistical Software, to appear.Google Scholar
- Zhang D, et al. JMFit: A SAS Macro for Assessing Model Fit in Joint Models of Longitudinal and Survival Data. J Stat Softw. 2015. in press.Google Scholar
- StataCorp. 2015. Stata Statistical Software: Release 14. College Station, TX: StataCorp LP.Google Scholar
- Lunn DJ, et al. WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility. Stat Comput. 2000;10(4):325–37.View ArticleGoogle Scholar
- Beal SL, S. L.B. NONMEM Users Guides. Icon Development Solutions. Ellicott City, Maryland, USA. 1989-98.Google Scholar
- ISO International Standard ISO/IEC 14882:2014(E) – Programming Language C++.Google Scholar
- Muthén L.K and Muthén B.O. Mplus User’s Guide. Sixth Edition. Los Angeles, CA: Muthén & Muthén; 1998-2010.Google Scholar