扩增子统计绘图1箱线图:Alpha多样性
来源:互联网 发布:农村网络服务站 编辑:程序博客网 时间:2024/06/03 11:16
本网对Markdown排版支持较差,对格式不满意的用户请跳转至 或“宏基因组”公众号阅读;
写在前面
优秀的作品都有三部分曲,如骇客帝国、教父、指环王等。
扩增子系列课程也分为三部曲:
第一部《扩增子图表解读》:加速大家对同行文章的解读能力。
第二部《扩增子分析解读》:学习数据分析的基本思路和流程。
第三部《扩增子统计绘图》:即是对结果进行可视和统计检验,达到出版级的图表结果。
《扩增子统计绘图》系列文章介绍
《扩增子统计绘图》是之前发布的《扩增子图表解读》和《扩增子分析解读》的进阶篇,是在大家可以看懂文献图表,并能开展标准扩增子分析的基础上,进行结果的统计与可视化。其章节设计与《扩增子图表解读》对应,为八节课八种常用图形(箱线图、散点图、热图、曼哈顿图、火山图、维恩图、三元图和网络图),基本满足文章常用的图片种类需求。
也适合对公司标准化分析返回结果的进一步统计、可视化及美化,达到出版级别,冲击高分文章。
本部分练习所需文件位于百度网盘,链接:http://pan.baidu.com/s/1hs1PXcw 密码:y33d。
绘制Alpha多样性线箱图
绘图和统计全部为R语言,建议复制代码,在Rstuido中运行,并设置工作目录为存储之前分析结果文件的result目录。
# 运行前,请在Rstudio中菜单栏选择“Session - Set work directory -- Choose directory”,弹窗选择之前分析目录中的result文件夹# 安装相关软件包,如果末安装改为TRUE运行即可安装if (FALSE){ source("https://bioconductor.org/biocLite.R") biocLite(c("ggplot2"))}# 加载相关软件包library("ggplot2") # load related packages# 读入实验设计和Alpha多样性值design = read.table("design.txt", header=T, row.names= 1, sep="\t") alpha = read.table("alpha.txt", header=T, row.names= 1, sep="\t")# 以Observed OTU为例进行可视化和统计分析,其它指数将observed_otus替换为shannon, chao1, PD_whole_tree即可计算# 合并Alpha指数与实验设计index = cbind(alpha, design[match(rownames(alpha), rownames(design)), ]) # 绘图代码、预览、保存PDFp = ggplot(index, aes(x=genotype, y=observed_otus, color=genotype))+ geom_boxplot(alpha=1, outlier.size=0, size=0.7, width=0.5, fill="transparent") + geom_jitter( position=position_jitter(0.17), size=1, alpha=0.7)+ labs(x="Groups", y="observed_otus index")pggsave(paste("alpha_observed_otus.pdf", sep=""), p, width = 5, height = 3)# 统计组间是否显著差异# anova对指数与分组统计observed_otus_stats <- aov(observed_otus ~ genotype, data = index)# 使用TukeyHSD对组间进行检验,效正pvalueTukey_HSD_observed_otus <- TukeyHSD(observed_otus_stats, ordered = FALSE, conf.level = 0.95)# 结果中提取需要的结果Tukey_HSD_observed_otus_table <- as.data.frame(Tukey_HSD_observed_otus$genotype)# 预览结果Tukey_HSD_observed_otus_table# 保存结果到文件,按Pvaule值由小到大排序write.table(Tukey_HSD_observed_otus_table[order(Tukey_HSD_observed_otus_table$p, decreasing=FALSE), ], file="alpha_observed_otus_stats.txt",append = FALSE, quote = FALSE, sep="\t",eol = "\n", na = "NA", dec = ".", row.names = TRUE,col.names = TRUE)
Observed OTU多样性箱线图
详细的图片讲解,可参考1箱线图:Alpha多样性,老板再也不操心我的文献阅读
各组间的统计结果如下:主要看最后一列p adj(Adjust P-value)是否显著,本文数据不显著
diff lwr upr p adjOE-KO -7.52380952380952 -24.480725165752 9.43310611813294 0.515429907536906WT-KO -6.11111111111111 -21.9728532782553 9.75063105603303 0.604309699204896WT-OE 1.4126984126984 -15.5442172292441 18.3696140546409 0.976169656924344
其它箱线图绘制教程
推荐阅读并实际操作陈同博士写的文章:
R语言学习 - 箱线图(小提琴图、抖动图、区域散点图)
R语言学习 - 箱线图一步法
Reference
- http://mp.weixin.qq.com/s/5PkLAQ0yOxyGfl9MAH7veQ
- http://mp.weixin.qq.com/s/Zvmht0kOyOf02P8jQNjaOw
- http://mp.weixin.qq.com/s/8w6jV9MtJZ4h3ATaPP_Rsw
- http://mp.weixin.qq.com/s/MJWyUzeBZ7WJQhpAGX6Fgg
想了解更多宏基因组、16S文献阅读和分析相关文章,快关注“宏基因组”公众号,干货第一时间推送。
系统学习生物信息,快关注“生信宝典”,那里有几千志同道合的小伙伴一起学习。
- 扩增子统计绘图1箱线图:Alpha多样性
- 扩增子图表解读1箱线图:Alpha多样性,老板再也不操心的我文献阅读了
- 扩增子统计绘图2散点图:Beta多样性
- 扩增子分析解读6进化树,Alpha,Beta多样性
- 宏基因组扩增子统计绘图大全:中文首发,最详系,零基础(箱线图、散点图、热图、曼哈顿图、火山图、韦恩图、三元图、网络图)
- 扩增子统计绘图7三元图
- 扩增子统计绘图8网络图-MENA
- 扩增子统计绘图3热图:样品相关分析,差异OTU
- 扩增子统计绘图4曼哈顿图:差异OTU和Taxonomy
- 扩增子统计绘图5火山图:差异OTU数量及变化规律
- 扩增子统计绘图6韦恩图:比较组间共有和特有OTU或分类单元
- 扩增子图表解读2散点图:组间整体差异分析(Beta多样性)
- 零基础测序分析图表解读大全:老板再也不愁我的文献阅读了! 宏基因组扩增子(箱线图、散点图、热图、曼哈顿图、火山图、韦恩图、三元图、网络图)
- 扩增子分析QIIME2. 7 实验设计和统计结果元数据Metadata
- 扩增子分析解读7物种分类统计,筛选进化树和其它
- python第三方库-pygal绘图功能介绍1-线图
- 扩增子分析流程QIIME. 1 使用Docker配置QIIME
- 扩增子分析QIIME2. 1简介和安装
- 【POJ2449】Remmarguts' Date (A* k短路)
- "could not create the java virtual machine"启动eclipse报错
- 原来这就是 ’面向对象‘, ‘编程’。。
- [转]hadoop,spark,storm,pig,hive,mahout等到底有什么区别和联系?
- 在windows上传一个新的项目到GitHub上
- 扩增子统计绘图1箱线图:Alpha多样性
- Building Serverless Architectures pdf
- ElasticSearch-深入理解系列8-高级查询
- 树莓派之SPI编程
- I2S波形记录
- socket
- HDU 4746-Mophues (莫比乌斯反演)
- Python模块包中__init__.py文件的作用
- 国产 Android 权限申请最佳适配方案 —— permissions4m