2016-03-28 2 views
1

Я управляю многоуровневой логистической регрессией, используя функцию glmer из пакета lme4 в R. Моя переменная биномиального исхода (или ответа) кодируется как и g.строка результата ответа glmer

Мой вопрос: как я могу узнать, какой уровень переменной результата принимается за исходный результат?

Кажется, ссылка по умолчанию является первой в алфавитном порядке (так c в моем случае), это правда?

Я предполагаю, что это после того, как с помощью levels и relevel функции:

levels(data$Outcome) 
# [1] "c" "g" 

test <- relevel(data$Outcome, ref = "g") 
# levels(test) 
# [1] "g" "c" 

Это, как бы сказать, что c была ссылка раньше, но это скорее косвенными.

Мои домашние задания:
Я не нашел ответ, используя ?glmer или онлайн руководство по эксплуатации в формате PDF lme4, или эти связанные сообщения:

Для случая 0/1 кодированной ответ:
Using glmer for logistic regression, how to verify response reference

Got наконечник с помощью функции «перекодировать» от:
Logistic regression - defining reference level in R

ответ

0

Вы должны изменить свой контрольный уровень. This post демонстрирует, как это сделать. В вашем случае написать это так:

data$Outcome <- factor(data$Outcome, levels = c("g", "c")) 

Редактировать основанный на OP-х comment- Чтобы ответить на ваш вопрос: Да, уровни фактора расположены в алфавитном порядке по умолчанию. Это post этого R-Blogger обсуждает его больше.

+0

Следуя вашему предложению, грубо изменился знак для всех оценок логарифмов (например, мой перехват -1.19 стал 1.19). Я также понял, что 'g' - мой самый редкий результат, означающий, что' c' действительно был закодирован как '0' (который я назвал дефолтом) и' g' как '1'. Другими словами: результат, который является первым в алфавитном порядке, действительно является «дефолтом»; ответ на мой вопрос: «ДА». – UncleVinnie

Смежные вопросы