用于分析许多文件的 for 循环

Question

我正在尝试使用 for 循环编写一个 R 脚本来执行统计分析并为每个文件生成图表。

我的文件以 .txt 结尾 这是我的脚本，但需要改进如下。

files <- list.files(path="/Users/MD/Desktop/Files/", pattern="*.txt")
for (i in files){
  i_automean <- mean(i[1:5651,10])
  i_xmean <- mean(i[5652:7977,10])
  barplot(i[,10]/i_automean, ylim = c(0,2.5))
  abline(h=1, col ="red", lwd =3)
  abline(h=22.06126/i_automean, col ="red", lwd =3)
  abline(v=2300*1.2, col="blue", lwd =3)
  abline(v=5651*1.2, col="blue", lwd=3)
  pdf(i.pdf, "/Users/MD/Desktop/Files/");
  dev.off()
} 

Error in i[1:5651, 10] : incorrect number of dimensions

当我为每个文件一个一个地运行每个命令时，我没有收到错误：

我一一运行的命令是：

File1 <- read.table(file = "/Users/MD/Desktop/Files/File1.txt")
File1
File1_automean <- mean(File1[1:5651,10])
File1_automean
File1_Xmean <- mean(File1[5652:7977,10])
File1_Xmean
barplot(File1[,10]/18.02876, ylim = c(0,2.5))
abline(h=1, col ="red", lwd =3)
abline(h=22.06126/18.02876, col ="red", lwd =3)
abline(v=2300*1.2, col="blue", lwd =3)
abline(v=5651*1.2, col="blue", lwd=3)

注意：18.02876 是 File1_automean 结果。

Answer 1

你觉得这个解决方案怎么样？

files <- list.files(path="/Users/MD/Desktop/Files/", pattern="*.txt")
setwd("/Users/MD/Desktop/Files")
for (j in 1:length(files)){
  i <- read.table(files[j])
  i_automean <- mean(i[1:5651,10])
  i_xmean <- mean(i[5652:7977,10])
  barplot(i[,10]/i_automean, ylim = c(0,2.5))
  abline(h=1, col ="red", lwd =3)
  abline(h=22.06126/i_automean, col ="red", lwd =3)
  abline(v=2300*1.2, col="blue", lwd =3)
  abline(v=5651*1.2, col="blue", lwd=3)
  pdf(paste0(files[j],".pdf"))
  graphics.off()
}