R语言嵌套循环遍历嵌套列表调用Seurat FeaturePlot问题

原代码存在几处语法和逻辑错误，逐一修正后即可实现需求：

• genotypes列表定义时第三个元素缺失赋值符号=，会直接触发语法报错
• 外层循环写法for (i in length(genotypes))逻辑错误：length(genotypes)仅返回列表长度的单个数值，循环只会执行最后一个索引的单次遍历，无法覆盖所有Seurat对象
• 内层循环未匹配与外层索引一一对应的基因子集：定义的top_genes列表顺序和genotypes顺序完全对应，第i个基因型对应的目标基因集就是top_genes[[i]]
• 额外注意：Seurat的FeaturePlot函数的features参数接收字符型向量，原代码存储基因名用的是嵌套list结构，传入时需要展平为字符向量避免参数类型报错。

修正后的可运行代码如下：

condition_1_top_genes = list("abc", "def", "efg")
condition_2_top_genes = list("a35", "2353", "rea3")
condition_3_top_genes = list("fae", "wai", "wtes")

top_genes = list(condition_1_top_genes, 
                 condition_2_top_genes,
                 condition_3_top_genes)

# 补全缺失的赋值符号
genotypes <- list("genotype1" = genotype1, 
                  "genotype2" = genotype2,
                  "genotype3" = genotype3)

# 用seq_along遍历所有列表索引，避免length写法的逻辑错误
for (i in seq_along(genotypes)){
  # 取当前索引i对应的基因子集，unlist将list结构转为字符向量适配FeaturePlot参数要求
  current_genes <- unlist(top_genes[[i]])
  for (j in current_genes){
    FeaturePlot(genotypes[[i]], features = j)
    # 若需要自动保存出图，可在此处添加ggsave等保存逻辑
  }
}

如果后续需要调整基因集和基因型的对应关系，只需要保证top_genes中子列表的顺序和genotypes中Seurat对象的顺序一致即可，不需要修改循环主体代码。

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