Patients who have undergone renal transplantation are monitored longitudinally at irregular time intervals over 10 years or more.This yields a set of biochemical and physiological markers containing valuable information to anticipate a failure of the graft.A general linear, generalized linear, or nonlinear mixed model is used to describe the longitudinal profile of each marker.To account for the correlation between markers, the univariate mixed models are combined into a multivariate mixed model (MMM) by specifying a joint distribution for the random effects.Due to the high number of markers, a pairwise modeling strategy, where all possible pairs of bivariate mixed models are fitted, is used to obtain parameter estimates for the MMM.These estimates are used in a Bayes rule to obtain, at each point in time, the prognosis for long-term success of the transplant.It is shown that allowing the markers to be correlated can improve this prognosis.