R dummy.coef 提取原始編碼中的係數

說明

這根據係數的原始級別而不是編碼變量來提取係數。

用法

dummy.coef(object, ...)

## S3 method for class 'lm'
dummy.coef(object, use.na = FALSE, ...)

## S3 method for class 'aovlist'
dummy.coef(object, use.na = FALSE, ...)

參數

object	線性模型擬合。
use.na	奇異模型中係數的邏輯標誌。如果use.na為真，則未確定的係數將丟失；如果為 false，他們將獲得一個可能的值。
...	傳遞給其他方法或從其他方法傳遞的參數。

細節

擬合線性模型具有因子項對比係數，其數量通常比級別數少 1。該函數重新表達了原始編碼中的係數；由於係數將在減少的基礎上進行擬合，因此將遵守任何隱含的約束(例如 contr.helmert 或 contr.sum 的零和)。使用 dummy.coef 進行 contr.treatment 對比沒有什麽意義，因為根據定義，缺失係數為零。

使用的方法有一些限製，並且對於 poly(x, 2) 等術語會給出不完整的結果。然而，它足以滿足其主要目的，aov 模型。

值

給出每項係數值的列表。對於多層 aov 模型，每個層都有這樣的列表。

警告

此函數用於人工檢查輸出：不應將其用於計算。使用編碼變量進行所有計算。

對於奇異值，結果與 S 不同，其中 S 可能不正確。

例子

options(contrasts = c("contr.helmert", "contr.poly"))
## From Venables and Ripley (2002) p.165.
npk.aov <- aov(yield ~ block + N*P*K, npk)
dummy.coef(npk.aov)

npk.aovE <- aov(yield ~  N*P*K + Error(block), npk)
dummy.coef(npk.aovE)

也可以看看

aov , model.tables