glmer 参考结果字符串
glmer reference outcome string
我是运行多级逻辑回归,使用R中包lme4
中的函数glmer
。我的二项式结果(或响应-)变量编码为c
和 g
.
我的问题是:如何知道结果变量的哪个水平作为参考结果?
似乎默认引用是按字母顺序排列的第一个(所以 c
在我的例子中),这是真的吗?
我在使用 levels
和 relevel
函数后猜测:
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
从以下位置获得了使用 'recode' 函数的提示:
Logistic regression - defining reference level in R
您需要更改参考级别的顺序。 This post 演示了如何操作。在你的情况下这样写:
data$Outcome <- factor(data$Outcome, levels = c("g", "c"))
根据 OP 的评论进行编辑- 回答您的问题:是的,因子水平默认按字母顺序排列。这位 R-Blogger 的 post 对此进行了更多讨论。
我是运行多级逻辑回归,使用R中包lme4
中的函数glmer
。我的二项式结果(或响应-)变量编码为c
和 g
.
我的问题是:如何知道结果变量的哪个水平作为参考结果?
似乎默认引用是按字母顺序排列的第一个(所以 c
在我的例子中),这是真的吗?
我在使用 levels
和 relevel
函数后猜测:
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
从以下位置获得了使用 'recode' 函数的提示:
Logistic regression - defining reference level in R
您需要更改参考级别的顺序。 This post 演示了如何操作。在你的情况下这样写:
data$Outcome <- factor(data$Outcome, levels = c("g", "c"))
根据 OP 的评论进行编辑- 回答您的问题:是的,因子水平默认按字母顺序排列。这位 R-Blogger 的 post 对此进行了更多讨论。