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

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等。

Answer 4

您可以使用 relevel 函数： dataframe$x1 <- relevel(datafrmae$x1, "type here the reference category")