尔云间生信代码|基于高通量数据的样本相似性分析

本代码基于高通量测序得到的原始count，提供样本相似性分析，主要基于vegan软件包来分析数据。用户需要输入基因表达矩阵文件和样本分组文件，设置相关参数，软件可自行根据用户选择的计算相似性矩阵的方法，从而对输入数据进行组间差异显著性的判别，以量化的数值来精确评判分组是否具有意义。本软件针对3组类别的样本数据，绘制出三组的降维分析PCA和NMDS图，并在图片上添加了用于评判组间差异大小及显著性的相关指标，从而使用户不仅直观感受样本分组情况，而且能够评估组间实验数据的差异性，以及差异性是否显著。

使用方法：

Rscript anosim_mrpp.r -e=  -g=  -m=  -ms=

参数说明：

USAGE:

anosim_mrpp.r -e=-g=-m=-ms=

PARAMETERS:

-e      exprmatrix input csv format.

-g      sample_group file,input txt format.

-m     the method of calculating the difference and significance

between groups,optional parameters :anosim,mrpp.

-ms  the method of calculating similarity index, optional parameters :manhattan,euclidean, canberra, clark, bray,kulczynski, jaccard, gower, altGower,morisita, horn, mountford, raup,binomial, chao, cao or mahalanobis.

操作步骤：

1、打开命令行界面，输入“Rscript anosim_mrpp.r”调阅帮助文档，确定该程序所需的输入文件。

2、用户根据帮助文档中的参数说明内容，对参数进行设置。这里，必须输入参数有4个，分别是-e，表示基因表达矩阵文件，-g表示样本分组文件，第二列的列名必须为group，-m表示计算组间差异显著性的方法，可选参数为anosim和mrpp，-ms表示计算样本相似性矩阵，或者说样本间距离的方法，可选参数有euclidean, bray等。

3、完成参数提交后，按下回车键，整个程序即正式开始进入执行。每步执行内容都会给出提示。程序执行完毕后，界面会显示”Program execution is completed"结束语。

结果展示：

本软件输出*.pdf 图片

特别说明：本代码经申请软件著作权，仅转让使用权，不转让所有权

如需代码及示例数据等文件，请扫码聊天框回复 “代码”领取！ 

写在文末：

如果您近期想做生信方面的文章而苦于没有思路，或者不知道如何来入手生信分析，或者兑具体的某一个图有作图需求，都可以扫码咨询小云，我们有专业的技术团队，生信热点思路设计、生信分析、热点方向生信挖掘等，如有需要，可扫码下方二维码了解详情 !