我正在尝试为 probitmfx 函数的边际效应和 p 值生成一个很好的回归表,其中 p 值在每个协变量的边际效应下报告。我希望它看起来像的图片示例是来自 Stata 的类似输出。我按照这里的建议尝试了观星功能,但如果我没有 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 中设计这样的表格吗?