我正在使用 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",...)）