Re: [R] HLM Model

2011-01-31 Thread Belle
Hi Silvano: Could you tell me what correlation=corSymm(form = ~ 1 |id) represents? In our case, team is random effect, trt, pairs, grade, school are fixed effect, and each team is within school. I still got the different results from both SAS and R. unstruct -

Re: [R] HLM Model

2011-01-28 Thread Silvano
-- - Original Message - From: Belle ping...@gmail.com To: r-help@r-project.org Sent: Thursday, January 27, 2011 5:43 PM 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

[R] HLM Model

2011-01-27 Thread Belle
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

Re: [R] HLM Model

2011-01-27 Thread Doran, Harold
] On Behalf Of Belle Sent: Thursday, January 27, 2011 2:44 PM To: r-help@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

Re: [R] HLM Model

2011-01-27 Thread Belle
Hi Harold: Yes, this was the R code that I tried, and got different result from SAS. Is that mean I cannot actually use R to run unstructured covariance matrix? How can I solve this problem if I need an unstructured covariance matrix method? Thanks for the help. -- View this message in

Re: [R] HLM Model

2011-01-27 Thread Doran, Harold
Sent: Thursday, January 27, 2011 3:54 PM To: r-help@r-project.org Subject: Re: [R] HLM Model Hi Harold: Yes, this was the R code that I tried, and got different result from SAS. Is that mean I cannot actually use R to run unstructured covariance matrix? How can I solve this problem

Re: [R] HLM Model

2011-01-27 Thread Belle
Hi Harold: I know the outputs are different between SAS and R, but the results that I got have big difference. Here is part of the result based on the SAS code I provided earlier: Cov Parm SubjectEstimate Error Value Pr Z

Re: [R] HLM Model

2011-01-27 Thread Jeremy Miles
The empirical statement on the proc mixed line gives you robust standard errors, I don't think you get them in R. In SAS you specify that the predictors are to be dummy coded using the class . Are they factors in R? I can't tell from the SAS output, because the formatting has been lost.