Question

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

请帮助我，谢谢。

Answer 1

类似于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属性是原始模型这一事实。

Answer 2

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

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

乳胶的Probit边际效应输出

2 个答案: