Joint models for longitudinal and time-to-event data have recently attracted a lot of attention in statistics and biostatistics. Even though these models enjoy a wide range of applications in many different statistical fields, they have not yet found their rightful place in the toolbox of modern applied statisticians mainly due to the fact that they are rather computationally intensive to fit. The main difficulty arises from the requirement for numerical integration with respect to the random effects. This integration is typically performed using Gaussian quadrature rules whose computational complexity increases exponentially with the dimension of the random-effects vector. A solution to overcome this problem is proposed using a pseudo-adaptive GaussHermite quadrature rule. The idea behind this rule is to use information for the shape of the integrand by separately fitting a mixed model for the longitudinal outcome. Simulation studies show that the pseudo-adaptive rule performs excellently in practice, and is considerably faster than the standard GaussHermite rule.

GaussHermite rule, Numerical integration, Random effects, Survival analysis, Time-dependent covariates
dx.doi.org/10.1016/j.csda.2011.09.007, hdl.handle.net/1765/65530
Computational Statistics & Data Analysis
Department of Biostatistics

Rizopoulos, D. (2012). Fast fitting of joint models for longitudinal and event time data using a pseudo-adaptive Gaussian quadrature rule. Computational Statistics & Data Analysis, 56(3), 491–501. doi:10.1016/j.csda.2011.09.007