Displaying 1 result from an estimated 1 matches for "treatpla".
Did you mean:
treat_a
2002 May 27
1
nlme cross-over and fixed nested
...quot;Pla"))*rnorm(NSubj,2
,3)
# Sequence as fixed effect
glu.lme<-lme(Glu~Treat*Seq,data=glu,random=~1|Patient)
print(summary(glu.lme))
# result looks perfectly valid for me
# Value Std.Error DF t-value p-value
#(Intercept) 24.694617 1.602160 18 15.413324 <.0001
#TreatPla 5.106889 1.096139 18 4.658980 0.0002
#Seq2 0.466712 2.265797 18 0.205981 0.8391
#TreatPla:Seq2 3.525737 1.550174 18 2.274413 0.0354
# Sequence as random ???????
glu.lme<-lme(Glu~Treat,data=glu,random=~1|Seq/Patient)
print(summary(glu.lme))
Dieter Menne
-----------...