多个数据集整合神器RobustRankAggreg包的示例分析

今天就跟大家聊聊有关多个数据集整合神器RobustRankAggreg包的示例分析，可能很多人都不太了解，为了让大家更加了解，小编给大家总结了以下内容，希望大家根据这篇文章可以有所收获。

成都创新互联专注为客户提供全方位的互联网综合服务，包含不限于成都做网站、成都网站设计、江城网络推广、微信小程序定制开发、江城网络营销、江城企业策划、江城品牌公关、搜索引擎seo、人物专访、企业宣传片、企业代运营等，从售前售中售后，我们都将竭诚为您服务，您的肯定，是我们最大的嘉奖；成都创新互联为所有大学生创业者提供江城建站搭建服务，24小时服务热线：13518219792，官方网址：www.cdcxhl.com

RobustRankAggreg包在各种数据挖掘文章里面亮相的频次之高，无需我多言，大家可以去查看一下引用它的文献，基本上都是GEO数据库挖掘文章：

RobustRankAggreg文章

比如发表在peerJ的BIOINFORMATICS AND GENOMICS的文章Identification of key candidate genes and biological pathways in bladder cancer 里面的：

多个数据集整合神器RobustRankAggreg包的示例分析

4个GEO数据集

作者把这4个数据集，分别独立走差异分析，火山图，热图等等标准流程，基本上读一下我在生信技能树的表达芯片的公共数据库挖掘系列推文就明白了；

解读GEO数据存放规律及下载，一文就够
解读SRA数据库规律一文就够
从GEO数据库下载得到表达矩阵一文就够
GSEA分析一文就够（单机版+R语言版）
根据分组信息做差异分析- 这个一文不够的
差异分析得到的结果注释一文就够

你也可以很轻松的分析这几个数据集：GSE7476, GSE13507, GSE37815 and GSE65635 ，然后作者就使用了RobustRankAggreg包对这4个数据集的差异分析结果进行整合，如下：

The integrated DEGs were screened using the RRA package (corrected P < 0.05, logFC > 1 or −logFC < −1).
The RRA method is based on the assumption that each gene in each dataset is randomly arranged.
If the gene ranks high in all datasets, the associated P-value is lower, the possibility of differential gene expression is greater.
Through rank analysis, 343 integrated DEGs, consisting of 111 upregulated genes and 232 downregulated genes, were identified by the RRA method

并且把top20的上调基因和下调基因的差异倍数进行热图可视化，如下：

多个数据集整合神器RobustRankAggreg包的示例分析

top20的上调基因和下调基因的差异倍数进行热图可视化

当然了，不仅仅是mRNA的表达芯片，其它，比如circRNA芯片也是如此，同样是发表于2018的文章：A circRNA–miRNA–mRNA network identification for exploring underlying pathogenesis and therapy strategy of hepatocellular carcinoma

就是下载了3个GEO数据集，走差异分析，并且使用RobustRankAggreg包进行整合，最后仅仅是确定了6个circRNA。

多个数据集整合神器RobustRankAggreg包的示例分析

circRNA芯片整合

几百篇文章我们就不用一一解读啦，反正都是独立的数据集自己做自己的差异分析，然后把多个数据集的差异基因拿去使用RobustRankAggreg包进行整合。

RobustRankAggreg包说明书

这个RobustRankAggreg包超级简单，有意思的是居然并不在bioconductor列表哦，可能是因为它最开始并不是为生物信息学领域的数据分析而创造的吧！因为不在bioconductor，所以它的示例教程一塌糊涂，需要一点背景才能理解。其重点就是aggregateRanks函数而已：

options(BioC_mirror="https://mirrors.tuna.tsinghua.edu.cn/bioconductor/")
options("repos" = c(CRAN="https://mirrors.tuna.tsinghua.edu.cn/CRAN/"))
# https://bioconductor.org/packages/release/bioc/html/GEOquery.html
if (!requireNamespace("BiocManager", quietly = TRUE))
  install.packages("BiocManager")
BiocManager::install("RobustRankAggreg",ask = F,update = F) 
library(RobustRankAggreg)
?aggregateRanks

一般来说，正常R包的函数，都是可以通过问号来调取其帮助文档的，aggregateRanks函数也不例外。我们直接看一下示例代码：

set.seed(1234567)
glist <- list(sample(letters, 4), sample(letters, 10), sample(letters, 12))
freq=as.data.frame(table(unlist(glist)))
# Aggregate the inputs
ag=aggregateRanks(glist = glist, N = length(letters))
ag$Freq=freq[match(ag$Name,freq$Var1),2]

的确是超级简单，可以看到，我们有26个字母，假设是26个基因，然后做了3次随机抽样，假设是3个数据集的差异分析，拿到的上调基因，列表如下：

多个数据集整合神器RobustRankAggreg包的示例分析