【问题标题】:glmer reference outcome stringglmer 参考结果字符串
【发布时间】:2016-03-28 15:04:22
【问题描述】:

我正在运行多级逻辑回归,使用 R 中 lme4 包中的函数 glmer。我的二项式结果(或响应)变量编码为 cg

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

似乎默认引用是按字母顺序排列的第一个(在我的例子中是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获得了使用“重新编码”功能的提示

【问题讨论】:

    标签: r logistic-regression lme4 mixed-models


    【解决方案1】:

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

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

    根据 OP 的评论进行编辑- 回答您的问题:是的,因子水平默认按字母顺序排列。这位 R-Blogger 的 post 对此进行了更多讨论。

    【讨论】:

    • 根据您的建议,所有对数赔率估计值都成功反转了符号(例如,我的截距 -1.19 变为 1.19)。我还意识到g 是我最罕见的结果,这意味着c 确实被编码为0(我称之为默认值),g 被编码为1。换句话说:按字母顺序排列的结果确实是“默认”;我的问题的答案是“是”。
    猜你喜欢
    • 1970-01-01
    • 2019-06-19
    • 2020-10-04
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2016-09-27
    相关资源
    最近更新 更多