jjzjj

KEGG通路的从属/注释信息如何获取

TOP生物信息 2023-09-30 原文

富集分析经常做，之前只知道GO term有从上到下的层次关系，今天才知道KEGG pathway也有类似的分层关系。

起因是我准备更新一下自己的kegg富集结果展示图，之前一直画的这种图，略显朴素了。


正好搜到了这张图，还挺好看的（主要是配色很清新）

https://doi.org/10.1016/j.phymed.2021.153714

这张图并不难（等有空把代码整理出来发给大家），可是这个kegg pathway annotation我真没见过呀！去官网看了看，能找到，但没法下载

https://www.genome.jp/kegg/pathway.html

在朋友圈求助后，费了很大劲才把这个问题解决。中间好几个小伙伴都很热心地给了我帮助，感谢他们！

以下是我整理的KEGG pathway annotation文件的网盘链接，用到的代码和数据也在里面，想改善一下kegg富集图的朋友可以看看。

链接：https://pan.baidu.com/s/18pwYZGGZSuk2_LGW8_nQFQ
提取码：ihyn
(PS: 觉得有用可以给我本期推送的第1篇推文点个赞呀，谢谢了！)

最终能得到这个表格

整理这个文件的思路如下

  1. 浏览器打开这个网页(https://www.genome.jp/kegg/pathway.html)，然后查看网页源代码（一般是鼠标右键），就能看到这个：
先不要被它吓到，后面的信息提取都是基于这些字符串

  1. 然后复制粘贴到一个文本文件kegg_html.txt，删掉前188行左右，后17行左右（删掉的这些行明显不含有用信息），得到kegg_html_copy.txt。

  2. 然后运行我编写的代码pre.R就能得到最终表格了

library(tidyverse)
tmp1=readLines("kegg_html_copy.txt")

tmp2=c()
for (li in 1:length(tmp1)) {
  if(str_detect(tmp1[li],">[0-9]")) {
    tmp2=append(tmp2,tmp1[li])
  }
}

###
big_anno=""
small_anno=""
final_lines=c()
for (li in 1:length(tmp2)) {
  if (str_detect(tmp2[li],">[0-9]\\. ")) {
    tmp_anno=str_extract(tmp2[li],">.*<")
    tmp_anno=str_split(tmp_anno,"\\. ")[[1]][2]
    tmp_anno=str_split(tmp_anno,"<")[[1]][1]
    big_anno=tmp_anno
  } else if (str_detect(tmp2[li],">[0-9]\\.[0-9]{1,2} ")) {
    tmp_anno2=str_extract(tmp2[li],">.*<")
    tmp_anno2=str_replace(tmp_anno2,"^>[0-9]\\.[0-9]{1,2} ","")
    tmp_anno2=str_split(tmp_anno2,"<")[[1]][1]
    small_anno=tmp_anno2
  } else if (str_detect(tmp2[li],">[0-9]{5}")){
    element1=str_extract(tmp2[li],">[0-9]{5}")
    element1=str_split(element1,">")[[1]][2]
    
    if (!str_detect(tmp2[li],"hsa\\+pathogen")) {
      element2=str_extract(tmp2[li],"pathway\\/[a-zA-Z]{2,4}[0-9]{5}")
      element2=str_split(element2,"\\/")[[1]][2]
      
      element3=str_extract(tmp2[li],"pathway.*?<") #非贪婪匹配
      element3=str_extract(element3,">.*<")
      element3=str_split(element3,">")[[1]][2]
      element3=str_split(element3,"<")[[1]][1]
    } else {
      element2="organism:hsa+pathogen"
      element3=str_extract(tmp2[li],"hsa\\+pathogen.*?<") #非贪婪匹配
      element3=str_extract(element3,">.*<")
      element3=str_split(element3,">")[[1]][2]
      element3=str_split(element3,"<")[[1]][1]
    }
    
    tmp_line=paste(as.character(element1),element2,element3,big_anno,small_anno,sep = ";")
    final_lines=append(final_lines,tmp_line)
  }else{
    print(tmp2[li])
  }
}

###
final_df=as.data.frame(final_lines)
colnames(final_df)="V1"
final_df=final_df%>%apply(1, function(x){as.data.frame(str_split(x,";")[[1]])})
final_df=as.data.frame(final_df)
final_df=as.data.frame(t(final_df))
rownames(final_df)=NULL
colnames(final_df)=NULL
colnames(final_df)=c("ID","Pathway Identifier","Pathway","big annotion","small annotion")

###
library(xlsx)
write.xlsx(final_df,file = "kegg_info.xlsx",col.names = T,row.names = F)
从属KEGGsectionelementtmp

有关KEGG通路的从属/注释信息如何获取的更多相关文章

  1. ruby - 如何使用 Nokogiri 的 xpath 和 at_xpath 方法 - 2

    我正在学习如何使用Nokogiri，根据这段代码我遇到了一些问题:require'rubygems'require'mechanize'post_agent=WWW::Mechanize.newpost_page=post_agent.get('http://www.vbulletin.org/forum/showthread.php?t=230708')puts"\nabsolutepathwithtbodygivesnil"putspost_page.parser.xpath('/html/body/div/div/div/div/div/table/tbody/tr/td/div

  2. ruby - 如何从 ruby​​ 中的字符串运行任意对象方法？ - 2

    总的来说，我对ruby​​还比较陌生，我正在为我正在创建的对象编写一些rspec测试用例。许多测试用例都非常基础，我只是想确保正确填充和返回值。我想知道是否有办法使用循环结构来执行此操作。不必为我要测试的每个方法都设置一个assertEquals。例如:describeitem,"TestingtheItem"doit"willhaveanullvaluetostart"doitem=Item.new#HereIcoulddotheitem.name.shouldbe_nil#thenIcoulddoitem.category.shouldbe_nilendend但我想要一些方法来使用

  3. python - 如何使用 Ruby 或 Python 创建一系列高音调和低音调的蜂鸣声？ - 2

    关闭。这个问题是opinion-based.它目前不接受答案。想要改进这个问题？更新问题，以便editingthispost可以用事实和引用来回答它.关闭4年前。Improvethisquestion我想在固定时间创建一系列低音和高音调的哔哔声。例如:在150毫秒时发出高音调的蜂鸣声在151毫秒时发出低音调的蜂鸣声200毫秒时发出低音调的蜂鸣声250毫秒的高音调蜂鸣声有没有办法在Ruby或Python中做到这一点？我真的不在乎输出编码是什么(.wav、.mp3、.ogg等等)，但我确实想创建一个输出文件。

  4. ruby-on-rails - 如何验证 update_all 是否实际在 Rails 中更新 - 2

    给定这段代码defcreate@upgrades=User.update_all(["role=?","upgraded"],:id=>params[:upgrade])redirect_toadmin_upgrades_path,:notice=>"Successfullyupgradeduser."end我如何在该操作中实际验证它们是否已保存或未重定向到适当的页面和消息？ 最佳答案 在Rails3中，update_all不返回任何有意义的信息，除了已更新的记录数(这可能取决于您的DBMS是否返回该信息)。http://ar.ru

  5. ruby-on-rails - 'compass watch' 是如何工作的/它是如何与 rails 一起使用的 - 2

    我在我的项目目录中完成了compasscreate.和compassinitrails。几个问题:我已将我的.sass文件放在public/stylesheets中。这是放置它们的正确位置吗？当我运行compasswatch时，它不会自动编译这些.sass文件。我必须手动指定文件:compasswatchpublic/stylesheets/myfile.sass等。如何让它自动运行？文件ie.css、print.css和screen.css已放在stylesheets/compiled。如何在编译后不让它们重新出现的情况下删除它们？我自己编译的.sass文件编译成compiled/t

  6. ruby-on-rails - Rails 常用字符串(用于通知和错误信息等) - 2

    大约一年前，我决定确保每个包含非唯一文本的Flash通知都将从模块中的方法中获取文本。我这样做的最初原因是为了避免一遍又一遍地输入相同的字符串。如果我想更改措辞，我可以在一个地方轻松完成，而且一遍又一遍地重复同一件事而出现拼写错误的可能性也会降低。我最终得到的是这样的:moduleMessagesdefformat_error_messages(errors)errors.map{|attribute,message|"Error:#{attribute.to_s.titleize}#{message}."}enddeferror_message_could_not_find(obje

  7. ruby - 如何将脚本文件的末尾读取为数据文件(Perl 或任何其他语言) - 2

    我正在寻找执行以下操作的正确语法(在Perl、Shell或Ruby中):#variabletoaccessthedatalinesappendedasafileEND_OF_SCRIPT_MARKERrawdatastartshereanditcontinues. 最佳答案 Perl用__DATA__做这个:#!/usr/bin/perlusestrict;usewarnings;while(){print;}__DATA__Texttoprintgoeshere 关于ruby-如何将脚

  8. ruby - 解析 RDFa、微数据等的最佳方式是什么，使用统一的模式/词汇(例如 schema.org)存储和显示信息 - 2

    我主要使用Ruby来执行此操作，但到目前为止我的攻击计划如下:使用gemsrdf、rdf-rdfa和rdf-microdata或mida来解析给定任何URI的数据。我认为最好映射到像schema.org这样的统一模式，例如使用这个yaml文件，它试图描述数据词汇表和opengraph到schema.org之间的转换:#SchemaXtoschema.orgconversion#data-vocabularyDV:name:namestreet-address:streetAddressregion:addressRegionlocality:addressLocalityphoto:i

  9. ruby - 如何指定 Rack 处理程序 - 2

    Rackup通过Rack的默认处理程序成功运行任何Rack应用程序。例如:classRackAppdefcall(environment)['200',{'Content-Type'=>'text/html'},["Helloworld"]]endendrunRackApp.new但是当最后一行更改为使用Rack的内置CGI处理程序时，rackup给出“NoMethodErrorat/undefinedmethod`call'fornil:NilClass”:Rack::Handler::CGI.runRackApp.newRack的其他内置处理程序也提出了同样的反对意见。例如Rack

  10. ruby - 如何每月在 Heroku 运行一次 Scheduler 插件？ - 2

    在选择我想要运行操作的频率时，唯一的选项是“每天”、“每小时”和“每10分钟”。谢谢!我想为我的Rails3.1应用程序运行调度程序。 最佳答案 这不是一个优雅的解决方案，但您可以安排它每天运行，并在实际开始工作之前检查日期是否为当月的第一天。 关于ruby-如何每月在Heroku运行一次Scheduler插件？，我们在StackOverflow上找到一个类似的问题： https://stackoverflow.com/questions/8692687/

随机推荐