如何通过循环批量用列表数据框列子集化矩阵？

当然可以通过循环实现批量处理，这里给你两种实用的方法，推荐用列表统一管理结果，避免生成大量零散变量：

方法1：for循环（直观易理解） #

通过遍历列表metadata_river的每个元素，提取对应的sequence_id来子集化矩阵，最后把结果存入一个命名列表：

# 创建空列表存储所有子集结果
otu_subsets <- list()

# 按河流名称遍历metadata_river
for (river in names(metadata_river)) {
  # 提取当前河流的sequence_id列
  target_ids <- metadata_river[[river]]$sequence_id
  # 子集化otu_table并写入列表，用河流名作为元素名
  otu_subsets[[river]] <- otu_table[, target_ids]
}

之后你可以通过otu_subsets$Alabama调用对应河流的子集结果，和你之前单独生成的al_otu完全一致。

方法2：lapply函数（更简洁高效） #

R中处理列表类数据，lapply是更高效的选择，一行代码就能完成批量处理，结果同样是带命名的列表：

otu_subsets <- lapply(metadata_river, function(df) {
  otu_table[, df$sequence_id]
})

这个方法会自动沿用metadata_river里的元素名称作为结果列表的名称，效果和for循环完全相同。

额外说明（不推荐） #

如果一定要生成37个单独的变量（比如Alabama_otu这类），可以用list2env把列表元素转成全局变量，但这种方式会让环境里变量过多，不利于管理，谨慎使用：

list2env(otu_subsets, .GlobalEnv)

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