2023-01-17

¥69.90

购买后会显示你购买的服务器的账号和密码。并且发送服务器的账号和密码到你的邮箱作为备份。格式为xxxx:zzzz （以“：”分割，前边是账号，后边是密码），例如 xiaoyun:998899 那么你的服务器的账号和密码分别是xiaoyun和998899。

有了账号和密码后，不会登录服务器的参考链接：http://www.biocloudservice.com/wordpress/?p=292

有问题咨询客服，客服微信：18502195490

分类：代码

描述

描述

通路富集分析方法已经是生物学功能分析中非常有效的手段，通过利用富集算法如超几何分布，fisher精确检验等可以获得与表型显著相关的生物学通路。然而几乎所有的富集工具所利用统计学方法都是从通路基因与表型相关基因之间的重合率或关联度来评价通路的重要性的，通常我们通过显著性P值来进行评估。然而在生物学角度讲，尽管一些基因富集到了通路中，但是这些基因本身有着不同程度的差异表达，即有些基因异常变化程度很高，体现出与表型的高度关联性，而有些基因尽管差异表达，但是差异程度不明显，很可能是受到关键基因诱导发生共表达，或噪音干扰造成的。因此在评价一个通路与表型变化是否相关的过程中，除了基因数目的影响，我们也应该考虑基因自身表达情况的差异。本代码充分考虑了基因本身的差异表达程度（pvalue值），利用逆累计分布函数将其转化为Z值，通过算法得到通路在这些差异表达基因的差异表达下的偏差值，以挖掘显著偏差的通路以及对通路的显著偏差做出最大贡献的差异表达基因集合。

只需要基因表达谱、样本分组文件及通路富集结果，软件将自行计算出差异表达基因及相应的p.value值，同时挖掘出实验条件或疾病状态下显著偏差的通路以及对通路的显著偏差做出最大贡献的差异表达基因集合。

使用方法：Rscript pathway_score.r –exp= –group= –lfc= –KEGG_pathway=参数说明：

USAGE:

pathway_score.r –exp=–group=–lfc=–KEGG_pathway=

PARAMETERS:

–exp the matrix of gene expression ,input csv format.

–group the group of sampple ,input txt format,“group“ column must exsist.

–KEGG_pathway the KEGG pathway that DEGs were enriched,input txt format.

–lfc the cutoff of logFC.default:1

操作步骤：

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

2、用户根据帮助文档中的参数说明内容，对参数进行设置。这里，必须输入参数有3个，分别是–exp，表示基因表达矩阵文件，这里必须保存为csv格式；–group表示样本分组文件，输入格式为txt文本文件，并且必须包含group列且顺序一定要和基因表达矩阵列对应，相应的分组这里必须设置为“Tumor”，“Normal”，如需修改，需进入程序源代码进行相应修改。–KEGG_pathway表示通路富集结果文件，这里默认的是用clusterProfiler跑出来的结果，要用到p.adjust列。