单细胞-GSEA结果解读

最新推荐文章于 2024-03-07 00:04:39 发布

生信学徒

最新推荐文章于 2024-03-07 00:04:39 发布

阅读量6.5k

点赞数 3

文章标签：生物信息学

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/qq_45478665/article/details/119642140

版权

做单细胞富集分析的，大家都应该知道GSEA这个软件吧，但是GSEA分析的结果是一个很复杂的图，大家很难看懂，我在这里给大家说下这个结果图怎么看。

这个图分为三个部分

第一部分：为基因 Enrichment Score 的折线图，横轴为该基因下的每个基因，纵轴为对应的 Running ES，在折线图中有个峰值，该峰值就是这个基因集的 Enrichemnt score，峰值之前的基因就是该基因集下的核心基因，核心基因就是在单细胞很重要很关键的基因，如果基因在顶部富集，我们可以说，从总体上看，该基因集是上调趋势，反之，如果在底部富集，则是下调趋势；

第二部分：也就是hit值，也就是图中的黑线出现的位置，这个黑线可以看出第一部分每次折现点的基因的位置

第三部分：这是rank值图，采用了 Signal2Noise 算法，我对他的理解是，纵轴它反映了基因集中的每个基因的相对表达量，横轴可以看出这是这个基因在基因数据集中出现的位置。

我所做的所有分析与教程的代码都会在我的个人公众号中，请打开微信搜索“生信学徒”进行关注，欢迎生信的研究人员和同学前来讨论分析。

ps：公众号刚刚建立比较简陋，但是该有的内容都不会少。

————————————————

关注

3
点赞
踩
16

收藏

觉得还不错? 一键收藏
0
评论
单细胞-GSEA结果解读

做单细胞富集分析的，大家都应该知道GSEA这个软件吧，但是GSEA分析的结果是一个很复杂的图，大家很难看懂，我在这里给大家说下这个结果图怎么看。这个图分为三个部分第一部分：为基因 Enrichment Score 的折线图，横轴为该基因下的每个基因，纵轴为对应的 Running ES，在折线图中有个峰值，该峰值就是这个基因集的 Enrichemnt score，峰值之前的基因就是该基因集下的核心基因，核心基因就是在单细胞很重要很关键的基因，如果基因在顶部富集，我们可以说，从总体上看，该基
复制链接

扫一扫

生信学徒 CSDN认证博客专家 CSDN认证企业博客

码龄5年

26: 原创

13万+: 周排名

47万+: 总排名

22万+: 访问

: 等级

812: 积分

1069: 粉丝

179: 获赞

71: 评论

1162: 收藏

私信

关注

热门文章

最新评论

单细胞测序流程（十）单细胞的kegg富集分析和圈图（终章）
2401_84573921: use strict; use warnings; my %hash = (); open(RF, "id.txt") or die $!; while(my $line = <RF>) { chomp($line); my @arr = split(/\t/, $line); $hash{$arr[2]} = $arr[0]; } close(RF); open(KEGG, "keggId.txt") or die $!; open(WF, ">kegg.txt") or die $!; while(my $line = <KEGG>) { if ($. == 1) { print WF $line; next; } chomp($line); my @arr = split(/\t/, $line); my @idArr = split(/\//, $arr[$#arr - 1]); my @symbols = (); foreach my $id (@idArr) { if (exists $hash{$id}) { push(@symbols, $hash{$id}); } } $arr[$#arr - 1] = join("/", @symbols); print WF join("\t", @arr) . "\n"; } close(WF); close(KEGG);
单细胞测序流程（十）单细胞的kegg富集分析和圈图（终章）
2401_84573921: 可以解决吗？
单细胞测序流程（五）t-sne聚类分析和寻找marker基因
m0_51028803: 你试试把绘图框拉放大，我上次就是这样然后出来了
单细胞测序流程（七）单细胞的细胞类型轨迹分析
m0_63761949: 画图celltype2 和cluster一样的话，celltype赋值那儿改一下，直接view（pdata（cds））看一下整体结构，可以看出拿他的行名再在前面粘一个cluster就和clusterann的键都对应上了，不要用pdata（cds）＄cluster去对应，和clusterann的键不一样。
单细胞测序流程（七）单细胞的细胞类型轨迹分析
沙茶还是土笋冻: 引用「cds <- newCellDataSet(data, phenoData = pd, featur」 Error in warn_version(graph) : This graph was created by a now unsupported old igraph version. Call upgrade_version() before using igraph functions on that object. 请问这个怎么办

最新文章

目录

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

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