使用R包clusterProfiler 进行富集分析与可视化
clusterProfiler官方网址:
http://www.bioconductor.org/packages/release/bioc/html/clusterProfiler.html
ID 转换功能
函数:bitr;bitr_kegg
功能:实现GeneSymbol, GeneID, Ensembl,UniprotID, ncbi-proteinid等之间的相互准换。支持约20物种。
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/64f1366cf61b4d8241d377f5f2e76c4a.png)
图片说明:上图中x是一个list,里面放的是基因名,也就是SYMBOL这一列,也可以传入其他,只要在fromType一列写上对应的类别即可。
GO/KEGG 富集分析
一条命令即可完成分析
1. GO分析
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/14476e1788d09d10ae85b52481c27b96.png)
2. KEGG分析
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/4ad52c74bb88b3d1feccd13cf25b2f67.png)
3. 输出结果
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/d0d96f49e37209ddc83a0d41cee1dd5d.png)
图片说明:第一张图为GO分析,第二张图为KEGG分析,传入的gene 是一个基因列表,该列表可以是geneID GeneSymbol 甚至是ENSEMBL, 只要在keyType 一列说明具体是哪一种类型就可以了。第三张图为输出结果,p.adjust和qvalue都是对p值的校正,p.adjust 在分析时可以选用方法,比如图中选择的是BH, qvalue是FDR校正。
GO/KEGG 富集分析可视化
直接使用上述分析结果ego 和kk,一条命令即可完成可视化
柱状图
barplot(ego, showCategory=20)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/3c0e1574ccd6e42b7d649d8f7c86a122.png)
气泡图
dotplot(kk)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/8ffd0cf39749fc42e3d45006afe0976b.png)
通路-基因网络图
cnetplot(ego, categorySize=“pvalue”, foldChange=geneList)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/5e648fdf621414e5854197c95fe2e176.png)
cnetplot(ego, foldChange=geneList, circular = TRUE, colorEdge = TRUE)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/0abbc0b1e1d2cb305aa02c5c70035301.png)
热图
heatplot(ego, foldChange=geneList)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/6e349108dcd341ea3fe58f021bb98364.png)
pathway-pathway Network
emapplot(ego)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/3a43f84ac5098fee9f4a190b2e8eab29.png)
GO通路网络图
goplot(ego)
![在这里插入图片描述](https://i-blog.csdnimg.cn/blog_migrate/a8d056810ceba7e8111aea2825000ec4.png)
KEGG通路注释
![](https://i-blog.csdnimg.cn/blog_migrate/f5616696073b7f62c7fc6c2b6e27654c.png)
图片说明:在这里,geneList 是一个耗油基因名及对应的flodchange 的列表,图中绿色即为下调,红色为上调。
其他可视化如GSEA分析等,可参照官方说明。