如何在R中从线性模型对象获取全模型ANOVA表?
获取线性模型的全模型F检验ANOVA表
你可以通过对比全模型和仅含截距的零模型来得到全模型的F检验ANOVA表,直接用anova()函数就能实现:
# 拟合零模型(仅包含截距项) null_mod <- lm(y ~ 1) # 对比零模型与你的全模型,输出全模型F检验的ANOVA表 anova(null_mod, cube_mod)
或者更简洁的写法,无需单独存储零模型:
anova(lm(y ~ 1), cube_mod)
这个输出的结果就是检验x、x_2、x_3联合起来是否对y有显著影响的标准ANOVA表,和summary(cube_mod)里的F检验统计量、p值完全对应,但呈现为ANOVA表的格式。
如果只是需要提取全模型F检验的统计量和p值,也可以直接从模型的summary结果中获取:
# 获取F检验统计量 summary(cube_mod)$fstatistic # 计算对应的p值 pf(summary(cube_mod)$fstatistic[1], summary(cube_mod)$fstatistic[2], summary(cube_mod)$fstatistic[3], lower.tail = FALSE)
内容的提问来源于stack exchange,提问作者spencergw
相关产品推荐
相关产品推荐

