从 R 中的 lme4 对象获取单个斜率答案

【问题标题】：Obtaining individual slopes from an lme4 object in R从 R 中的 lme4 对象获取单个斜率
【发布时间】：2020-08-05 23:36:19
【问题描述】：

我是 R 中 lme4 包的新手。在下面的示例中，我想知道是否有可能在安装我的 @987654325 后为每个 dep 获得 gender 斜率（即差异） @模型？

dat <- data.frame(dep = rep(LETTERS[1:6],each=2), gender = rep(c("Ma","Fe"),6), 
       admit=c(512,89,353,17,120,202,138,131,53,94,22,24), 
       reject=c(313,19,207,8,205,391,279,244,138,299,351,317))

lme4::glmer(cbind(admit,reject) ~ gender+dep + (gender|dep), data=dat, family=binomial)

【问题讨论】：

对于具有单个斜坡的固定效果，您可能需要查看feisr 包。

标签： r regression logistic-regression lme4 mixed-models

【解决方案1】：

在 lme4 中，您可以从 ranef 获得估计的斜率，但在您的模型中，您需要对全局和单位特定项求和，如下例所示。

library(lme4)

dat <- data.frame(dep = rep(LETTERS[1:6],each=2), gender = rep(c("Ma","Fe"),6), 
                  admit=c(512,89,353,17,120,202,138,131,53,94,22,24), 
                  reject=c(313,19,207,8,205,391,279,244,138,299,351,317))

mod1 <- glmer(cbind(admit,reject) ~ gender+dep + (gender|dep), data=dat, family=binomial)
summary(mod1)

ran_gender <- ranef(mod1)$dep

fe_mod1 <- fixef(mod1)

slopes <- fe_mod1[[2]] + ran_gender[,2]
slopes

【讨论】：