Dear Nicholas, The sas code may estimate a different variance for each study_id. I'm not sure, but mixed in spss would do this if repeated is used. In that case, add something like "weights = varIdent(form = ~1|study_id)" to your lme code in R. Regards, Ben. [[alternative HTML version deleted]]