我想对存储在我的工作目录中的多个数据集执行anova。到目前为止，我想出了：
files <- list.files(pattern = ".csv")
for (i in seq_along(files)) {
mydataset.i <- files[i]
AnovaModel.1 <- aov(DES ~ DOSE, data=mydataset.i)
summary(AnovaModel.1)
}
如您所见，我对循环非常陌生，无法完成这项工作。我也知道我需要添加一个代码以将所有摘要输出附加到一个文件中。我将不胜感激您可以提供任何帮助来指导工作循环，该循环可以在目录中的多个 .csv 文件（相同的标题）上执行 anovas 并生成记录输出。
【问题讨论】：
您可能希望将list.files 与full.names = TRUE 一起使用，以防您不在同一条路径上。
files <- list.files("path_to_my_dir", pattern="*.csv", full.names = T)
# use lapply to loop over all files
out <- lapply(1:length(files), function(idx) {
# read the file
this.data <- read.csv(files[idx], header = TRUE) # choose TRUE/FALSE accordingly
aov.mod <- aov(DES ~ DOSE, data = this.data)
# if you want just the summary as object of summary.aov class
summary(aov.mod)
# if you require it as a matrix, comment the previous line and uncomment the one below
# as.matrix(summary(aov.mod)[[1]])
})
head(out)
这应该给你一个list，列表中的每个条目都有一个summary matrix，与输入文件列表相同的顺序。
【讨论】：
你的错误是你的循环没有加载你的数据。您的文件名列表在“文件”中，然后您开始浏览该列表并将 mydataset.i 设置为与您的迭代器 i 匹配的文件的名称...但是您尝试在文件名上运行 aov存储在 mydataset.i 中！
您正在寻找将输出重定向到文件的命令是 sink。考虑以下几点：
sink("FileOfResults.txt") #starting the redirect to the file
files <- list.files("path_to_my_dir", pattern="*.csv", full.names = T) #using the fuller code from Arun
for (i in seq_along(files)){
mydataset.i <- files[i]
mydataset.d <- read.csv(mydataset.i) #this line is new
AnovaModel.1 <- aov(DES ~ DOSE, data=mydataset.d) #this line is modified
print(summary(AnovaModel.1))
}
sink() #ending the redirect to the file
我更喜欢 Arun 的这种方法，因为结果直接存储到文件中，无需跳过列表，然后必须弄清楚如何以可读的方式将列表存储到文件中。
【讨论】：