Next: Generalized linear models, Previous: Analysis of variance and model comparison, Up: Statistical models in R
函數 update() 是一個非常便利的函數。它允許
擬合一個比原先模型增加或減少一個項的模型。
它的形式是
> new.model <- update(old.model, new.formula)
在 new.formula 中,公式中包含的句點,
.,
僅僅表示“舊的公式模型中
的對應部分”。例如
> fm05 <- lm(y ~ x1 + x2 + x3 + x4 + x5, data = production)
> fm6 <- update(fm05, . ~ . + x6)
> smf6 <- update(fm6, sqrt(.) ~ .)
這將擬合從資料框 production 中得到的五個變數的
多重回歸模型,擬合額外增加一個變數的
六個回歸量的模型,和進一步對回應值進行
平方根轉換後的模型擬合。
注意參數 data= 在最開始調用模型擬合函數的時候指定。
這個資訊將會通過擬合模型物件傳遞給
函數 update() 及其相關者。
符號 . 同樣可以用在其他情況下,不過含義 有點不同。如
> fmfull <- lm(y ~ . , data = production)
它將會擬合回應量 y 對
資料框 production 中所有變數回歸的模型。
其他研究逐步回歸的函數是
add1(), drop1() 和 step()。
從字面上就可以看出這些函數的意義,更細節的內容可以參考
説明文檔。