我正在尝试为边际效应制作一个很好的回归表&来自probitmfx函数的p值,其中p值在每个协变量的边际效应下报告。我希望它看起来是什么样子的一个图片示例是在这里Similar Output from Stata。我尝试了观星功能,如建议的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中设计这样一个表吗?
发布于 2021-02-23 18:20:46
您可能可以使用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表:
https://stackoverflow.com/questions/66337650
复制相似问题