Advertisement
Guest User

Untitled

a guest
Jun 18th, 2019
65
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.23 KB | None | 0 0
  1. library(KEGGREST)
  2. library(png)
  3.  
  4. listDatabases() #list avail databases
  5.  
  6. keggInfo("kegg") #kegg release information
  7. keggInfo("brite") #kegg release information
  8. keggInfo("module") #kegg release information
  9.  
  10.  
  11. #list kegg maps
  12. #keggList("pathway")
  13. write.table(data.frame(keggList("pathway")), "kegg-maps.tsv", col.names=F, row.names=T, quote=F, sep="\t")
  14.  
  15. #list ECs
  16. #keggList("enzyme")
  17. write.table(data.frame(keggList('enzyme')), "kegg-ECs.tsv", col.names=F, row.names=T, quote=F, sep="\t")
  18.  
  19. #list modules
  20. #keggList("module")
  21. write.table(data.frame(keggList('module')), "kegg-modules.tsv", col.names=F, row.names=T, quote=F, sep="\t")
  22.  
  23. #list kegg orthologies
  24. #keggList("ko")
  25. write.table(data.frame(keggList('ko')), "kegg-orthologies.tsv", col.names=F, row.names=T, quote=F, sep="\t")
  26.  
  27. #list kegg orthologies
  28. #keggList("brite")
  29. write.table(data.frame(keggList('brite')), "kegg-BRITE.tsv", col.names=F, row.names=T, quote=F, sep="\t")
  30.  
  31. #list kegg orthologies
  32. #keggList("genome")
  33. write.table(data.frame(keggList('genome')), "kegg-genome.tsv", col.names=F, row.names=T, quote=F, sep="\t")
  34.  
  35.  
  36.  
  37.  
  38.  
  39.  
  40.  
  41. temp2<-data.frame(keggConv("afu", "ncbi-proteinid"))
  42.  
  43. keggFind("module", "map00680")
  44.  
  45. keggFind("module", "methan")
  46.  
  47.  
  48. png<-keggGet("path:map01200", "aaseq")
  49. write(png,"test.png")
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement