Converting SAS code into R

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.