Displaying 1 result from an estimated 1 matches for "adstest".
Did you mean:
bdstest
2012 Nov 07
2
LMER vs PROC MIXED estimates
...0 METHOD=REML MAXITER=1000;
CLASS GEOGRAPHY;
MODEL y = X1 X2 X3/SOLUTION;
RANDOM INTERCEPT X2/SOLUTION SUBJECT = GEOGRAPHY;
ODS OUTPUT SOLUTIONR=RANDOM_EFFECT;
ODS OUTPUT SOLUTIONF= FIXED_EFFECT;
RUN;
the equivalent code 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 t...