Displaying 1 result from an estimated 1 matches for "inv_2".
Did you mean:
in8_2
2008 Nov 19
2
GAMM and anove.lme question
...Hi R user,
I am using the gamm() function of the mgcv-package. Now I would like to
decide on the random effects to include in the model. Within a GAMM
framework, is it allowed to compare the following two models
inv_1<-gamm(y~te(sat,inv),data=daten_final, random=list(proband=~1))
inv_2<-gamm(y~te(sat,inv),data=daten_final, random=list(proband=~sat))
with a likelihood ratio test for a traditional GLMM, like this:
anova(inv_1$lme, inv_2$lme)
The output is as follows:
Model df AIC BIC logLik Test L.Ratio p-value
inv_2$lme 1 10 21495.90 2155...