Latex 的 Probit 边际效应输出

Question

我正在计算 R mfx 包的概率边际效应。我想为边际效应输出生成 Latex 代码。我尝试了 OLS 和 probit 系数的 stargazer 包，它对两者都适用，但是对于 probit 边际效应（通过使用 probitmfx 命令）它不起作用。

请在这方面帮助我，谢谢。

Answer 1

我还没有找到一个简单的解决方案，但这是我使用的技巧：

ols <- lm(your_formula, data=your_data, family="binomial")
probit <- glm(your_formula, data=your_data, family="binomial")
probit_margins <- probitmfx(your_formula, your_data, atmean=FALSE)$mfxest
your_table <- stargazer(ols, probit, coef = list(NULL, probit_margins[,1]), 
                                     se = list(NULL, probit_margins[,2]))

它与 logit 的工作原理相同，您不一定需要 atmean=FALSE 选项。我在示例中包含了一个 OLS，以向您展示如何包含 stargazer 可以处理的其他规范，而无需指定系数和标准误差。

Answer 2

类比bixiou的回答，可能稍微优雅一点（至少我们不适合两次）：

ols <- lm(your_formula, data=your_data, family="binomial")
probit_margins <- probitmfx(your_formula, your_data, atmean=FALSE)

your_table <- stargazer(ols, probit_margins$fit,
coef = list(NULL, probit_margins$mfxest[,1]),
se = list(NULL, probit_margins$mfxest[,2]))

这利用了 probitmfx 对象的 fit 属性是原始模型这一事实。

Latex 的 Probit 边际效应输出

Probit Marginal Effects output for Latex

r

marginal-effects