[R] HLM Model

[Doran, Harold]

I think it should be

fm <- lmer(score ~ trt + pair + grade + school + (1|team), test)

The unstructured covariance matrix you use in proc mixed is not available in Rs function for mixed models

> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On
> Behalf Of Belle
> Sent: Thursday, January 27, 2011 2:44 PM
> To: r-help at r-project.org
> Subject: [R] HLM Model
> 
> 
> Hi
> 
> I am trying to convert SAS codes to R, but some of the result are quite
> different from SAS.
> 
> When I ran proc mixed, I have an option ddfm=bw followed by the model. How
> can I show this method in R　(I am thinking that this maybe the reason that I
> can't get the similar results)
> 
> below is my SAS codes:
> 
> proc mixed data=test covtest empirical;
> class pair grade team school;
> 	model score = trt pair grade school/ solution covb ddfm=bw ;
> 	random int /  sub=team solution type=un;
> run;
> 
> I have tried both lmer and hglm, but non of them works.
> 
> Could anyone tell me how can I covert this SAS codes to R? Thanks
> --
> View this message in context: http://r.789695.n4.nabble.com/HLM-Model-
> tp3242999p3242999.html
> Sent from the R help mailing list archive at Nabble.com.
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.