We apply a linear mixed-effects model to multivariate failure time data. Co
mputation of the regression parameters involves the Buckley-James method in
an iterated Monte Carlo expectation-maximization algorithm, wherein the Mo
nte Carlo E-step is implemented using the Metropolis-Hastings algorithm. Fr
om simulation studies, this approach compares favorably with the marginal i
ndependence approach, especially when there is a strong within-cluster corr
elation.