Skip to main content
medRxiv
  • Home
  • About
  • Submit
  • ALERTS / RSS
Advanced Search

Applying time series analyses on continuous accelerometry data – a clinical example in older adults with and without cognitive impairment

View ORCID ProfileTorsten Rackoll, Konrad Neumann, Sven Passmann, Ulrike Grittner, Nadine Külzow, Julia Ladenbauer, Agnes Flöel
doi: https://doi.org/10.1101/2020.03.24.20042226
Torsten Rackoll
1Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt – Universität zu Berlin, and Berlin Institute of Health, NeuroCure Excellence Cluster (NCRC), Berlin, Germany
2Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt – Universität zu Berlin, and Berlin Institute of Health, Center for Stroke Research Berlin (CSB), Berlin, Germany
3Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt – Universität zu Berlin, and Berlin Institute of Health, QUEST – Center for Transforming Biomedical Research, Berlin, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Torsten Rackoll
  • For correspondence: torsten.rackoll@charite.de
Konrad Neumann
4Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt – Universität zu Berlin, and Berlin Institute of Health, Institute of Biometry and Clinical Epidemiology, Berlin, Germany
5Berlin Institute of Health (BIH), Berlin, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Sven Passmann
1Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt – Universität zu Berlin, and Berlin Institute of Health, NeuroCure Excellence Cluster (NCRC), Berlin, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ulrike Grittner
4Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt – Universität zu Berlin, and Berlin Institute of Health, Institute of Biometry and Clinical Epidemiology, Berlin, Germany
5Berlin Institute of Health (BIH), Berlin, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Nadine Külzow
1Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt – Universität zu Berlin, and Berlin Institute of Health, NeuroCure Excellence Cluster (NCRC), Berlin, Germany
6Kliniken Beelitz GmbH, Beelitz-Heilstätten, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Julia Ladenbauer
7Department of Neurology, University Medicine Greifswald, Greifswald, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Agnes Flöel
2Charité – Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt – Universität zu Berlin, and Berlin Institute of Health, Center for Stroke Research Berlin (CSB), Berlin, Germany
7Department of Neurology, University Medicine Greifswald, Greifswald, Germany
8German Center for Neurodegenerative Diseases (DZNE), partner site Rostock/Greifswald, Germany
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Abstract
  • Full Text
  • Info/History
  • Metrics
  • Supplementary material
  • Data/Code
  • Preview PDF
Loading

Abstract

Introduction Current analysis approaches of accelerometry data use sum score measures which do not provide insight in activity patterns over 24 hours, and thus do not adequately depict circadian activity patterns. Here, we used a functional approach to analyze accelerometer data that models activity pattern and circadian rhythm. As a test case, we demonstrated its application in patients with mild cognitive impairment (MCI) and age-matched healthy older volunteers (HOV). Moreover, we assessed the impact of chronotype on distribution of activity data.

Methods Data of two studies were pooled for this analysis. Following baseline cognitive assessment participants were provided with accelerometers for seven consecutive days. A function on scalar regression (FoSR) approach was used to analyze 24 hours accelerometer data. In a second step, analyses were controlled for chronotype using the German version of the morningness-eveningness questionnaire (d-MEQ).

Results Information on 47 HOV (mean age 66 SD 6 years) and 13 patients with MCI (mean age 69, SD 8 years) were available for this analysis. MCI patients displayed slightly higher activity in the morning hours as compared to HOV (maximum relative activity at 7:35 am: 75.6%, 95% CI 2.6 to 200.4%, p = 0.031). After controlling for d-MEQ, disturbed activity patterns were found in MCI of intermediate or evening chronotype, compared to HOV, i.e., MCI presented with higher activities in the morning hours (peak at 8:40 am: 357.6%, 95% CI 92.9 to 985.1, p < 0.001) and early afternoon hours (peak at 1:40 pm: 401.8%, 95% CI 63.9 to 1436.4, p < 0.001).

Discussion Using a novel approach of FoSR, we found timeframes with higher activity levels in MCI patients compared to HOV which were not evident if sum scores of amount of activity were used. In addition, we found that previously described activity patterns as a function of chronotype swere altered in MCI patients, possibly indicating that changes in circadian rhythmicity in neurodegenerative disease are detectable using easy-to-administer accelerometry.

Clinical Trials Effects of Brain Stimulation During Nocturnal Sleep on Memory Consolidation in Patients With Mild Cognitive Impairments, https://clinicaltrials.gov/ct2/show/NCT01782391?term=NCT01782391&rank=1,

ClinicalTrial.gov identifier: NCT01782391

Effects of Brain Stimulation During a Daytime Nap on Memory Consolidation in Patients With Mild Cognitive Impairment,

https://clinicaltrials.gov/ct2/show/NCT01782365?term=NCT01782365&rank=1,

ClinicalTrial.gov identifier: NCT01782365

Introduction

Activity assessment in natural environment is a widely used approach to determine habitual differences between cohorts (1) and has been used in different disease populations (2,3) or as a predictor for overall mortality (4). Devices to detect rest-activity patterns of body movements include small inertial measurement units (IMUs). These IMUs assess electrical current caused, for example, by acceleration, which is a proxy measure of a body’s movement in space. These so-called accelerometers or actigraphs capture bodily movement with a high precision in time. Still, common accelerometer assessments report aggregated data on number of steps per day or on amount of time spent in high, low or no activity as well as number of bouts without activity (5). The problem of aggregated data is that it reduces acquired information, thereby losing the high dimensionality of continuous 24 hours recordings. Of note, a more holistic view on activity distributions in daily rhythms has been postulated as an emerging topic lately (6). Thus, more detailed assessments of biological rhythms are needed (7).

Measures of circadian organization of activity patterns also have severe drawbacks: They are mostly limited to the amplitude, mesor and acrophase (time of the peak within the fitted 24 h rhythm) of activity which highlights the extremes of daily activity but have a poor resolution in time (8). Changes in circadian activity responses might be subtle so more precise tools are needed to analyze biological rhythms data. Researchers from Johns-Hopkins University (9,10) and Columbia University (11) introduced the concept of functional data analysis (FDA) to accelerometry data. In a first methodological study Goldsmith et al (11) demonstrated that while standard procedures showed inconclusive results, FDA was able to demonstrate differences in distribution of activity pattern over time. This approach thus might be able to detect subtle changes in circadian organization of activity patterns in a variety of disease states, including incipient neurodegenerative disease.

The pathology of Alzheimer’s disease (AD) serves here as a promising clinical example as AD starts years or even decades before first clinical symptoms become apparent (12). Thus, early diagnosis is mandatory for early treatments in order to delay the progression of cognitive decline. Several markers, including disrupted sleep behavior have been described before as a potential screening target (13). Likewise a disruption of circadian rhythm has been proposed (14). Daily rest-activity patterns are hypothesized to be controlled by circadian functions (15,16) which might thus be well-suited to detect early changes in AD pathology. Moreover, such assessments with small IMUs would constitute non-invasive and low-cost approaches that may allow for screening of large cohorts. In addition, such tools could be employed for repeated follow-up in long-term preventive or therapeutic interventions like lifestyle changes or dietary supplementation.

Here, as a clinical example we applied a functional data analysis algorithm on healthy older volunteers (HOV) and patients with mild cognitive impairment (MCI) in a previously acquired data set to describe the circadian rhythm in both groups, and to test if their activity patterns differed throughout a 24 hours cycle. We hypothesized that the distribution of activity would differ with regard to timing and magnitude between HOV and MCI. In a second step, based on previous reports showing that circadian activity is highly influenced by chronotypes (17), we asked if and how chronotypes (morningness vs. eveningness-types) would modify activity patterns between HOV and MCI.

Methods

Data reported here were taken from baseline measurements of two intervention studies that assessed the effects of oscillatory tDCS during sleep (daytime nap vs night sleep) on cognitive performance (Effects of Brain Stimulation During a Daytime Nap on Memory Consolidation (study 1) & Effects of Brain Stimulation During Nocturnal Sleep on Memory Consolidation (study 2) in young and older healthy subjects and subjects with mild cognitive impairment (MCI)). These studies are described in detail elsewhere (18,19) and included comprehensive assessment of subjective (using sleep diary and sleep questionnaires) and objective sleep-wake behavior among others, as well as neuropsychological testing and structural imaging of the brain using magnetic resonance imaging (MRI). Objective sleep assessment included seven days of accelerometry. In the current analysis we used all available baseline data of patients with MCI and age-matched HOV.

Study approval

The studies were approved by the institutional review board of the Charité Universitätsmedizin Berlin, Germany (EA1_295_12 & EA1_028_12), and were conducted in accordance with the declaration of Helsinki (Version 2008). All participants gave written informed consent prior to participation, and received a small reimbursement for their time.

Participants

HOV (50 – 90 years) were recruited via advertisements in the local database of the Charité Universitätsmedizin Berlin, Germany. MCI patients (50-90 years) were referred to the study from the memory clinic of the Charité University Hospital. The original studies included also young healthy participants; for the current analyses we focused on HOV and MCI, given our aim to compare circadian organization of activity in older adults. Participants underwent a structured telephone interview to exclude the presence of manifest sleep disturbances, contraindications for MRI, and non-native German speakers. Out of 242 potential participants that had entered the study, 66 HOV (mean age 66, SD 6) and 31 MCI (mean age 70, SD 8) were then invited to the laboratory to determine study eligibility which included a clinical interview, neurological examination, structural MRI, and standardized cognitive testing. Inclusion and exclusion criteria for both groups are detailed in the supplementary. In short HOV participants had to show no signs of dementia, and no present episode of depression (monitored with the Beck’s Depression Inventory (BDI (20))) while MCI patients had to fulfill core clinical criteria for the diagnosis of MCI outlined by Petersen and others (21). Further, clinical assessment and structural MRI revealed no systemic or brain diseases accounting for declined cognition. Patients diagnosed with amnestic or amnestic plus MCI were included.

Neuropsychological Assessment

Both groups were assessed using a neuropsychological test battery addressing various cognitive functions to ensure that HOV performed within age and gender matched normal range. A detailed description of all domains and tests can be found in the supplementary material (Table S2).

In addition we determined the affective state using the Positive and Negative Affect Schedule (PANAS (22)), anxiety state using the State Trait Anxiety Inventory (STAI-G X1 (23)), self-reported well-being using the World Health Organization Quality of Life-questionnaire (WHOQoL (24)), and the ability to cope with stressful events with a 120-items questionnaire (Stressverarbeitungsfragebogen, SVF (25)) which measured the individual tendency of using specific strategies of stress coping in general (trait). The SVF differentiates 10 strategies aiming on reduction of stress summarized to the following three positive strategies: Positive 1 (Pos 1; including cognitive strategies like reevaluation), Positive 2 (Pos 2, involving distraction-like strategies) and Positive 3 (Pos 3, focusing on control of stress and stress reactions). Moreover, the SVF cover seven strategies which were merged to “negative strategies” (Neg) and are assumed to cause an increase in stress (e.g., rumination, resignation). Lastly, questionnaires were administered to assess subjective sleep habits with the German version of the Pittsburgh Sleep Quality Index (PSQI (26)), daytime sleepiness with the German version of the Epworth Sleepiness Scale (ESS (27)), and individual chronotypes with the German version of Morningness–Eveningness-Questionnaire (d-MEQ (28)).

Accelerometry

The ActiGraph GT3X+ (ActiGraph, Pensacola, FL, USA) was used in this study. It is able to assess acceleration in the vertical, antero-posterior and medio-lateral axes. It has shown high inter-instrument reliability (Intraclass correlation 0.97 (29)) and intra-instrument reliability within frequencies that are common in human activities, and is described as a reasonable tool for longitudinally measuring sleep (30). During the week following baseline cognitive assessment each participant wore a GT3X+ on the hip (figure S2) for continuous seven-day recordings to fully capture daily physical activity. Subjects were asked to wear the device at all times and just to remove it for any activity involving water (showering, swimming, etc). The devices were pre-programmed with default settings (30 Hz in three axis, with a fixed start and stop time). Data were downloaded in 60 sec epochs. Accelerometry data download and sum score descriptives were performed using ActiLife Software 6.8.2 (ActiGraph, Pensacola, FL, USA).

After download, data were first checked for its actual wear time. The default algorithm according to Choi (31) was used to identify non-wear time. Days of wear time of < 60 % were excluded from further analysis. Raw values of accelerometry assessment are changes in current which are conversed to activity counts (32). Activity counts are an abstract value, usually used in literature. Even though they constitute an abstract value derived from company owned algorithms, higher values represent more intense activity compared to lower values. Different cut-offs to classify activity are available. We picked 0 – 99 activity counts per minute as sedentary activity, 100 – 2019 as light activity, 2020 – 5998 as moderate and everything above 5999 as vigorous activity (33) as these cut-points derived from a large NHANES cohort of 6329 participants and thus were assumed to correctly reflect levels of activity.

Statistical Analyses

We applied a Function-on-Scalar-Regression (FoSR) approach that is in line with recent publications by Goldsmith et al (11) and the group of Biostatistics from Johns Hopkins University (9,10,34). Function on scalar in this study means that activity counts as main outcome are functions of time (response), while covariates, here MCI patients/HOV etc. are scalars. However, for our analysis we made two important modifications:

  • – The set of basis functions used in the regression model are discrete wavelet functions instead of cubic splines (7).

  • – In addition to the fixed effects, the model equation contains a random intercept to account for clustering of data within individuals.

Pre-processing steps for FoSR as well as all statistical analyses were performed using statistical software R (Version 3.2.1) and the R package ‘wavelets’ (version 0.3-0) (35).

For FDA activity count data needed to be aggregated to at least 5 min epochs according to Xiao et al (10) resulting in 288 count data sets per participant per day (24h=1440 minutes=288*5 minutes). Furthermore, we log-transformed activity counts as distribution of count data were skewed. For high resolution accelerometer data robustness of data is a problem as changes between zero activity and high activity peaks might be present within small time periods. Furthermore, rapid changes in activity might not reflect circadian rhythm activity patterns. Therefore, a smoothing algorithm is needed to distinguish between signal (circadian activity pattern) and noise (changes between zero activity and high activity peaks within small time periods). Cubic splines or wavelet transformation are both suggested for smoothing of this type of data. Discrete wavelet functions and cubic splines have both a localized support, meaning that each wavelet function / each piece of the cubic spline model is able to fit the activity data well in a particular time range and reduces noise in this time range. In contrast to cubic splines, wavelets are periodic basis functions (with a period of 24 hours) and are thus better suited for modelling circadian patterns. Here discrete wavelet transformation was used for smoothing activity counts.

Discrete Wavelet Analysis

The discrete wavelet transformation is a well-established method in the analysis of time series (c.f. (36) and (37) pp. 174-179). Different series of wavelets are in use for the analysis of time-series (Coiflets, Least Asymetric, Best Localized, Daubechies wavelets, all implemented in the R package ‘wavelets’). The full discrete wavelet transformation is a regular linear (matrix) transformation from a finite dimensional real vector space into a space of the same dimension. Here, we used incomplete wavelet transformation with a non-trivial kernel that maps the time series from a high dimensional space (dim=288) into a real space with far less dimensions (dim=18, 18 basis functions). The purpose of this approach is twofold:

  • – Smoothing the activity data (removing noise).

  • – Defining a set of basis functions which can be used in the function-on-scalar regression analysis for testing group/subgroup specific activity patterns.

The incomplete discrete wavelet transformation behaves like a low-pass filter, meaning that rapid changes between non-activity and high activity will be smoothed out. For our analysis, we chose Daubechies wavelets of length 10 (d10, which corresponds to a time range of 50 minutes for the first level wavelet functions), meaning that an activity pattern of 50 minutes is modelled by one function.

The degree of smoothing depends on the number of basis functions used, which results from the length of the functions. Shorter wavelet filters tend to give poor results in smoothing the data (not enough smoothing, too much noise) whereas wavelets with long filter length are not well localized (too much smoothing, no sufficient model fit with regard to pattern relevant changes in activity). The use of 18 basis functions seemed to be a good compromise between good smoothing results and preserving a sufficiently high resolution for detecting activity patterns within circadian rhythm. The incomplete wavelet transform in concise matrix notation is Embedded Image

x denotes the vector containing the raw data whereas xsm is the vector of smoothed data (the subscript sm stands for “smoothed”). The columns of the matrix Ylp are the 18 basis function of the incomplete wavelet transform. Since each basis function is a vector of length 288, the matrix Ylp has 288 rows and 18 columns. Finally, the subscript lp refers to the low pass property of the incomplete wavelet transform.

Function-on-Scalar-Regression model

The presentation of the technical details of Function-on-Scalar-Regression closely follows the presentation in the appendix of (36). The model equation of FoSR is largely analogous to the model equation of a linear mixed effects model. We included a random intercept in the model equation since some participants contributed several correlated 24 hours activity records to the study.

Let xi,j be the vector of a log(1+count) transformed 24 hour activity record for day j = 1,…, ji of participant i =1,…,N aggregated in 5 minutes time blocks. Since each volunteer wore the actigraph device between one and eight days, ji may range from 1 to 8. Hence, the model equation is Embedded Image

In equation [1] ξ ik enotes the value of the kth covariate of the ith participant and εi,j (t) are the normally distributed residual terms with expectation zero and common variance σ2. The random intercept terms bi for the participants are normally distributed with expectation zero and variance τ2. We expanded each (row) coefficient vectors βk(t) using the wavelet basis functions that are the columns of Ylp Embedded Image

Omitting the index t=0,…, 287, the vector β k has length 288 and γk is a row vector of length 18 (number of wavelet basis functions). Details of the estimation of model parameters are presented in the appendix.

The FoSR model provides for each covariate (e.g. a group indicator variable) absolute and relative differences of activity counts with 95% confidence intervals and p-value for each five minutes époque. We report periods where the covariate (e.g. an indicator variable defining group differences) is significant in all five minutes époques and report only the largest relative difference in percent with two-sided 95% CI and p-value. All analyses were controlled only for age in order not to overfit the model.

Addressing the problem of multiple testing

All p-values were Bonferroni adjusted by a factor of 18. The Bonferroni factor was chosen equal to the number of basis functions used in the FoSR model. Accordingly, the confidence level of the confidence bands was risen from 1-α to 1-α /18.

Secondary analysis

For our secondary analysis we focused on the impact of chronotypes on the circadian organization of activity within both groups. We therefore stratified the analysis by chronotype groups: morning (moderate and definite morning types) and intermediate/evening chronotypes (intermediate, moderate, definite evening types) according to d-MEQ. We compared activity patterns of MCI and HOV within each strata.

Results

97 HOV and MCI patients were available for inclusion (patient flow is displayed in Figure S1 of the appendix). Of these 97 participants, 37 had to be excluded either due to missing actigraph data, abnormal MRI findings that suggested a differential diagnosis (e. g., stroke), or due to insufficient wear time of the actigraph (19 HOV, 18 MCI patients). An additional four MCI patients had to be excluded due to incoherent zero values distributed over a 24 hours cycle. Thus 47 HOV (mean age 66, SD 6) and 13 MCIs (mean age 69, SD 8) remained for final analysis with 228 valid accelerometer days in HOV and 60 days in MCI respectively. We compared basic characteristics of included and excluded subjects but could not identify major differences.

Baseline characteristics

Information on baseline characteristics including sociodemographic characteristics and cognition can be found in Table 1. Overall, MCI patients differed from HOV with regard to age, (mean 69 vs. 65 years, SMD = 0.51), depressive symptoms (median 7 vs. 3 in BDI, SMD=0.7), cognitive functioning (mean MMSE 28 vs. 29, SMD=1.07) and in positive and negative stress coping strategies (lower values in MCI with SMD from 0.99 to 1.39). The d-MEQ displayed higher rates of morning types in MCI compared to OLD (moderate morning type: n=17 (37.0%) vs. n=6 (46.2%); definite morning type: n=6 (13.0%) vs. n=3 (23.1), SMD=0.53) while subjective sleep ratings did not show substantial differences (median 4 vs. 3 in PSQI, SMD=0.24).

View this table:
  • View inline
  • View popup
Table 1: Baseline characteristics and aggregated descriptive accelerometer data for HOV and MCI

Participants spent most of their day in a sedentary lifestyle or in light activity with only small amounts of time in higher activity levels (Table 1), lacking substantial differences between groups (SMD for all measures <0.5). Aggregated analyses of sleep assessed via accelerometry exhibited comparable characteristics between the two groups (most SMD 0.2 or lower) only for wake after sleep onset and average awakening MCIs had somewhat higher values (SMD=0.42/0.40).

Function-on-Scalar-Regression

We compared the distribution of activity in MCI patients with HOV over the averaged 24 hours day-night cycle using FoSR. Variability of activity levels was higher in MCI patients compared to HOV. Fig 1a displays the smoothed average absolute activity over the course of a 24 hours cycle. MCI showed higher activity levels throughout the waking hours with a slight decrease below levels of HOV in the evening hours. For direct comparison between groups we analyzed relative activity between MCI and HOV adjusted for age (Fig 1b). We observed higher relative activity in MCI patients compared to HOV in the morning hours between 7:25 and 8:00 am peaking at 7:35 am (75.6%, 95% CI 2.6 to 200.4%, p = 0.031). In addition, during the day between 2:35 and 5:05 pm slightly higher level of activity with a peak at 4:10 pm (48.5%, 95% CI −11.7 to 149.6%, p = 0.412) was detected in MCI patients compared to HOV. In the evening hours between 8:35 and 10:05 pm MCI patients showed a slightly relative lower activity with its lowest value at 9:30 pm compared to HOV (−34.1%, 95% CI −60.8 to 10.7, p = 0.294). Over the course of the sleeping hours (mean time to bed and time out of bed hours in each group are displayed in Table 1) activity patterns of MCI patients and HOV were comparable. In sum MCI patients were more active in the morning hours and activity in MCI had higher variance compared to HOV.

Figure 1:
  • Download figure
  • Open in new tab
Figure 1: Comparison of absolute and relative activity distribution between MCI and HOV

Distribution of a) mean absolute activity between MCI and HOV over time and b) mean relative activity between MCI and HOV over time with 95% Confidence Interval (CI) adjusted for age.

Chronotypes

According to the d-MEQ, the majority of participants in the HOV group were intermediate types (44%) while in MCI patients the majority was of moderate morning type (46%) (Table 1).

The HOV group showed the expected pattern of activity with regard to questionnaire-based chronotypes. Compared to moderate and definite morning types combined, the group of intermediate, moderate evening and definite evening types combined displayed higher activity levels in the evening hours between 8:25 pm and 0:20 am with a peak at 11:20 pm (144.4%, 95% CI 48.2 to 303.0, p < 0.001) and lower activity in the morning hours between 6:20 and 9:30 am with the lowest activity at 7:50 am (−60.1%, 95% CI −75.2 to −35.8, p < 0.001) (Fig 2a).

Figure 2:
  • Download figure
  • Open in new tab
Figure 2: Comparison of activity distribution for different chronotypes

Distribution of mean relative activity over time clustered for the different chronotypes (assessed with the D-MEQ) a) definite evening, moderate evening, and intermediate type vs moderate morning, and definite morning types in MCI, b) definite evening, moderate evening, and intermediate type vs moderate morning, and definite morning types in HOV, c) mean relative activity over time between MCI and HOV for clustered intermediate and evening types, and d) mean relative activity over time between MCI and HOV for morning types. Data are shown with 95% confidence intervals. All analyses are adjusted for age.

In the MCI group definite evening chronotypes were not present. We opposed moderate and definite morning types to clustered intermediate, and moderate evening types as done before (38). Over the course of the day only early afternoon hours between 1:05 and 2:25 pm displayed differences with higher levels of activity in evening and intermediate types peaking at 1:40 pm (401.8%, 95% CI 63.9 to 1436.4, p < 0.001) compared to morning types (Fig 2b).

Under the assumption that activity patterns are linked to chronotypes, we would expect no intergroup differences (between MCI and HOV) within specific chronotypes. However, since circadian rhythm might be disrupted in MCI participants, we hypothesized to find group differences within chronotype strata. Subsequently we evaluated intergroup comparison between MCI and HOV for morning chronotypes as well as for evening and intermediate chronotypes separately. While activity patterns in morning types were comparable between groups (Fig 2c) combined groups of intermediate and evening types showed disrupted activity patterns in the morning hours between 7:45 and 10:00 am with a higher peak at 8:40 am (357.6%, 95% CI 92.9 to 985.1, p < 0.001) in MCI vs. HOV and early afternoon hours between 1:00 and 2:25 pm with a peak at 2 pm (345.2%, 95% CI 87.7 to 955.8, p < 0.001) with higher activity levels in the MCI group (Fig 2d).

To conclude we identified disrupted activity patterns in MCI compared to HOV in the combined group of intermediate and evening chronotypes that were not discernable before this stratified analysis by chronotype.

Discussion

In this study our goal was to apply the previously published Function-on-Scalar-Regression (FoSR) algorithm for the analysis of accelerometry data to a clinical population to provide an example for the algorithm’s potential in observational or intervention studies. We analyzed data from two pooled, previously published trials with the commonly applied approach in accelerometry, i.e., accelerometer derived sum scores (5,39). Here, we did not find substantial differences in proportions of time spent in specific activity levels or in total volume of activity. Next, we administered FoSR to investigate the circadian organization of activity patterns in patients with MCI and HOV using time series analysis and found higher relative activity in the morning and early afternoon hours and lower activity levels in the evening hours in MCI patients compared to HOV. While the net difference in amount of activity is likely to be negligible, our approach allowed us to demonstrate different time-based distributions of activity between the two populations.

Addressing circadian organization of activity, several authors have used five or ten hours spent in lowest or highest activity to assess phase shifts in activity in patients with MCI or dementia (16,40,41). Overall, findings from these studies have not been unequivocal, which indicates that especially in the early stages of cognitive impairment changes are subtle. Musiek et al. (16) described in their study differential activity patterns distributed throughout the day with lower amplitude between lowest and highest amount of activity and an increased variability over the course of 24 hours in preclinical AD patients. The authors compared different methods to analyze accelerometer data but did not apply time sensitive functional analysis. In contrast, time sensitive functional analysis was central in our study to identify time zones during which MCI patients were more or less active compared to HOV.

In a second step we controlled our regression model for circadian predisposition. We found that MCI patients with either intermediate or evening chronotype showed higher amounts of activity in the early afternoon hours compared to morning types. Moreover, we demonstrated that distribution of activity during a 24-hours cycle was different in intermediate and evening types from activity patterns found in HOV in the present study, or previous reports (38). To the best of our knowledge altered activity patterns in patients with MCI mediated by chronotype have not been described before. Individuals with the eveningness type have been reported to show higher risk of depressive symptoms (42) or insomnia (43), both independent risk factors for cognitive decline. Moreover, a disruption of circadian activity synchrony has been demonstrated in mice models of AD (44,45) and been described in human AD cohorts (46,47). Still contribution of diurnal predisposition on health outcomes is a relatively understudied field especially in AD. Even though our study cannot determine if the differences found between HOV and patients with MCI would lead to earlier detection of MCI in screening of individuals in the community, but this issue can now be addressed in future large-scale studies.

Limitations

Several limitations should be considered when interpreting our findings. First, groups of HOV and MCI patients differed slightly with regard to age, which might introduce a bias due to altered sleep behavior in the process of aging. However, we believe that these differences did not substantially impact our main results, given that the analysis was controlled for age. Second, we did not adjust for sex in our regression model to avoid overfitting in our small sample. Future studies with larger cohorts should adjust for sex differences to address the issue of potential sex differences. Third, the sample of our MCI group was small. However, we used our data as a test case to demonstrate an analysis approach with potential in a wide range of future applications. The validity of this approach has to be confirmed or refuted in future larger trials. Fourth, we cannot determine from the present data if disturbed activity patterns in MCI patients with the eveningness chronotypes are causally related to the cognitive impairment. Also, morningness types did not show altered activity patterns. However, our results may now be used as a starting point to generate new hypotheses, giving the robustness of the results.

Future research

First, to better differentiate changes in activity patterns due to chronotypes replication in large longitudinal studies are needed to corroborate our results and confirm the FoSR algorithm as a clinical measure. Also assessment in different pathological and healthy conditions are warranted to extract specific alterations in activity patterns for a certain pathology like AD. Second, results derived from accelerometers with FoSR approach should be compared to wrist-worn accelerometers, to facilitate studies of larger cohorts with commercially available activity monitors. Third, to inform inferences on the underlying mechanisms, biological markers such as melatonin and cortisol secretion or body temperature should be determined to more accurately define circadian rhythms. Finally, it should be assessed if parameters described here are sensitive to changes over time in response to therapies.

Conclusion

Using a novel approach of FoSR, we were able to demonstrate improved precision in analysis of time based high dimensional activity data. We found that MCI patients exhibited increased relative activity in the morning hours compared to HOV. In addition, we found that previously described activity patterns as a function of chronotype were altered in MCI patients of the eveningness chronotype. The approach of FoSR to analyze activity patterns may constitute an important screening tool for incipient neurodegeneration in the population at large, particularly if accelerometers of commercial smart watches are found to provide with valid measures for FoSR.

Data Availability

Data and analysis scripts are going to be available with publication at https://doi.org/10.5281/zenodo.3718578

https://doi.org/10.5281/zenodo.3718578

Author contributions

Design of the original trials and securing of funding: AF, NK

Design of the analysis: TR, KN, UG

Data collection: SP, JL Data analysis: KN, TR

Statistical analysis: KN, TR

Paper writing: TR, KN, SP

All authors read and approved the final manuscript.

Conflict of interest

None declared.

Data sharing

Data and analysis scripts are available upon request at https://doi.org/10.5281/zenodo.3718578

Supporting information

S1 Table. Inclusion and exclusion criteria

S2 Table. Neuropsychological test battery

S1 Appendix. Sleep and activity assessment

S2 Appendix. Statistical appendix

S1 Figure. Flow Chart

S2 Figure. Positioning of actigraph device

Acknowledgements

This work was supported by grants from the Deutsche Forschungsgemeinschaft (Fl 379-10/1; Fl 379-11/1, and DFG-Exc 257). We thank Lena Reich for their help with data acquisition. Special thanks to Patrizia Müller and Madlee Einsiedler for data preparation and actigraphy sleep ratings. Additionally, we want to thank Luo Xiao, Jeff Goldsmith and Ciprian Crainiceanu for provision of their R-scripts and helpful comments on the implementation of functional data analysis in accelerometry data. Furthermore, we would like to thank all subjects for their participation.

List of Abbreviations

BDI
Becks Depression Inventory
CI
Confidence interval
d-MEQ
Morningness-Eveningness-Questionnaire (german version)
ESS
Epworth Sleepiness Scale
FDA
Functional data analysis
FoSR
Function-on-scalar-regression
MCI
Mild cognitive impairment
MMSE
Mini-Mental-State-Examination
MRI
Magnetic Resonance Imaging
HOV
Healthy older volunteers
PANAS
Positive and Negative Affect Schedule
PSQI
Pittsburgh Sleep Quality Index
SMD
Standardized mean difference
STAIG-X1
State-Trait Anxiety Inventory (german version)
SVF
Stress-Verarbeitungs Fragebogen (translated: stress coping questionnaire)
tDCS
transcranial direct current stimulation
WHOQoL
World Health Organization Quality of Life

References

  1. 1.↵
    Perry B, Herrington W, Goldsack JC, Grandinetti CA, Vasisht KP, Landray MJ, et al. Use of Mobile Devices to Measure Outcomes in Clinical Research, 2010–2016: A Systematic Literature Review. Digit Biomarkers. 2018;2(1):11–30.
    OpenUrl
  2. 2.↵
    García-Hermoso A, Notario-Pacheco B, Recio-Rodríguez JI, Martínez-Vizcaíno V, Rodrigo de Pablo E, Magdalena Belio JF, et al. Sedentary behaviour patterns and arterial stiffness in a Spanish adult population - The EVIDENT trial. Atherosclerosis. 2015;243(2):516–22.
    OpenUrl
  3. 3.↵
    Hamer M, Coombs N, Stamatakis E. Associations between objectively assessed and self-reported sedentary time with mental health in adults: An analysis of data from the health survey for England. BMJ Open. 2014;4(3):1–7.
    OpenUrlCrossRef
  4. 4.↵
    Diaz KM, Howard VJ, Hutto B, Colabianchi N, Vena JE, Safford MM, et al. Patterns of sedentary behavior and mortality in U.S. middle-aged and older adults a national cohort study. Ann Intern Med. 2017;167(7):465–75.
    OpenUrlCrossRefPubMed
  5. 5.↵
    Quante M, Kaplan ER, Rueschman M, Cailler M, Buxton OM, Redline S. Practical considerations in using accelerometers to assess physical activity, sedentary behavior, and sleep. Sleep Heal. 2015;1(4):275–84.
    OpenUrl
  6. 6.↵
    Rosenberger ME, Fulton JE, Buman MP, Troiano RP, Grandner MA, Buchner DM, et al. The 24-Hour Activity Cycle: A New Paradigm for Physical Activity. Vol. 51, Medicine and Science in Sports and Exercise. 2019. 454–464 p.
    OpenUrl
  7. 7.↵
    Leise TL. Analysis of Nonstationary Time Series for Biological Rhythms Research. J Biol Rhythms. 2017;32(3):187–94.
    OpenUrl
  8. 8.↵
    Refinetti R, Cornélissen G, Halberg F. Procedures for numerical analysis of circadian rhythms. Vol. 38, Biological Rhythm Research. 2007. 275–325 p.
    OpenUrlCrossRefPubMedWeb of Science
  9. 9.↵
    Bai J, Di C, Xiao L, Evenson KR, LaCroix AZ, Crainiceanu CM, et al. An activity index for raw accelerometry data and its comparison with other activity metrics. PLoS One. 2016;11(8):1–14.
    OpenUrlCrossRefPubMed
  10. 10.↵
    Xiao L, Huang L, Schrack JA, Ferrucci L, Zipunnikov V, Crainiceanu CM. Quantifying the lifetime circadian rhythm of physical activity: A covariate-dependent functional approach. Biostatistics. 2015;16(2):352–67.
    OpenUrlCrossRefPubMed
  11. 11.↵
    Goldsmith J, Liu X, Jacobson JS, Rundle A. New Insights into Activity Patterns in Children, Found Using Functional Data Analyses. Med Sci Sports Exerc. 2016;48(9):1723–9.
    OpenUrl
  12. 12.↵
    Reisa A. Sperling, Paul S. Aisen, Laurel A. Beckett, David A. Bennet SC. Toward defining the preclinical stages of alzheimer’s disease: Recommendations from the national institute on aging. Alzheimer’s Dement. 2011;7(3):280–92.
    OpenUrlCrossRefPubMedWeb of Science
  13. 13.↵
    Mander BA, Winer J, Jagust WJ, Walker MP. Sleep: A Novel Mechanistic Pathway, Biomarker, and Treatment Target in the Pathology of Alzheimer’s Disease? Trends Neurosci [Internet]. 2016;xx:1–15. Available from:http://dx.doi.org/10.1016/j.tins.2016.05.002
    OpenUrl
  14. 14.↵
    Leng Y, Musiek ES, Hu K, Cappuccio FP, Yaffe K. Association between circadian rhythms and neurodegenerative diseases. Lancet Neurol. 2019;18(3):307–18.
    OpenUrlCrossRef
  15. 15.↵
    Lim ASP, Yu L, Costa MD, Leurgans SE, Buchman AS, Bennett DA, et al. Increased fragmentation of rest-activity patterns is associated with a characteristic pattern of cognitive impairment in older individuals. Sleep. 2012;35(5):633–40B.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    Musiek ES, Bhimasani M, Zangrilli MA, Morris JC, Holtzman DM, Ju Y-ES. Circadian Rest-Activity Pattern Changes in Aging and Preclinical Alzheimer Disease. JAMA Neurol. 2018 May;75(5):582–90.
    OpenUrl
  17. 17.↵
    Gaggioni G, Maquet P, Schmidt C, Dijk DJ, Vandewalle G. Neuroimaging, cognition, light and circadian rhythms. Vol. 8, Frontiers in Systems Neuroscience. 2014. p. 1–12.
    OpenUrl
  18. 18.↵
    Ladenbauer J, Külzow N, Passmann S, Antonenko D, Grittner U, Tamm S, et al. Brain stimulation during an afternoon nap boosts slow oscillatory activity and memory consolidation in older adults. Neuroimage. 2016 Nov 15;142:311–23.
    OpenUrlCrossRefPubMed
  19. 19.↵
    Paßmann S, Külzow N, Ladenbauer J, Antonenko D, Grittner U, Tamm S, et al. Boosting Slow Oscillatory Activity Using tDCS during Early Nocturnal Slow Wave Sleep Does Not Improve Memory Consolidation in Healthy Older Adults. Brain Stimul. 2016;accepted.
  20. 20.↵
    Kühner C, Bürger C, Keller F, Hautzinger M. Reliabilität und Validität des revidierten Beck-Depressionsinventars (BDI-II). Nervenarzt. 2007 Jun 7;78(6):651–6.
    OpenUrlCrossRefPubMedWeb of Science
  21. 21.↵
    Petersen RC. Mild cognitive impairment as a clinical entity and treatment target. Arch Neurol. 2004;62(7):1160–3; discussion 1167.
    OpenUrl
  22. 22.↵
    Watson D, Clark LA, Tellegen A. Development and validation of brief measures of positive and negative affect: the PANAS scales. J Pers Soc Psychol. 1988;54(6):1063–70.
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.↵
    Spielberger C, Gorsuch R, Lushene R. Manual for State-Trait Anxiety Inventory. Palo Alto, CA: Consulting Psychologist Press; 1970.
  24. 24.↵
    Angermeyer MC, Reinhold K, Matschinger H. WHOQOL-100 und WHOQOL-BREF: Handbuch für die deutschsprachige Version der WHO Instrumente zur Erfassung von Lebensqualität. Göttingen: Hogrefe;
  25. 25.↵
    Erdmann G, Janke W. Stress, Stressverarbeitung und ihre Erfassung durch ein mehrdimensionales Testsystem. 4th ed. Göttingen (Germany): Hogrefe; 2008.
  26. 26.↵
    Buysse DJ, Reynolds CF, Monk TH, Berman SR, Kupfer DJ. The Pittsburgh sleep quality index: A new instrument for psychiatric practice and research. Psychiatry Res. 1989 Apr;28(2):193–213.
    OpenUrlCrossRefPubMedWeb of Science
  27. 27.↵
    Johns M. A new method for measuring daytime sleepiness: the Epworth sleepiness scale. Sleep. 1991;14(6):540–5.
    OpenUrlCrossRefPubMedWeb of Science
  28. 28.↵
    Griefahn B, Künemund P, Bröde P, Mehnert P. Zur Validität der deutschen Übersetzung des Morningness-Eveningness-Questionnaires von Horne und Östberg. Somnologie. 2001;5:71–80.
    OpenUrlCrossRef
  29. 29.↵
    Jarret H, Fitzgerald L, Routen AC. Inter-instrument reliability of the Actigraph GT3X+ Ambulatory Activity Monitor during free-living conditions in adults. J Phys Act Health. 2015;12(3):382–7.
    OpenUrlCrossRefPubMed
  30. 30.↵
    Marino M, Li Y, Rueschman MN, Winkelman JW, Ellenbogen JM, Solet JM, et al. Measuring Sleep: Accuracy, Sensitivity, and Specificity of Wrist Actigraphy Compared to Polysomnography. Sleep. 2013;36(11):1747–55.
    OpenUrlCrossRefPubMed
  31. 31.↵
    Choi L, Liu Z, Matthews, Charles E, Buchowski MS. Validation of Accelerometer Wear and Nonwear Time Classification Algorithm. Med Sci Sport Exerc. 2011;43(2):357–64.
    OpenUrlCrossRefPubMedWeb of Science
  32. 32.↵
    Fulk GD, Sazonov ES. Using Sensors to Measure Activity in People with Stroke. Top Stroke Rehabil. 2011;18(6):746–57.
    OpenUrlCrossRefPubMedWeb of Science
  33. 33.↵
    Troiano RP, Berrigan D, Dodd KW, Masse LC, Tilert T, Mcdowell M. Physical Activity in the United States measured by Accelerometer. Med Sci Sport Exerc. 2008;40(1):181–8.
    OpenUrlCrossRefPubMedWeb of Science
  34. 34.↵
    Crainiceanu CM, Staicu AM, Ray S, Punjabi N. Bootstrap-based inference on the difference in the means of two correlated functional processes. Stat Med. 2012;31(26):3223–40.
    OpenUrlCrossRefPubMed
  35. 35.↵
    Aldrich E. wavelets: Functions for Computing Wavelet Filters, Wavelet Transforms and Multiresolution Analyses. 2019; Available from: https://cran.r-project.org/package=wavelets
  36. 36.↵
    Percival DB, Walden AT. Wavelet Methods for Time Series Analysis [Internet]. Cambridge: Cambridge University Press; 2000. Available from: http://ebooks.cambridge.org/ref/id/CBO9780511841040
  37. 37.↵
    Hastie T, Tibshirani R, Friedman JH (Jerome H. The elements of statistical learningL: data mining, inference, and prediction. Second Edition. Springer Series in Statistics; 2009. 745 p.
  38. 38.↵
    Quante M, Mariani S, Weng J, Marinac CR, Kaplan ER, Rueschman M, et al. Zeitgebers and their association with rest-activity patterns. Chronobiol Int. 2019;36(2):203–13.
    OpenUrl
  39. 39.↵
    Matthews CE, Hagströmer M, M PD, Bowles HR. Best Practices for Using Physical Activity Monitors. Med Sci Sports Exerc. 2013;44(18):1–17.
    OpenUrl
  40. 40.↵
    Ortiz-Tudela E, Martinez-Nicolas A, Díaz-Mardomingo C, García-Herranz S, Pereda-Pérez I, Valencia A, et al. The Characterization of Biological Rhythms in Mild Cognitive Impairment. Biomed Res Int. 2014;2014:1–7.
  41. 41.↵
    Cochrane A, Robertson IH, Coogan AN. Association between circadian rhythms, sleep and cognitive impairment in healthy older adults: An actigraphic study. J Neural Transm. 2012 Oct 10;119(10):1233–9.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    Coleman MY, Cain SW. Eveningness is associated with greater subjective cognitive impairment in individuals with self-reported symptoms of unipolar depression. J Affect Disord. 2019 Sep 28;256:404–15.
    OpenUrl
  43. 43.↵
    Chan JWY, Lam SP, Li SX, Yu MWM, Chan NY, Zhang J, et al. Eveningness and Insomnia: Independent Risk Factors of Nonremission in Major Depressive Disorder. Sleep. 2014 May 1;37(5):911–7.
    OpenUrl
  44. 44.↵
    Stevanovic K, Yunus A, Joly-Amado A, Gordon M, Morgan D, Gulick D, et al. Disruption of normal circadian clock function in a mouse model of tauopathy. Exp Neurol. 2017;294:58–67.
    OpenUrlCrossRef
  45. 45.↵
    Oyegbami O, Collins HM, Pardon M-C, Ebling FJP, Heery DM, Moran PM. Abnormal Clock Gene Expression and Locomotor Activity Rhythms in Two Month-Old Female APPSwe/PS1dE9 Mice. Curr Alzheimer Res. 2017;14(8):850–60.
    OpenUrl
  46. 46.↵
    Lu Z, Harris TB, Shiroma EJ, Leung J, Kwok T. Patterns of Physical Activity and Sedentary Behavior for Older Adults with Alzheimer’s Disease, Mild Cognitive Impairment, and Cognitively Normal in Hong Kong. J Alzheimer’s Dis. 2018;66(4):1453–62.
    OpenUrl
  47. 47.↵
    1. Cermakian N
    Weissová K, Bartoš A, Sládek M, Nováková M, Sumová A. Moderate Changes in the Circadian System of Alzheimer’s Disease Patients Detected in Their Home Environment. Cermakian N, editor. PLoS One. 2016;11(1):e0146200.
    OpenUrl
Back to top
PreviousNext
Posted March 26, 2020.
Download PDF

Supplementary Material

Data/Code
Email

Thank you for your interest in spreading the word about medRxiv.

NOTE: Your email address is requested solely to identify you as the sender of this article.

Enter multiple addresses on separate lines or separate them with commas.
Applying time series analyses on continuous accelerometry data – a clinical example in older adults with and without cognitive impairment
(Your Name) has forwarded a page to you from medRxiv
(Your Name) thought you would like to see this page from the medRxiv website.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Applying time series analyses on continuous accelerometry data – a clinical example in older adults with and without cognitive impairment
Torsten Rackoll, Konrad Neumann, Sven Passmann, Ulrike Grittner, Nadine Külzow, Julia Ladenbauer, Agnes Flöel
medRxiv 2020.03.24.20042226; doi: https://doi.org/10.1101/2020.03.24.20042226
Reddit logo Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools
Applying time series analyses on continuous accelerometry data – a clinical example in older adults with and without cognitive impairment
Torsten Rackoll, Konrad Neumann, Sven Passmann, Ulrike Grittner, Nadine Külzow, Julia Ladenbauer, Agnes Flöel
medRxiv 2020.03.24.20042226; doi: https://doi.org/10.1101/2020.03.24.20042226

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Subject Area

  • Public and Global Health
Subject Areas
All Articles
  • Addiction Medicine (280)
  • Allergy and Immunology (579)
  • Anesthesia (139)
  • Cardiovascular Medicine (1943)
  • Dentistry and Oral Medicine (252)
  • Dermatology (184)
  • Emergency Medicine (333)
  • Endocrinology (including Diabetes Mellitus and Metabolic Disease) (697)
  • Epidemiology (11097)
  • Forensic Medicine (8)
  • Gastroenterology (622)
  • Genetic and Genomic Medicine (3164)
  • Geriatric Medicine (308)
  • Health Economics (561)
  • Health Informatics (2041)
  • Health Policy (862)
  • Health Systems and Quality Improvement (782)
  • Hematology (310)
  • HIV/AIDS (682)
  • Infectious Diseases (except HIV/AIDS) (12713)
  • Intensive Care and Critical Care Medicine (707)
  • Medical Education (317)
  • Medical Ethics (92)
  • Nephrology (334)
  • Neurology (2981)
  • Nursing (164)
  • Nutrition (463)
  • Obstetrics and Gynecology (589)
  • Occupational and Environmental Health (614)
  • Oncology (1549)
  • Ophthalmology (476)
  • Orthopedics (185)
  • Otolaryngology (266)
  • Pain Medicine (202)
  • Palliative Medicine (57)
  • Pathology (402)
  • Pediatrics (912)
  • Pharmacology and Therapeutics (381)
  • Primary Care Research (355)
  • Psychiatry and Clinical Psychology (2781)
  • Public and Global Health (5588)
  • Radiology and Imaging (1091)
  • Rehabilitation Medicine and Physical Therapy (631)
  • Respiratory Medicine (759)
  • Rheumatology (338)
  • Sexual and Reproductive Health (311)
  • Sports Medicine (289)
  • Surgery (343)
  • Toxicology (48)
  • Transplantation (159)
  • Urology (132)