Displaying 1 result from an estimated 1 matches for "fm1nidl".
2008 May 22
1
mixed model resuts from SAS and R
...nR = out1.randomnidltest2;
proc mixed data = a1 ;
class sub ;
model nidl = time / solution ;
random int time / sub = sub solution;
run;
ods output close;
2. mixed model in R:
====================
a1<-read.table(file="c:\\aldi\\a1.txt",sep=",",header=T)
library(nlme)
fm1nidl.lme<-lme(nidl~time,data=a1,random=~time | sub)
plot(coef(fm1nidl.lme))
3. SAS output:
==============
Plot of nidl*time. Symbol used is '*'.
40 ?
?
?*
?
?*
?* *
?* *
?* * *
20 ?* * *
?*...