應用場景:比如說我們有50個marker,需要輸出50張UMAP圖,每張圖上顯示一個marker的表達情況。
pdf("./name.pdf")
for( i in 1:length(myMarker)){
p1<-FeaturePlot(object = pbmc, features =myMarker[i],pt.size = 1.5)+ggtitle(myMarker_geneName[i])
print(p1)}
dev.off()
其中,myMarker是一個向量,為很多基因的Ensembl向量,而myMarker_geneName是相對應的Gene Symbol,放在title位置比較美觀。