LMER vs PROC MIXED estimates

Hi David,

Thanks for the reply. sorry for typo in the R script ...what I meant to
write was
testdata <- read.csv("adstest.csv",header=TRUE,sep=",") 
attach(testdata) 
library(lme4) 
options(contrasts = c(factor = "contr.SAS",ordered = "contr.poly")) 
lmm.2=lmer(y~X1+X2+X3 + (X2|Geography),REML=TRUE,data=testdata)

not sure if attaching the dataset causes any problem.
I'm not _sure_ either, but I would definitely advise detach()-ing it. You should not need to attach when you are also passing testdata to the data argument.
Also, I would take
your suggestion and post this in Mixed Models SIG.

Regards,
Sandip

--
View this message in context: http://r.789695.n4.nabble.com/LMER-vs-PROC-MIXED-estimates-tp4648653p4648671.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.
David Winsemius, MD
Alameda, CA, USA

LMER vs PROC MIXED estimates

Thread (5 messages)