probitmfx的回归输出

4xy9mtcn  于 2023-06-19  发布在  其他
关注(0)|答案(2)|浏览(82)

我试图为边际效应和probitmfx函数的p值生成一个很好的回归表,其中p值在每个协变量的边际效应下报告。我希望它看起来像的图片示例在这里Similar Output from Stata。我尝试了stargazer功能,如建议的here,但这似乎不工作,如果我没有一个OLS / probit。

data_T1 <- read_dta("xxx")
#specification (1)
T1_1 <- probitmfx(y ~ x1 + x2 + x3, data=data_T1)
#specification (1)
T1_2 <- probitmfx(y ~ x1 + x2 + x3 + x4 + x5, data=data_T1)
#this is what I tried but does not work
table1 <- stargazer(coef=list(T1_1$mfxest[,1], T1_2$mfxest[,1]),
p=list(T1_2$mfxest[,4],T1_2$mfxest[,4]), type="text")

有什么建议我如何设计这样的表在R?

rhfm7lfc

rhfm7lfc1#

你可以使用parameters包来生成一个漂亮的表:
代码:

library(mfx)
library(parameters)

# simulate some data
set.seed(12345)
n <- 1000
x <- rnorm(n)

# binary outcome
y <- ifelse(pnorm(1 + 0.5 * x + rnorm(n)) > 0.5, 1, 0)

data <- data.frame(y, x)
mod <- probitmfx(formula = y ~ x, data = data)

print_html(model_parameters(mod))

在Rmarkdown中使用的HTML表:

wfsdck30

wfsdck302#

您可以使用texregdocumentation):

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

相关问题