Next: , Previous: Generic functions for extracting model information, Up: Statistical models in R



11.4 方差分析和模型比較

模型擬合函數 aov(formula, data=data.frame) 和函數 lm() 非常的相似,在Generic functions for extracting model information 部分列出的泛型函數同樣適用。

需要注意的是 aov() 還允許分析 多方差層次的模型,如裂區實驗設計, 利用區組內資訊進行的平衡不完全區組設計(balanced incomplete block design)等。 模型公式

     response ~ mean.formula + Error(strata.formula)

strata.formula 定義了一個多層次實驗設計中 多個方差因素。最簡單的情況是,strata.formula 是單因素的。 它定義了一個雙層次的實驗,也就是研究在這些因數的水準內 或者水準間的實驗回應。

例如,考慮到所有的決定變數因數,模型公式可以 設計如下:

     > fm <- aov(yield ~ v + n*p*k + Error(farms/blocks), data=farm.data)

這常常用來描述一個同時含有均值模型 v + n*p*k 和三個誤差層次(“農田之間”, “農田內但區組之間”和“區組內”)的實驗。