8

我正在尝试为一组回归创建一个观星表,我在其中对我的数据子集运行每个回归。我认为,这样做的自然方法是使用split从我的数据创建一个 data.frames 列表,通过lapply在 data.frames 列表上使用创建一个 lm 对象列表,然后将该列表提供给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.

如果我手动将它们分开,则效果很好:

# 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")

但是对于我的真实数据,我要拆分的东西有几个值,我宁愿不手动拆分它们。为什么我得到“%错误:无法识别的对象类型”的任何想法。错误?

4

1 回答 1

8

有一个简单的解决方法,由 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)
于 2015-09-20T20:12:26.727 回答