r - 相当于R中的对比度(Stata)
问题描述
我想知道R中是否有与Stata的contrast
命令等效的东西。
具体来说,我感兴趣sex@agegroup#group2
我知道margins
计算模型对象的边际(或部分)影响的包。
从Stata网页我提取了一个例子: https ://www.stata.com/features/overview/contrasts/
regress y agegroup##sex
contrast r.sex@agegroup
Contrasts of marginal linear predictions
Margins : asbalanced
df F P>F
sex@agegroup
(female vs male) 10-19 1 1.13 0.2945
(female vs male) 20-29 1 3.36 0.0743
(female vs male) 30-39 1 5.00 0.0310
(female vs male) 40-59 1 0.41 0.5279
(female vs male) 60-79 1 2.71 0.1076
Joint 5 2.52 0.0448
Denominator 40
Contrast Std. Err. [95% Conf. Interval]
sex@agegroup
(female vs male) 10-19 6.841855 6.441542 -6.176987 19.8607
(female vs male) 20-29 -11.80631 6.441542 -24.82515 1.212534
(female vs male) 30-39 -14.40607 6.441542 -27.42491 -1.387228
(female vs male) 40-59 -4.101691 6.441542 -17.12053 8.917151
(female vs male) 60-79 -10.60137 6.441542 -23.62022 2.417469
解决方案
感谢 Ben Bolker,我发现emmeans
并且可以复制 Stata 的对比命令结果:
library(emmeans)
lm.model= lm(y ~ agegroup*sex,data= df)
emmeans(lm.model, pairwise ~ sex*agegroup,simple = c("sex"))
推荐阅读
- mysql - 防止两个外键同时为 NULL
- python - 有没有一种更简洁的方法可以在 Python 中编写这个布尔比较?
- mern - 部署到heroku时找不到.scss文件
- python - Google Places API - 按非标准类别/类型搜索
- javascript - 无法从我的 create-react-app 中的“redux-saga”导入 createSagaMiddleware
- php - 使用 python 将数组打包为 16 位字符串
- reactjs - 自定义材质 UI 组件
- aws-cli - 如何使用 pyenv 使安装在一个环境中的工具在其他环境中可用
- azure - Azure 应用程序网关 DNS 返回 307 到后端池
- php - phpmyadmin 正在导出文本文件而不是 sql 文件