如何解决我正在使用 R 中的 ... 作为我正在制作的函数中的控制变量,如何将控制变量转换为列表?
我有一个进行地理加权回归的函数,我有我的形状文件,x
用于 x 变量,y
用于 y 变量,...
用于控制变量.例如,我想先调用回归但有多个控制变量 - 我该怎么做?
GWR.function <- function(shape1,x,y,...) {
model <- lm(shape1[[x]] ~ shape1[[y]] + shape1[[...]])
return(summary(model))
}
我认为它需要使用列表,但我不知道如何做到这一点。
解决方法
通常,do.call()
是将列表作为函数参数引入的方式。在这里,我将 y
与 ...
结合,并将结果 character
向量传递给 reformulate
。如果您想处理 ...
的单个元素,list(...)
会将其转换为“常规”列表。
## set up example
dd <- data.frame(x=1:10,y=rnorm(10),z=rnorm(10))
GWR.function <- function(shape1,x,y,...) {
predvars <- do.call("c",c(list(y),...))
form <- reformulate(predvars,response=x)
model <- lm(form,data=shape1)
return(summary(model))
}
GWR.function(dd,"x","y","z")
您可以将函数重新定义为 function(shape1,...)
以稍微简化代码(然后是 predvars <- do.call("c",...)
)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。