Displaying 1 result from an estimated 1 matches for "phqmodel2".
Did you mean:
phqmodel1
2012 May 31
1
anova of lme objects (model1, model2) gives different results depending on order of models
...tter
which unfortunately only returns results about ANOVA which isn't much
help. I'm using the latest version of R and nlme, just checked both.
Here is the code and output:
> PHQmodel1=lme(PHQ~Age+Gender+Date*Treatment, data=compfinal,
random=~1|Case, na.action=na.omit)
>
> PHQmodel2=lme(PHQ~Age+Gender+Date*Treatment, data=compfinal,
random=~1|Case, na.action=na.omit,
+ correlation=corAR1(form=~Date|Case))
> anova(PHQmodel1, PHQmodel2) # accept model 2
Model df AIC BIC logLik Test
L.Ratio p-value
PHQmodel1 1 8 48784....