在 R 中复制 Stata 边距答案

【问题标题】：Replicate Stata margins in R在 R 中复制 Stata 边距
【发布时间】：2016-08-04 23:20:15
【问题描述】：

我有一个回归模型，其中包含一堆与回归变量交互的指标变量。在这种情况下，P 值可能适用于模型选择。我知道 Stata 有一个名为 margins 的命令，在这种情况下它真的很有帮助。示例：

margins, dydx(*)

来自relevant Stata manual

说明

边距是根据先前拟合模型在某些固定值下的预测计算得出的统计数据协变量和平均或以其他方式整合剩余的协变量。 margins 命令估计协变量和指定值的响应边际将结果显示为表格。能力包括估计边际均值、最小二乘均值、平均和条件边际效应和部分效应（可以报告为衍生品或弹性）、平均和条件调整预测和预测余量。

R中有没有类似的功能？

【问题讨论】：

没有理由不引用最新的手册版本（我写的是 14 版）。
首先查看effects 和lsmeans 包

标签： r replication regression stata margins

【解决方案1】：

我制作了一个 R 端口的 Stata 边距命令，直观地称为 margins。行为是您可能期望的：

library("margins")
x <- lm(mpg ~ cyl * hp + wt, data = mtcars)
summary(margins(x))
##  factor     AME     SE       z      p   lower   upper
##     cyl  0.0381 0.5999  0.0636 0.9493 -1.1376  1.2139
##      hp -0.0463 0.0145 -3.1909 0.0014 -0.0748 -0.0179
##      wt -3.1198 0.6613 -4.7176 0.0000 -4.4160 -1.8236

与Stata相比：

quietly reg mpg c.cyl##c.hp wt
margins, dydx(*)
------------------------------------------------------------------------------
             |            Delta-method
             |      dy/dx   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         cyl |   .0381376   .5998897     0.06   0.950    -1.192735     1.26901
          hp |  -.0463187    .014516    -3.19   0.004     -.076103   -.0165343
          wt |  -3.119815    .661322    -4.72   0.000    -4.476736   -1.762894
------------------------------------------------------------------------------

【讨论】：