基因集变异分析(GSVA)是一种非参数、无监督的方法,用于通过表达数据集的样本估计基因集富集的变异。GSVA在坐标系中执行更改,将数据从逐个样本矩阵的基因转换为逐个样本矩阵的基因集,从而允许对每个样本的路径富集进行评估。这种新的GSVA富集分数矩阵有助于以途径为中心应用标准分析方法,如功能富集、生存分析、聚类、CNV途径分析或跨组织途径分析。

# if (!requireNamespace("BiocManager", quietly = TRUE))
#   install.packages("BiocManager")
#
# BiocManager::install("GSVA")library(GSEABase)
library(GSVA) p <- 10000 ## number of genes
n <- 30    ## number of samples
## simulate expression values from a standard Gaussian distribution
X <- matrix(rnorm(p*n), nrow=p,dimnames=list(paste0("g", 1:p), paste0("s", 1:n)))
X[1:5, 1:5]## sample gene set sizes
gs <- as.list(sample(10:100, size=100, replace=TRUE))
## sample gene setsgs <- lapply(gs, function(n, p)paste0("g", sample(1:p, size=n, replace=FALSE)), p)names(gs) <- paste0("gs", 1:length(gs))# Estimates GSVA enrichment scores.
gsva.es <- gsva(X, gs, verbose=FALSE)
class(gsva.es) # "matrix" "array"
dim(gsva.es)
gsva.es[1:5, 1:5]## gsva.es 列为样本品,行为不同基因集的富集分数,
## gsva.es矩阵可以用limma,DESeq2等进行差异化分析

GSVA包进行基因集变异分析相关推荐

  1. 一个R包完成单细胞基因集富集分析 (全代码)

    singleseqgset | 单细胞RNA-Seq基因集富集分析 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (Ch ...

  2. ClusterProfiler在线基因集富集分析,支持自定义基因集、任意物种

    为什么pathway富集分析结果没有我感兴趣的通路? GO和KEGG富集分析使用差异基因(上调基因,下调基因,或者上下调合起来的基因)作为输入,使用超几何分布等算法计算显著富集的GO term或者通路 ...

  3. 单个基因集富集分析泡泡图绘制

    富集分析是生物信息分析中快速了解目标基因或目标区域功能倾向性的最重要方法之一.其中代表性的计算方式有两种: 一是基于筛选的差异基因,采用超几何检验判断上调或下调基因在哪些GO或KEGG或其它定义的通路 ...

  4. 基因集富集分析(GSEA)

    GSEA简介 首先简单介绍一下GSEA,它是2005年在PNAS上发扬光大的方法,沿用至今,目的是看差异表达的基因在哪些基因集中富集.相比于Over-representation只关注显著差异表达的基 ...

  5. 单基因gsea_GSEA:基因集富集分析和ssGSEA:单样本基因集富集分析

    传统富集分析(基于超几何分布或者Fisher精确检验):关注一列差异基因是否是随机分布在某一感兴趣的基因集中(某通路的基因) 得到通路富集的结果时: (1).一条通路中既有上调基因又有下调基因,无法确 ...

  6. 评分addmodule 绝对值评分 8种方法可视化你的单细胞基因集打分gsea 缺氧评分

    基因集合打分不考虑基因高低表达 基因集合评分不考虑高低表达 绝对值 基因集合里面的正负表达影响addmodulescore的结果 构建的基因集进行评分时 基因的高低正负表达基因集合评分不考虑高低表达_ ...

  7. 基因集泛癌分析,强势来袭!!高分必选!!

    本团队最近筛选到一些比较好的基因,需要对单基因泛癌做深入分析的请与本团队联系,先到先得. 1.快来看看7+的单基因-泛癌-纯生信怎么做,一看就会! 2.单基因泛癌发7+!可复现性极强!你值得拥有! 3 ...

  8. gsva gsea ssgsea gaochao 使用GSVA方法计算某基因集在各个样本的表现

    傻傻分不清!GSEA & GSVA有啥差别?史上最全教程来了! - 知乎 (zhihu.com) 文章发表于2013年,GSVA: gene set variation analysis fo ...

  9. 如何利用clusterProfiler进行基因集的KEGG富集分析?

    NGS 测序项目,不管是基因组测序,还是转录组测序,通常会得到一个基因列表,记录了基因突变,或者高/低表达量. 对成百上千甚至上万个基因进行解读,往往是困难的,对基因进行分组以帮助对数据的理解就非常有 ...

最新文章

  1. 第十五章 异常检测-机器学习老师板书-斯坦福吴恩达教授
  2. 和cnn结合_写给小白的R-CNN介绍
  3. Java动态解析域名
  4. Django Cookie于Session
  5. 【Java】java中this$0 this$1 this$2
  6. 通讯(transport)
  7. unity用代码生成的物体或line renderer让其在确定的某一个层里面的方法
  8. linux mysql配置_Linux下MySQL安装配置 MySQL配置参数详解
  9. 统计学习导论(ISLR) 第四章分类算法课后习题
  10. 显示器尺寸对照表_电脑显示器尺寸对照表(电脑显示器尺寸一览表)
  11. 计算机十二种常用密码的破解方法
  12. 什么软件可以更换手机ip地址
  13. 圆圈1计算机语言,LOGO编程新理念  画圆和圆弧1★
  14. 大数据技术之_16_Scala学习_11_客户信息管理系统+并发编程模型 Akka+Akka 网络编程-小黄鸡客服案例+Akka 网络编程-Spark Master Worker 进程通讯项目
  15. 微信支付分免密支付接口服务如何开通?
  16. 结合python使用ffmpeg将批量ts文件合成mp4。
  17. 关注灾情,心系九寨18
  18. python+opencv+图像特效(图像灰度处理、颜色翻转、图片融合,边缘检测,浮雕效果,颜色映射)
  19. ccf等会议排行参考
  20. 高级语言(C语言)、汇编语言、机器语言区别?编译器如何将高级语言编译成机器语言?

热门文章

  1. C#用Tesseract进行OCR识别,可识别中英日韩所有语言
  2. 博士生在没有导师指导的情况下,该如何自己选题发 CVPR ?
  3. 目标检测一卷到底之后,终于有人为它挖了个新坑|CVPR2021 Oral
  4. oracle 统计一年中每个月数据总和_excel表格有每月数据 怎样统计全年的-用excel公式怎样计算每年每个月的数据总和?...
  5. SpringBoot 多个src文件夹,入口类的位置
  6. RGB+Depth规则箱体测量方案
  7. 初探Django2.1:Win10+Python3.6.6环境下安装Django
  8. 第五课.高斯判别分析
  9. c语言判断s1是否大于s2,C语言程序设计_复习资料一.doc
  10. JGG:口腔微生物组助力妊娠糖尿病识别