如何在R中批量处理同一文件夹内的79个转移矩阵文件并生成平均流图？

问题背景 #

已完成单个转移矩阵文件的R分析流程，需批量处理同一文件夹内78个同格式文件，最终生成所有个体的平均流图，需用循环实现批量操作。

原单文件分析代码 #

library(RcmdrMisc)
# 读取单个文件
BHFF_14h31 <- read.table("BHFF_14h31.txt", header=TRUE, row.names=1, sep="\t", na.string="", dec=",")
BHFF_14h31_m <- as.matrix(BHFF_14h31)

# 计算列百分比
BHFF_14h31 <- colPercents(BHFF_14h31_m, 2)
# 计算行列和与香农多样性
transc <- colSums(BHFF_14h31)
transr <- rowSums(BHFF_14h31)
diversity(transc, index = "shannon")

# 绘制流图
plotmat(BHFF_14h31)

批量处理实现步骤 #

1. 环境准备 #

设置工作目录为包含所有转移矩阵文件的文件夹，并加载依赖包：

library(RcmdrMisc)

# 替换为你的实际文件夹路径
setwd("your/file/directory/path")

2. 获取目标文件列表 #

提取文件夹内所有.txt格式的转移矩阵文件：

# 获取所有txt文件的完整路径
file_list <- list.files(pattern = "\\.txt$", full.names = TRUE)

3. 初始化存储对象 #

创建空矩阵用于累积所有文件的处理结果，同时创建向量存储每个文件的香农多样性值：

# 读取第一个文件获取矩阵维度，初始化累积矩阵
first_data <- read.table(file_list[1], header=TRUE, row.names=1, sep="\t", na.string="", dec=",")
accumulated_matrix <- matrix(0, nrow = nrow(first_data), ncol = ncol(first_data), 
                            dimnames = dimnames(first_data))

# 初始化存储香农多样性的向量
shannon_values <- numeric(length(file_list))
names(shannon_values) <- basename(file_list)

4. 循环处理所有文件 #

遍历每个文件，执行与单文件相同的分析流程，并累积结果：

for (i in seq_along(file_list)) {
  # 读取当前文件
  current_data <- read.table(file_list[i], header=TRUE, row.names=1, sep="\t", na.string="", dec=",")
  current_matrix <- as.matrix(current_data)
  
  # 计算列百分比
  current_percent <- colPercents(current_matrix, 2)
  
  # 计算香农多样性并存储
  transc <- colSums(current_percent)
  shannon_values[i] <- diversity(transc, index = "shannon")
  
  # 累积当前百分比矩阵
  accumulated_matrix <- accumulated_matrix + current_percent
  
  # 可选：绘制单个文件的流图（按需开启）
  # plotmat(current_percent, main = basename(file_list[i]))
}

5. 计算平均矩阵并绘制平均流图 #

# 计算所有文件的平均百分比矩阵
average_matrix <- accumulated_matrix / length(file_list)

# 绘制平均流图
plotmat(average_matrix, main = "所有个体的平均转移流图")

6. 保存结果（可选） #

# 保存平均矩阵到文件
write.table(average_matrix, "average_transition_matrix.txt", sep="\t", quote=FALSE)

# 保存香农多样性结果
write.table(shannon_values, "shannon_diversity_results.txt", sep="\t", quote=FALSE)

内容的提问来源于stack exchange，提问作者Shyam Sharma