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) fro

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 o

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

[R] Extracting overdispersion estimates from lmer amd glm objects

2006-09-11 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 unable