我当前的数据集如下所示:
wage <- rnorm(10000)
educ <- rnorm(10000)
age <- rnorm(10000)
tce<- rnorm(10000)
work <- rbinom(n=10000,size=1,prob=0.05)
manu <- rbinom(n=10000,size=1,prob=0.05)
id <- rnorm(10000)
df <- data.frame(wage, educ, age, tce, work, manu, id)
wage
、work
和manu
是因变量,其余变量是自变量.
目前,我正在重复语法,但只是更改结果变量:
library(fixest)
model1 <- feols(work ~ educ + age + tce | id, data = df)
model2 <- feols(manu ~ educ + age + tce | id, data = df)
model2 <- feols(wage~ educ + age + tce | id, data = df)
有没有一种方法可以使用for循环来运行这种回归?
此外,在运行回归后,我还想绘制回归系数:
library(modelsummary)
modelplot(
list(model1, model2, model3)
)
然而,由于for循环不会创建新对象,我如何绘制系数?
非常感谢.