将qiime2学习进行到底2

发现自己还是比较菜的,不能在qiime2的结果文件里进行较好的操作来达到获得各种属相对丰度的表,于是决定把数据导出到biom和txt格式来进行后续操作。这样就可以利用qiime2的新算法的优势,和之前已有的宏基因组公众号上介绍的流程来进行后续操作。毕竟,对于16S扩增子分析最重要的就是获得物种分类丰度信息。

继续阅读 →

将QIIME2学习进行到底

qiime2-2019.1已经发布,程序稳定性越来越好,鉴于官方已经停止支持qiime1,有必要把qiime2的所有细节都理清,学好,这样才能对自己的数据进行实战分析,并将结果运用于实验和生产过程中。发现文档更新也相当快,感谢公众号宏基因组翻译的文档,让我在看许多专业术语方面扫清不少障碍,但是你介于翻译过来的命令却已经过时,还是对照着看最新版的,基本上很少改动,当然,专业英语好的除外。

继续阅读 →

qiime2-2019.1更新学习笔记

QIIME2 2019.1发布啦,继续紧跟步伐,看看更新了哪些内容,一并备忘。 这个发布生版本主要针对更新依赖环境,升级到了Python 3.6,一个新的 r-vegan Adonis visualizer,修复了一些小bugs,虚拟机镜像和新的预训练分类器立马可得。

继续阅读 →

PCA画图学习

昨天,想用R语言画个PCA图,来确认一批样本重复测序的吻合度,虽然已经单个数据审视和检查过了,还是用PCA看一下整体的情况,于是决定学画PCA图。使用搜狗微信搜索,找到一篇百迈客微信公众号的视频讲解,讲的不错,跟着以1.5倍速听完就可以画出自己的图了,挺开心的。

继续阅读 →

R包decontam学习笔记

这两天的工作都是从早上打开热心肠日报开始的,不得不说基本上每天干货满满呀,每天一篇给力的综述的节奏。ps,发现ubiome的科学家和志愿者也在维护着一个国外版的热心肠日报,地址放在这:https://microbiomedigest.com/ 值得点赞,英语好的同鞋可以参考参考,不知道国内的灵感是不是来自这个,但是平心而论,国内的做得已经水平不次于国外的了。 很感谢这些人,为大众做着科普的事情,有许多人把今天的生命科学比做几十年前的计算机产业,当年,计算机还是庞然大特,像供神一样地供在一间特殊的房间,多么像今天的测序仪呀。而牛津纳米孔,直接把测序仪缩小到了手持设备的级别,简直从当年的大型机飞跃到了今天的智能手机呀!相信未来,测序仪会像手机一样人手一部来检测自己的各项基因、蛋白和代谢指标,这才是精准医学嘛。

继续阅读 →

探索消费级基因检测–祖源

祖源是消费级基因检测比较感兴趣的问题,今天我也来学习一下那如何做一个简陋的祖源分析。使用的是SPA(Spatial Ancestry analysis)这个软件,在google中搜索基因型填充无意中搜到了这个软件,由几个作者写成,第一作者还是华人,有点骄傲,可惜没有针对中国人或者东亚人的模型,只有世界和欧洲人的模型,瞬间又觉得有点失落,好在许多中国公司在积累了几万以上的数据后也开始进行中国人自己的模型构建,而且做得还不错。

继续阅读 →

qiime2和美国肠道计划学习笔记

qiime2正在从一个单纯的marker gene分析软件向代谢组学和基因组学发展,并且,正逐步形成一个生态系统,方便任何开发者为之开发插件。.qza和.qzv其实是标准的zip格式,可以方便解压得到数据,和之前一直以为的qiime2是封闭的,不让人把数据解开的印象完全不一样嘛。

继续阅读 →