R C 设置因子的对比

说明

设置因子的 "contrasts" 属性。

用法

C(object, contr, how.many, ...)

参数

object	因子或有序因子
contr	这与使用形成对比。可以是一个矩阵，每个因子水平一行，或者一个合适的函数，如 contr.poly 或给出函数名称的字符串
how.many	要设置的对比度数量，默认情况下比 nlevels(object) 少 1。
...	函数 contr 的附加参数。

细节

为了与 S 兼容，contr 可以是 treatment 、 helmert 、sum 或 poly (不带引号)作为 contr.treatment 的简写，依此类推。

值

设置了"contrasts" 属性的因子object。

例子

## reset contrasts to defaults
options(contrasts = c("contr.treatment", "contr.poly"))
tens <- with(warpbreaks, C(tension, poly, 1))
attributes(tens)
## tension SHOULD be an ordered factor, but as it is not we can use
aov(breaks ~ wool + tens + tension, data = warpbreaks)

## show the use of ...  The default contrast is contr.treatment here
summary(lm(breaks ~ wool + C(tension, base = 2), data = warpbreaks))


# following on from help(esoph)
model3 <- glm(cbind(ncases, ncontrols) ~ agegp + C(tobgp, , 1) +
     C(alcgp, , 1), data = esoph, family = binomial())
summary(model3)

参考

Chambers, J. M. and Hastie, T. J. (1992) Statistical models. Chapter 2 of Statistical Models in S eds J. M. Chambers and T. J. Hastie, Wadsworth & Brooks/Cole.

也可以看看

contrasts 、contr.sum 等