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

问题描述

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

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`?

解决方法

类似这样的东西

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