Hello --
I am comparing two
GLMs (binomial dependent variable)
, the results are the following:
> m1<-glm(symptoms ~ phq_index, data=data2)
> m2<-glm(symptoms ~ 1, data=data2)
Trying to compare these models using
> anova (m1, m2)
I do not obtain chi-square values or a chi-square difference test; instead,
I get loglikelihood ratios:
> Likelihood ratio tests of cumulative link models:
> formula: link: threshold:
> m2 sym_bin ~ 1 logit flexible
> m1 sym_bin ~ phq_index logit flexible
> no.par AIC logLik LR.stat df Pr(>Chisq)
> m2 1 10947 -5472.5
> m1 9 9711 -4846.5 1252 8 < 2.2e-16 ***
Since reviewers would like me to report chi-square values: how to I obtain
them when comparing GLMs? I'm looking for an output similar to the output
of the GLMER function in LME4, e.g.:
> anova(m3,m4)
...
> Df AIC BIC logLik Chisq Chi Df Pr(>Chisq)
> m3 13 11288 11393 -5630.9
> m4 21 11212 11382 -5584.9 92.02 8 < 2.2e-16 ***
Thank you!
Eiko
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
