如何在R中以pdf格式批量输出图

应用场景:比如说我们有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位置比较美观。

你可能感兴趣的:(如何在R中以pdf格式批量输出图)