 Research article
 Open access
 Published:
Using joint models to disentangle intervention effect types and baseline confounding: an application within an intervention study in prodromal Alzheimer’s disease with Fortasyn Connect
BMC Medical Research Methodology volume 19, Article number: 163 (2019)
Abstract
Background
Many prodromal Alzheimer’s disease trials collect two types of data: the time until clinical diagnosis of dementia and longitudinal patient information. These data are often analysed separately, although they are strongly associated. By combining the longitudinal and survival data into a single statistical model, joint models can account for the dependencies between the two types of data.
Methods
We illustrate the major steps in a joint modelling approach, motivated by data from a prodromal Alzheimer’s disease study: the LipiDiDiet trial.
Results
By using joint models we are able to disentangle baseline confounding from the intervention effect and moreover, to investigate the association between longitudinal patient information and the time until clinical dementia diagnosis.
Conclusions
Joint models provide a valuable tool in the statistical analysis of clinical studies with longitudinal and survival data, such as in prodromal Alzheimer’s disease trials, and have several added values compared to separate analyses.
Background
Alzheimer’s disease (AD) is a neurodegenerative disorder characterised by a slow progressive deterioration of cognitive capacity. The pathophysiological changes begin long before clinical manifestations of the disorder, and the disease spectrum spans from clinically asymptomatic to severely impaired [1]. The terminology of prodromal AD designates the initial mild state of cognitive impairment, whereas the dementia state represents the subsequent clinically manifest severe cognitive impairment. The specific transition between prodromal AD and the clinical diagnosis of AD dementia can be challenging [2] as AD should not be viewed with discrete and defined clinical stages, but as a multifaceted process moving along a biological and clinical continuum [1]. Given this underlying continuum, the moment of receiving the dementia diagnosis does not represent a discrete biological event. Nonetheless, having received the dementia diagnosis does indicate a certain level of disease progression. As such, the event ‘AD dementia diagnosis’ has been used in many studies that focus on risk factors, see for example: [3–5], and has obvious impact on patient care.
Prodromal AD trials frequently collect the time until clinical dementia diagnosis in combination with longitudinal patient information. These longitudinal patient information include clinical biomarkers or performance of patients in psychometric tests and can help to describe or understand disease progression. Yet, most studies dealing with longitudinal and survival (i.e., timetoevent) data analyse the data separately, mostly by relying on wellestablished statistical methods such as linear mixed models for longitudinal data and Cox proportional hazard models for survival data. However, a method that allows the simultaneous modelling of longitudinal measurements with a survival outcome is the joint model for longitudinal and survival data, see for example: Wulfsohn & Tsiatis (1997) [6], Henderson et al. (2000) [7], Tsiatis & Davidian (2004) [8] and Rizopoulos (2012) [9]. By combining the longitudinal and survival data into a single statistical model, joint models can account for or infer the dependencies between the two types of data. In certain situations, e.g., when it is of interest to study the association between a clinical biomarker or cognitive measure over time and the time until clinical diagnosis, a joint modelling approach is even required. More specifically, when it is of interest to study the association between a survival outcome and an endogenous timevarying covariate, such as a biomarker or another covariate measured on patients during the study, the traditional Cox model is not appropriate [10, 11]. First approaches to fit joint models have focused on the socalled twostage methods, in which as a first step, a model is fit to the longitudinal data, and as a second step, the fitted longitudinal values are inserted in the Cox model. Many authors, such as Dafni & Tsiatis (1998) [12], Tsiatis & Davidian (2001) [13] and Sweeting & Thompson (2011) [10], have shown that the twostage method still provides potentially biased and inefficient estimates. In comparison, the joint model simultaneously estimates the parameters in the longitudinal and survival parts of the model, for example by relying on maximum likelihood estimation.
Joint modelling is an active area in biostatistics with numerous methodological papers (within AD research, see for example: [14–17]) and has already been adopted in several clinical research fields such as cancer [18, 19] and cardiovascular disease [20, 21]. However handson introductions for clinicians are still limited. This paper aims to provide an introduction into the application of joint models, motivated by data from a prodromal AD trial: the LipiDiDiet trial [22].
The LipiDiDiet trial is a randomised controlled trial, with the objective of assessing the effect of medical nutrition (Souvenaid) on cognitive functioning in patients with prodromal AD. The active component of Souvenaid is Fortasyn Connect, a specific nutrient combination designed to address nutritional requirements in the presence of AD pathology [23]. In the paper on the LipiDiDiet trial’s main results, longitudinally measured variables of cognition and time to dementia diagnosis were analysed separately. In the LipiDiDiet trial the effect on the longitudinally measured primary endpoint related to cognition did not reach significance in the primary model, while in secondary models significance was reached. In addition, benefits were seen on longitudinal measures of cognition and function, and brain atrophy measures, which were secondary outcome measures in the trial [22]. A worsening of cognition is among the criteria for AD dementia diagnosis [24]. One could hypothesise that an intervention that is effective in decreasing or preventing cognitive decline would also prevent or delay the clinical diagnosis. In this paper we show how we can use joint models to optimally utilise the relationship between the longitudinal information and the event times in order to gain understanding into the process of how an intervention affects disease progression. In doing so, the application of joint models reveals relevant information about the strength and the type of the associations between the longitudinal measures of cognition and the risk of an event. Moreover, we investigate the effect of differences in baseline characteristics on study outcome. Using a joint model, we can disentangle baseline confounding from the intervention effect. Throughout the analysis of the data at hand, we aim to introduce and illustrate the major steps in a joint modelling approach for the nonstatistical reader.
Methods
LipiDiDiet trial
The LipiDiDiet trial is a 24month randomised, controlled, doubleblind, multicentre trial, performed at 11 study sites across different countries. The goal of the LipiDiDiet trial was to investigate the effects of Fortasyn Connect on cognition and related measures in prodromal AD patients. For this purpose, several longitudinal measures of cognitive functioning were recorded. In this paper we include two of them: the Clinical Dementia Rating sum of boxes (CDRSB) and memory domain from a neuropsychological test battery (NTB memory domain).
The CDRSB score reflects global clinical impression and ranges from a score of 0 to 18, with a higher score indicating a worse status. It is obtained through a semistructured interview of patients and informants, summing scores of cognitive functioning on each of the following domain box scores: memory, orientation, judgement and problem solving, community affairs, home and hobbies, and personal care.
NTB memory domain is a composite zscore based on Consortium to Establish a Registry for AD (CERAD) 10word list learning immediate recall, CERAD 10word delayed recall, and CERAD 10word recognition. A higher zscore indicates a better memory.
Individual patients’ scores were measured at baseline, where randomisation to either the test or control group took place, as well as around months 12 and 24 with an additional visit around 6 months for NTB memory domain. At each visit it was recorded whether patients had received the diagnosis of dementia. Progression to dementia was diagnosed according to criteria defined by DSMIV, the National Institute of Neurological and Communicative Disorders and Stroke, and the AD and Related Disorders Association criteria for AD.
In this article, we focus on AD dementia as a specific form of dementia. The study sample consisted of 311 patients (modified intentiontotreat population in the LipiDiDiet main paper [22]), of whom 57 (36%) patients in the control group and 62 (41%) in the test group had received the AD dementia diagnosis. The median followup times were respectively 1.96 years in the control, and 1.94 years in the test group. Despite the randomisation procedure, a statistically significant difference between the intervention groups was found in baseline Mini–Mental State Examination (baseline MMSE, p=0.039, twosided ttest), reflecting baseline cognitive performance. The higher baseline MMSE score in the control group denotes better performance and suggests a lower risk of receiving the dementia diagnosis in this group at baseline. Figure 1 displays the histograms of baseline MMSE scores in the test and control group. For further information regarding the LipiDidiet trial, including information on the randomisation procedure, we refer to the LipiDiDiet main paper [22].
Methodology for the standard joint model
As the name suggests, a joint model for longitudinal and survival data consists of a longitudinal submodel and a survival submodel. The longitudinal submodel is typically a mixed effects model aiming to describe the shapes of the patientspecific longitudinal profiles. For continuous longitudinal data, linear mixed models can take into account that repeated measurements from the same patient may be more correlated than measurements from other patients, by including not only fixed effects but also patientspecific random effects. For background information on mixed models, we refer to Verbeke & Molenberghs (1997) [25] and Fitzmaurice et al. (2008) [26].
In order to formulate our longitudinal submodel for the longitudinal trajectories, as a first step we investigated the observed longitudinal profiles for six randomly selected patients. Figures 2 and 3 show the longitudinal profiles for respectively their CDRSB and NTB memory domain observations; these figures show that there is a lot of variation between patients. Therefore, we allowed each patient to have its own trajectory, by incorporating patientspecific intercepts and slopes. For the average CDRSB and NTB memory domain trajectories we used linear effects of time (β_{1}) but more complicated functions of time such as quadratic or higher order polynomials, e.g., using splines, are also possible [27]. We also tried trajectory functions for CDRSB and NTB memory domain using quadratic time effects, but these were found to give similar results (results not shown). To model the effect of Fortasyn Connect, we included both a main effect of the intervention (β_{2}) and an interaction of intervention by time (β_{3}) in order to allow the trajectories of the intervention groups to be different over time. This is necessary, because the intervention is expected to have a gradual effect, possibly resulting in CDRSB and NTBmemory domain levels for the test group that are worsening more slowly. Further, we included and intercept (β_{0}) and main effects for baseline MMSE (β_{4}) and site (β_{5}). This resulted in the following longitudinal submodel for the CDRSB observations, and similarly defined for NTB memory domain
where CDR_{i}(t) are the observed values of CDRSB for patient i at actual time points t, and the time points at which measurements take place may vary between patients. Further b_{i0} and b_{i1} denote respectively the patientspecific intercept and slope. The longitudinal profile of observed values CDR_{i}(t) is broken down in a trajectory function \(\tilde {\text {CDR}}_{i}(t)\) and a random error term ε_{i}(t), which is assumed to be normally distributed. The trajectory function is assumed to describe the true but unobserved trajectory of the longitudinal marker, and as will be seen later, is used in the survival submodel ‘joining’ the two submodels. The main effect of the intervention, β_{2}, denotes the difference between the intervention groups at baseline, while the interaction effect β_{3}, describes the intervention effect over time.
A common choice for the survival submodel is a Cox model, which is used to model the hazard of experiencing the event, i.e., in this case receiving the dementia diagnosis. For background information on Cox models, see Cox (1972) [28], Klein & Moeschberger (1997) [29] and Therneau & Grambsch (2013) [30]. In case the proportional hazard assumption of the Cox model is violated, alternative modelling frameworks for the survival submodel exist, such as the accelerated failure time model [31]. In this paper we formulated our joint model using a Cox model. Additionally, we fitted a joint model using an accelerated failure time model which gave similar findings (results not shown). In the Cox model we included the intervention as a timeindependent effect, and the estimated true trajectory of the longitudinal marker as a timevarying effect. Since there can be variation across sites in how early a patient is diagnosed, we also corrected for site in the survival submodel. The hazard λ_{i}(t) of dementia diagnosis at time t for patient i is therefore modeled using the following survival submodel,
where the parameter α links the longitudinal process, i.e., the trajectory function \(\tilde {\text {CDR}}_{i}(t)\), or similarly \(\tilde {\text {NTB}}_{i}(t)\), to the survival process. More specifically, the quantity exp(α) denotes the hazard ratio at time t for a oneunit increase in the trajectory of the longitudinal marker at the same time point. Further, λ_{0}(t) is the baseline hazard and γ_{1} denotes a direct effect on the survival outcome. To gain a better understanding of how the intervention affects the risk of receiving the dementia diagnosis, and to explain what we mean by a ‘direct effect’, we distinguish three types of coefficients. These are schematically illustrated in Fig. 4. β describes the intervention effect on the longitudinal marker. As indicated before, there are two types of β’s here; β_{2} denoting the difference in the longitudinal outcome between the intervention groups at baseline, and β_{3} describing the intervention effect on the longitudinal outcome over time. Secondly, since the parameter α measures the effect of the longitudinal process on the survival outcome, together, β_{3} and α, quantify the timevarying intervention effect on the risk of receiving the dementia diagnosis manifesting through the longitudinal marker. The third type of parameter involving the intervention is γ_{1} and is directly related to the risk of receiving the dementia diagnosis. Within the joint model we can therefore distinguish the direct process (Fig. 4, bottom arm), capturing the direct effect on the survival outcome, and the indirect process (Fig. 4, upper arm), in which the coefficients quantify the indirect intervention effect on the survival outcome through the longitudinal marker.
As is the case for the Cox model, the intervention effect in the joint model is the hazard ratio of the test versus the control group. In particular, the total intervention effect is the hazard ratio between two generic patients, i in the test group (fortasyn_{i}=1) and i^{′} in the control group (\(\text {\texttt {fortasyn}}_{i'} = 0\phantom {\dot {i}\!}\)) who do not further differ concerning other covariates. In the joint model, this hazard ratio is a combination of the indirect and direct process. For our formulated joint model the hazard ratio for the total intervention effect denotes exp{γ_{1} +α(β_{2}+β_{3}×t)}, with the first part (i.e., γ_{1}), for the direct process and the latter (timevarying) for the indirect process.
Methodology for investigating the baseline confounding
The two processes of the joint model differ in how they handle the aspect of time. The indirect process can model how an intervention effect varies over time by modelling the intervention effect in the mixed model as a divergence of trajectories. In the direct process however, we are dealing with the proportional hazard assumption of the Cox regression model meaning that the direct effect on the survival outcome is assumed to be constant over the whole period. This arm is therefore likely to capture the effects already present right at the start of the intervention period. In a situation such as this one, where the effect of the intervention on the survival outcome  manifesting through the longitudinal marker  is expected to increase gradually over time, but an effect of any possible baseline confounding on the survival outcome is expected to be immediate, the baseline confounding will for a large extent end up in the direct arm of the model. This is a very appealing property of the joint model that makes it a very effective tool to investigate and control for the effect of potential baseline confounding.
MMSE, found to be significantly different at baseline, is noted to be an important predictor for outcome parameters [22]. This suggests that, before the start of the intervention, the test group might on average have been more likely to receive the dementia diagnosis than the control group due to an imbalance of baseline characteristics. This hampers the interpretability of the results as post baseline outcomes are a combination of the intervention effect and the effect of differences already present at baseline. To investigate this, we examined whether the lower baseline MMSE scores in the test group were related to higher risks of receiving the dementia diagnosis at the start of the trial, therefore possibly counteracting the intervention effect. This required fitting an additional joint model in which we corrected for the effect of the baseline MMSE score on dementia diagnosis by including its value in the survival submodel, given by
We will illustrate how the coefficients of this extended joint model can be an effective tool to investigate and control for the effect of potential baseline confounding using the CDRSB data from the LipiDiDiet trial.
Naturally, apart from being a useful property in investigating possible baseline confounding, the combination of an immediate (direct) and a progressive (indirect) effect helps us to understand the process by which the intervention affects the risk of dementia diagnosis.
Methodology for investigating the association between the longitudinal and survival process
Another aspect of the process by which the intervention affects the timing of dementia diagnosis is determined by the type of the association between the longitudinal and the survival process. The joint model defined in the previous section is the standard joint model and assumes that the value of the longitudinal outcome at any time t is related to the risk of an event at the same time point. However, the underlying relationship between the two processes could have a more complex nature. Examples of longitudinal characteristics possibly related to dementia diagnosis, are the current value, the stability at the current moment, the history of the longitudinal profile up to now or combinations of these characteristics [9]. For demonstration purposes we compared joint models that vary with respect to the type of association that is assumed between the longitudinal data i.e, NTB memory domain, and the survival process i.e., timing of dementia diagnosis. We investigated whether, given the current value of NTB memory domain, the rate of change (i.e., the slope) contains any additional information on the risk of receiving a dementia diagnosis. More specifically, the slope indicates by how much the NTB memory domain for a particular patient is increasing or decreasing at a specific time point. This required fitting a joint model in which we included the slope of NTB memory domain as an additional term in the survival submodel, given by
where the slope of NTB memory domain is obtained by taking the derivative of the trajectory function, consisting of the fixed and random effects, that is,
The parameter α_{1} has the same interpretation as the parameter α before, and the parameter α_{2} measures the association between the slope of the NTB memory domain trajectory and the risk of an event at the same time point, holding \(\tilde {\text {NTB}}_{i}(t)\) constant. Using this joint model, two patients with the same level of NTB memory domain at the current moment do not necessarily have to be at equal risk of receiving the dementia diagnosis. For example, if one patient’s NTB memory domain level is decreasing very rapidly while another patient’s NTB memory domain level is remaining constant, it might be more realistic to assume that the first patient has a higher risk of receiving the dementia diagnosis than the latter  although they have the same value at the current moment.
A similarity between this joint model and the standard joint model, is that the risk of an event at the current moment is related to characteristics of the trajectory at that same time point only. However, the risk of receiving the dementia diagnosis may not depend solely on the level of NTB memory domain or its rate of change at the current moment, but it might also be related to the history of the NTB memory domain levels. That is, two patients with the same characteristics at the current moment are not necessarily at the same risk of receiving a dementia diagnosis if their history of NTB memory domain levels were very different. One approach to take the history of NTB memory domain levels into account is by summarising its cumulative effect i.e., the area under the curve (AUC). The area under the curve indicates the cumulative effect of NTB memory domain values for a particular patient up to the current time point. We also investigated this type of association by fitting a joint model with the following survival submodel
where α_{3} measures how strongly the risk of an event at time t is related to the cumulative effect of NTB memory domain for patient i by time point t. A possible limitation of this joint model is that it gives all past values of NTB memory domain the same weight in terms of their impact on the risk of receiving the dementia diagnosis at the current time point. This may not always be a reasonable assumption. As an alternative, a weight function can be used that places different weights at different time points, for example to give more weight to more recent values of the longitudinal marker. For information on how to use this weight function we refer to [9].
Figure 5 gives a graphical representation of different ways of modelling the association, respectively using the current value, the current value plus the rate of change and the cumulative effect of the longitudinal trajectory.
All the statistical analyses in this paper were performed with statistical software package R, using Rpackage JM [32]. The package uses maximum likelihood for the parameter estimation and assumes rightcensoring. The R code to fit the joint models can be found in the web appendix (Additional file 1).
Results
Note that results in this paper can to some extent differ from results in the LipiDiDiet main paper [22], since different types of models are used. In the main paper, mixed models were used that included the outcome baseline value as a covariate, according to a prespecified statistical analysis plan. In the results presented below, the mixed model approach is part of the joint models and in this mixed model approach, the outcome baseline values are included in the longitudinal trajectory. Modelling outcome baseline values as part of the trajectory is preferred in the joint model context as it maximises the amount of information that is used to estimate the association between the longitudinal data and the survival data.
Results of the standard joint model
Parameter estimates, standard errors, and associated pvalues for the standard joint model are presented in Tables 1 and 2a, respectively for CDRSB and NTB memory domain. Not surprisingly, from the longitudinal submodels we observe that the CDRSB and NTB memory domain scores significantly worsen over time, reflected by an increase of on average 0.61 (95% CI: 0.520.70) per year for CDRSB and a decrease of on average 0.10 (95% CI: 0.040.16) per year for NTB memory domain. For the CDRSB score, however, we see that there is significantly less worsening over time in the test group than in the control group, with the average increase being 0.23 (95% CI: 0.100.37) per year less in the test group compared to the control group.
Further, we observe that both scores have strong associations with the risk of receiving the dementia diagnosis. In particular, a unit increase in CDRSB corresponds to a exp(α) = 2.0fold increase (95% CI: 1.72.3), and a 0.2 unit decrease in NTB memory domain corresponds to a exp(−α×0.2) = 1.3fold increase (95% CI: 1.21.4) in the risk of receiving the dementia diagnosis. Thus, as expected, high values for CDRSB and low values for NTB memory domain are associated with higher risks of receiving the dementia diagnosis. Note that the association for NTB memory domain (zscore) is reported per 0.2unit increase, instead of per 1 unit, since the former denotes a more realistic increase.
Results investigating the baseline confounding
We notice from the results in Table 1a that the coefficients β_{3} and α, which together quantify the indirect intervention effect, are both significant. These results suggest that the intervention decreases the risk of receiving the dementia diagnosis through its effect on CDRSB. Simultaneously, not surprisingly given the baseline imbalance, we observe a nearly significant direct effect with the test group being exp(γ_{1})=1.5fold more likely to receive the dementia diagnosis than the control group. As explained above, the direct effect measures a constant effect over time, due to the proportional hazard assumption of the survival submodel, and is therefore likely to capture possible effects of the baseline confounding. The significant direct effect in favour of the control group is therefore an indication of baseline confounding, also supported by the baseline difference in MMSE.
Comparing the results for γ_{1} of the model with (Table 1b) versus the model without baseline MMSE correction (Table 1a), we observe that by correcting for baseline MMSE in the survival submodel, the direct effect shrinks. This is also illustrated in Fig. 6 where the effects of the coefficients for the separate components of the joint models are displayed over time. Comparing the effects of exp(γ_{1}) (the dashed lines) for 6a versus 6b shows that including the baseline MMSE correction, made the estimate for the direct effect shift towards a hazard ratio of 1, meaning no difference. Table 1 and Fig. 6 also show that the estimates for the indirect effect components (β_{2}, β_{3} and α) are hardly affected by the in  or exclusion  of baseline MMSE in the survival submodel. Based on these results, we hypothesise that the baseline confounding in MMSE is indeed directly related to dementia diagnosis and that it masks the total intervention effect, being a combination of the indirect and direct processes. The latter is graphically illustrated in Fig. 7, in which the total intervention effect from the joint model  that is, the combination of the separate components of Fig. 6  is displayed as a solid line.
Figure 7 also shows the hazard ratios for the intervention effect on dementia diagnosis as estimated from a separately run Cox model (dashed lines). We observe that by using a joint model, and more specifically by incorporating the increasing intervention effect on the longitudinal marker, we can model an increasing intervention effect over time on the risk of dementia diagnosis. While by using the (standard) Cox model, with the underlying proportional hazard assumption, the intervention effect is assumed to be constant over time from baseline onward.
Results for the association between the longitudinal and survival process
Table 2 presents the results of joint models using the current value plus slope (b), and the cumulative effect (c) of the NTB memory domain trajectory for the link between the two processes. From the two types of joint models we observe similar results on the longitudinal process. From the association parameters, we observe that, as expected, decreasing trajectories and small cumulative values for NTB memory domain are associated with higher risks of receiving the dementia diagnosis. Both the rate of increase and the cumulative effect are strongly associated with the risk for dementia diagnosis. For example, if a patient’s NTB memory domain score decreases by 0.2 units faster per year, or 1/60 units faster per month, then the risk of dementia diagnosis is associated with a exp(−α_{2}×0.2) = 3.9fold (95% 1.98.0) increase in the hazard. In the same way, if the cumulative effect of the history of the NTB memory domain levels (i.e., AUC) decreases with one unit, then this corresponds to a exp(−α) = 2.0fold increase (95% CI: 1.62.5) in the risk of dementia diagnosis.
We compared the two alternative types of joint models with the standard joint model based on measures for the model fit (information criteria AIC and BIC). Both measures indicated that the joint model using the current value plus slope is the best fitting joint model, suggesting that inclusion of the slope of the NTB memory domain trajectory improves the fit of the model compared to the standard joint model. The model using a cumulative effect was not found to have a better fit to the data than the standard joint model.
Discussion
Scientists within the (prodromal) AD research field have much to gain from joint models for longitudinal and survival data. When estimating the time until clinical dementia diagnosis, while accounting for the effect of a longitudinal biomarker or cognitive measure, joint models can not only provide estimates for their association, but they can also further investigate the type of association.
This paper aimed to provide an introduction into the application of joint models with special interest in the relationship between the longitudinal information and the event times, using data from a prodromal AD trial. First of all, we reanalysed the data, combining the longitudinal data on cognitive functioning with the survival data on dementia diagnosis in order to account for their dependencies. Both longitudinal outcomes, CDRSB and NTB memory domain, were strongly associated with the risk of dementia diagnosis. For CDRSB we observed a statistically significant intervention effect on the longitudinal trajectory. Secondly, for NTB memory domain we investigated the type of association between the longitudinal profiles and the risk of dementia diagnosis. Specifically, we investigated three association types: the current value, the current value in combination with the rate of increase and the cumulative effect. We concluded that it was the current value in combination with the rate of increase of the longitudinal trajectory, that best captures the association with dementia diagnosis.
Additionally, this paper demonstrated the added value of a typical characteristic of joint models, namely the combination of the direct and indirect processes, both with different possibilities in modelling the effect of time. The joint model suggested an increased hazard ratio for the test versus the control group at the beginning of the trial. Given that there was no intervention before or at baseline, this increased hazard ratio was hypothesised to be caused by an imbalance between the intervention groups in characteristics at baseline. The groups were found to have a statistically significant imbalance at baseline in MMSE. MMSE is known to reflect cognitive performance and having an imbalance in MMSE between the intervention groups at baseline suggested that the groups  despite the randomisation process  might have on average differed in where they were in the disease continuum at baseline. Including baseline MMSE in the joint model markedly decreased the hazard ratio at the beginning of the trial, which fits into the hypothesis that the increased hazard ratio at the beginning of the trial was caused by baseline imbalance.
The imbalance between the intervention groups in characteristics at baseline might have been composed of several factors for which baseline MMSE was only a proxy. However, including the baseline MMSE in the joint models provided a tool to disentangle baseline confounding from the intervention effect.
Further, this paper illustrated another positive feature of joint models which is to model intervention effects on the hazard ratio that are changing over time. In the standard Cox models, the intervention effect is assumed to be constant during the entire followup, an assumption that is often not biologically meaningful. Using a joint model, it is possible to model a timevarying intervention effect on the survival outcome by incorporating a timevarying intervention effect on the longitudinal marker. In this prodromal AD trial, the joint model revealed an indication of an increasing intervention effect over time, suggesting a decreased hazard ratio for the test group at the end of the 24month trial.
Using time to dementia diagnosis as an outcome measure within the limited timeframe of a clinical trial has practical issues which complicate its use. First, a large part of the diagnoses cluster around the study visits when cognitive testing is performed and progression to dementia is thus detected. As a consequence, a part of the observed event times is intervalcensored, although the statistical software used for the analyses, did not cover this type of censoring. Another aspect is that, the diagnosis represents a single time point when the disease is thought of as a process moving along a continuum. Time to dementia diagnosis provided therefore only a rough measure of disease progression. However, using the information on time to dementia diagnosis was found to have an added value, by applying a statistical approach that combines every patient’s moment of diagnosis with his or her longitudinal trajectory.
Conclusion
Joint models provide a valuable tool in the statistical analysis of clinical studies with longitudinal and survival data, such as in prodromal Alzheimer’s disease trials, and have several added values compared to separate analyses.
Availability of data and materials
The data are proprietary information of the LipiDiDiet clinical study group.
Abbreviations
 AD:

Alzheimer’s disease
 AUC:

Area under the curve
 CDRSB:

Clinical dementia rating sum of boxes
 CERAD:

Consortium to establish a registry for AD
 MMSE:

Mini–mental state examination
 NTB:

Neuropsychological test battery
References
Aisen PS, Cummings J, Jack CR, Morris JC, Sperling R, Frölich L, Jones RW, Dowsett SA, Matthews BR, Raskin J, et al.On the path to 2025: understanding the Alzheimer’s disease continuum. Alzheimer’s Res Ther. 2017; 9(1):60.
Petersen RC. Mild cognitive impairment as a diagnostic entity. J Intern Med. 2004; 256(3):183–94.
Vaughan RM, Coen RF, Kenny R, Lawlor BA. Semantic and phonemic verbal fluency discrepancy in mild cognitive impairment: Potential predictor of progression to alzheimer’s disease. J Am Geriatr Soc. 2018; 66(4):755–9.
Munro CE, Donovan NJ, Amariglio RE, Papp KV, Marshall GA, Rentz DM, PascualLeone A, Sperling RA, Locascio JJ, Vannini P. The impact of awareness of and concern about memory performance on the prediction of progression from mild cognitive impairment to alzheimer disease dementia. Am J Geriatr Psychiatr. 2018; 26(8):896–904.
Baldeiras I, Santana I, Leitão MJ, Gens H, Pascoal R, TábuasPereira M, BeatoCoelho J, Duro D, Almeida MR, Oliveira CR. Addition of the A β42/40 ratio to the cerebrospinal fluid biomarker profile increases the predictive value for underlying alzheimer’s disease dementia in mild cognitive impairment. Alzheimer’s Res Ther. 2018; 10(1):33.
Wulfsohn MS, Tsiatis AA. A joint model for survival and longitudinal data measured with error. Biometrics. 1997; 53(1):330–9.
Henderson R, Diggle P, Dobson A. Joint modelling of longitudinal measurements and event time data. Biostatistics. 2000; 1(4):465–80.
Tsiatis AA, Davidian M. Joint modeling of longitudinal and timetoevent data: an overview. Stat Sin. 2004; 14(3):809–34.
Rizopoulos D. Joint Models for Longitudinal and Timetoevent Data: With Applications in R. Boca Raton: Chapman and Hall/CRC; 2012.
Sweeting MJ, Thompson SG. Joint modelling of longitudinal and timetoevent data with application to predicting abdominal aortic aneurysm growth and rupture. Biom J. 2011; 53(5):750–63.
Prentice R. Covariate measurement errors and parameter estimation in a failure time regression model. Biometrika. 1982; 69(2):331–42.
Dafni UG, Tsiatis AA. Evaluating surrogate markers of clinical outcome when measured with error. Biometrics. 1998; 54(4):1445–62.
Tsiatis AA, Davidian M. A semiparametric estimator for the proportional hazards model with longitudinal covariates measured with error. Biometrika. 2001; 88(2):447–58.
Li K, Luo S. Functional joint model for longitudinal and timetoevent data: an application to Alzheimer’s disease. Stat Med. 2017; 36(22):3560–72.
Yu B, Ghosh P. Joint modeling for cognitive trajectory and risk of dementia in the presence of death. Biometrics. 2010; 66(1):294–300.
Dantan E, Joly P, Dartigues JF, JacqminGadda H. Joint model with latent state for longitudinal and multistate data. Biostatistics. 2011; 12(4):723–36.
ProustLima C, Dartigues JF, JacqminGadda H. Joint modeling of repeated multivariate cognitive measures and competing risks of dementia and death: a latent process and latent class approach. Stat Med. 2016; 35(3):382–98.
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.
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.
Crowther MJ, Lambert PC, Abrams KR. Adjusting for measurement error in baseline prognostic biomarkers included in a timetoevent analysis: a joint modelling approach. BMC Med Res Methodol. 2013; 13(1):146.
Andrinopoulou ER, Rizopoulos D, Jin R, Bogers AJ, Lesaffre E, Takkenberg JJ. An introduction to mixed models and joint modeling: analysis of valve function over time. Ann Thorac Surg. 2012; 93(6):1765–72.
Soininen H, Solomon A, Visser PJ, Hendrix SB, Blennow K, Kivipelto M, Hartmann T, Hallikainen I, Hallikainen M, Helisalmi S, et al.24month intervention with a specific multinutrient in people with prodromal alzheimer’s disease (lipiDiDiet): a randomised, doubleblind, controlled trial. Lancet Neurol. 2017; 16(12):965–75.
de Wilde MC, Vellas B, Girault E, Yavuz AC, Sijben JW. Lower brain and blood nutrient status in alzheimer’s disease: results from metaanalyses. Alzheimer’s Dement: Transl Res Clin Interv. 2017; 3(3):416–31.
McKhann GM, Knopman DS, Chertkow H, Hyman BT, Jack Jr CR, Kawas CH, Klunk WE, Koroshetz WJ, Manly JJ, Mayeux R, et al.The diagnosis of dementia due to alzheimer’s disease: Recommendations from the national institute on agingalzheimer’s association workgroups on diagnostic guidelines for Alzheimer’s disease. Alzheimer’s Dement. 2011; 7(3):263–9.
Verbeke G, Molenberghs G. Linear Mixed Models for Longitudinal Data. New York: Springer; 2000.
Fitzmaurice G, Davidian M, Verbeke G, Molenberghs G. Longitudinal Data Analysis. Boca Raton: CRC press; 2008.
Brown ER, Ibrahim JG, DeGruttola V. A flexible Bspline model for multiple longitudinal biomarkers and survival. Biometrics. 2005; 61(1):64–73.
Cox DR. Regression models and lifetables. J R Stat Soc Ser B Methodol. 1972; 34(2):187–202.
Klein JP, Moeschberger ML. Survival Analysis: Techniques for Censored and Truncated Data. New York: Springer; 2006.
Therneau TM, Grambsch PM. Modeling Survival Data: Extending the Cox Model. New York: Springer; 2013.
Tseng YK, Hsieh F, Wang JL. Joint modelling of accelerated failure time and longitudinal data. Biometrika. 2005; 92(3):587–603.
Rizopoulos D. JM: An R package for the joint modelling of longitudinal and timetoevent data. J Stat Softw Online. 2010; 35(9):1–33.
Acknowledgements
The authors thank the LipiDiDiet clinical study group for the valuable data that was used in this paper. The last author would like to acknowledge support by the Netherlands Organization for Scientific Research (VIDI grant number 016.146.301).
Funding
The LipiDiDiet study project was funded by the European Commission under the 7th framework programme of the European Union (grant agreement number 211696). The funding body had no role in the study design, data collection, data analysis, data interpretation, or in writing the manuscript.
Author information
Authors and Affiliations
Contributions
Authors TH and HS designed the trial. Author FO drafted the manuscript and analysed the results under supervision of SS and DR. Authors SS, DR, TH and AH provided critical input to the manuscript. All authors read and approved the final manuscript.
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
FO, SS and AH are employees of Danone Nutricia Research. HS and TH were supported by a grant from the European Commission for the LipiDiDiet study (FP7211696 LipiDiDiet). HS has served as advisory board member for ACImmune and MERCK. Institution, UEF, has received funding from Nutricia for extension studies of LipiDiDiet Trial (no personal payment). Author DR declares to have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file
Additional file 1
R code to fit joint models. (PDF 35 kb)
Rights and permissions
Open Access This 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.
About this article
Cite this article
van Oudenhoven, F.M., Swinkels, S.H., Hartmann, T. et al. Using joint models to disentangle intervention effect types and baseline confounding: an application within an intervention study in prodromal Alzheimer’s disease with Fortasyn Connect. BMC Med Res Methodol 19, 163 (2019). https://doi.org/10.1186/s128740190791z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s128740190791z