java8阻止gsea运行_过了个年,你连GSEA都不认得了!

基因集富集分析(Gene Set Enrichment Analysis,GSEA)是麻省理工学院和哈佛大学的broad institute研究团队开发的一个针对全基因组表达谱芯片数据进行分析的工具。基本思想是使用预定义的基因集,将基因按照在两类样本中的差异表达程度排序,然后检验预先设定的基因集合是否在这个排序表的顶端或者底端富集。基因集合富集分析检测基因集合而不是单个基因的表达变化,因此可以包含这些细微的表达变化,预期得到更为理想的结果。

GSEA富集过程包括三步骤:

1.计算富集分数(Enrichment Score);

2.估计富集分数的显著性水平;

3.矫正多重假设检验。

131494504_1_20180428110301222

那今天,我们就来主要讲讲GSEA富集分析的操作,和结果解释吧~

首先注册GSEA

(http://software.broadinstitute.org/gsea/register.jsp)

131494504_2_20180428110301285

直接通过邮箱进行注册,十分方便快捷。

然后,再下载一个GSEA软件

(http://software.broadinstitute.org/gsea/login.jsp;jsessionid=F71E86DCBFB4F9CBED84C0495E5B02A7)

131494504_3_20180428110301394

用刚才注册的箱邮登录,我们就可以直接选择符合自已电脑的版本进行下载

131494504_4_20180428110301488

好了,有了GSEA,还得保证你有个够格java环境(http://www.oracle.com/technetwork/java/javase/downloads/jdk8-downloads-2133151.html)

131494504_5_20180428110301660

下载java,选择适合自己电脑的版本131494504_6_20180428110301785

做好一切准备功夫,我们就开始使用GSEA!

制作输入表格

一个gct文件,后缀后是gct,可先在excel做好,然后复制到记事本中(推荐sublime或notepat ),另存为gct后缀的文件

131494504_7_20180428110301910

一个分组文件cls

131494504_8_20180428110301988

“9 3 1”代表着9个样品,分为3组,1是默认值,以空格或tab分开

第二行是分组信息

第三行与gct文件中的样品对应,分属哪个组

做好这两个文件后,开始使用GESA了

131494504_9_2018042811030266

选择你自已的GSEA版本,如果java环境装好的话,就可以运行了

131494504_10_20180428110302128

加载之前做好的两个文件

131494504_11_20180428110302347

出现no errors 就行了

131494504_12_20180428110302488

1.选择表达矩阵

2.富集背景,这里只选了KEGG和GO

131494504_13_20180428110302753

3. Number of permutations这里最多可以选择1000次(次数越多结果相对越可靠,但是占用CPU越多,我这里选的是 1000,如果你电脑支持不了就调小一些,别勉强喔~)

4.选择比较组

5.每组样品都有多于7个样品,则Permutation type选择phenotype,结果理论上更好;否则选择gene_set(这里选的是gene_set.)

6. chip平台,这里用的Ensemb mmu的数据

131494504_14_20180428110302894

7.选择输出结果的路径

填好所有信息之后,结果好了在右边会显示success,点击Success可以直接显示结果

131494504_15_20180428110302988

这样富集分析就成了~

131494504_16_2018042811030350

对于富集分析的结果是如何分析的,大家可以关注明天的内容喔~

编辑:咸鱼君

无损读文献

无伤看干货

这里有逗比的科研老司机

车票不要998,也不要888,只要666

谈笑间,带你走遍科研的雄山峻岭

赛哲生物︱成就创业与创新梦想

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值