I don't know if this belongs here or in StackExchange, it is a mixed but probably pretty simple question. How do I normally report a Likelihood Ratio Test? I would love a good reference in your answer, I have searched but could not find any good answers.
> glmm0 <- glmer(yngel ~ (1|nest), data, family=poisson(link="log"))
> glmm <- glmer(yngel ~ age.level + (1|nest) + 0, data, family=poisson(link="log"))
> anova(glmm0, glmm)
Data: data
Models:
glmm0: yngel ~ (1 | nest)
glmm: yngel ~ age.level + (1 | nest) + 0
Df AIC BIC logLik deviance Chisq Chi Df Pr(>Chisq)
glmm0 2 682.33 689.38 -339.16 678.33
glmm 3 672.37 682.95 -333.18 666.37 11.959 1 0.000544 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
My best guess so far is: I used likelihood ratio test to compare the model with the fixed effect to a model without it. The model including the fixed effect (age-level) was a better fit ($\chi^2(df=?)=11.96$, $p=0.00054$).
And I cannot actually figure out how many df to report from that. The there is 2 for one model and 3 for the other, and 1 between them.
Thank you for your help.