Estimating multivariate longitudinal trajectories using mixed-effects models with crossed random effects