在 lme4 包中获得两个回归斜率差异的输出答案

【问题标题】：getting an output for difference of two regression slopes in lme4 package在 lme4 包中获得两个回归斜率差异的输出
【发布时间】：2021-02-11 13:07:42
【问题描述】：

在下面的混合回归模型中，我得到了meanses 的输出和ses 估计值的输出。是否也可以在此模型中为meanses - ses (3.675037 - 2.191165) 获得输出（例如Std. Error）（我愿意使用任何包）？

library(lme4)

hsb <- read.csv('https://raw.githubusercontent.com/rnorouzian/e/master/hsb.csv')

fit <- lmer(math ~ ses + meanses + (1|sch.id), data = hsb)

coef(summary(fit))

             Estimate Std. Error   t value
(Intercept) 12.661262  0.1493726 84.762956
ses          2.191165  0.1086673 20.163983
meanses      3.675037  0.3776607  9.731055 # can we have `Std. Error` for `meanses - ses`?

【问题讨论】：

检查this。 multcomp 包可用于线性混合效果模型，如包说明中所述。

标签： r statistics regression lme4 multi-level

【解决方案1】：

类似的东西

> summary(multcomp::glht(fit, "ses - meanses = 0"))

     Simultaneous Tests for General Linear Hypotheses

Fit: lmer(formula = math ~ ses + meanses + (1 | sch.id), data = hsb)

Linear Hypotheses:
                   Estimate Std. Error z value Pr(>|z|)    
ses - meanses == 0   -1.484      0.422  -3.517 0.000437 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Adjusted p values reported -- single-step method)

【讨论】：