Latex 的 Probit Marginal Effects 输出答案

【问题标题】：Probit Marginal Effects output for LatexLatex 的 Probit Marginal Effects 输出
【发布时间】：2017-07-12 05:52:25
【问题描述】：

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

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

【问题讨论】：

xtable 和 Hmisc 包都具有输出 Latex 表格的功能，但您可能需要将矩阵放在一起以手动输出。
谢谢@Imo，你能举个小例子，这样我就可以继续了。
只需滚动浏览these posts。你应该会找到很多有用的例子。
谢谢，我浏览了所有链接，但仍然没有成功。
我设法在这里找到了答案 [link] (r-bloggers.com/porting-stata-like-marginal-effects-to-latex)

标签： r marginal-effects

【解决方案1】：

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

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 可以处理的其他规范，而无需指定系数和标准误差。

【讨论】：

如果您有两个概率模型并希望将它们的边际效应包括在表中，您会怎么做？谢谢！

【解决方案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 属性是原始模型这一事实。

【讨论】：