我只是以一种良好、灵活的方式运行了一系列模型,强制执行数据代码分离。我的配置部分有一个很好的公式和模型列表,我可以使用它 lapply获取模型对象的列表。现在我想将它们显示在 stargazer, ,但它不需要列表对象。如何在不必键入每个列表元素的情况下执行此操作?

可重现的例子:

require(stargazer)
l <- list()
l$lm1 <- lm(rating ~ complaints + privileges + learning + raises + critical,
data=attitude)
l$lm2 <- lm(rating ~ complaints + privileges + learning, data=attitude)
## create an indicator dependent variable, and run a probit model
attitude$high.rating <- (attitude$rating > 70)
l$prbt <- glm(high.rating ~ learning + critical + advance, data=attitude,
family = binomial(link = "probit"))
stargazer( l[[1]], l[[2]], l[[3]], title="Results", align=TRUE, type="text")
有帮助吗?

解决方案

请确保您使用的是该软件包的最新版本。从版本 4.5.3 开始(自 2013 年 11 月起在 CRAN 上可用), stargazer 已经能够以您期望的方式接受对象列表:

stargazer(l, title="Results", align=TRUE, type="text")

其他提示

使用 do.call:

do.call( stargazer, l ) 

但是,这妨碍了以通常的方式传递参数:

> do.call( stargazer, l, type="text" )
Error in do.call(stargazer, l, type = "text") : 
  unused argument (type = "text")

因此,您必须将命名参数添加到列表中:

l$type <- "text"
l$align <- TRUE
l$title <- "Results"
do.call( stargazer, l )

另一种方法是柯里化 stargazer 函数:

require(functional)
sgCurried <- Curry( stargazer, type="text" ) # all arguments to stargazer go in here
do.call( sgCurried, l )
许可以下: CC-BY-SA归因
不隶属于 StackOverflow
scroll top