2016-03-28 123 views
1

我使用lme4中的函数glmer运行多级逻辑回归。我的二项结果(或响应)变量编码为cgglmer参考结果字符串

我的问题是:我怎么知道结果变量的哪个级别作为参考结果?

它似乎默认的引用是第一个字母(所以c在我的情况),这是真的吗?

我使用的levelsrelevel功能后猜测这:

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

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

这似乎是说,c是之前的参考,但是这是比较间接的。

我的家庭作业:
我没有找到使用?glmer答案,或lme4在线PDF手册,或者这些相关文章:

对于0/1编码的响应的情况下:
Using glmer for logistic regression, how to verify response reference

得到使用“重新编码”功能,从尖:
Logistic regression - defining reference level in R

回答

0

您需要更改参考级别的顺序。 This post演示如何做。在你的情况写这样说:

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

编辑基于OP的comment-要回答你的问题:是的,因子水平在默认情况下按字母顺序排列。此R-Blogger的post对此进行了更多讨论。

+0

根据您的建议,所有对数赔率估计值都会完全颠倒符号(例如,我的拦截-1.19变为1.19)。我也意识到'g'是我最罕见的结果,这意味着'c'确实被编码为'0'(我将其称为默认值),'g'编码为'1'。换句话说,首先按字母顺序排列的结果确实是'默认';我的问题的答案是'是'。 – UncleVinnie