A fast EM algorithm for fitting joint models of a binary response and multiple longitudinal covariates subject to detection limits

Computational Statistics & Data Analysis, 85, 37–53.