代码来源:
(1条消息) R语言完美重现STAMP结果图_刘永鑫Adam的博客-CSDN博客
先看一下STAMP软件中的丰度差异比图的参数有那些
看到上面的P值过滤条件是选择0.05,最重要的是上面有一个多次矫正方法选择,软件里面是没有选择的,代码是有选择的
唯一要做调整的就是平均丰度低于1%,STAMP软件没有做这样的调整。
后面根据代码得到的结果图,不难看出结果存在一定的偏差
其实我最开始用STAMP软件,没有选择显著性检验方式,结果显示的是九个,选择后就是七个,但是后面重新导入文件的时候,结果就一直是九个,后面发现是没有进行校正,选择一样的P值校正方式,结果是一样的,不过软件的精确度保存不同。
仔细看这里要选择跟代码一样的校正方式
代码编写的程序要方便很多,尤其是在做多组分析的时候,不用像软件要一个一个保存。