Изменение ссылочной группы для категориальной предикторной переменной в логистической регрессии

Я выполняю логистическую регрессию с категориальной предикторной переменной с уровнями 0-6. По умолчанию R рассматривает уровень 0 в качестве контрольной группы.

Как я могу сказать R использовать, например, уровень 3 вместо уровня 0 как ссылку, не переименовывая уровни?

Ответ 1

Используйте функцию C, чтобы определить контрасты в кадре данных.

Если ваш dataframe DF, а фактор-переменная fct, то

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

(непроверенные).

Ответ 2

Это очень легко с пакетами Фрэнка Харрелла (что позволяет использовать многие дополнительные функции).
Например, для dataframe с именем '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 для вывода небольшого набора данных)

Ответ 3

Функция C уже предложена, а также посмотрите на contrasts, relevel и reorder, среди прочих.

Ответ 4

Вы можете использовать функцию release: dataframe $x1 < - release (datafrmae $x1, "введите здесь ссылочную категорию" )