r - 在 R 的 for 循环中使用字符串 - 将字符串作为 for 循环中的方差分析中的变量包含在内时出错

Question

我在 R 中有以下方差分析，效果很好:

fit <- aov(dependent1 ~ X + Z + X*Z, data=dataset)
drop1(fit,~.,test="F")

“dependent1”、“X”和“Z”是列名称。

我想创建一个 for 循环，在其中循环一定数量的因变量，我尝试了以下方法:

dependent_variables <- c("dependent1", "dependent2", "dependent3")

for (i in dependent_variables) {
  fit <- aov(i ~ X + Z + X*Z, data=dataset)
  drop1(fit,~.,test="F")
}

如果我运行此命令，我会收到一条错误消息:

Error in model.frame.default(formula = i ~ X + Z + X *  : 
  variable lengths differ (found for 'X')

知道这里出了什么问题吗？

Answer 1

示例数据(可能满足也可能不满足方差分析的标准)

X <- rnorm(100)
Z <- rnorm(100)
dependent1 <- rnorm(100)
dependent2 <- rnorm(100)
dependent3 <- rnorm(100)

dataset <- cbind(data.frame(X, Z, dependent1, dependent2, dependent3))

以下脚本可以工作，您需要输入因变量的行列号:

for (i in 3:5) {
fit <- aov(dataset[ , i] ~ X + Z + X*Z, data=dataset)
drop <- drop1(fit,~.,test="F")
print(fit)
print(drop)
}