Question

我正在尝试编写一个可以应用于字符串向量或列表而不是编写循环的函数。我的目标是对不同的内生变量运行回归并保存结果表。由于有经验的R用户告诉我们应该学习应用函数，我想试一试。这是我的尝试：

破碎的例子：

library(ExtremeBounds)
Data <- data.frame(var1=rbinom(30,1,0.2),var2=rbinom(30,1,0.2),var3=rnorm(30),var4=rnorm(30),var5=rnorm(30),var6=rnorm(30))

spec1 <- list(y=c("var1"),freevars=("var3"),doubtvars=c("var4","var5"))
spec2 <- list(y=c("var2"),freevars=("var4"),doubtvars=c("var3","var5","var6"))
specs <- c("spec1","spec2")

myfunction <- function(x){
eba <- eba(data=Data, y=x$y, 
       free=x$freevars, 
       doubtful=x$doubtvars, 
       reg.fun=glm, k=1, vif=7, draws=50, se.fun = se.robust, weights = "lri", family = binomial(logit))

output <- eba$bounds
output <- output[,-(3:7)]
}

lapply(specs,myfunction)

这给了我一个错误，让我猜出R在x应该"spec1"或"spec2"时无法理解。此外，我不太明白lapply会在这里尝试收集什么。你能否向我提供一些最佳实践/提示如何将这些事情传达给R？

错误：Error in x$y : $ operator is invalid for atomic vectors

工作示例：

这是spec1的一个工作示例，不使用应用程序，显示我正在尝试做的事情。我想通过7个规范循环这个例子，但我正试图摆脱循环。输出不必保存为csv，所有输出或任何其他集合的列表都会很棒！

eba <- eba(data=Data, y=spec1$y, 
       free=spec1$freevars, 
       doubtful=spec1$doubtvars, 
       reg.fun=glm, k=1, vif=7, draws=50, se.fun = se.robust, weights = "lri", family = binomial(logit))
output <- eba$bounds
output <- output[,-(3:7)]
write.csv(output, "./Results/eba_pmr.csv")

Answer 1

根据@ user20650的评论，解决方案非常简单：

在lapply命令中，使用lapply(mget(specs),myfunction)获取规范列表元素的名称而不是列表本身。

或者，可以将specs定义为列表：specs <- list(spec1,spec2)但是有一个缺点，即lapply命令将返回不同规范编号的列表。第一个版本保留了规范的名称（spec1和spec2），这使得生成列表的工作变得更加容易。

R：关于函数何时理解其论点的一些提示？

1 个答案: