三行搞定GSVA分析

果不其然，昨天的推送阅读完成率不到50%：文献阅读系列（九）、一文了解GSVA原理。

大篇幅的原理介绍确实容易让人脑溢血，那我必然知道大家喜欢的是什么，上一次的三行搞定还是：三行命令搞定热图。三行奉上（两行都花在读数据上了）：

myexpr <- read.csv('myexpr.csv')
my.gene.set <- readRDS('my.gene.set.rds')
myresult <- GSVA::gsva(myexpr,my.gene.set)

得到的结果如下：

如果你也想自己试试的话，测试文件在这里，如果过期了可以添加我的qq领取。

链接：https://pan.baidu.com/s/1BHbrBgElbqmPrRd4xi3vFQ?pwd=a5ec

提取码：a5ec

下一篇推送主要来教大家如何详细的通过代码实现GSVA，以及GSVA得到的结果又可以进行如何的下游分析，0点准时更。

最近发现后台中有一些消息我没能及时看到并答复，微信后台中超过48h后便不允许回复读者消息，这里还是再给大家留一下答疑的扣扣号，方便大家随时交流：1913507043。