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 型”方差分析，其中每個術語一次刪除一個，尊重其層次結構。