如何解决在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)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。