Journal
STATISTICS IN MEDICINE
Volume 35, Issue 7, Pages 1193-1209Publisher
WILEY
DOI: 10.1002/sim.6779
Keywords
joint modelling; delayed entry; left truncation; survival analysis; mixed effects; adaptive Gauss-Hermite quadrature
Categories
Funding
- National Institute for Health Research [DRF-2012-05-409, NF-SI-0512-10159] Funding Source: researchfish
- National Institutes of Health Research (NIHR) [DRF-2012-05-409] Funding Source: National Institutes of Health Research (NIHR)
- Department of Health [DRF-2012-05-409] Funding Source: Medline
Ask authors/readers for more resources
A now common goal in medical research is to investigate the inter-relationships between a repeatedly measured biomarker, measured with error, and the time to an event of interest. This form of question can be tackled with a joint longitudinal-survival model, with the most common approach combining a longitudinal mixed effects model with a proportional hazards survival model, where the models are linked through shared random effects. In this article, we look at incorporating delayed entry (left truncation), which has received relatively little attention. The extension to delayed entry requires a second set of numerical integration, beyond that required in a standard joint model. We therefore implement two sets of fully adaptive Gauss-Hermite quadrature with nested Gauss-Kronrod quadrature (to allow time-dependent association structures), conducted simultaneously, to evaluate the likelihood. We evaluate fully adaptive quadrature compared with previously proposed non-adaptive quadrature through a simulation study, showing substantial improvements, both in terms of minimising bias and reducing computation time. We further investigate, through simulation, the consequences of misspecifying the longitudinal trajectory and its impact on estimates of association. Our scenarios showed the current value association structure to be very robust, compared with the rate of change that we found to be highly sensitive showing that assuming a simpler trend when the truth is more complex can lead to substantial bias. With emphasis on flexible parametric approaches, we generalise previous models by proposing the use of polynomials or splines to capture the longitudinal trend and restricted cubic splines to model the baseline log hazard function. The methods are illustrated on a dataset of breast cancer patients, modelling mammographic density jointly with survival, where we show how to incorporate density measurements prior to the at-risk period, to make use of all the available information. User-friendly Stata software is provided. (C) 2015 The Authors. Statistics in Medicine Published by JohnWiley & Sons Ltd.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available