如何在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位置比较美观。

©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容