A general model for two-level multivariate data, with responses possibly missing at random, is described. The model combines regressions on fixed explanatory variables with structured residual covariance matrices. The likelihood function is reduced to a form enabling computational methods for estimating the model to be devised. (Author)