Я выполняю логистическую регрессию с категориальной предикторной переменной с уровнями 0-6. По умолчанию R рассматривает уровень 0 в качестве контрольной группы.
Как я могу сказать R использовать, например, уровень 3 вместо уровня 0 как ссылку, не переименовывая уровни?
Используйте функцию C, чтобы определить контрасты в кадре данных.
Если ваш dataframe DF, а фактор-переменная fct, то
DF$fct <- C(DF$fct, contr.treatment, base=3)
(непроверенные).
Это очень легко с пакетами Фрэнка Харрелла (что позволяет использовать многие дополнительные функции).
Например, для 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 для вывода небольшого набора данных)
Функция C уже предложена, а также посмотрите на contrasts, relevel и reorder, среди прочих.
Вы можете использовать функцию release: dataframe $x1 < - release (datafrmae $x1, "введите здесь ссылочную категорию" )