0

私は Seurat で単一セル rna-seq に取り組んでおり、Seurat オブジェクトに対して for() ループを作成して、平均遺伝子発現のヒートマップをいくつか描画しようとしています。

for(i in c(seuratobject1, seuratobject2, seuratobject3)){
  cluster.averages <- data.frame(AverageExpression(i, features = genelist))
  cluster.averages$rowmeans <- rowMeans(cluster.averages)
  genelist.new <- as.list(rownames(cluster.averages))
  cluster.averages <- cluster.averages[order(cluster.averages$rowmeans),]
  HMP.ordered <- DoHeatmap(i, features = genelist.new, size = 3, draw.lines = T)
  ggsave(HMP.ordered, file=paste0(i, ".HMP.ordered.png"), width=7, height=30)

i を seurat オブジェクトとして使用するため、ggsave 行は機能しません。したがって、私の質問: "i" に格納されている私の seurat オブジェクトの名前を使用するように ggsave() を取得するにはどうすればよいですか?

私はsubstitute(i)とdeparse(substitute(i))を試しましたが成功しませんでした。

4

1 に答える 1