Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2023 Mar;32(3):474-492.
doi: 10.1177/09622802221146307. Epub 2022 Dec 26.

Bivariate joint models for survival and change of cognitive function

Shengning Pan  1 , Ardo van den Hout  1
Affiliations

Bivariate joint models for survival and change of cognitive function

Shengning Pan et al. Stat Methods Med Res. 2023 Mar.

Abstract

Changes in cognitive function over time are of interest in ageing research. A joint model is constructed to investigate. Generally, cognitive function is measured through more than one test, and the test scores are integers. The aim is to investigate two test scores and use an extension of a bivariate binomial distribution to define a new joint model. This bivariate distribution model the correlation between the two test scores. To deal with attrition due to death, the Weibull hazard model and the Gompertz hazard model are used. A shared random-effects model is constructed, and the random effects are assumed to follow a bivariate normal distribution. It is shown how to incorporate random effects that link the bivariate longitudinal model and the survival model. The joint model is applied to the English Longitudinal Study of Ageing data.

Keywords: Joint model; bivariate binomial distribution; cognitive function; shared random-effects model; survival analysis.

PubMed Disclaimer

Conflict of interest statement

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Figures

Figure 1.
Figure 1.
Bias and corresponding Monte Carlo 95% confidence interval. Circles represent biases, and geometry bars represent Monte Carlo 95% confidence interval.
Figure 2.
Figure 2.
Frequency distribution of recalled words at baseline interview. Y1 represents immediate-recall words; Y2 represents delayed-recall words.
Figure 3.
Figure 3.
Recalled words trajectories. Left-hand side: immediate-recall words Y1 trajectories. Right-hand side: delayed-recall words Y2 trajectories. Individuals are represented in the same colour in both plots.
Figure 4.
Figure 4.
Fitted bivariate binomial distribution for immediate recall ( Y1) and delayed recall ( Y2). Left hand side: conditional on age 50; Right hand side: conditional on age 80.
Figure 5.
Figure 5.
Comparison of K-M survival curves: predicted survival curves (grey lines), the mean of those survival curves (blue line) and the median of survival curves (red line).
See this image and copyright information in PMC

References

    1. Diggle PJ, Heagerty P, Liang KY. et al. Analysis of Longitudinal Data. Oxford: Oxford University Press, 2002.
    1. Banks J, Breeze E, Lessof C. et al. Retirement, health and relationships of the older population in england: The 2004 english longitudinal study of ageing (wave 2), 2006.
    1. Hickey GL, Philipson P, Jorgensen A. et al. Joint modelling of time-to-event and multivariate longitudinal outcomes: Recent developments and issues. BMC Med Res Methodol 2016; 16: 1–15. - PMC - PubMed
    1. Van Den Hout A, Muniz-Terrera G. Joint models for discrete longitudinal outcomes in aging research. J R Stat Soc: Ser C (Applied Statistics) 2016; 65: 167–186.
    1. Laird NM, Ware JH. Random-effects models for longitudinal data. Biometrics 1982; 38: 963–974. - PubMed

Publication types