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