# 扩增子统计绘图5火山图：差异OTU数量及变化规律

《扩增子统计绘图》系列文章介绍

《扩增子统计绘图》是之前发布的《扩增子图表解读》和《扩增子分析解读》的进阶篇，是在大家可以看懂文献图表，并能开展标准扩增子分析的基础上，进行结果的统计与可视化。其章节设计与《扩增子图表解读》对应，为八节课八种常用图形(箱线图、散点图、热图、曼哈顿图、火山图、维恩图、三元图和网络图)，基本满足文章常用的图片种类需求。

5火山图：差异OTU数量及变化规律# 绘制火山图if (max(x$logFC)>4){x[x$logFC>4,]$logFC = 4} # norm x axisif (min(x$logFC)< -4){x[x$logFC< -4,]$logFC = -4} # norm x axisx$level = as.factor(ifelse(x$sig==1, "enriched",ifelse(x\$sig==-1, "depleted","nosig")))# Volcanol plot of fold change vs abundance plotp = ggplot(x, aes(x=logFC, y=logCPM, color=level)) + geom_point()  +   scale_colour_manual(values=c("green","red","grey"))+ xlim(-4, 4)+  labs(x="log2(fold change)",y="log2(count per million)", title=paste("OE vs WT", sep=" "))pggsave(file=paste("vol_otu.pdf", sep=""), p, width = 8, height = 5)ggsave(file=paste("vol_otu.png", sep=""), p, width = 8, height = 5)

# 绘制有物种信息和相对丰度的火山图p = ggplot(x, aes(x=logFC, y=logCPM, color=level, size=logCPM, shape=tax)) + geom_point()  +   scale_colour_manual(values=c("red","green","grey"))+ xlim(-4, 4)+  labs(x="log2(fold change)",y="log2(count per million)", title=paste("OE vs WT", sep=" "))p

https://m.sciencenet.cn/blog-3334560-1072641.html

## 全部精选博文导读

GMT+8, 2023-1-27 06:55