如何在R中以pdf格式批量输出图
2020-03-17 本文已影响0人
台式计算机
应用场景:比如说我们有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位置比较美观。