q2cli 1.在查看插件的详细信息时清…
最近觉得半瓶醋的16S数据分析可以继续进…
这两天在用mac下载ena的fastq数据,发现aspera connect,windows和mac版本的死活不会用呢,于是就学习了一下,发现一样是可以使用命令行的,这里是mac的,后面再整理下win的。
最近工作用到了一些回归分析方面的知识,整理一下,放在这里备忘。主要是简单的多元线性回归和分位数回归两种,基本上都只需要几个R语言命令就能解决了。
发现自己还是比较菜的,不能在qiime2的结果文件里进行较好的操作来达到获得各种属相对丰度的表,于是决定把数据导出到biom和txt格式来进行后续操作。这样就可以利用qiime2的新算法的优势,和之前已有的宏基因组公众号上介绍的流程来进行后续操作。毕竟,对于16S扩增子分析最重要的就是获得物种分类丰度信息。
qiime2-2019.1已经发布,程序稳定性越来越好,鉴于官方已经停止支持qiime1,有必要把qiime2的所有细节都理清,学好,这样才能对自己的数据进行实战分析,并将结果运用于实验和生产过程中。发现文档更新也相当快,感谢公众号宏基因组翻译的文档,让我在看许多专业术语方面扫清不少障碍,但是你介于翻译过来的命令却已经过时,还是对照着看最新版的,基本上很少改动,当然,专业英语好的除外。
QIIME2 2019.1发布啦,继续紧跟步伐,看看更新了哪些内容,一并备忘。
这个发布生版本主要针对更新依赖环境,升级到了Python 3.6,一个新的 r-vegan Adonis visualizer,修复了一些小bugs,虚拟机镜像和新的预训练分类器立马可得。
昨天,想用R语言画个PCA图,来确认一批样本重复测序的吻合度,虽然已经单个数据审视和检查过了,还是用PCA看一下整体的情况,于是决定学画PCA图。使用搜狗微信搜索,找到一篇百迈客微信公众号的视频讲解,讲的不错,跟着以1.5倍速听完就可以画出自己的图了,挺开心的。
这两天的工作都是从早上打开热心肠日报开始的,不得不说基本上每天干货满满呀,每天一篇给力的综述的节奏。ps,发现ubiome的科学家和志愿者也在维护着一个国外版的热心肠日报,地址放在这:https://microbiomedigest.com/ 值得点赞,英语好的同鞋可以参考参考,不知道国内的灵感是不是来自这个,但是平心而论,国内的做得已经水平不次于国外的了。
很感谢这些人,为大众做着科普的事情,有许多人把今天的生命科学比做几十年前的计算机产业,当年,计算机还是庞然大特,像供神一样地供在一间特殊的房间,多么像今天的测序仪呀。而牛津纳米孔,直接把测序仪缩小到了手持设备的级别,简直从当年的大型机飞跃到了今天的智能手机呀!相信未来,测序仪会像手机一样人手一部来检测自己的各项基因、蛋白和代谢指标,这才是精准医学嘛。
祖源是消费级基因检测比较感兴趣的问题,今天我也来学习一下那如何做一个简陋的祖源分析。使用的是SPA(Spatial Ancestry analysis)这个软件,在google中搜索基因型填充无意中搜到了这个软件,由几个作者写成,第一作者还是华人,有点骄傲,可惜没有针对中国人或者东亚人的模型,只有世界和欧洲人的模型,瞬间又觉得有点失落,好在许多中国公司在积累了几万以上的数据后也开始进行中国人自己的模型构建,而且做得还不错。