[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