将变量输入回归函数

Question

我有这个 feature_list，其中包含几个可能的值，例如“A”、“B”、“C”等。并且 time_list 中有时间。

因此，我将有一个循环，我将在其中遍历每个不同的值并将其放入公式中。

类似 for(i in ...) ，然后 my_feature <- feature_list[i] 和 my_time <- time_list[i] 。

然后我将时间和所选特征放入数据帧中以用于回归

feature_list<- c("GPRS")
time_list<-c("time")
calc<-0

feature_dim <- length(feature_list)
time_dim <- length(time_list)

data <- read.csv("data.csv", header = TRUE, sep = ";")
result <- matrix(nrow=0, ncol=5)
errors<-matrix(nrow=0, ncol=3)

for(i in 1:feature_dim) {
    my_feature <- feature_list[i]
    my_time <- time_list[i]

    fitdata <- data.frame(data[my_feature], data[my_time])

    for(j in 1:60) {

        my_b <- 0.0001 * (2^j)

        for(k in 1:60) {
            my_c <- 0.0001 * (2^k)
            cat("Feature: ", my_feature, "\t")
            cat("b: ", my_b, "\t")
            cat("c: ", my_c, "\n")

            err <- try(nlsfit <- nls(GPRS ~ 53E5*exp(-1*b*exp(-1*c*time)), data=fitdata, start=list(b=my_b, c=my_c)), silent=TRUE)
            calc<-calc+1

            if(class(err) == "try-error") {
                next
            }

            else {
                coefs<-coef(nlsfit)
                ess<-deviance(nlsfit)
                result<-rbind(result, c(coefs[1], coefs[2], ess, my_b, my_c))
            }
    }
} 
}

现在在 nls() 调用中我希望能够调用 my_feature 而不仅仅是“A”或“B”或某件事，然后到列表中的下一个。但我在那里遇到错误。我做错了什么？

Answer 1

您可以使用粘贴创建公式的字符串版本，包括所需的变量名称，然后使用 as.formula 或 formula 函数将其转换为要传递的公式到 nls。

as.formula(paste(my_feature, "~ 53E5*exp(-1*b*exp(-1*c*time))"))

另一种选择是使用 bquote 函数将变量名称插入到函数调用中，然后评估函数调用。

Answer 2

我不久前使用过 R，也许你可以尝试一下：

你想要的是创建一个包含变量列表的公式，对吧？

因此，如果响应变量是列表的第一个元素，其他元素是解释变量，您可以这样创建公式：

my_feature[0] ~ reduce("+",my_feature[1:]) 。这可能有效。

这样您就可以创建取决于 my_features 中的变量的公式。