kegg富集分析之:KEGGREST包(9大功能)
这个包依赖极有可能是这个:https://www.kegg.jp/kegg/docs/keggapi.html ,如果可以看懂会很好理解
由于KEGG数据库分享数据的策略改变,因此KEGG.db包不在能用,推荐KEGGREST包
But a number of years ago,KEGG changed their policy about sharing their data and so the KEGG.db package is no longer allowed to be current. Users who are interested in a more current pathway data are encouraged to look at the KEGGREST or reactome.db packages.
1、安装
if("KEGGREST" %in% rownames(installed.packages()) == FALSE) {source("http://bioconductor.org/biocLite.R");biocLite("KEGGREST")}
suppressMessages(library(KEGGREST))
ls('package:KEGGREST')
2、所有的对象及功能
keggConv:Convert KEGG identifiers to/from outside identifiers(ID转换功能)
keggFind:Finds entries with matching query keywords or other query data in a given database(即搜索功能)
keggGet:Retrieves given database entries(获取序列,图片等功能)
keggInfo:Displays the current statistics of a given database(即统计功能)
keggLink:Find related entries by using database cross-references(交互功能).
keggList:Returns a list of entry identifiers and associated definition for a given database or a given set of database entries.
listDatabases:Lists the KEGG databases which may be searched.
mark.pathway.by.objects:Client-side interface to obtain an url for a KEGG pathway diagram with a given set of genes marked(标记功能)
3、每个对象的功能简单使用
3.1、keggConv(Convert KEGG identifiers to/from outside identifiers,即ID转换功能)
语法:keggConv(target, source, querySize = 100)
例子:
## conversion from NCBI GeneID to KEGG ID for E. coli genes##
head(keggConv("eco", "ncbi-geneid"),2)
head(keggConv("ncbi-geneid", "eco"),2)
########conversion from KEGG ID to NCBI GI##########
head(keggConv("ncbi-proteinid", c("hsa:10458", "ece:Z5100")),2)
3.2)keggFind(Finds entries with matching query keywords or other query data in a given database,即检索功能)
语法:keggFind(database, query, option = c("formula", "exact_mass", "mol_weight"))
例子:
keggFind("genes", c("shiga", "toxin")) ## for keywords "shiga" and "toxin"
keggFind("genes", "shiga toxin") ## for keywords "shiga toxin"
keggFind("compound", "C7H10O5", "formula") ## for chemical formula "C7H10O5"
keggFind("compound", 174.05, "exact_mass") ## for 174.045 =< exact mass < 174.055
keggFind("compound", 300:310, "mol_weight") ## molecular weight 300 =< 310
3.3)keggGet (Retrieves given database entries,获取序列,图片等功能)
语法:keggGet(dbentries, option = c("aaseq", "ntseq", "mol", "kcf","image", "kgml"))
keggGet(c("cpd:C01290", "gl:G00092")) ## retrieves a compound entry and a glycan entry
keggGet(c("C01290", "G00092")) ## same as above, without prefixes
keggGet(c("hsa:10458", "ece:Z5100")) ## 检索 a human gene entry and an E.coli O157 gene entry
keggGet(c("hsa:10458", "ece:Z5100"), "aaseq") ## 检索 amino acid sequences of a human gene and an E.coli O157 gene
png <- keggGet("hsa05130", "image") ## retrieves the image file of a pathway map
t <- 'hsa05130.png'
library(png)
writePNG(png, t)
keggGet("hsa05130", "kgml")
hsa05130.png
3.4)keggInfo(Displays the current statistics of a given database,即统计功能)
语法:keggInfo(database)
head(keggInfo("kegg") ## displays the current statistics of the KEGG database
keggInfo("pathway") ## displays the number pathway entries including both the reference and organism-specific pathways
keggInfo("hsa") ##### displays the number of gene entries for the KEGG organism Homo sapiens
3.5)keggLink(keggLink Find related entries by using database cross-references.交互功能)
语法:keggLink(target, source)
head(keggLink("pathway", "hsa"),5) ### KEGG pathways linked from each of the human genes
head(keggLink("hsa", "pathway"),5) ## human genes linked from each of the KEGG pathways
head(keggLink("pathway", c("hsa:10458", "ece:Z5100")),5) ## KEGG pathways linked from a human gene and an E. coli O157 gene
head(keggLink("hsa:126"),5)
3.6) keggList (Returns a list of entry identifiers and associated definition for a given database or set of database entries,返回信息表)
用法:keggList(database, organism)
keggList("pathway") ## returns the list of reference pathways
keggList("pathway", "hsa") ## returns the list of human pathways
keggList("organism") ## returns the list of KEGG organisms with taxonomic classification
keggList("hsa") ## returns the entire list of human genes
keggList("T01001") ## same as above
keggList(c("hsa:10458", "ece:Z5100")) ## returns the list of a human gene and an E.coli O157 gene
keggList(c("cpd:C01290","gl:G00092")) ## returns the list of a compound entry and a glycan entry
keggList(c("C01290+G00092")) ## same as above (prefixes are not necessary)
3.7) listDatabases(Lists the KEGG databases which may be searched.查看有数据库)
用法:listDatabases()
例子:
listDatabases() ##查看有哪些数据库
3.8)mark.pathway.by.objects:Client-side interface to obtain an url for a KEGG pathway diagram with a given set of genes marked(标记功能,例如上调y用红框,下调用绿框)
语法:
mark.pathway.by.objects(pathway.id, object.id.list)
color.pathway.by.objects(pathway.id, object.id.list,fg.color.list, bg.color.list)
例子:
url <- mark.pathway.by.objects("path:eco00260",c("eco:b0002", "eco:c00263"))
if(interactive()){browseURL(url)}
url <- color.pathway.by.objects("path:eco00260",c("eco:b0002", "eco:c00263"),c("#ff0000", "#00ff00"), c("#ffff00", "yellow"))
kegg富集分析之:KEGGREST包(9大功能)的更多相关文章
- 【R】clusterProfiler的GO/KEGG富集分析用法小结
前言 关于clusterProfiler这个R包就不介绍了,网红教授宣传得很成功,功能也比较强大,主要是做GO和KEGG的功能富集及其可视化.简单总结下用法,以后用时可直接找来用. 首先考虑一个问题: ...
- KEGG富集分析散点图.md
输入数据格式 pathway = read.table("kegg.result",header=T,sep="\t") pp = ggplot(pathway ...
- DAVID 进行 GO/KEGG 功能富集分析
何为功能富集分析? 功能富集分析是将基因或者蛋白列表分成多个部分,即将一堆基因进行分类,而这里的分类标准往往是按照基因的功能来限定的.换句话说,就是把一个基因列表中,具有相似功能的基因放到一起,并和生 ...
- python scipy包进行GO富集分析p值计算
最近总是有需要单独对某一个类型的通路进行超几何分布的p值计算,这里记录一下python包的计算方法 使用scipy的stat里面的hypergeom.sf方法进行富集分析的p值计算 hsaxxxxx ...
- 富集分析DAVID、Metascape、Enrichr、ClueGO
前言 一般我们挑出一堆感兴趣的基因想临时看看它们的功能,需要做个富集分析.虽然公司买了最新版的数据库,如KEGG,但在集群跑下来嫌麻烦.这时网页在线或者本地化工具派上用场了. DAVID DAVID地 ...
- R: 修改镜像、bioconductor安装及go基因富集分析
1.安装bioconductor及go分析涉及的相关包 source("http://bioconductor.org/biocLite.R") options(BioC_mirr ...
- GO富集分析示例【华为云技术分享】
版权声明:本文为博主原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接和本声明. 本文链接:https://blog.csdn.net/devcloud/article/detai ...
- DEPICT实现基因优化(gene prioritization)、gene set富集分析(geneset enrichment)、组织富集分析(tissue enrichment)
全基因组关联分析除了找到显著的关联位点,我们还可以做基因优化.geneset富集分析.组织富集分析,下面具体讲一讲怎么利用GWAS的summary数据做这个分析. summary数据就是关联分析的结果 ...
- GO富集分析 信号通路
基因富集分析是分析基因表达信息的一种方法,富集是指将基因按照先验知识,也就是基因组注释信息进行分类. 信号通路是指能将细胞外的分子信号经细胞膜传入细胞内发挥效应的一系列酶促反应通路.这些细胞外的分子信 ...
随机推荐
- SqlBulkCopy(批量复制)使用方法
SqlBulkCopy提供了一种将数据复制到Sql Server数据库表中高性能的方法.SqlBulkCopy 包含一个方法 WriteToServer,它用来从数据的源复制数据到数据的目的地. Wr ...
- 创建mysql 用户并限定其操作主机 和 数据库
参考链接 http://www.cnblogs.com/top5/archive/2010/09/14/1825571.html ******** GRANT ALL ON push.* TO pus ...
- php webservice实例(转载)
首先大家要简单了解了何谓webservice,接下来就做两个非常简单的例子,webservice还是逃不开server端(服务器端) 与client端(客户端) . 做这个测试之前,要确认你的php配 ...
- JAVA课程设计——多源教学数据管理系统
团队简介 团队名称: 419圣斗士 团队成员 姓名 成员介绍 任务分配 周炳辉(组长) 来自网络的一个大佬,穿女装很合适 poi与servlet 徐宏伟 网络中一个具有强大隐藏实力的大哥 css,部分 ...
- Spring+hibernate+mysql事物不回滚的原因以及处理
最近项目突然出了点问题,然后发现用Service层下面的一个类的一个方法里的事务居然没有回滚.然后自己写了一个测试方法经过了N次测试都是不回滚.以下是测试方法的一部分: @Transactional( ...
- WDA编译失败问题
1.放假回来,wda编译失败,报错如下 2018-09-25 10:03:09.020964+0800 WebDriverAgentRunner-Runner[335:33309] +[CATrans ...
- ios ideviceintaller安装
1.安装brew 打开终端输入: /usr/bin/ruby -e "$(curl -fsSL https://raw.githubusercontent.com/Homebrew/inst ...
- 接口测试maven管理
接口测试框架选择 界面化工具,针对不会编码的测试人员: 1.Jmeter性能测试工具,不具备完备的接口测试框架功能 2.Robotframerwork 3.PostMan 推荐框架: ResrAssu ...
- Leetcode 之Simplify Path @ python
Given an absolute path for a file (Unix-style), simplify it. For example,path = "/home/", ...
- Ps操作技巧(快捷键大全)
一.工具箱(多种工具共用一个快捷键的可同时按[Shift]加此快捷键选取) 矩形.椭圆选框工具 [M] 移动工具 [V] 套索.多边形套索.磁性套索 [L] 魔棒工具 [W] 裁剪工具 [C] 切片工 ...