SUMMARY
Background Dementia is defined by cognitive decline that affects functional status. Longitudinal ageing surveys often lack a clinical diagnosis of dementia though measure cognitive and function over time. We used unsupervised machine learning and longitudinal data to identify transition to probable dementia.
Methods Multiple Factor Analysis was applied to longitudinal function and cognitive data of 15,278 baseline participants (aged 50 years and more) from the Survey of Health, Ageing, and Retirement in Europe (SHARE) (waves 1, 2 and 4–7, between 2004 and 2017). Hierarchical Clustering on Principal Components discriminated three clusters at each wave. We estimated probable or “Likely Dementia” prevalence by sex and age, and assessed whether dementia risk factors increased the risk of being assigned probable dementia status using multistate models. Next, we compared the “Likely Dementia” cluster with self-reported dementia status and replicated our findings in the English Longitudinal Study of Ageing (ELSA) cohort (waves 1–9, between 2002 and 2019, 7,840 participants at baseline).
Findings Our algorithm identified a higher number of probable dementia cases compared with self-reported cases and showed good discriminative power across all waves (AUC ranged from 0.754 [0.722–0.787] to 0.830 [0.800–0.861]). “Likely Dementia” status was more prevalent in older people, displayed a 2:1 female/male ratio and was associated with nine factors that increased risk of transition to dementia: low education, hearing loss, hypertension, drinking, smoking, depression, social isolation, physical inactivity, diabetes, and obesity. Results were replicated in ELSA cohort with good accuracy.
Interpretation Machine learning clustering can be used to study dementia determinants and outcomes in longitudinal population ageing surveys in which dementia clinical diagnosis is lacking.
Funding French Institute for Public Health Research (IReSP), French National Institute for Health and Medical Research (Inserm), NeurATRIS Grant (ANR-11-INBS-0011), and Front-Cog University Research School (ANR-17-EUR-0017).
INTRODUCTION
Major neurocognitive disorder (MND), commonly known as dementia, is a clinical syndrome characterised by a decline in cognitive performance that compromises patient’s independence1. Repeated clinical visits and assessments reveal the progression from a healthy state to dementia. International diagnostic criteria are available to identify dementia cases. Yet, more than half of the cases in high income countries (HIC)2 and up to 90% in low and middle income countries (LMIC)3 remain undetected. For such, new methods are needed to identify dementia cases and to study dementia determinants at the level of countries or continents.
Several population-based surveys modelled on the United-States Health and Retirement Study (HRS) are conducted in multiple countries to study the impact of the transition between late-life work and retirement in aging people4. The “HRS family” studies offer the opportunity to compare ageing outcomes internationally5. Yet, in these and in many other surveys, clinical dementia status is either not available or only self-reported by participants or their families, which underestimates the number of cases.
In the absence of clinical diagnosis in population ageing surveys, unsupervised machine learning, generally used to discover clusters or patterns within datasets6, can identify probable dementia cases. In a previous work, we applied an unsupervised clustering method to cross-sectional data from HRS and Survey of Health, Ageing and Retirement in Europe (SHARE) to identify high likelihood of dementia7 based on variables related to demographics, comorbidities, functional status, mobility, cognition, and neuropsychiatric symptoms. However, applying this clustering method to cross-sectional data did not allow us to investigate longitudinal transition from normal to impaired functional status, or to assess risk factors associated with transition to dementia status.
Herein, we built a clustering analysis for identifying transition to high likelihood of dementia in population ageing surveys using repeated measurements of cognition and functional status with a modified unsupervised machine learning algorithm. Our objectives were to demonstrate that this method can identify probable dementia in population aging surveys where dementia is either poorly or non-diagnosed, and that this method is also efficient to study dementia risk factors. Three analyses were used to ascertain the internal validity of “Likely Dementia” status: (1) we compared “Likely Dementia” identification with self-reported dementia, (2) we studied the prevalence of “Likely Dementia” status according to sex and age, (3) we tested whether traditional dementia risk factors were associated with a higher risk of transition to “Likely Dementia” cluster. To demonstrate replicability, we conducted our study using SHARE survey and replicated it in the English Longitudinal Study of Ageing (ELSA).
METHODS
Study design and participants
We used the harmonised dataset provided by the Gateway to Global Aging5 of SHARE, a longitudinal panel study across multiple countries in Europe and Israel8. This population survey takes place every two years and follows a representative sample of individuals aged 50 years or older from each participating country. The harmonised version of SHARE consists of seven waves so far (the third being retrospective) conducted between 2004 and 2017. We included subjects from countries who have participated in SHARE since the first wave (ie, Austria, Belgium, Denmark, France, Germany, Greece, Israel, Italy, The Netherlands, Spain, Sweden, and Switzerland), aged 50 years or older with consecutive follow-ups.
Selected variables
We used variables related to cognition and function to remain close to the DSM-5 criteria of MND. The selected variables are listed in the supplementary tables (Supplementary Table 1 & 2). All variables with more than 30% missing values were discarded and the remaining data were imputed using the imputeMFA command of the missMDA R package9.
Clustering
We ran Multiple Factor Analysis (MFA) followed by Hierarchical Clustering on Principal Components (HCPC) using FactoMineR R package10 and longitudinal data from all waves at the same time. MFA is a principal component method that balances for differences in the number of active variables per domain by forming active groups (procedure details are shown in supplementary data). For the clustering, we retained only active groups that represented participants’ function or cognition (see supplementary tables). Each participant, at each wave, was assigned to one of the three possible clusters (ie each participant could transition from one cluster to another, from one wave to another longitudinally). The number of clusters was set at three based on previous work for identification of high likelihood of dementia7. At the first wave, we singled out a cluster with a high probability of dementia (named “Likely dementia” cluster) based on the impaired cognition and function detected in its participants. Any participant classified in “Likely dementia” cluster was permanently assigned to it (ie, making any incident case a prevalent one).
We took into account the attrition due to study dropout and death across waves. We applied Inverse Probability Weighting (IPW) using the ipw R package11. For each wave, a logistic regression model was built based on the participants’ age, sex, and country of origin characteristics collected at the previous wave. Weights were derived by inverting the product of the predicted probabilities computed by the model. We integrated those weights into both imputation and clustering methods.
Self-reported diagnosis of dementia
We evaluated the discrimination power of our clustering method counting on its identification of “Likely dementia” status compared with the self-reported dementia status, data collected from the second wave of SHARE, using Sensitivity, Specificity and Area Under the Curve (AUC) metrics.
Effect of age, sex, and risk factors for dementia
We computed the prevalence of “Likely dementia” status of each wave by sex and by age. Participants were divided into six age groups (under 65 years, 65–69 years, 70–74 years, 75–79 years, 80–85 years, and more than 85 years).
We examined the role of several modifiable risk factors12 in transitioning to “Likely dementia” cluster: low education, hearing loss, hypertension, excessive alcohol drinking, current smoking, depression, social isolation, physical inactivity, diabetes, obesity, and air pollution. Past history of traumatic brain injury was not available in the database and could not be tested. All risk factors were measured at baseline and were imputed if necessary.
We dichotomised all ordinal risk factors variables. Education level was categorised as high (upper secondary and vocational training or tertiary education) or low (less than upper secondary). For hearing loss, we used self-reported hearing capacity as a proxy considering it either being normal (excellent, very good, and good) or bad (fair or poor). Moderate and vigorous physical activity were merged into physically active (frequency: more than once per week, once per week, one to three times a month) or inactive (hardly ever or never). The remaining risk factors were treated as dichotomous as they were in the database: hypertension (ever had high blood pressure vs never had high blood pressure), drinking (21 units or more of alcohol per week vs less than 21 units of alcohol per week), smoking (current smoker vs non-current smoker), depression (Centre for Epidemiologic Studies Depression [CES-D] scale score greater than or equal to five vs CES-D scale score less than five), social isolation (participating in social activities weekly vs non-participating in social activities weekly), diabetes (ever had diabetes vs never had diabetes), obesity (Body Mass Index [BMI] ≥ 30 kg/m2 vs BMI < 30 kg/m2), air pollution (living in urban area vs living in rural area).
Multistate models
In each wave, a participant could be classified in one of the three clusters (Cluster 1, Cluster 2 or Cluster 3; see above). Data being interval-censored, we applied multistate models using MSM package13 to study the impact of dementia risk factors on the risk of transition to “Likely Dementia” cluster.
We used age as the time-scale. It was calculated as the difference between birth date and interview date in years. In the multistate models, age was divided by 10 to facilitate the computational process without altering the Hazard Ratios (HR) results. Sex was considered as binary (male or female). All transitions were adjusted for sex, and transition towards “Likely Dementia” cluster was further adjusted for age. All covariates were set at baseline. For each risk factor, we computed its corresponding HR.
We checked the robustness of the multistate models in two steps. First, we considered death as a competing risk and added it as an absorbing state in the models. This was investigated in SHARE where vital status was reported consistently. Second, we excluded early prevalent and incident dementia cases by excluding participants categorised with a likelihood of dementia at first and second waves and ran multistate analyses again.
Replication cohort
In order to confirm our results, we chose the harmonised version of ELSA14 as a replication cohort. It is a representative longitudinal panel study of people aged 50 years and over in England. It comprises nine waves ranging from 2002 to 2019.
Ethical approval and guidelines
All participants provided informed consent and both studies obtained ethical approvals from local research committees. We followed both STROBE (STrengthening the Reporting of OBservational studies in Epidemiology) and MELODEM (The MEthods in LOngitudinal research on DEMentia) guidelines15,16 for the reporting of this study.
Role of the funding source
Sponsors of the study had no role in study design, data collection, data analysis, data interpretation, or writing of the report.
RESULTS
Identification of probable dementia
Of the initial sample of SHARE (n=30,419), we restricted our analyses to participants aged 50 years and over at baseline (n=29,102), who had consecutive follow-ups (n=15,278) (Figure 1). After running the clustering, the distribution between the clusters was uneven. At baseline, the first cluster (n=11,369) and the second (n=1,294) encompassed the majority of the sample, leaving a small part for the third cluster (n=535) (Table 1). Participants of the first and second clusters had similar baseline characteristics evoking healthy ageing. Participants of the third cluster were older (mean age 76.5 years [SD 11]), often female (n=368 [68.6%]), had lower education level (n=426 [79.6%] attained less than upper secondary education), more mobility impairment (mean mobility impairment score 4.9 [SD 1.5]), more functional impairment (mean Activities of Daily Living [ADL] score 3.1 [SD 1.7] and mean Instrumental Activities of Daily Living [IADL] score 4.2 [SD 1.8]), and more impaired cognition (mean immediate word recall test 2.6 [SD 1.9] and mean verbal fluency 10.4 [SD 6]) than participants of the first and second clusters at baseline. These characteristics corroborated that the third cluster was the one reflecting a high likelihood of dementia, thus we named it “Likely Dementia” cluster. Conversely, the first and second clusters were composed of participants considered dementia-free.
(*) Values were imputed using MissMDA package.
Discrimination power
We compared our algorithm identification with the self-reported dementia diagnosis in the SHARE dataset, which was available from wave 2 (Table 2). Our clustering algorithm allowed the identification of a higher number of “Likely Dementia” cases compared with self-reported dementia cases. The AUC metric ranged from 0.754 (0.722– 0.787) to 0.830 (0.800–0.861), suggesting good discrimination power. Sensitivity peaked at wave 4 reaching 0.714 (0.659–0.770) then slowly decreased after. Specificity remained high (> 0.9) in all waves. Results by country are given in supplementary data (Supplementary Table 1).
Abbreviations: AUC, Area Under the Curve; CI, confidence interval; NA, not available
Effect of age and sex
Older age and being female were both associated with an increased risk of entering “Likely Dementia” cluster. The prevalence of “Likely Dementia” was higher in women with approximatively a 2:1 female to male ratio across all waves (Fig 2.A). The number of “Likely Dementia” cases increased with age (Fig 2.B). For instance, at wave 2, the prevalence of “Likely Dementia” cases gradually rose with age: 1.8% in those under 65 years, 3.1% in 65–69 years, 5.9% in 70–74 years, 10.2% in 75–79 years, 18.9% in 80–85 years, and 37.4% in more than 85 years old participants.
Multistate models
To assess the associations of dementia risk factors with the risk of transitioning to ‘‘Likely Dementia’’ cluster (Table 3), we computed a multistate model (Figure 3.A). Nine of the eleven dementia risk factors chosen a priori were associated with an increased risk of transition from cluster 1 to “Likely Dementia” cluster: low education level (Hazard Ratio [HR] 1. 92% CI [1.58−2.33]), poor hearing (1.74 [1.45−2.09]), hypertension (1.35 [1.14−1.16]), smoking (1.45 [1.13−1.87]), depression (2.51 [1.06−3.07]), social isolation (1.66 [1.39−1.98]), physical inactivity (3.66 [2.97−4.51]), diabetes (2.4 [1.94−2.96]), and obesity (1.7 [1.39−2.07]). Some of these associations were also significant for transition from cluster 2 to “Likely Dementia” cluster: depression (2.39 [1.62−3.53]), social isolation (2.31 [1.51−3.53]), physical inactivity (3.21 [2.12−4.87]), and obesity (1.58 [1.08−2.32]).
Analyses using age as time-scale. All transitions were adjusted for sex. Transition towards the third cluster (“Likely dementia”) was further adjusted for age and each risk factor individually. All risk factors were taken at baseline. Main analysis was based on a multistate model (Model 1). Sensitivity analyses were based on a multistate survival model with death as an absorbing state. First, 105 participants were removed because of inconsistencies of dates (Model 2). Second, cases identified either at the first or the second wave were removed (Model 3). Abbreviations: HR, hazard ratio; CI, confidence interval.
In the first sensitivity analysis which took into account death (Figure 3.B), we excluded 105 participants due to inconsistencies between interview and death dates. All of the above-described associations between dementia risk factors and transition to “Likely dementia” cluster remained significant albeit with lower HR, except for hypertension. Of more, smoking became significantly associated with the risk of transition from cluster 2 to “Likely Dementia” cluster (2.23 [1.57−3.16]). In the second sensitivity analysis, where prevalent and incident cases at wave 1 (2004−05) and wave 2 (2006−07) (n=1,000) were further removed, HRs of transition from cluster 1 to ‘‘Likely Dementia’’ cluster did not change, but excessive alcohol drinking became a significant risk factor (1.34 [1.17−1.53]). As for transitions from cluster 2 to ‘‘Likely Dementia’’ cluster, only smoking (2.68 [1.79−4.03]) and depression (1.65 [1.07−2.54]) remained significant.
Replication in ELSA
Of the initial sample of ELSA (n=12,099), we restricted our analyses to participants over 50 years at baseline (n=11,522) and further restricted to participants who had consecutive follow-ups (n=7,840) (Supplementary Figure 2). Overall, results obtained with ELSA participants were similar to those found in the SHARE cohort (Supplementary Table 4).
At baseline, participants of the third cluster (n=659) were more likely older (mean age 69.8 [SD 11]), more likely female (n=401 [60.8%]), of lower education level (n=423 [64.2%] attained less than upper secondary education), had more mobility impairments (mean mobility impairment score 4.9 [SD 1.5]), more functional impairment (mean ADL score 2.7 [SD 1.5] and mean IADL score 2.5 [SD 1.4]) and worse cognition (mean immediate word recall test 4.6 [SD 1.9] and mean verbal fluency 16.5 [SD 6]) than the other clusters.
Our clustering algorithm identified a higher number of “Likely Dementia” cases compared with self-reported dementia cases. Except for wave 1 (2002−03) in which the number of self-reported dementia cases was small (n=24), the algorithm identification AUC metric values were similar to those found with SHARE (Supplementary Table 5). Sensitivity and specificity were balanced.
Women were more likely to be in the “Likely Dementia” group, and prevalence of “Likely Dementia” status rose with age (Supplementary Figure 3).
Ten dementia risk factors were tested (not air pollution due to missing urbanicity data). Their associations with transition to “Likely Dementia” cluster remained similar to those found with the SHARE dataset (Supplementary Table 6) except for excessive alcohol drinking which was protective for the transition from cluster 1 to “Likely Dementia” cluster (HR 0.6 [0.43−0.83]). Four risk factors were associated with an increased risk of transition from cluster 2 to “Likely Dementia” cluster: hypertension (1.64 [1.13−2.38]), depression (2 [1.26−3.17]), physical inactivity (2.69 [1.73−4.18]), and diabetes (2.23 [1.26−3.95]). We did not take death into account in the multistate models as vital status data were not available for each wave in sensitivity analysis.
Removing prevalent and incident cases at wave 1 (2002−03) and wave 2 (2004−05) in sensitivity analysis led to similar results with few exceptions. Excessive alcohol drinking was no longer significant for the transition from cluster 1 to ‘‘Likely Dementia’’ cluster (0.79 [0.58−1.08]). Only physical inactivity remained significant for the risk of transition from cluster 2 to ‘‘Likely Dementia’’ cluster (2.02 [1.1−3.69]).
DISCUSSION
Unsupervised clustering applied to two longitudinal population-based surveys of ageing (SHARE and ELSA) identified participants with high likelihood of dementia using longitudinal data related to functional and cognitive measures. In both surveys, this method had a good discrimination performance when compared with self-reported diagnosis of dementia. “Likely Dementia” status was more common in older participants and in women with a 2:1 sex ratio. Low education, hearing loss, hypertension, smoking, depression, social isolation, physical inactivity, diabetes, and obesity were associated with a higher risk of subsequent transition to “Likely dementia” cluster. Results for excessive alcohol drinking and air pollution were inconclusive. Applying clustering to longitudinal cohorts for the identification of high likelihood of dementia paves the way for researchers to conduct future secondary analyses on population ageing surveys worldwide.
Although supervised machine learning algorithms have already been used in population surveys to identify persons with dementia17, they have their limitations, eg, they require a subsample of data to be labelled “diagnosis of dementia”, and their external validity remains variable. Conversely, unsupervised machine learning may overcome such limitations as suggested in a previous cross-sectional study7. Here, we used an improved clustering method combining longitudinal data and a limited number of variables related to participants’ cognition and daily functions. Our clustering algorithm identified a greater number of people with a high likelihood of dementia in both SHARE and ELSA compared with self-reported dementia cases. Identifying a higher number of probable dementia cases in population ageing surveys might give a better statistical power to future studies of dementia determinants and outcomes. Moreover, this clustering method relies on cognitive and functional status data, largely available in HRS family studies and in several population ageing surveys, which makes it very suitable to apply to other ageing surveys including those in LMIC. Noteworthy, our study took into account many biases inherent to longitudinal studies, in particular attrition18 due to loss to follow-up or death. Internal validity was assessed using different approaches: comparison with self-reported diagnosis of dementia, impact of age and sex on dementia prevalence, and impact of known dementia risk factor on the risk of being classified as a “Likely Dementia” case. Results were obtained using data of 12 countries participating in SHARE, and then replicated in ELSA.
Nonetheless, our results should be carefully examined. Our algorithm detects a “Likely Dementia” status which cannot, by any stretch, be taken as a diagnosis of the disease without clinical validation. Future studies that compare our identification method with the recently developed cognitive assessment in HRS family cohorts using the Harmonized Cognitive Assessment Protocol (HCAP)19 are warranted. Our method cannot distinguish the aetiology of dementia, whether Alzheimer’s disease (AD) or other origins. Contrary to the results of our prior cross-sectional study, Cluster 1 and Cluster 2 participants were similar in terms of daily function, cognition, and mobility, yet they differed in their risk of transition to Cluster 3. However, we cannot rule out the possibility that the non-significant HRs observed for the transition from Cluster 2 to “Likely Dementia” cluster resulted from a lack of statistical power. Although this three-cluster partition remains consistent with our earlier work7, future investigation will test the interest of further simplification by merging the first two clusters together. The lack of biological or imaging biomarkers in this study could be seen as a limitation. Yet, biomarkers are often costly, human expert-dependent and rarely available in large population ageing studies. As for genetics, Apolipoprotein E (APOE)20 and polygenic scores21 are associated with a higher risk of AD, but the role of genetic factors in explaining future risk of dementia remains modest21,22. Although most dementia risk factors were associated with a higher risk of being assigned a “Likely Dementia” status, results for excessive alcohol drinking were ambiguous. We observed a deleterious drinking effect in SHARE, whereas it was protective in ELSA. Excessive drinking has been entangled for the brain damage it causes23, yet its exact relationship with dementia risk is debated since alcohol thresholds and time of exposure differ between studies24,25. The association between air pollution and dementia was inconclusive in SHARE and could not be explored in ELSA. Urbanicity (ie, geographical variation between urban and rural areas) was used as a proxy for air pollution as proposed recently12. Yet, people living in rural areas have shown higher rates of dementia compared with their urban counterparts 26,27. Switching to quantifiable pollution markers (fine particulate matter or ozone) that have been linked to an increased risk of dementia28 is more than desirable.
Unsupervised clustering is an efficient method to detect people with probable dementia in population ageing surveys using their cognitive and functional characteristics in a longitudinal setting. This approach opens new perspectives for the analyses of population data sets already available worldwide in HIC and LMIC to better compare and understand dementia determinants and outcomes.
Data Availability
ELSA raw data can be downloaded via the UK Data Service (https://ukdataservice.ac.uk/find-data/). SHARE raw data can be accessed on their website (http://www.share-project.org/data-access.html).
Contributors
AGM and LCL contributed to the conceptualisation and design of the study. AGM analysed the data. FH granted statistical methods expertise. HV helped with the computational analyses. AGM and LCL wrote the first draft. All authors reviewed and approved the final version of the manuscript. AGM decided to submit this manuscript for publication.
Declaration of interests
EB is an advisory board member of SafelyYou. EB has vested stock options in SafelyYou stock. KF is a participant on a data safety monitoring board Eli Lilly and NIH-sponsored trials. KF is a board member of Alector. ACBL has received academic grants from the French National Research Agency (ANR-17-EURE-0017) and NeurATRIS (ANR-11-INBS-0011). ACBL has received fees for consultancy work from Roche. All other authors declare no competing interests.
Data sharing
All the data we used are publicly available. We used the harmonised versions of SHARE and ELSA provided by the Gateway to Global Aging Data (https://g2aging.org/). ELSA raw data can be downloaded via the UK Data Service (https://ukdataservice.ac.uk/find-data/). SHARE raw data can be accessed on their website (http://www.share-project.org/data-access.html).
Supplementary documents
Multiple Factor Analysis (MFA) is a principal component method similar to Principal Component Analysis (PCA) but, compared to PCA, it balances the influence of groups of variables. In our case, it balances the influence of two components, participants’ function and cognition, which have equal importance in the clinical definition of dementia (data structure for SHARE is given in Supplementary Figure 1). Then, as PCA, it gives a representation of individuals in such a way that individuals are close on the representation if they have close values from the point of view of all the variables of all the groups. Subsequently, clustering can be performed on the MFA results. Full description of the variables used for SHARE and ELSA studies are available in Supplementary Table 1 and 2, respectively.
Each participant (from 1 to N) was seen multiple times (at least twice). Variables informative of functional status (from 1 to K) and variables informative of cognition (from 1 to J) formed the two active groups in the Multiple Factor Analysis (MFA).
*Values obtained using bootstrapping
Abbreviations: AUC, Area under the Curve; CI, confidence interval.
Analyses using age as time-scale. All transitions were adjusted for sex. Transition towards the third cluster (“Likely dementia”) was further adjusted for age and each risk factor individually. All risk factors were taken at baseline. Main analysis was based on a multistate model (Model 1). In sensitivity analysis, cases identified either at the first or the second wave were removed (Model 2). Abbreviations: HR, hazard ratio; CI, confidence interval; NA, not available
Acknowledgements
The Gateway to Global Aging Data project is developed by Centre for Economic and Social Research (CESR) at University of Southern California (https://cesr.usc.edu/). The project is funded by National Institute on Aging, National Institutes of Health (R01 AG030153, RC2 AG036619, R03 AG043052, R24 AG048024)
SHARE was funded by the European Commission, through FP5 (QLK6-CT-2001-00360), FP6 (SHARE-I3: RII-CT-2006-062193, COMPARE: CIT5-CT-2005-028857, SHARELIFE: CIT4-CT-2006-028812), FP7 (SHARE-PREP: GA 211909, SHARE-LEAP: GA 227822, SHARE M4: GA 261982, DASISH: GA 283646), and Horizon 2020 (SHARE-DEV3: GA 676536, SHARECOHESION: GA 870628, SERISS: GA 654221, SSHOC: GA 823782), and by DG Employment, Social Affairs and Inclusion (VS 2015/0195, VS 2016/0135, VS 2018/0285, VS 2019/0332 and VS 2020/0313). Additional funding from the German Ministry of Education and Research, the Max Planck Society for the Advancement of Science, the US National Institute on Aging (U01_AG09740-13S2, P01_AG005842, P01_AG08291, P30_AG12815, R21_AG025169, Y1-AG-4553-01, IAG_BSR06-11, OGHA_04-064, HHSN271201300071C, and RAG052527A) and from various national funding sources is gratefully acknowledged.
ELSA was funded by the National Institute on Aging (R01AG017644) and by a consortium of UK government departments: Department for Health and Social Care; Department for Transport; Department for Work and Pensions, which is coordinated by the National Institute for Health Research (NIHR, 198-1074). Funding has also been provided by the Economic and Social Research Council (ESRC).
First and foremost, we would like to express our gratitude to Julie Josse for all the methodological recommendations she gave for this project.
We would like to thank both Bioinformatics services of Institut Mondor de Recherche Biomédicale (IMRB) and Ecole Normale Supérieure (ENS) for their help in employing computer clusters for our analysis.