On Fri, 9 Apr 2010, Sven Garbade wrote:
Hi list,
can anybody point me to the trick how glm is computing the dispersion
parameter in quasi-poisson regression, eg.
glm(...,family="quasipoisson")?
It's the sum of squared Pearson residuals divided by the residual degrees
of freedom. For example:
example("glm")
fm <- glm(counts ~ outcome + treatment, family = quasipoisson)
summary(fm)
sum(residuals(fm, type = "pearson")^2)/df.residual(fm)
hth,
Z
Thanks ®ards, Sven
______________________________________________
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.
______________________________________________
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.