在逻辑回归中更改分类预测变量的参考组

Question

我正在使用具有0-6级别的分类预测变量运行逻辑回归.默认情况下,R将0级视为参考组.

如何在不重命名级别的情况下告诉R使用级别3而不是级别0作为参考？

Answer 1

使用此C功能可在数据框中定义对比度.

如果您的数据框是DF,因子变量是fct,那么

DF$fct <- C(DF$fct, contr.treatment, base=3)

(未经测试).

Answer 2

使用Frank Harrell的软件包(允许许多其他功能)非常容易.
例如,对于名为"df"的数据帧

library(Hmisc) 
library(rms)
dd=datadist(df)
options(datadist='dd')
(m=lrm(y ~ catvar, data=df)) #uses the mode as the reference group
summary(m, catvar=3) #using level=3 as the reference group
summary(m, catvar=0) #using level=0 as the reference group

(并且,一些示例数据本来不错 - 您总是可以dput用来输出一个小数据集)

Answer 3

该C功能已经被建议，还可以查看contrasts、relevel、 和reorder等。