在 lme4 包中获取两个回归斜率之差的输出

Question

在下面的混合回归模型中，我得到了 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`?

Answer 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)

在 lme4 包中获取两个回归斜率之差的输出

getting an output for difference of two regression slopes in lme4 package

r

statistics

regression

lme4

multi-level