我正在尝试为一组回归创建观星表,在该表中,我对数据的一个子集运行了每个回归。我认为,执行此操作的自然方法是使用split从我的数据创建data.frames列表,使用data.frames列表创建lm对象lapply列表,然后将该列表提供给stargazer。例如,
library(MASS)
library(stargazer)
data(Boston)
# This doesn't work
by.river <- split(Boston, Boston$chas)
fit <- lapply(by.river, lm, formula = crim ~ indus)
stargazer(fit, type = "text")
# % Error: Unrecognized object type.
# % Error: Unrecognized object type.
Run Code Online (Sandbox Code Playgroud)
如果我将它们手动分割,效果很好:
# This works
fit2 <- vector(mode = "list", length = 2)
fit2[[1]] <- lm(crim ~ indus, data = Boston, subset = (chas == 0))
fit2[[2]] <- lm(crim ~ indus, data = Boston, subset = (chas == 1))
stargazer(fit2, type = "text")
Run Code Online (Sandbox Code Playgroud)
但是根据我的真实数据,我要拆分的东西具有多个值,我宁愿不手工拆分它们。为什么我得到“%错误:无法识别的对象类型”的任何想法。错误?
有一个简单的解决方法,由BondedDust暗示,并通过仔细阅读的帮助来建议lapply。
fit <- lapply(by.river, function(dd)lm(crim ~ indus,data=dd))
stargazer(fit, type = "text")
fit[[1]]$call
#lm(formula = crim ~ indus, data = dd)
Run Code Online (Sandbox Code Playgroud)