R anova.lm 线性模型拟合的方差分析

说明

计算一个或多个线性模型拟合的方差分析表。

用法

## S3 method for class 'lm'
anova(object, ...)

## S3 method for class 'lmlist'
anova(object, ..., scale = 0, test = "F")

参数

object, ...	类 lm 的对象，通常是调用 lm 的结果。
test	指定要使用的检验统计量的字符串。可以是 "F" 、 "Chisq" 或 "Cp" 之一，允许部分匹配，或者 NULL 不进行测试。
scale	数字。噪声方差的估计\sigma^2。如果为零，则将从所考虑的最大模型中进行估计。

细节

指定单个对象会给出该拟合的方差表的顺序分析。也就是说，随着公式的每一项依次相加，残差平方和的减少量以表的行加上残差平方和的形式给出。

该表将包含将行均方与残差均方进行比较的 F 统计数据(和 P 值)。

如果指定了多个对象，则表中的一行用于显示每个模型的剩余自由度和平方和。对于除第一个模型之外的所有模型，还给出了自由度和平方和的变化。 (这仅在模型嵌套时才具有统计意义。)按照惯例，从最小到最大列出模型，但这取决于用户。

可选地，该表可以包括测试统计数据。通常，F 统计量是最合适的，它将一行的均方与所考虑的最大模型的残差平方和进行比较。如果指定scale，则可以使用卡方检验。 Mallows 的 C_p 统计量是残差平方和加上 \sigma^2 估计值的两倍乘以残差自由度。

值

类 "anova" 的对象继承自类 "data.frame" 。

警告

仅当两个或多个模型适合同一数据集时，它们之间的比较才有效。如果存在缺失值，这可能会出现问题R的默认值是na.action = na.omit被使用，并且anova.lmlist将检测到此错误。

例子

## sequential table
fit <- lm(sr ~ ., data = LifeCycleSavings)
anova(fit)

## same effect via separate models
fit0 <- lm(sr ~ 1, data = LifeCycleSavings)
fit1 <- update(fit0, . ~ . + pop15)
fit2 <- update(fit1, . ~ . + pop75)
fit3 <- update(fit2, . ~ . + dpi)
fit4 <- update(fit3, . ~ . + ddpi)
anova(fit0, fit1, fit2, fit3, fit4, test = "F")

anova(fit4, fit2, fit0, test = "F") # unconventional order

参考

Chambers, J. M. (1992) Linear models. Chapter 4 of Statistical Models in S eds J. M. Chambers and T. J. Hastie, Wadsworth & Brooks/Cole.

也可以看看

模型拟合函数lm、anova。

drop1 用于所谓的“II 型”方差分析，其中每个术语一次删除一个，尊重其层次结构。