Skip to main content

Dynamical indicators in time series of healthcare expenditures predict mortality risk of older adults following spousal bereavement



The process of aging renders older people susceptible for adverse outcomes upon stress. Various indicators derived from complex systems theory have been proposed for quantifying resilience in living organisms, including humans. We investigated the ability of system-based indicators in capturing the dynamics of resilience in humans who suffer the adversity of spousal bereavement and tested their predictive power in mortality as a finite health transition.


Using longitudinal register data on weekly healthcare consumption of all Danish citizens over the age of 65 from January 1st, 2011, throughout December 31st, 2016, we performed statistical comparisons of the indicators ‘average’, ‘slope’, ‘mean squared error’, and ‘lag-1 autocorrelation’ one year before and after spousal bereavement, stratified for age and sex. The relation between levels of these indicators before bereavement and mortality hazards thereafter was determined by time to event analysis. We assessed the added value for mortality prediction via the time dependent area (AUC) under the receiver operating characteristic curve.


The study included 934,003 citizens of whom 51,890 experienced spousal bereavement and 2862 died in the first year thereafter. Healthcare consumption is increased, more volatile and accelerating with aging and in men compared to women (all p-values < 0.001). All dynamic indicators before bereavement were positively related with mortality hazards thereafter (all p-values < 0.001). The average discriminative performance for the 1-year mortality risk of the model with only age as a predictor (AUC: 68.9% and 70.2%) was significantly increased with the addition of dynamical indicators (78.5% and 82.4%) for males and females, respectively.


Dynamic indicators in time series of health care expenditures are strong predictors of mortality risk and could be part of predictive models for prognosis after life stressors, such as bereavement.

Peer Review reports


Resilience as a term has been used in a variety of different fields to describe a certain system’s ability to preserve specific characteristics or functions when being perturbed by different stressors [1]. Originally developed within the physical sciences, insights from the theory of complex dynamical systems have been extrapolated to understand and predict failure of organisms in the life sciences [2,3,4]. In human health, resilience to stressors can be interpreted as the potential of a person to bounce back into normal functioning after facing mental or physical disturbances [5,6,7]. This adaptive capacity of humans diminishes across time due to the process of aging, a phenomenon called frailty, describing the functional decline of various physiological systems rendering older people at an increased vulnerability to stressors [8].

A more formal application of complex dynamical systems theory to human aging has led to the development of indicators for quantifying and estimating resilience [9]. According to theory some systems may exhibit a phenomenon called ‘Critical Slowing Down’ (CSD) during which the system recovers in a slower rate from disturbances [9,10,11,12]. This decline of resilience can be manifested as the system approaches a tipping point transitioning from a stable steady state, e.g. health towards a contrasting one, e.g. disease [11, 13]. Using statistical terms, slower recovery rate can be reflected in high temporal autocorrelation, which measures the similarity of system’s dynamics in consecutive intervals (lags), accompanied with high variability, hence fluctuations, from the initial steady state. The aforestated statistical metrics can be found in the literature as Dynamical Indicators of Resilience (DIORs), since they are used to observe subtle changes in the dynamics of a given continuous (physiological) variable [9]. In general, the list of available, suggested DIORs is not restrictive, but ranges from plain univariate indicators such as autocorrelation, standard deviation, skewness and kurtosis [11, 14,15,16] to multivariate ones such as cross-correlation [17]. Empirical testing of such dynamical indicators in time-series of self-rated health, postural balance, heart rate and physical activity data has provided evidence for their usefulness in the gerontology field, by their ability to quantify physical resilience of older adults and predict their respective recovery potential from stressors in the future [6, 18,19,20,21,22]. The present study will empirically test the epidemiological potential of dynamical indicators in time-series of individuals’ healthcare expenditures available from the Danish (nationwide) registers, with the latter reflecting a person’s health status and its dynamics.

Spousal bereavement, i.e. death of a spouse, is a major stressor of health, and according to numerous scientific findings appears to be associated with increased hospitalization and mortality hazard [23,24,25,26,27,28,29,30,31]. The phenomenon in which bereaved individuals have higher mortality risk compared with their non-bereaved counterparts is termed as ‘widowhood effect’[32]. In addition to that, the effect of bereavement on mortality does not appear to behave similarly across different age groups and sexes, also being dependent on the time since the realization of the event [29,30,31, 33,34,35]. Apart from the factors above, there is further indication that resilience mechanisms to the stressor is not homogeneous across individuals [36,37,38,39].

Since spousal bereavement is one of the most stressful and inevitable encounters during one’s life course [40], there is a need for the creation of tools and computational models which quantify resilience at an individual level. Such tools can be then consequently utilized to extract personalized risks of adverse outcomes, i.e. death, after the stressor’s realization [41, 42]. The aforementioned addresses the gap of the existing literature to address who specifically is at high risk after experiencing bereavement. It needs to be determined whether such a DIORs setting can not only provide useful insights on how human resilience changes in response to bereavement, but also provide individualized risks of approaching a tipping point, i.e., death, after being perturbed by the loss of a loved one, given the resilience status before the event happening.

The aim of this paper is to calculate various metric-based indicators, i.e., DIORs, to capture changes in resilience in response to loss of one’s spouse, as well as function as predictors of individualized mortality risks after. To this end, we studied the association of spousal bereavement—as a standardized stressor—with dynamics of health status, as manifested in the healthcare consumption of Danish elderly over the age of 65. The idea is that healthcare consumption can serve as an indicator variable of a person’s health and resilience. We first aim to analyse the behaviour of these indicators, pre- and post- bereavement reflecting loss of resilience. Our second aim is to determine whether these indicators are predictors of failure by testing the hypothesis of them being associated with death, an ultimate ‘tipping point’ upon perturbation. The aforementioned would be fulfilled by observing significantly higher mortality in the groups of people with the highest levels of these indicators before bereavement.


Data processing

Data were collected from national population registers of Statistics Denmark ( containing information on the healthcare consumption of Danish Elderly over the age of 65, who had been resident in Denmark for at least 5 years before the date of death, from January 1st, 2011 to December 31st, 2016 [43, 44]. Healthcare consumption was divided into weekly measurements and was initially grouped into different kinds of costs, while being measured in thousands of Danish Kroners (1 unit = 1000 DKKR). Information on prescription drug, hospital (inpatient and outpatient), home care, residential care, and primary care costs was available and a new variable, aggregating all the types of healthcare costs, was constructed for each individual, Variables describing sex, date of birth, date of death and bereavement date (date when spouse died) were virtually complete.

Availability of data & materials

According to Danish law, scientific organizations can be authorized to work with data within Statistics Denmark and can provide access to individual scientists inside and outside of Denmark. Data are available via the Research Service Department at Statistics Denmark: ( for researchers who meet the criteria for access to confidential data.

Definition of DIORs

We computed for each individual the average (AVERAGE) healthcare consumption one year before and after bereavement as a measure of health status. Initial data exploration indicated that many individuals gradually increase their healthcare costs approximately linearly with time. Hence, we decided to compute the slope (SLOPE), i.e. the average linear temporal change of a person’s healthcare costs. Thus, the AVERAGE and SLOPE, serve as indicators which parameterize the overall trend of individuals, being reflected by an average level of expenditures and their respective linear rate of change. We also calculated Dynamical Indicators of Resilience (DIORs) or early warning signals which capture the individuals’ dispersion dynamics from their overall, health expenditures trend. For that purpose, we further computed the mean squared error (MSE) of the linear models regressing health costs on weeks, measuring the magnitude of dispersion and variability around that linear trajectory of costs. Higher values of MSE indicate larger fluctuations from that linear trend of healthcare patterns and may reflect large shifts of health status for the individual, indicative of critical slowing down. Finally, we evaluated the memory of those fluctuations across consecutive weeks, by calculating the lag-1 autocorrelation of the detrended expenditures (residuals) between subsequent moments (lagging by 1 unit of time). High levels of lag-1 autocorrelation are reflective of a stable fluctuating behavior, whereas low ones of a rapid changing one. We note that while there are many summary statistics which describe other aspects of the time-series distribution and are subject to empirical testing, in the current research we restrict to those depicting overall centrality and trend of costs (AVERAGE, SLOPE), as well as those reflecting the magnitude of variability (MSE) and its memory (lag-1 autocorrelation), in similar fashion with previous works on the literature [6, 18, 21, 22], for comparability purposes (See Appendix for the mathematical formulas of DIORs).

Data analysis

Comparisons of the DIORs one year before and after the spousal bereavement were performed for the whole sample, but also stratified based on sex and various age categories. The distributions of the weekly healthcare costs and the metrics describing them, were positively skewed (50% of the total cost variable had the value 0) and non-normally distributed. Hence, we assessed the significance of the statistical differences in DIORs between the two bereavement periods via non-parametric two-tailed Wilcoxon signed-rank tests.

Comparison analysis of the DIORs one year before bereavement, between individuals who survived within one year after the bereavement event with those who did not, was performed to evaluate the hypothesis of these metrics signaling early loss of health resilience. Statistically significant results were assessed with the use of two-tailed Mann–Whitney U tests.

We deployed quantile-median regression models [45] using each DIOR at a time as response (Y) variable and the age of the individuals as predictor (X), while adjusting for sex, to predict the change in the median of the DIORs in response to an increase in age. Therefore, we evaluated if those dynamic indicators follow the theory of aging, by observing the models’ slope-coefficient, with older individuals becoming frailer, thus exhibiting higher DIORs values. Ιt should be noted that the slopes of the median regression models are different from the SLOPE DIOR, with the former indicating the between-individuals variability and the latter the within-individual one. Sex-based statistical comparisons (two-tailed Mann–Whitney U tests), both pre- and post-bereavement, were furtherly performed to examine if loss of health resilience, thus frailty, is more present in one of the sexes.

To investigate if extreme DIORs of healthcare expenditures before bereavement are predictive of death after the stressor, we performed a time to event (mortality) analysis using Kaplan–Meier survival estimates and Cox regression stratified on sex, with a follow-up period up to one year after the bereavement with the outcome being death from any cause. Each DIOR was categorized into tertiles, with the reference category being the bottom tertile for each respective DIOR. Autocorrelation for those individuals with constantly zero healthcare costs across the follow-up could not be computed and were classified in the bottom tertile of the respective metric as this corresponds to our intuitive expectation that persons with zero healthcare costs are more resilient than persons with varying healthcare costs. We note that some individuals may have healthcare trajectories with many consecutive zeros (i.e., zero healthcare expenses for and extended time-period) leading to high autocorrelation measures, and that for these cases lag-1 autocorrelation is not an appropriate indicator of critical slowing down, as the observed “memory” in the data is not due to a slowing down per se, but rather due to lack of signal. This issue may lead to some noise in the intended use of the lag-1 autocorrelation measure (the intention being that high autocorrelation corresponds to slowing down phenomena indicating higher frailty), and we anticipate that this may cause decrease in the significance of statistical associations between the lag-1 autocorrelation and dependent variables.

To assess the added value of DIORs in the discrimination ability of a Cox regression model including only age to predict one-year mortality after spousal bereavement, we calculated the time-dependent area under the one-year receiver-operating-characteristic curve (AUC) [46,47,48]. Estimation of prediction error was performed using the leave-one-out-bootstrap method, on 50 bootstrap (drawn with replacement) samples [49]. Within a 1-year timespan, the AUC is the probability of an individual receiving a higher mortality risk, given that he/she died within that year, compared to an individual who survived in the same period. All the analyses were conducted using R version 4.0.3 [50] (See Appendix for the mathematical formulas of the models).


Population statistics

The general population consisted of 934,003 Danish citizens, with their age ranging from 65 to 105 years old. Weekly healthcare expenditures for each individual were available from January 1st, 2011, to December 31st, 2016. For the purposes of calculating their respective DIORs one year before and after bereavement, our analysis was restricted to 51,890 individuals who experienced a spousal bereavement within the 2012–2015 time period (5.55% of the initial sample). The final sample consisted of 17,813 males (34.3%) and 34,077 females (65.7%) with their median age (interquartile range) at start of the study being 76 (11) and 74 (10) years respectively. From the sample of 51,890 bereaved individuals, 2862 died in the first year thereafter, with more than 90% being followed for one year after bereavement. An illustration of the sample processing is depicted in Fig. 1.

Fig. 1
figure 1

Flowchart of the analysis of Danish healthcare register data. Upper box presents the initial sample of Danes over 65 years for the general population. Middle box explains the process from the original sample towards the final one containing only the spousal bereaved individuals suitable for the analysis. The box at the bottom illustrates the process of time to event analysis

The population average healthcare costs were investigated by a stratified-on sex analysis for the weekly average expenditures and residual time-series and portrayed in Fig. 2. Increasing trends are observed on the weekly average, manifesting the effect of aging in healthcare costs. Both the weekly average costs and residuals are increased after bereavement, a pattern that can be also observed in the sex and age stratified analysis (see Additional File 1; Supplementary Fig. 1 & 2).

Fig. 2
figure 2

Sex stratified weekly average and residuals of healthcare expenditures pre- and post-bereavement. Week 0 depicts the date of the standardized stressor of bereavement, splitting the time interval one year before and after bereavement. Linear regression lines are fitted both before and after the stressor, indicating the differences in the healthcare consumption pattern between the two periods. Average healthcare consumption after bereavement is higher than before, for both males and females. The average of weekly-based residuals also seems to be following the same pattern, that being increased after bereavement. Females seem to exhibit a more resilient behavior, with smaller deviations when compared with males, manifested in the weekly average and the residuals of their respective healthcare costs

DIORs across bereavement states

A comparison of the dynamical indicators of resilience (DIORs) pre- and post-bereavement is listed in Table 1. The average of healthcare expenditures (DKK * 1000) post-bereavement (median (IQR) = 0.299 (1.377)) was significantly higher (p-value < 0.001) than pre-bereavement (median (IQR) = 0.208 (0.697)). A comparison of slope and MSE of healthcare costs pre- and post-bereavement followed a similar pattern. Lag-1 autocorrelation of linearly detrended healthcare costs exhibited a different behaviour than the rest of the metrics, that being not significantly different (p-value = 0.065) after bereavement (median (IQR) = 0.019 (0.378)) when compared with before (median (IQR) = 0.017 (0.391).

Table 1 Median (IQR) of DIORs one year before and after spousal bereavement in subgroups

Statistical comparisons indicated that males manifest significantly higher (p-value < 0.001) DIORs than females (p-value < 0.001), both one year before and after bereavement (Additional File 1; Supplementary Table S1).

We applied quantile-median regression to verify the association between DIORs and aging. When the DIORs were regressed on age, all the extracted slope coefficients were positive and statistically significant (p-value < 0.001) (Additional File 1; Supplementary Table S2). We observed similar trends on the DIORs in both males and females, as well as within the strata of age (see Table 1).

Predictive value of DIORs before bereavement for death thereafter

In order to assess the hypothesis of DIORs before bereavement being quantitative indicators of the level of resilience of a system, we compared them across those bereaved individuals who survived during the first year after bereavement, i.e. ‘Survivors’ and those who did not, i.e. ‘Deceased’. The average was significantly higher (p-value < 0.001) for the latter group (median (IQR) = 1.825 (5.326)) compared with the former (median (IQR) = 0.191 (0.588)). The slope coefficient for the Deceased group (median (IQR) = 0.004 (0.062)) was also significantly higher (p-value < 0.001) than the Survivors group (median (IQR) = 0 (0.006)). The MSE exhibited a similar pattern, that being significantly higher (p-value < 0.001) for the Deceased (median (IQR) = 8.005 (31.389)) when compared with the Survivors (median (IQR) = 0.456 (3.434)). Lag-1 autocorrelation of detrended healthcare expenses was also found to be significantly higher (p-value < 0.001) for those who died earlier (median (IQR) = 0.356 (0.600)) compared with those who survived (median (IQR) = 0.008 (0.360)). Similar differences were found for both sexes separately, as well as in all age restricted strata (statistical-test results available in Additional File 1; Supplementary Table S3, S4 and S5).

Time to event analysis was performed to investigate the value of tertile-grouped DIORs before bereavement to predict death in the first year thereafter. The Kaplan–Meier curves of the analysis are presented in Fig. 3, depicting the cumulative percentage of deaths belonging in each tertile group of all the DIORs for the period of 52 weeks after bereavement. The cumulative percentage of deaths within one year after spousal bereavement was the highest for the group of people belonging in the top tertile of each DIOR. The pattern was observed both in males and females, with the latter exhibiting smaller percentage of deaths than the former. The middle and bottom tertile seem to have moderate to weak differences in the survival curves.

Fig. 3
figure 3

Sex stratified Kaplan–Meier curves of DIORs on cumulative deaths (%) up to one-year after bereavement. Time is measured in weeks with each respective number at risk in the tables below each curve. P-values indicate the significance of log-rank test for difference in survival between the tertiles. Shaded areas represent 95% CIs. First row depicts the results for males, whereas second row is for females. The colors depict the tertiles of each respective DIOR (red: top-3rd tertile, blue: middle -2nd tertile, green: bottom-1st tertile)

The estimated hazard ratios (HRs) of DIORs on mortality risk, along with their respective 95% confidence intervals (CIs), of the Cox regression adjusted for age, are reported in Table 2. When compared with the reference category (bottom tertile), the group of males in the top tertile of average appeared to have a 4.81-fold higher mortality risk [95% confidence interval [CI], 4.04–5.72]. The corresponding mortality risk was found to be 7.28-fold higher [95% CI, 6.03–8.78] for females. Similar risk differences were found for the slope coefficient and MSE. Although we did not necessarily expect autocorrelation to exhibit the same predictive power, a similar trend was observed in the top tertile group for lag-1 autocorrelation. Males exhibited a 2.86-fold higher mortality risk [95% CI, 2.44–3.34] while females’ risk being 4.65-fold higher [95% CI, 3.95–5.48], when compared with the bottom tertile. We further assessed the joint effect of extreme DIORs in mortality risk by creating a Score variable, counting the number of DIORs for which an individual belongs in the top tertile category. The 9% of individuals belonging in the top tertile category of all the DIORs had 14.59-fold [95% CI, 12.33–17.26] higher risk compared with those not belonging into any top tertile DIORs category (Additional File 1; Supplementary Table S6).

Table 2 Hazard ratios (95% CI) of all-cause mortality after bereavement dependent on various DIORs in tertiles

Discrimination ability of DIORs

We assessed the added value of DIORs to the discrimination ability, next to age, in predicting mortality within one year after spousal bereavement among all the individuals (N = 51,890) aged 65 years and older. Table 3 summarises the results of the AUC analysis stratified on sex. Age seems to be the main factor of discriminative value with an AUC of 68.9% [95% CI, 67.4%—70.3%] and 70.2% [95% CI, 68.8%—71.6%] for males and females, respectively. The addition of average to the model with age, increased AUC by 7.1% [95% CI, 6%—8.3%] and 8.9% [95% CI, 7.7%—10.1%] for males and females, respectively. The discriminative value, based on the AUC, was increased with the inclusion of slope by an extra 2.5% [95% CI, 1.8%—3.2%] for males and by 4.1% [95% CI, 3.2%—4.9%]. The MSE of healthcare consumption improved the discriminatory power (AUC) to predict all-cause mortality by an extra 8.1% [95% CI, 6.9%—9.2%] for males and 10% [95% CI, 8.8%—11.2%] for females. Lag-1 autocorrelation increased the discriminatory value of the model with age, by a further 3.4% [95% CI, 2.6% – 4.3%] in males and 6.5% [95% CI, 5.4% – 7.5%] in females.

Table 3 Discrimination ability to predict one-year mortality based on DIORs stratified by sex

The final model’s predictor variables consisted of age, average, slope coefficient, MSE and lag-1 autocorrelation. The discriminatory power according to the AUC was 78.5% [95% CI, 77.4%—79.6%] for males and 82.4% [95% CI, 81.5%—83.3%] for females. In regards to the model’s covariates, the average indicator was highly correlated with MSE (ρ = 0.79, Spearman’s p-value < 0.001), while the cross-correlations among the rest of the indicators were ranging from 0.07 to 0.33.


The current research focused on the perturbation of spousal bereavement and the ability of DIORs to capture transitions in human health in response to this perturbation. To this end, we analyzed the healthcare consumption of Danish citizens over the age of 65 by applying the logic of systems failure, as expressed in the dynamical indicators of resilience. Aggregation of all the different types of costs for each person merges all the discrete layers of healthcare together (from treatment to care), thus serving as a read out of their health status and level of resilience.

The analysis of bereaved individuals produced the following observations. First, all the indicators describing the healthcare expenditures data were significantly increasing with age and higher in males, suggesting lower resilience. In addition, the average, slope coefficient and MSE were significantly higher the year after bereavement when compared with the year before. The aforementioned is a strong indication that spousal bereavement seems to be a life-event which apart from causing intense physical and mental distress to peoples’ lives, poses the risk of increased, more volatile healthcare consumption and a rate of increase across time. Last, all the dynamical indicators were found to be signals of early loss of resilience before critical transitions in human health, i.e. -death. We inferred the latter by observing that the group of individuals with the highest indicators of resilience, before bereavement, exhibited a significantly higher mortality risk after bereavement when compared to those with low DIORs. Ultimately, we found evidence that these indicators are able to capture the dynamics of human resilience and predict mortality, indicating that with ageing comes frailty, with males being the frail sex.

Time to event analysis evaluated the hypothesis of DIORs predictive value. Individuals with high DIORs one year before bereavement were found to have significantly increased mortality risk for the year thereafter when compared with low DIORs group. These outcomes were consistent for each respective indicator, both for males and females. Therefore, this finding further enhances the hypothesis of these metrics being indicators of resilience in the life sciences, with higher DIORs indicating less resilient systems.

The AUC analysis indicated that all of the DIORs seem to increase the AUC when included in a model with age. We observed the highest increase in the AUC by adding the MSE, followed by the average, then lag-1 autocorrelation and last, the slope. Since our goal was to evaluate the predictive performance of DIORs on mortality, we implemented both the average and MSE in the final model, even though their correlation was high, following the recommendations of recently published work [51]. The final predictive model with average, MSE, lag-1 autocorrelation, slope, and age as predictors reached an AUC of 78.5% and 82.4% for males and females, respectively. Therefore, it appears that the various signals extracted from the dynamic behavior of health-related time-series are capable of discriminating mortality risks amongst older persons well, both for males and females. The AUC extracted from models of previous studies investigating the added value of various health indicators to mortality predictions in older adults ranged from 0.71 to 0.82 [48, 52]. Thereby, from the abundance of predictive models using various clinical and physiological hallmarks, ours using four dynamical indicators reached a score amongst the highest performing ones, adding increased value in the literature related to prediction of mortality risk at old age.

Analysis of time series data consisting of longitudinal measurements provides the opportunity to study the trajectory of an individual’s health across time, which can be reflected in several dynamical indicators, i.e., variability, rate of change and correlation, moving beyond the mean or just a snapshot value at a given time point. Therefore, we propose a methodological approach that utilizes not only static but also dynamical factors of human health, emphasizing the need for fine-grained, longitudinal measurements in datasets of epidemiological interest. That approach allows for the acquisition of a various, distinct signals which can be further researched for their respective effect, either individually or combined, on specific clinical endpoints. We showed that implementing indicators (average, slope coefficient, MSE, and lag-1 autocorrelation) as covariates which measure dynamical properties of the human health trajectory can increase the predictive performance of models on crucial clinical endpoints, as shown in the values of the hazard ratios and the AUC of the analysis. Therefore, we propose the use of such method on future epidemiological studies of data consisting of longitudinal measurements of health-related variables.

We consider the current study as an extension of the ideas of previous work on complex systems and DIORs, in the field of life sciences [6, 18, 21, 22], aiming to make the logic applicable in register-based longitudinal healthcare data. For that purpose, indicators were chosen primarily based on their ability to capture distinct dynamics of healthcare consumption, without being restricted to match the paradigm of critical slowing down. Interestingly, using DIORs as explanatory variables, along with age and sex, reached an AUC of over 80% for a multifactorial event such as death. We believe that by combining dynamical metrics with more static variables, such as socioeconomic status, education, marital status, Body Mass Index, or comorbidities could greatly enhance the overall performance while also contributing to an even better understanding of resilience and thus death. In addition, the recent advances of wearable and sensor data can provide us with time series of plenty of data, which when analyzed properly by extracting dynamical indicators from them, can depict the dynamics of human health across time. Of course, the applications of these DIORs and their research can be extended beyond healthcare costs, investigating physiological responses across time, such as blood pressure, temperature, glycose levels, electric activity of the heart etc.

The strengths of this research are the sample size, the completeness of healthcare costs, which is based on national register data and is representative of the bereaved population, as well as the time to event analysis to study the behavior of these dynamical indicators in the field of Health sciences. To the best of our knowledge, no other studies under the DIORs framework have used health-expense related variables to quantify human resilience and predict mortality risk after the appearance of a stressor.

However, the study also has its limitations. Our current research uses healthcare expenditures as a measure of health. While in general, the expenditures reflect the health status of individuals and their need for cure and care, there might be potential bias since the null-costs behavior might be explained by the unwillingness of individuals to utilize the medical services rather than their good state of health. Nevertheless, the tax-based healthcare system of Denmark might be one of the factors attenuating the specific bias since it restricts possible financial burden and out of pocket expenses. With regards to the metrics’ computation, while there were over 50 thousand individuals included in the study, in order to study an abrupt stressor such as bereavement, DIORs were measured for a maximum of 52 consecutive weeks before and after bereavement. We are aware that the aforementioned is marginally enough to study the behavior of these indicators and their accompanied dynamics in rolling windows. Regarding the analysis, the results cannot be generalized for the whole population of Denmark since the data include only individuals with age 65 or older who have experienced spousal bereavement. Last, further research is needed to investigate if the constructed tertiles and overall analysis can be applied to bereaved individuals over the age of 65 in countries other than Denmark. The reason for that is that healthcare expenditures are bonded with the healthcare system and social welfare of each country, a factor which may lead to different patterns in the data.

In conclusion, the current research investigated the aggregated time series of healthcare expenditures in Danish bereaved older adults available from the nationwide registers. The tax-based Danish healthcare system, which manages to provide free access to services in its residents, renders those time-series of expenditures proximal data on peoples’ level of health and its dynamics due to the presence of longitudinal measurements. We proposed a novel methodological framework, which via several indicators, investigated the signals of such valuable data and their potential in quantifying health resilience. With regards to the DIORs, we focused on the metrics of average, slope, MSE and lag-1 autocorrelation and evaluated their ability to signal loss of resilience within the health sciences. Spousal bereavement is a stressor which is known to shift the dynamics of human health. Older people who suffer bereavement seem to have increased average consumption, rate of increase, and variability in the healthcare consumption patterns. We showed that these indicators which signal loss of stability in the physical sciences, seem to also be strong predictors of loss of health resilience as well. Indeed, individuals exhibiting high values of the respective DIORs before bereavement, were the ones who had increased mortality risk after the stressor. Hence, the higher these indicators are, the less resilient an individual seems to be. We believe that the implementation of a dynamical approach, which not only considers static predictors but also these ones that capture physiological and dynamic responses (variability, drift), can add significant value to the understanding and assessment of human’s health status.

Availability of data and materials

According to Danish law, scientific organizations can be authorized to work with data within Statistics Denmark and can provide access to individual scientists inside and outside of Denmark. Data are available via the Research Service Department at Statistics Denmark: ( for researchers who meet the criteria for access to confidential data.


  1. Baggio JA, Brown K, Hellebrandt D. Boundary object or bridging concept? A citation network analysis of resilience. Ecol Soc. 2015;20(2).

  2. Veraart AJ, Faassen EJ, Dakos V, van Nes EH, Lürling M, Scheffer M. Recovery rates reflect distance to a tipping point in a living system. Nature. 2012;481(7381):357–9.

    Article  CAS  Google Scholar 

  3. Dai L, Korolev KS, Gore J. Slower recovery in space before collapse of connected populations. Nature. 2013;496(7445):355–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Romero-Ortuno R, Cogan L, O’Shea D, Lawlor BA, Kenny RA. Orthostatic haemodynamics may be impaired in frailty†. Age Ageing. 2011;40(5):576–83.

    Article  PubMed  Google Scholar 

  5. Gijzel SMW, Whitson HE, van de Leemput IA, Scheffer M, van Asselt D, Rector JL, et al. Resilience in Clinical Care: Getting a Grip on the Recovery Potential of Older Adults. J Am Geriatr Soc. 2019;67(12):2650–7.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Gijzel SMW, van de Leemput IA, Scheffer M, van Bon GEA, Weerdesteyn V, Eijsvogels TMH, et al. Dynamical Indicators of Resilience in Postural Balance Time Series Are Related to Successful Aging in High-Functioning Older Adults. Journals Gerontol Ser A. 2019;74(7):1119–26.

    Article  Google Scholar 

  7. Whitson HE, Duan-Porter W, Schmader KE, Morey MC, Cohen HJ, Colón-Emeric CS. Physical resilience in older adults: Systematic review and development of an emerging construct. Journals Gerontol - Ser A Biol Sci Med Sci. 2016;71(4):489–95.

    Article  Google Scholar 

  8. Hoogendijk EO, Afilalo J, Ensrud KE, Kowal P, Onder G, Fried LP. Frailty: implications for clinical practice and public health. Lancet. 2019;394(10206):1365–75.

    Article  PubMed  Google Scholar 

  9. Scheffer M, Elizabeth Bolhuis J, Borsboom D, Buchman TG, Gijzel SMW, Goulson D, et al. Quantifying resilience of humans and other animals. Proc Natl Acad Sci U S A. 2018;115(47):11883–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Scheffer M, Carpenter SR, Lenton TM, Bascompte J, Brock W, Dakos V, et al. Anticipating critical transitions. Science. 2012;338(6105):344–8.

    Article  CAS  PubMed  Google Scholar 

  11. Scheffer M, Bascompte J, Brock WA, Brovkin V, Carpenter SR, Dakos V, et al. Early-warning signals for critical transitions. Nature. 2009;461(7260):53–9.

    Article  CAS  PubMed  Google Scholar 

  12. Thom R. Topological models in biology. Topology. 1969;8(3):313–35.

    Article  Google Scholar 

  13. van Nes EH, Arani BMS, Staal A, van der Bolt B, Flores BM, Bathiany S, et al. What Do You Mean, ‘Tipping Point’? Trends Ecol Evol. 2016;31(12):902–4.

    Article  PubMed  Google Scholar 

  14. Guttal V, Jayaprakash C. Changing skewness: an early warning signal of regime shifts in ecosystems. Ecol Lett. 2008;11(5):450–60.

    Article  PubMed  Google Scholar 

  15. Dakos V, Carpenter SR, Brock WA, Ellison AM, Guttal V, Ives AR, et al. Methods for detecting early warnings of critical transitions in time series illustrated using simulated ecological data. PLoS One. 2012;7(7).

  16. Carpenter SR, Brock WA. Rising variance: a leading indicator of ecological transition. Ecol Lett. 2006;9(3):311–8.

    Article  CAS  PubMed  Google Scholar 

  17. Dakos V, van Nes EH, Donangelo R, Fort H, Scheffer M. Spatial correlation as leading indicator of catastrophic shifts. Theor Ecol. 2010;3(3):163–74.

    Article  Google Scholar 

  18. Gijzel SMW, van de Leemput IA, Scheffer M, Roppolo M, Olde Rikkert MGM, Melis RJF. Dynamical Resilience Indicators in Time Series of Self-Rated Health Correspond to Frailty Levels in Older Adults. Journals Gerontol Ser A. 2017;72(7):991–6.

    Article  Google Scholar 

  19. Olde Rikkert MGM, Dakos V, Buchman TG, de Boer R, Glass L, Cramer AOJ, et al. Slowing Down of Recovery as Generic Risk Marker for Acute Severity Transitions in Chronic Diseases. Crit Care Med. 2016;44(3):601–6.

    Article  CAS  PubMed  Google Scholar 

  20. Nelson B, McGorry PD, Wichers M, Wigman JTW, Hartmann JA. Moving From Static to Dynamic Models of the Onset of Mental Disorder. JAMA Psychiatry. 2017;74(5):528.

    Article  PubMed  Google Scholar 

  21. Gijzel SMW, Rector J, van Meulen FB, van der Loeff RS, van de Leemput IA, Scheffer M, et al. Measurement of Dynamical Resilience Indicators Improves the Prediction of Recovery Following Hospitalization in Older Adults. J Am Med Dir Assoc. 2020;21(4):525-530.e4.

    Article  PubMed  Google Scholar 

  22. Rector JL, Gijzel SMW, van de Leemput IA, van Meulen FB, Olde Rikkert MGM, Melis RJF. Dynamical indicators of resilience from physiological time series in geriatric inpatients: Lessons learned. Exp Gerontol. 2021;149:111341.

    Article  PubMed  Google Scholar 

  23. Das A. Spousal Loss and Health in Late Life. J Aging Health. 2013;25(2):221–42.

    Article  PubMed  Google Scholar 

  24. Holm AL, Severinsson E, Berland AK. The Meaning of Bereavement Following Spousal Loss: A Qualitative Study of the Experiences of Older Adults. SAGE Open. 2019;9(4):215824401989427.

    Article  Google Scholar 

  25. Stroebe M, Schut H, Stroebe W. Health outcomes of bereavement. Lancet. 2007;370(9603):1960–73.

    Article  PubMed  Google Scholar 

  26. Tseng F-M, Petrie D, Wang S, Macduff C, Stephen AI. The impact of spousal bereavement on hospitalisations: Evidence from the Scottish Longitudinal Study. Health Econ. 2018;27(2):e120-38.

    Article  PubMed  Google Scholar 

  27. Seiler A, von Känel R, Slavich GM. The Psychobiology of Bereavement and Health: A Conceptual Review From the Perspective of Social Signal Transduction Theory of Depression. Front Psychiatry [Internet]. 2020 Dec 3;11. Available from:

  28. Parkes CM. Coping with loss: Bereavement in adult life. BMJ. 1998;316(7134):856–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Blanner C, Mejldal A, Prina AM, Munk-Jørgensen P, Ersbøll AK, Andersen K. Widowhood and mortality: a Danish nationwide register-based cohort study. Epidemiol Psychiatr Sci. 2020;29:e149.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Hart CL, Hole DJ, Lawlor DA, Davey Smith G, Lever TF. Effect of conjugal bereavement on mortality of the bereaved spouse in participants of the Renfrew/Paisley Study. J Epidemiol Community Heal. 2007;61(5):455–60.

    Article  Google Scholar 

  31. Shor E, Roelfs DJ, Curreli M, Clemow L, Burg MM, Schwartz JE. Widowhood and Mortality: A Meta-Analysis and Meta-Regression. Demography. 2012;49(2):575–606.

    Article  PubMed  Google Scholar 

  32. Sullivan AR, Fenelon A. Patterns of Widowhood Mortality. Journals Gerontol Ser B Psychol Sci Soc Sci. 2014;69B(1):53–62.

    Google Scholar 

  33. Stroebe M, Stroebe W, Schut H. Gender Differences in Adjustment to Bereavement: An Empirical and Theoretical Review. Rev Gen Psychol. 2001;5(1):62–83.

    Article  Google Scholar 

  34. Manor O, Eisenbach Z. Mortality after spousal loss: are there socio-demographic differences? Soc Sci Med. 2003;56(2):405–13.

    Article  PubMed  Google Scholar 

  35. Moon JR, Kondo N, Glymour MM, Subramanian S V. Widowhood and Mortality: A Meta-Analysis. Hartling L, editor. PLoS One [Internet]. 2011 Aug 17;6(8):e23465. Available from:

  36. Bonanno GA, Wortman CB, Lehman DR, Tweed RG, Haring M, Sonnega J, et al. Resilience to loss and chronic grief: A prospective study from preloss to 18-months postloss. J Pers Soc Psychol. 2002;83(5):1150–64.

    Article  PubMed  Google Scholar 

  37. Bonanno GA, Kaltman S. Toward an integrative perspective on bereavement. Psychol Bull. 1999;125(6):760–76.

    Article  CAS  PubMed  Google Scholar 

  38. Bonanno GA, Kaltman S. The varieties of grief experience. Clin Psychol Rev [Internet]. 2001 Jul;21(5):705–34. Available from:

  39. Galatzer-Levy IR, Bonanno GA. Beyond normality in the study of bereavement: Heterogeneity in depression outcomes following loss in older adults. Soc Sci Med. 2012;74(12):1987–94.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Holmes TH, Rahe RH. The social readjustment rating scale. J Psychosom Res. 1967;11(2):213–8.

    Article  CAS  PubMed  Google Scholar 

  41. Malgaroli M, Maccallum F, Bonanno GA. Machine yearning: How advances in computational methods lead to new insights about reactions to loss. Curr Opin Psychol. 2022;43:13–7.

    Article  PubMed  Google Scholar 

  42. Hadley EC, Kuchel GA, Newman AB, Allore HG, Bartley JM, Bergeman CS, et al. Report: NIA Workshop on Measures of Physiologic Resiliencies in Human Aging. Journals Gerontol Ser A. 2017;72(7):980–90.

    Article  Google Scholar 

  43. Hansen AV, Mortensen LH, Trompet S, Westendorp R. Health care expenditure in the last five years of life is driven by morbidity, not age: A national study of spending trajectories in Danish decedents over age 65. Ramagopalan S V., editor. PLoS One [Internet]. 2020 Dec 18;15(12):e0244061. Available from:

  44. Christensen BJ, Gørtz M, Kallestrup-Lamb M. Medical Spending in Denmark. Fisc Stud. 2016;37(3–4):461–97.

    Article  Google Scholar 

  45. Huang Q, Zhang H, Chen J, He M. Quantile Regression Models and Their Applications: A Review. J Biom Biostat [Internet]. 2017;08(03). Available from:

  46. Steyerberg EW, Vickers AJ, Cook NR, Gerds T, Gonen M, Obuchowski N, et al. Assessing the performance of prediction models: A framework for traditional and novel measures. Epidemiology. 2010;21(1):128–38.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Heagerty PJ, Zheng Y. Survival model predictive accuracy and ROC curves. Biometrics. 2005;61(1):92–105.

    Article  PubMed  Google Scholar 

  48. Kusumastuti S, Gerds TA, Lund R, Mortensen EL, Westendorp RGJ. Discrimination ability of comorbidity, frailty, and subjective health to predict mortality in community-dwelling older people: Population based prospective cohort study. Eur J Intern Med. 2017;42:29–38.

    Article  PubMed  Google Scholar 

  49. Efron B. Estimating the Error Rate of a Prediction Rule: Improvement on Cross-Validation. J Am Stat Assoc. 1983;78(382):316.

    Article  Google Scholar 

  50. R Core Team (2013). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL

  51. Leeuwenberg AM, van Smeden M, Langendijk JA, van der Schaaf A, Mauer ME, Moons KGM, et al. Performance of binary prediction models in high-correlation low-dimensional settings: a comparison of methods. Diagnostic Progn Res. 2022;6(1):1.

    Article  Google Scholar 

  52. Kusumastuti S, Rozing MP, Lund R, Mortensen EL, Westendorp RGJ. The added value of health indicators to mortality predictions in old age: A systematic review. Eur J Intern Med. 2018;57:7–18.

    Article  PubMed  Google Scholar 

Download references


The authors kindly thank Statistics Denmark for their precious help with data access.


This work is funded by grants from Novo Nordisk Fonden [NNF17OC0027812].

Author information

Authors and Affiliations



AK, PN and RGJW were involved in the conceptualization of the research question, identified and appraised relevant literature. MJ acted as fourth reviewer, providing critical input in the analysis, and assisted in the interpretation of the results. AK drafted the manuscript and performed the computational analysis with critical input from all other authors. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Alexandros Katsiferis.

Ethics declarations

Ethics approval and consent to participate

All methods were performed in accordance with the relevant guidelines and regulations. Danish legislation allows for register-based research of this type to be conducted without the consent of participants and without ethical committee approval. The study was conducted according to the rules of the Danish Data Protection Agency. All data was held at Statistics Denmark, which is the Danish national statistical institution.

Consent for publication

Not applicable.

Competing interests

The authors declare that they do not have any competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file1: Supplementary Figure 1.

Weekly average of healthcare expenditures pre- and post-bereavement between groups. Supplementary Figure 2. Weekly average residuals of healthcare expenditures pre- and post-bereavement between groups. Supplementary Table S1. Median (IQR) of DIORs between males and females across bereavement periods. Supplementary Table S2. Quantile-median regression coefficients of DIORs regressed on age adjusted for sex. Supplementary Table S3. Median (IQR) of DIORs between ‘Deceased’ and ‘Survivors’ one year before bereavement for males. Supplementary Table S4. Median (IQR) of DIORs between ‘Deceased’ and ‘Survivors’ one year before bereavement for females. Supplementary Table S5. Median (IQR) of DIORs before bereavement between ‘Deceased’ and ‘Survivors’ stratified on age categories. Supplementary Table S6. Hazard ratios (95% CI) of all-cause mortality after bereavement dependent on DIORs Score

Additional file2: 

Mathematical Formulas of Dynamical Indicators of Resilience (DIORs) and statistical models 

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Katsiferis, A., Nielsen, P.Y., Jensen, M.K. et al. Dynamical indicators in time series of healthcare expenditures predict mortality risk of older adults following spousal bereavement. BMC Geriatr 22, 301 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Resilience
  • Bereavement
  • Complex systems
  • Health
  • Aging