You can probably get there with some programming steps to model the covariance structure, but I don't think it will be easy. Using the mvn distribution in the RANDOM statement looks to me like the equivalent of an unstructured matrix in PROC MIXED, but programming to set parameters in the covariance matrix sigma to fit the specified error structure would be needed, and I think will take a lot of trial and error. Google doesn't seem to have much in the way of help. Perhaps there is an R package (probably nmle) that has the coding that could be adapted. Good luck on this one. If you find a good answer, I really, really, really want to know, as I have over thirty studies where this could be applied. Steve Denham
... View more