search for: bigad

...de that I was writting in R is as follows: 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=bigads) I am not sure if I have got the R script/options correct...but I seem to be getting different estimates from the same dataset.... any help on this would be highly appreciated!!!! -- View this message in context: http://r.789695.n4.nabble.com/LMER-vs-PROC-MIXED-estimates-tp4648653.html Sent...