![](index_files/figure-html/fig-1-1.png)
Dong, Milholland, and Vijg (2016) claimed in a Nature paper that the maximum reported age at death has been decreasing, and that this result is significant. Their main finding is summarized in the following figure:
The latter shows (a) linear regressions of the maximum reported age at death over time, with a changepoint the year after Jeanne Calment died, (b) cubic smoothing splines fitted separately to the \(r=5\) largest order statistics in each year and (c) sample average of the age at death of supercentenarians, smoothing using a cubic smoothing spline.
While the problems with the statistical analyses have been reported elsewhere Rootzén and Zholud (2017), the purpose is this post to use recent International Database on Longevity (IDL) (French Institute for Demographic Studies 2023) to revisit some of the points of contention and highlight their impact on inference.
When it was added online in 2010, the IDL included records of supercentenarians (i.e., people who died above the age of 110 years) who died during a specified period of time, denoted \(c_1\) and \(c_2\). These calendar dates define the sampling frame, a form of interval truncation. Indeed, people who were born in year 1900, say, who survived beyond 2010 would not be included. Similarly, for earlier cohorts, only the oldest individuals born in year \(l\) would be alive if \(l-c_1 > 110\). If the rate at which people reach old age was the same, these two effects would counterbalance one another, but this is not the case here since mortality has decreased over the course of the century and the world population has increased, so many ghosts are found near \(c_2\).
Dong, Milholland, and Vijg (2016) computed the sample age at death of people and the maximum reported age at death per birth year, sometimes referred to (incorrectly) as a birth cohort. As the IDL was augmented with more recent data, the pattern arising from the truncation disappears: this can be seen from Figure 2.3 from Jdanov, Shkolnikov, and Gellers-Barkmann (2021).
Changepoints
Dong, Milholland, and Vijg (2016) used two linear regression, splitting the data into two segments: before 1997 and after 1997. This year is notable because it is the age at which Jeanne Calment died. If data were independent and identically distributed (they aren’t), then we could indeed fit a linear regression and perform a Chow test (Chow 1960) to check. Choosing the maximum value is a form of selection mechanism, as it seems visually to be the point that maximizes the likelihood ratio statistic and the distribution is non-standard (Kim and Siegmund 1989). Brown, Albers, and Ritchie (2017) show that the downward trend would be removed altogether if we picked an earlier changepoint.
This approach is flawed from the start since the data are not identically distributed (due to truncation), the maximum of larger samples is stochastically larger than that of smaller samples (given the number of people over which we compute the maximum is different for each year) and the test procedure needs to account for the changepoint estimation. The truncation effects disappear if we consider longer collection period; for example, Figure 1 shows the data for France, England and Wales, USA and Japan from a more recent of the IDL (2023), showing indeed that much of the downward trend is gone.
![](index_files/figure-html/fig-1-1.png)
If we were to fit (incorrectly because it does not account for the truncation, nonstationarity of varying sample size for each year) a changepoint model with different intercept and slopes for year at death, we get:
<- chngpt::chngptm(
changepoint formula.1 = MRAD ~ 1,
formula.2 = MRAD ~ DYEAR,
data = idl_sub |> mutate(DYEAR = DEATH_YEAR - 1980),
type = "stegmented",
family = "gaussian",
var.type = "bootstrap",
bootstrap.type = "wild",
ci.bootstrap.size = 100L)
summary(changepoint)
Change point model threshold.type: stegmented
Coefficients:
est Std. Error* (lower upper) p.value*
(Intercept) 111.0443293 0.23299674 110.68618013 111.5995274 0.0000000000
DYEAR 0.1323626 0.03895209 0.02327091 0.1759631 0.0006785841
DYEAR>chngpt 1.6593439 1.81340901 -3.84508591 3.2634774 0.3601700308
(DYEAR-chngpt)+ -0.2545795 0.12671462 -0.33566793 0.1610534 0.0445288691
Threshold:
est Std. Error (lower upper)
15.000000 4.336735 11.000000 28.000000
While these results still fail to account for truncation and sample size, the ‘statistical significance’ is now gone.
Knots placement
The downward trend is mostly gone with the new data, but is less noticeable when we remove Jeanne Calment, would also be compensated if the data included Italy or more recent data from Japan (with reported people dying age 118 and 119). With limited data, the cubic smoothing splines shows noticeable artefacts of oversmoothing in Dong, Milholland, and Vijg (2016). The placement of the knots affects the conclusions drawn, as evidenced by automatic tuning in the right-hand panel of Figure 2, to be compared with equispaced knots in the left-hand panel.
![](index_files/figure-html/fig-2-1.png)
We notice in Figure 2 that, after a dip early 2000s (attributable to the heatwave?), the trend is going back up. There are thus hardly evidence for the claimed decrease.
Further effects of truncation
For positive random variables, we can estimate the expectation by computing the area under the survival curve. Since the sampling frame induces interval truncation and we seek to estimate the mean of the (untruncated) underlying distribution of excess lifetimes, we can instead use an EM algorithm to compute the nonparametric maximum likelihood estimator of the age at death and extract a mean estimator (Turnbull 1976). Figure 3 shows that the points are rather different from the sample mean and much less different, seeming to hover around a constant mean. We can also use the sample size as weighting in regression to draw a trend.
![](index_files/figure-html/fig-3-1.png)
Conclusion
This post tried to demonstrate that much of the claim that there is a downward trend in the maximum age at death are based on flawed statistical analyses and that small tweaks and newer data can lead to much different conclusions, highlighting their lack of robustness. While there may have been a decrease and Jeanne Calment’s record is still very much an outlier today, the evidence does not hold to scrutinity.
References
Citation
@online{belzile2023,
author = {Belzile, Léo},
title = {Trends in Longevity Records},
date = {2023-02-28},
url = {https://lbelzile.bitbucket.io/posts/trend-MRAD/},
langid = {en}
}