[R] Extracting overdispersion estimates from lmer amd glm objects

2006-09-12 Thread Toby Gardner
Dear list, I am needing to extract the estimate of overdispersion (deviance / residual degrees of freedom or c-hat) from multiple model objects - so they can then be used to compare the extent of overdispersion among alternative models as well as calculate qausi-AIC values. I have been

Re: [R] Extracting overdispersion estimates from lmer amd glm objects

2006-09-12 Thread Simon Blomberg
summary(modeltest)@sigma Toby Gardner wrote: Dear list, I am needing to extract the estimate of overdispersion (deviance / residual degrees of freedom or c-hat) from multiple model objects - so they can then be used to compare the extent of overdispersion among alternative models as

Re: [R] Extracting overdispersion estimates from lmer amd glm objects

2006-09-12 Thread David Barron
You can extract it from a summary.lmer object, which has a slot names sigma that contains this scale parameter. For example, sum.modeltest - summary(modeltest) [EMAIL PROTECTED] On 11/09/06, Toby Gardner [EMAIL PROTECTED] wrote: Dear list, I am needing to extract the estimate of

Re: [R] Extracting overdispersion estimates from lmer amd glm objects

2006-09-12 Thread David Barron
PS, the equivalent for a glm is: sum.modeltest - summary(modeltest) sum.modeltest$dispersion Hope this helps. On 11/09/06, Toby Gardner [EMAIL PROTECTED] wrote: Dear list, I am needing to extract the estimate of overdispersion (deviance / residual degrees of freedom or c-hat) from