Re: [R] Estimating QAIC using glm with the quasibinomial family

2008-05-09 Thread Darren Gillis
I am glad to have joined such an active online community! I assume from your R code that the deviance reported by R for quasi... families is not quasi-deviance and that the following function summarizes the calculation (note that I have added a parameter to k to account for the estimation of

Re: [R] Estimating QAIC using glm with the quasibinomial family

2008-05-07 Thread Darren Gillis
PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Kunio takezawa Sent: Tuesday, May 06, 2008 8:41 PM To: r-help@r-project.org Subject: Re: [R] Estimating QAIC using glm with the quasibinomial family R-users E-mail: r-help@r-project.org My question is: Will this calculation be valid with the residual

Re: [R] Estimating QAIC using glm with the quasibinomial family

2008-05-07 Thread Kunio takezawa
R-users E-mail: r-help@r-project.org but I was specificly interested in calculating QAIC and QAICc from a glm fitted with the family=quasibinomial option. If you use family=quasibinomial(link = logit) in glm(), the program will be: function () { xx - c(1,2,3,4,5,6,7,8,9,10) yy -

Re: [R] Estimating QAIC using glm with the quasibinomial family

2008-05-06 Thread Kunio takezawa
R-users E-mail: r-help@r-project.org My question is: Will this calculation be valid with the residual deviance returned by the glm() function using the quasibinomial family as reported in R? Let me show you a simple example, assuming c=2.5: function () { xx - c(1,2,3,4,5,6,7,8,9,10) yy