qiime2类型文件:读入的feature不能太多,亲测47000多个feature消耗内存110G,导入好几个小时依然失败。而一万多的feature就大概一分钟就能正常读入。
且看文件大小看不出来。
reads count数比较少的,例如只有
几十条
的这种样本需要去除,不然出来的alpha稀释曲线会以这个最少的为底限,最终出图x轴readsNums就只到几十。画图支持多分组:
.group = c("class1", "class2")。
alpha多样性的图可以自己设置比较组:
comparison = list(c("g1", "g2"), c("g1", "g3"))
。可以用step_increase = 0.1
调整比较组的连线距离,避免重叠。可以搭配使用ggsci包进行颜色修改。
MicrobiotaProcess注意事项
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 原文:https://www.sohu.com/a/270075112_99971433[https://www....
- 给粉丝朋友们带来了很多理解上的挑战,所以我们开辟专栏慢慢介绍其中的一些概念性的问题,上一期: 箱线图的生物学含义 ...
- 背景:qiime artifact 是用于存储qiime2的输入输出以及相关的元数据,并提供结果是如何产生的信息,...