Probit Marginal Effects вывод для латекса

Question

Я вычисляю предельные эффекты пробита от R mfx пакет. Я хочу сгенерировать латексный код для вывода предельных эффектов. Я старался stargazer пакет для OLS и пробитовых коэффициентов, он отлично работает для обоих, однако для пробитовых предельных эффектов (используя 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 Атрибут оригинальной модели.

Answer 3

Вы можете использоватьtexreg(документация ):

      model <- probitmfx(formula = y ~ x, data = your_data)
texreg(model)