Recent Publications

More Publications

Selected Publications

JAMA Intern Med

Recent & Upcoming Talks

More Talks

Recent Posts

I have submitted updates to both joineR (now version 1.2.1) and joineRML (now version 0.3.0) to CRAN this week. In addition, we now have fantastic looking hex sticker badges from the remarkable hexSticker package. As always, we value user feedback on these packages and will continue to update them.

joineRML v0.3.0


Description As of version 1.2.0, joineR fits the joint model proposed by Williamson et al. (2008) for joint models of longitudinal data and competing risks. Here, the longitudinal data submodel remains as per that analyzed in Henderson et al. (2000), namely \[Y_i(t) = X_i(t)^\top \beta + Z_i(t)^\top U_i + \epsilon_i(t),\] where \(Y_i(t)\) is a repeated measurement on subject \(i\) at time \(t\), \(U_i\) is a latent vector that follows a zero-mean multivariate normal distribution, \(X_i(t)\) and \(Z_i(t)\) are vectors of explanatory variables that may be time-constant or time-varying, and the \(\epsilon_i(t)\) are mutually independent errors, \(\epsilon_i(t) \sim N(0, \tau^2)\).


The joineRML package implements methods for analyzing data from multiple longitudinal studies in which the responses from each subject consists of time-sequences of repeated measurements and a possibly censored time-to-event outcome. The modelling framework for the repeated measurements is the multivariate linear mixed effects model. The model for the time-to-event outcome is a Cox proportional hazards model with log-Gaussian frailty. Stochastic dependence is captured by allowing the Gaussian random effects of the linear model to be correlated with the frailty term of the Cox proportional hazards model.



Cardiovascular Surgery Outcomes

Leveraging clinical registry data for research into outcomes following cardiovascular surgery.

JoineR-M: Joint Modelling of Repeated Measurements and Time-to-Event Data

The development of methodology and software for fitting joint models to multivariate longitudinal and time-to-event data.

Medical Statistics Education

Medical statistics teaching for clinicians.

Statistical Ecotoxicological Modelling

Ecotoxicological risk assessment for chemicals using species sensitivity distributions.

joineR: Joint Modelling of Repeated Measurements and Time-to-Event Data

An R package to fit joint models to a single repeated measure and a time-to-event outcome (single or competing risks) using an EM algorithm.

joineRML: Joint Modelling of Multivariate Longitudinal Data and Time-to-Event Outcomes

An R package to fit joint models to multivariate repeated measures data and a time-to-event outcome using an MCEM algorithm.


I was a teaching instructor/tutor for the following courses/workshops at Liverpool University and Durham University: