我一直在网上看到很多关于如何不建议循环的数据科学家的评论.但是,我最近发现自己处于使用一个有用的情况.我想知道以下过程是否有更好的替代方案(为什么替代方案会更好):
我需要运行一系列重复测量ANOVA,并且与您在下面看到的可重现示例类似地解决问题.
[我知道有关于运行多个ANOVA模型的其他问题,并且还有其他选项可用于这些类型的分析,但是现在我只想听听for循环的使用]
作为一个例子,四个重复测量ANOVA模型 - 四个因变量,每个变量分别在三个时间测量:
set.seed(1976)
code <- seq(1:60)
time <- rep(c(0,1,2), each = 20)
DV1 <- c(rnorm(20, 10, 2), rnorm(20, 10, 2), rnorm(20, 14, 2))
DV2 <- c(rnorm(20, 10, 2), rnorm(20, 10, 2), rnorm(20, 10, 2))
DV3 <- c(rnorm(20, 10, 2), rnorm(20, 10, 2), rnorm(20, 8, 2))
DV4 <- c(rnorm(20, 10, 2), rnorm(20, 10, 2), rnorm(20, 10, 2))
dat <- data.frame(code, time, DV1, DV2, DV3, DV4)
outANOVA <- list()
for (i in names(dat)) {
y <- dat[[i]]
outANOVA[i] …