GLM Anova Statistics
This is a utility function, used in
glm methods for
anova(…, test != NULL)
and should not be used by the average user.
stat.anova(table, test = c("Rao","LRT", "Chisq", "F", "Cp"), scale, df.scale, n)
- numeric matrix as results from
anova.glm(…, test = NULL).
- a character string, partially matching one of
- a residual mean square or other scale estimate to be used as the denominator in an F test.
- degrees of freedom corresponding to
- number of observations.
A matrix which is the original
table, augmented by a column
of test statistics, depending on the
Hastie, T. J. and Pregibon, D. (1992) Generalized linear models. Chapter 6 of Statistical Models in S eds J. M. Chambers and T. J. Hastie, Wadsworth & Brooks/Cole.
##-- Continued from '?glm': print(ag <- anova(glm.D93)) stat.anova(ag$table, test = "Cp", scale = sum(resid(glm.D93, "pearson")^2)/4, df.scale = 4, n = 9)