R stat.anova GLM 方差分析统计

说明

这是一个实用函数，用于 anova(..., test != NULL) 的 lm 和 glm 方法，普通用户不应使用。

用法

stat.anova(table, test = c("Rao","LRT", "Chisq", "F", "Cp"),
           scale, df.scale, n)

参数

table	数字矩阵作为 anova.glm(..., test = NULL) 的结果。
test	一个字符串，部分匹配 "Rao" 、 "LRT" 、 "Chisq" 、 "F" 或 "Cp" 之一。
scale	用作 F 检验中分母的残差均方或其他尺度估计。
df.scale	对应于 scale 的自由度。
n	观察次数。

值

一个矩阵，它是原始的 table ，根据 test 参数添加一列测试统计数据。

例子

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

参考

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.

也可以看看

anova.lm、anova.glm。