[R] quasibinomial glm
Stephane Dray
dray at biomserv.univ-lyon1.fr
Thu Jan 10 14:52:45 CET 2002
Hello list,
i have a glm with family=binomial, link=logit but there is
over-dispersion. So, in order to take into account for this problem i
choose to do a glm with family=quasibinomial(). I'm not an expert on
this subject and i ask if someone could validate my approach (i'm not
sure for the tests) :
quasi_glm(myformula,quasibinomial(),start=mystart)
summary(quasi) # test t for coefficients
anova(quasi,test="F") # F test (quasi glm) instead of Chisq test
(classical glm)
Thanks in advance.
--
Stéphane DRAY
---------------------------------------------------------------
Biométrie et Biologie évolutive - Equipe "Écologie Statistique"
Universite Lyon 1 - Bat 711 - 69622 Villeurbanne CEDEX - France
Tel : 04 72 43 27 56 Fax : 04 78 89 27 19
04 72 43 27 57 E-mail : dray at biomserv.univ-lyon1.fr
---------------------------------------------------------------
ADE-4 http://pbil.univ-lyon1.fr/ADE-4/ADE-4F.html
---------------------------------------------------------------
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list