我只是以一种很好的,灵活的方式运行了一系列模型,强制执行数据代码分离。我在配置部分中有一个很好的公式和模型列表,我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")
答案 0 :(得分:6)
请确保您使用的是最新版本的软件包。从版本4.5.3开始(自2013年11月起在CRAN上可用),stargazer
已经能够以您期望的方式接受对象列表:
stargazer(l, title="Results", align=TRUE, type="text")
答案 1 :(得分:4)
使用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 )
另一种方法是讨论观星者的功能:
require(functional)
sgCurried <- Curry( stargazer, type="text" ) # all arguments to stargazer go in here
do.call( sgCurried, l )