北京大学生物信息平台论坛

 找回密码
 立即注册
搜索
热搜: 通知 活动

差异表达基因富集分析容易使用的包 enrich

[复制链接]
wolf 发表于 2016-8-11 16:27:10 | 显示全部楼层 |阅读模式
本帖最后由 wolf 于 2016-8-11 16:32 编辑

enrich

enrich is a R package used for gene enrichment analysis.

Installation

This package can be installed from Github.
  1. devtools::install_github("wolfsonliu/enrich")
复制代码
Usage

For easy use, the most usefull functions are enrichKEGGPathway and enrichMSigDB. enrichKEGGPathway is used for pathway enrichment analysis. Input gene lists and pathway lists you care, and the function will output p value of enrichment. Make sure you have the internet connection, and function enrichKEGGPathway will automatically download KEGG pathway information from KEGG website.
  1. data("deg")

  2. de.entrez <- as.character(na.omit(deg$entrezgene))

  3. enrich.result <- enrichKEGGPathway(
  4.     entrez.gene.list  = de.entrez,
  5.     pathway.list      = NA,
  6.     whole.gene.number = 20000,
  7.     threshold         = 0.05,
  8.     adjust.p.method   = "BH",
  9.     organism          = "hsa"
  10. )

  11. enrich.result
复制代码
enrichMSigDB is used for enrichment of gene sets in MSigDB. Just like enrichKEGGPathway, input gene lists, and you will get the result. You can imput the msigdb data with appropriate data structure or use the data in the package. If you do not input msigdb data, the function will use the data from the package automatically. Gene sets in MSigDB include: hallmark gene sets, positional gene sets, curated gene sets, motif gene sets, computational gene sets, GO gene sets, oncogenic signatures, and immunologic signatures.
  1. data("deg")

  2. de.entrez <- as.character(na.omit(deg$entrezgene))

  3. enrich.result <- enrichMSigDB(
  4.     entrez.gene.list  = de.entrez,
  5.     msigdb            = NA,
  6.     whole.gene.number = 20000,
  7.     threshold         = 0.05,
  8.     adjust.p.method   = "BH"
  9. )

  10. enrich.result
复制代码
There are also other functions, you can find the useage in R with help function.

Information

Author: Zhiheng Liu
Version: 0.1
License: LGPL-3
Please contact me if you find bug in the enrich R package. Thanks!
回复

使用道具 举报

北京大学生物信息平台论坛

GMT+8, 2017-11-19 22:01 , Processed in 0.092196 second(s), 23 queries .

Powered by Discuz! X3

© 2001-2013 Comsenz Inc.

快速回复 返回顶部 返回列表