12

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

例如,我如何告诉 R 使用级别 3 而不是级别 0 作为参考而不重命名级别?

4

4 回答 4

9

使用该C函数在数据框中定义您的对比。

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

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

(未经测试)。

于 2012-06-12T19:20:15.823 回答
5

使用 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用来输出一个小数据集)

于 2012-06-12T19:20:05.290 回答
4

C功能已被建议,还请查看contrastsrelevelreorder等。

于 2012-06-12T20:10:04.397 回答
2

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

于 2015-05-31T06:44:50.177 回答