wwood-bioruby 1.2.11

Sign up to get free protection for your applications and to get access to all the features.
Files changed (386) hide show
  1. data/README.rdoc +205 -0
  2. data/README_DEV.rdoc +285 -0
  3. data/VERSION.yml +4 -0
  4. data/bin/bioruby +44 -0
  5. data/bin/br_biofetch.rb +47 -0
  6. data/bin/br_bioflat.rb +293 -0
  7. data/bin/br_biogetseq.rb +45 -0
  8. data/bin/br_pmfetch.rb +421 -0
  9. data/lib/bio.rb +306 -0
  10. data/lib/bio/alignment.rb +2518 -0
  11. data/lib/bio/appl/bl2seq/report.rb +334 -0
  12. data/lib/bio/appl/blast.rb +505 -0
  13. data/lib/bio/appl/blast/ddbj.rb +142 -0
  14. data/lib/bio/appl/blast/format0.rb +1438 -0
  15. data/lib/bio/appl/blast/format8.rb +83 -0
  16. data/lib/bio/appl/blast/genomenet.rb +263 -0
  17. data/lib/bio/appl/blast/ncbioptions.rb +220 -0
  18. data/lib/bio/appl/blast/remote.rb +105 -0
  19. data/lib/bio/appl/blast/report.rb +767 -0
  20. data/lib/bio/appl/blast/rexml.rb +144 -0
  21. data/lib/bio/appl/blast/rpsblast.rb +277 -0
  22. data/lib/bio/appl/blast/wublast.rb +635 -0
  23. data/lib/bio/appl/blast/xmlparser.rb +236 -0
  24. data/lib/bio/appl/blat/report.rb +530 -0
  25. data/lib/bio/appl/clustalw.rb +219 -0
  26. data/lib/bio/appl/clustalw/report.rb +152 -0
  27. data/lib/bio/appl/emboss.rb +203 -0
  28. data/lib/bio/appl/fasta.rb +235 -0
  29. data/lib/bio/appl/fasta/format10.rb +325 -0
  30. data/lib/bio/appl/gcg/msf.rb +212 -0
  31. data/lib/bio/appl/gcg/seq.rb +195 -0
  32. data/lib/bio/appl/genscan/report.rb +552 -0
  33. data/lib/bio/appl/hmmer.rb +126 -0
  34. data/lib/bio/appl/hmmer/report.rb +683 -0
  35. data/lib/bio/appl/iprscan/report.rb +374 -0
  36. data/lib/bio/appl/mafft.rb +259 -0
  37. data/lib/bio/appl/mafft/report.rb +226 -0
  38. data/lib/bio/appl/muscle.rb +52 -0
  39. data/lib/bio/appl/paml/baseml.rb +95 -0
  40. data/lib/bio/appl/paml/baseml/report.rb +32 -0
  41. data/lib/bio/appl/paml/codeml.rb +242 -0
  42. data/lib/bio/appl/paml/codeml/rates.rb +67 -0
  43. data/lib/bio/appl/paml/codeml/report.rb +67 -0
  44. data/lib/bio/appl/paml/common.rb +348 -0
  45. data/lib/bio/appl/paml/common_report.rb +38 -0
  46. data/lib/bio/appl/paml/yn00.rb +103 -0
  47. data/lib/bio/appl/paml/yn00/report.rb +32 -0
  48. data/lib/bio/appl/phylip/alignment.rb +133 -0
  49. data/lib/bio/appl/phylip/distance_matrix.rb +96 -0
  50. data/lib/bio/appl/probcons.rb +41 -0
  51. data/lib/bio/appl/psort.rb +548 -0
  52. data/lib/bio/appl/psort/report.rb +542 -0
  53. data/lib/bio/appl/pts1.rb +263 -0
  54. data/lib/bio/appl/sim4.rb +124 -0
  55. data/lib/bio/appl/sim4/report.rb +485 -0
  56. data/lib/bio/appl/sosui/report.rb +151 -0
  57. data/lib/bio/appl/spidey/report.rb +593 -0
  58. data/lib/bio/appl/targetp/report.rb +267 -0
  59. data/lib/bio/appl/tcoffee.rb +55 -0
  60. data/lib/bio/appl/tmhmm/report.rb +231 -0
  61. data/lib/bio/command.rb +593 -0
  62. data/lib/bio/compat/features.rb +157 -0
  63. data/lib/bio/compat/references.rb +128 -0
  64. data/lib/bio/data/aa.rb +353 -0
  65. data/lib/bio/data/codontable.rb +722 -0
  66. data/lib/bio/data/na.rb +223 -0
  67. data/lib/bio/db.rb +329 -0
  68. data/lib/bio/db/aaindex.rb +357 -0
  69. data/lib/bio/db/biosql/biosql_to_biosequence.rb +67 -0
  70. data/lib/bio/db/biosql/sequence.rb +508 -0
  71. data/lib/bio/db/embl/common.rb +352 -0
  72. data/lib/bio/db/embl/embl.rb +500 -0
  73. data/lib/bio/db/embl/embl_to_biosequence.rb +85 -0
  74. data/lib/bio/db/embl/format_embl.rb +190 -0
  75. data/lib/bio/db/embl/sptr.rb +1283 -0
  76. data/lib/bio/db/embl/swissprot.rb +42 -0
  77. data/lib/bio/db/embl/trembl.rb +41 -0
  78. data/lib/bio/db/embl/uniprot.rb +42 -0
  79. data/lib/bio/db/fantom.rb +597 -0
  80. data/lib/bio/db/fasta.rb +410 -0
  81. data/lib/bio/db/fasta/defline.rb +532 -0
  82. data/lib/bio/db/fasta/fasta_to_biosequence.rb +63 -0
  83. data/lib/bio/db/fasta/format_fasta.rb +97 -0
  84. data/lib/bio/db/genbank/common.rb +307 -0
  85. data/lib/bio/db/genbank/ddbj.rb +22 -0
  86. data/lib/bio/db/genbank/format_genbank.rb +187 -0
  87. data/lib/bio/db/genbank/genbank.rb +250 -0
  88. data/lib/bio/db/genbank/genbank_to_biosequence.rb +86 -0
  89. data/lib/bio/db/genbank/genpept.rb +60 -0
  90. data/lib/bio/db/genbank/refseq.rb +18 -0
  91. data/lib/bio/db/gff.rb +1846 -0
  92. data/lib/bio/db/go.rb +481 -0
  93. data/lib/bio/db/kegg/brite.rb +41 -0
  94. data/lib/bio/db/kegg/compound.rb +131 -0
  95. data/lib/bio/db/kegg/drug.rb +98 -0
  96. data/lib/bio/db/kegg/enzyme.rb +148 -0
  97. data/lib/bio/db/kegg/expression.rb +155 -0
  98. data/lib/bio/db/kegg/genes.rb +263 -0
  99. data/lib/bio/db/kegg/genome.rb +241 -0
  100. data/lib/bio/db/kegg/glycan.rb +166 -0
  101. data/lib/bio/db/kegg/keggtab.rb +357 -0
  102. data/lib/bio/db/kegg/kgml.rb +256 -0
  103. data/lib/bio/db/kegg/orthology.rb +136 -0
  104. data/lib/bio/db/kegg/reaction.rb +82 -0
  105. data/lib/bio/db/kegg/taxonomy.rb +331 -0
  106. data/lib/bio/db/lasergene.rb +209 -0
  107. data/lib/bio/db/litdb.rb +107 -0
  108. data/lib/bio/db/medline.rb +326 -0
  109. data/lib/bio/db/nbrf.rb +191 -0
  110. data/lib/bio/db/newick.rb +658 -0
  111. data/lib/bio/db/nexus.rb +1854 -0
  112. data/lib/bio/db/pdb.rb +29 -0
  113. data/lib/bio/db/pdb/atom.rb +77 -0
  114. data/lib/bio/db/pdb/chain.rb +210 -0
  115. data/lib/bio/db/pdb/chemicalcomponent.rb +224 -0
  116. data/lib/bio/db/pdb/model.rb +148 -0
  117. data/lib/bio/db/pdb/pdb.rb +1911 -0
  118. data/lib/bio/db/pdb/residue.rb +176 -0
  119. data/lib/bio/db/pdb/utils.rb +399 -0
  120. data/lib/bio/db/prosite.rb +597 -0
  121. data/lib/bio/db/rebase.rb +456 -0
  122. data/lib/bio/db/soft.rb +404 -0
  123. data/lib/bio/db/transfac.rb +375 -0
  124. data/lib/bio/db/url.rb +42 -0
  125. data/lib/bio/feature.rb +139 -0
  126. data/lib/bio/io/biosql/biodatabase.rb +64 -0
  127. data/lib/bio/io/biosql/bioentry.rb +29 -0
  128. data/lib/bio/io/biosql/bioentry_dbxref.rb +11 -0
  129. data/lib/bio/io/biosql/bioentry_path.rb +12 -0
  130. data/lib/bio/io/biosql/bioentry_qualifier_value.rb +10 -0
  131. data/lib/bio/io/biosql/bioentry_reference.rb +10 -0
  132. data/lib/bio/io/biosql/bioentry_relationship.rb +10 -0
  133. data/lib/bio/io/biosql/biosequence.rb +11 -0
  134. data/lib/bio/io/biosql/comment.rb +7 -0
  135. data/lib/bio/io/biosql/config/database.yml +20 -0
  136. data/lib/bio/io/biosql/dbxref.rb +13 -0
  137. data/lib/bio/io/biosql/dbxref_qualifier_value.rb +12 -0
  138. data/lib/bio/io/biosql/location.rb +32 -0
  139. data/lib/bio/io/biosql/location_qualifier_value.rb +11 -0
  140. data/lib/bio/io/biosql/ontology.rb +10 -0
  141. data/lib/bio/io/biosql/reference.rb +9 -0
  142. data/lib/bio/io/biosql/seqfeature.rb +32 -0
  143. data/lib/bio/io/biosql/seqfeature_dbxref.rb +11 -0
  144. data/lib/bio/io/biosql/seqfeature_path.rb +11 -0
  145. data/lib/bio/io/biosql/seqfeature_qualifier_value.rb +20 -0
  146. data/lib/bio/io/biosql/seqfeature_relationship.rb +11 -0
  147. data/lib/bio/io/biosql/taxon.rb +12 -0
  148. data/lib/bio/io/biosql/taxon_name.rb +9 -0
  149. data/lib/bio/io/biosql/term.rb +27 -0
  150. data/lib/bio/io/biosql/term_dbxref.rb +11 -0
  151. data/lib/bio/io/biosql/term_path.rb +12 -0
  152. data/lib/bio/io/biosql/term_relationship.rb +13 -0
  153. data/lib/bio/io/biosql/term_relationship_term.rb +11 -0
  154. data/lib/bio/io/biosql/term_synonym.rb +10 -0
  155. data/lib/bio/io/das.rb +461 -0
  156. data/lib/bio/io/dbget.rb +194 -0
  157. data/lib/bio/io/ddbjxml.rb +638 -0
  158. data/lib/bio/io/ebisoap.rb +158 -0
  159. data/lib/bio/io/ensembl.rb +229 -0
  160. data/lib/bio/io/fastacmd.rb +163 -0
  161. data/lib/bio/io/fetch.rb +195 -0
  162. data/lib/bio/io/flatfile.rb +482 -0
  163. data/lib/bio/io/flatfile/autodetection.rb +545 -0
  164. data/lib/bio/io/flatfile/bdb.rb +253 -0
  165. data/lib/bio/io/flatfile/buffer.rb +237 -0
  166. data/lib/bio/io/flatfile/index.rb +1381 -0
  167. data/lib/bio/io/flatfile/indexer.rb +805 -0
  168. data/lib/bio/io/flatfile/splitter.rb +297 -0
  169. data/lib/bio/io/higet.rb +73 -0
  170. data/lib/bio/io/hinv.rb +442 -0
  171. data/lib/bio/io/keggapi.rb +805 -0
  172. data/lib/bio/io/ncbirest.rb +733 -0
  173. data/lib/bio/io/ncbisoap.rb +155 -0
  174. data/lib/bio/io/pubmed.rb +307 -0
  175. data/lib/bio/io/registry.rb +292 -0
  176. data/lib/bio/io/soapwsdl.rb +119 -0
  177. data/lib/bio/io/sql.rb +186 -0
  178. data/lib/bio/location.rb +867 -0
  179. data/lib/bio/map.rb +410 -0
  180. data/lib/bio/pathway.rb +960 -0
  181. data/lib/bio/reference.rb +602 -0
  182. data/lib/bio/sequence.rb +456 -0
  183. data/lib/bio/sequence/aa.rb +152 -0
  184. data/lib/bio/sequence/adapter.rb +108 -0
  185. data/lib/bio/sequence/common.rb +310 -0
  186. data/lib/bio/sequence/compat.rb +123 -0
  187. data/lib/bio/sequence/dblink.rb +54 -0
  188. data/lib/bio/sequence/format.rb +358 -0
  189. data/lib/bio/sequence/format_raw.rb +23 -0
  190. data/lib/bio/sequence/generic.rb +24 -0
  191. data/lib/bio/sequence/na.rb +491 -0
  192. data/lib/bio/shell.rb +44 -0
  193. data/lib/bio/shell/core.rb +578 -0
  194. data/lib/bio/shell/demo.rb +146 -0
  195. data/lib/bio/shell/interface.rb +218 -0
  196. data/lib/bio/shell/irb.rb +95 -0
  197. data/lib/bio/shell/object.rb +71 -0
  198. data/lib/bio/shell/plugin/blast.rb +42 -0
  199. data/lib/bio/shell/plugin/codon.rb +218 -0
  200. data/lib/bio/shell/plugin/das.rb +58 -0
  201. data/lib/bio/shell/plugin/emboss.rb +23 -0
  202. data/lib/bio/shell/plugin/entry.rb +105 -0
  203. data/lib/bio/shell/plugin/flatfile.rb +101 -0
  204. data/lib/bio/shell/plugin/keggapi.rb +181 -0
  205. data/lib/bio/shell/plugin/midi.rb +430 -0
  206. data/lib/bio/shell/plugin/obda.rb +45 -0
  207. data/lib/bio/shell/plugin/psort.rb +56 -0
  208. data/lib/bio/shell/plugin/seq.rb +247 -0
  209. data/lib/bio/shell/plugin/soap.rb +87 -0
  210. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/bioruby_generator.rb +29 -0
  211. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/_classes.rhtml +4 -0
  212. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/_log.rhtml +27 -0
  213. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/_methods.rhtml +11 -0
  214. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/_modules.rhtml +4 -0
  215. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/_variables.rhtml +7 -0
  216. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/bioruby-bg.gif +0 -0
  217. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/bioruby-gem.png +0 -0
  218. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/bioruby-link.gif +0 -0
  219. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/bioruby.css +368 -0
  220. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/bioruby.rhtml +47 -0
  221. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/bioruby_controller.rb +144 -0
  222. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/bioruby_helper.rb +47 -0
  223. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/commands.rhtml +8 -0
  224. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/history.rhtml +10 -0
  225. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/index.rhtml +26 -0
  226. data/lib/bio/shell/rails/vendor/plugins/bioruby/generators/bioruby/templates/spinner.gif +0 -0
  227. data/lib/bio/shell/script.rb +25 -0
  228. data/lib/bio/shell/setup.rb +109 -0
  229. data/lib/bio/shell/web.rb +102 -0
  230. data/lib/bio/tree.rb +852 -0
  231. data/lib/bio/util/color_scheme.rb +191 -0
  232. data/lib/bio/util/color_scheme/buried.rb +59 -0
  233. data/lib/bio/util/color_scheme/helix.rb +59 -0
  234. data/lib/bio/util/color_scheme/hydropathy.rb +64 -0
  235. data/lib/bio/util/color_scheme/nucleotide.rb +31 -0
  236. data/lib/bio/util/color_scheme/strand.rb +59 -0
  237. data/lib/bio/util/color_scheme/taylor.rb +50 -0
  238. data/lib/bio/util/color_scheme/turn.rb +59 -0
  239. data/lib/bio/util/color_scheme/zappo.rb +50 -0
  240. data/lib/bio/util/contingency_table.rb +370 -0
  241. data/lib/bio/util/restriction_enzyme.rb +228 -0
  242. data/lib/bio/util/restriction_enzyme/analysis.rb +249 -0
  243. data/lib/bio/util/restriction_enzyme/analysis_basic.rb +217 -0
  244. data/lib/bio/util/restriction_enzyme/cut_symbol.rb +107 -0
  245. data/lib/bio/util/restriction_enzyme/double_stranded.rb +321 -0
  246. data/lib/bio/util/restriction_enzyme/double_stranded/aligned_strands.rb +130 -0
  247. data/lib/bio/util/restriction_enzyme/double_stranded/cut_location_pair.rb +103 -0
  248. data/lib/bio/util/restriction_enzyme/double_stranded/cut_location_pair_in_enzyme_notation.rb +38 -0
  249. data/lib/bio/util/restriction_enzyme/double_stranded/cut_locations.rb +76 -0
  250. data/lib/bio/util/restriction_enzyme/double_stranded/cut_locations_in_enzyme_notation.rb +107 -0
  251. data/lib/bio/util/restriction_enzyme/enzymes.yaml +7061 -0
  252. data/lib/bio/util/restriction_enzyme/range/cut_range.rb +24 -0
  253. data/lib/bio/util/restriction_enzyme/range/cut_ranges.rb +47 -0
  254. data/lib/bio/util/restriction_enzyme/range/horizontal_cut_range.rb +67 -0
  255. data/lib/bio/util/restriction_enzyme/range/sequence_range.rb +257 -0
  256. data/lib/bio/util/restriction_enzyme/range/sequence_range/calculated_cuts.rb +242 -0
  257. data/lib/bio/util/restriction_enzyme/range/sequence_range/fragment.rb +51 -0
  258. data/lib/bio/util/restriction_enzyme/range/sequence_range/fragments.rb +41 -0
  259. data/lib/bio/util/restriction_enzyme/range/vertical_cut_range.rb +77 -0
  260. data/lib/bio/util/restriction_enzyme/single_strand.rb +200 -0
  261. data/lib/bio/util/restriction_enzyme/single_strand/cut_locations_in_enzyme_notation.rb +135 -0
  262. data/lib/bio/util/restriction_enzyme/single_strand_complement.rb +23 -0
  263. data/lib/bio/util/restriction_enzyme/string_formatting.rb +111 -0
  264. data/lib/bio/util/sirna.rb +288 -0
  265. data/test/data/HMMER/hmmpfam.out +64 -0
  266. data/test/data/HMMER/hmmsearch.out +88 -0
  267. data/test/data/SOSUI/sample.report +11 -0
  268. data/test/data/TMHMM/sample.report +21 -0
  269. data/test/data/aaindex/DAYM780301 +30 -0
  270. data/test/data/aaindex/PRAM900102 +20 -0
  271. data/test/data/bl2seq/cd8a_cd8b_blastp.bl2seq +53 -0
  272. data/test/data/bl2seq/cd8a_p53_e-5blastp.bl2seq +37 -0
  273. data/test/data/blast/2.2.15.blastp.m7 +876 -0
  274. data/test/data/blast/b0002.faa +15 -0
  275. data/test/data/blast/b0002.faa.m0 +128 -0
  276. data/test/data/blast/b0002.faa.m7 +65 -0
  277. data/test/data/blast/b0002.faa.m8 +1 -0
  278. data/test/data/blast/blastp-multi.m7 +188 -0
  279. data/test/data/command/echoarg2.bat +1 -0
  280. data/test/data/embl/AB090716.embl +65 -0
  281. data/test/data/embl/AB090716.embl.rel89 +63 -0
  282. data/test/data/fasta/example1.txt +75 -0
  283. data/test/data/fasta/example2.txt +21 -0
  284. data/test/data/genscan/sample.report +63 -0
  285. data/test/data/iprscan/merged.raw +32 -0
  286. data/test/data/iprscan/merged.txt +74 -0
  287. data/test/data/paml/codeml/control_file.txt +30 -0
  288. data/test/data/paml/codeml/output.txt +78 -0
  289. data/test/data/paml/codeml/rates +217 -0
  290. data/test/data/prosite/prosite.dat +2233 -0
  291. data/test/data/refseq/nm_126355.entret +64 -0
  292. data/test/data/rpsblast/misc.rpsblast +193 -0
  293. data/test/data/soft/GDS100_partial.soft +92 -0
  294. data/test/data/soft/GSE3457_family_partial.soft +874 -0
  295. data/test/data/uniprot/p53_human.uniprot +1456 -0
  296. data/test/functional/bio/appl/test_pts1.rb +115 -0
  297. data/test/functional/bio/io/test_ensembl.rb +229 -0
  298. data/test/functional/bio/io/test_soapwsdl.rb +52 -0
  299. data/test/functional/bio/sequence/test_output_embl.rb +51 -0
  300. data/test/functional/bio/test_command.rb +301 -0
  301. data/test/runner.rb +14 -0
  302. data/test/unit/bio/appl/bl2seq/test_report.rb +134 -0
  303. data/test/unit/bio/appl/blast/test_ncbioptions.rb +112 -0
  304. data/test/unit/bio/appl/blast/test_report.rb +1135 -0
  305. data/test/unit/bio/appl/blast/test_rpsblast.rb +398 -0
  306. data/test/unit/bio/appl/genscan/test_report.rb +182 -0
  307. data/test/unit/bio/appl/hmmer/test_report.rb +342 -0
  308. data/test/unit/bio/appl/iprscan/test_report.rb +338 -0
  309. data/test/unit/bio/appl/mafft/test_report.rb +63 -0
  310. data/test/unit/bio/appl/paml/codeml/test_rates.rb +45 -0
  311. data/test/unit/bio/appl/paml/codeml/test_report.rb +45 -0
  312. data/test/unit/bio/appl/paml/test_codeml.rb +174 -0
  313. data/test/unit/bio/appl/sosui/test_report.rb +81 -0
  314. data/test/unit/bio/appl/targetp/test_report.rb +146 -0
  315. data/test/unit/bio/appl/test_blast.rb +277 -0
  316. data/test/unit/bio/appl/test_fasta.rb +130 -0
  317. data/test/unit/bio/appl/test_psort.rb +57 -0
  318. data/test/unit/bio/appl/test_pts1.rb +77 -0
  319. data/test/unit/bio/appl/tmhmm/test_report.rb +126 -0
  320. data/test/unit/bio/data/test_aa.rb +90 -0
  321. data/test/unit/bio/data/test_codontable.rb +107 -0
  322. data/test/unit/bio/data/test_na.rb +80 -0
  323. data/test/unit/bio/db/embl/test_common.rb +117 -0
  324. data/test/unit/bio/db/embl/test_embl.rb +214 -0
  325. data/test/unit/bio/db/embl/test_embl_rel89.rb +219 -0
  326. data/test/unit/bio/db/embl/test_embl_to_bioseq.rb +203 -0
  327. data/test/unit/bio/db/embl/test_sptr.rb +1812 -0
  328. data/test/unit/bio/db/embl/test_uniprot.rb +31 -0
  329. data/test/unit/bio/db/kegg/test_genes.rb +45 -0
  330. data/test/unit/bio/db/pdb/test_pdb.rb +152 -0
  331. data/test/unit/bio/db/test_aaindex.rb +197 -0
  332. data/test/unit/bio/db/test_fasta.rb +250 -0
  333. data/test/unit/bio/db/test_gff.rb +1190 -0
  334. data/test/unit/bio/db/test_lasergene.rb +95 -0
  335. data/test/unit/bio/db/test_medline.rb +127 -0
  336. data/test/unit/bio/db/test_newick.rb +293 -0
  337. data/test/unit/bio/db/test_nexus.rb +364 -0
  338. data/test/unit/bio/db/test_prosite.rb +1437 -0
  339. data/test/unit/bio/db/test_rebase.rb +101 -0
  340. data/test/unit/bio/db/test_soft.rb +138 -0
  341. data/test/unit/bio/db/test_url.rb +36 -0
  342. data/test/unit/bio/io/flatfile/test_autodetection.rb +375 -0
  343. data/test/unit/bio/io/flatfile/test_buffer.rb +251 -0
  344. data/test/unit/bio/io/flatfile/test_splitter.rb +369 -0
  345. data/test/unit/bio/io/test_ddbjxml.rb +80 -0
  346. data/test/unit/bio/io/test_ensembl.rb +109 -0
  347. data/test/unit/bio/io/test_fastacmd.rb +42 -0
  348. data/test/unit/bio/io/test_flatfile.rb +505 -0
  349. data/test/unit/bio/io/test_soapwsdl.rb +32 -0
  350. data/test/unit/bio/sequence/test_aa.rb +115 -0
  351. data/test/unit/bio/sequence/test_common.rb +373 -0
  352. data/test/unit/bio/sequence/test_compat.rb +69 -0
  353. data/test/unit/bio/sequence/test_dblink.rb +58 -0
  354. data/test/unit/bio/sequence/test_na.rb +330 -0
  355. data/test/unit/bio/shell/plugin/test_seq.rb +185 -0
  356. data/test/unit/bio/test_alignment.rb +1025 -0
  357. data/test/unit/bio/test_command.rb +349 -0
  358. data/test/unit/bio/test_db.rb +96 -0
  359. data/test/unit/bio/test_feature.rb +144 -0
  360. data/test/unit/bio/test_location.rb +599 -0
  361. data/test/unit/bio/test_map.rb +230 -0
  362. data/test/unit/bio/test_pathway.rb +499 -0
  363. data/test/unit/bio/test_reference.rb +252 -0
  364. data/test/unit/bio/test_sequence.rb +329 -0
  365. data/test/unit/bio/test_shell.rb +18 -0
  366. data/test/unit/bio/test_tree.rb +593 -0
  367. data/test/unit/bio/util/restriction_enzyme/analysis/test_calculated_cuts.rb +299 -0
  368. data/test/unit/bio/util/restriction_enzyme/analysis/test_cut_ranges.rb +103 -0
  369. data/test/unit/bio/util/restriction_enzyme/analysis/test_sequence_range.rb +240 -0
  370. data/test/unit/bio/util/restriction_enzyme/double_stranded/test_aligned_strands.rb +101 -0
  371. data/test/unit/bio/util/restriction_enzyme/double_stranded/test_cut_location_pair.rb +75 -0
  372. data/test/unit/bio/util/restriction_enzyme/double_stranded/test_cut_location_pair_in_enzyme_notation.rb +73 -0
  373. data/test/unit/bio/util/restriction_enzyme/double_stranded/test_cut_locations.rb +53 -0
  374. data/test/unit/bio/util/restriction_enzyme/double_stranded/test_cut_locations_in_enzyme_notation.rb +104 -0
  375. data/test/unit/bio/util/restriction_enzyme/single_strand/test_cut_locations_in_enzyme_notation.rb +83 -0
  376. data/test/unit/bio/util/restriction_enzyme/test_analysis.rb +246 -0
  377. data/test/unit/bio/util/restriction_enzyme/test_cut_symbol.rb +44 -0
  378. data/test/unit/bio/util/restriction_enzyme/test_double_stranded.rb +115 -0
  379. data/test/unit/bio/util/restriction_enzyme/test_single_strand.rb +147 -0
  380. data/test/unit/bio/util/restriction_enzyme/test_single_strand_complement.rb +147 -0
  381. data/test/unit/bio/util/restriction_enzyme/test_string_formatting.rb +60 -0
  382. data/test/unit/bio/util/test_color_scheme.rb +33 -0
  383. data/test/unit/bio/util/test_contingency_table.rb +94 -0
  384. data/test/unit/bio/util/test_restriction_enzyme.rb +42 -0
  385. data/test/unit/bio/util/test_sirna.rb +245 -0
  386. metadata +543 -0
@@ -0,0 +1,297 @@
1
+ #
2
+ # = bio/io/flatfile/splitter.rb - input data splitter for FlatFile
3
+ #
4
+ # Copyright (C) 2001-2008 Naohisa Goto <ng@bioruby.org>
5
+ #
6
+ # License:: The Ruby License
7
+ #
8
+ # $Id:$
9
+ #
10
+ #
11
+ # See documents for Bio::FlatFile::Splitter and Bio::FlatFile.
12
+ #
13
+
14
+ require 'bio/io/flatfile'
15
+
16
+ module Bio
17
+
18
+ class FlatFile
19
+
20
+ # The Bio::FlatFile::Splitter is a namespace for flatfile splitters.
21
+ # Each splitter is a class to get entries from a buffered input stream.
22
+ #
23
+ # It is internally called in Bio::FlatFile.
24
+ # Normally, users do not need to use it directly.
25
+ module Splitter
26
+
27
+ # This is a template of splitter.
28
+ class Template
29
+ # Creates a new splitter.
30
+ def initialize(klass, bstream)
31
+ @dbclass = klass
32
+ @stream = bstream
33
+ @entry_pos_flag = nil
34
+ end
35
+
36
+ # skips leader of the entry.
37
+ def skip_leader
38
+ raise NotImplementedError
39
+ end
40
+
41
+ # rewind the stream
42
+ def rewind
43
+ @stream.rewind
44
+ end
45
+
46
+ # Gets entry as a string. (String)
47
+ def get_entry
48
+ raise NotImplementedError
49
+ end
50
+
51
+ # Gets entry as a data class's object
52
+ def get_parsed_entry
53
+ ent = get_entry
54
+ if ent then
55
+ self.parsed_entry = dbclass.new(ent)
56
+ else
57
+ self.parsed_entry = ent
58
+ end
59
+ parsed_entry
60
+ end
61
+
62
+ # the last entry string read from the stream (String)
63
+ attr_reader :entry
64
+
65
+ # The last parsed entry read from the stream (entry data class).
66
+ # Note that it is valid only after get_parsed_entry is called,
67
+ # and the get_entry may not affect the parsed_entry attribute.
68
+ attr_reader :parsed_entry
69
+
70
+ # a flag to write down entry start and end positions
71
+ attr_accessor :entry_pos_flag
72
+
73
+ # start position of the entry
74
+ attr_reader :entry_start_pos
75
+
76
+ # (end position of the entry) + 1
77
+ attr_reader :entry_ended_pos
78
+
79
+ #--
80
+ #private
81
+ #
82
+ ## to prevent warning message "warning: private attribute?",
83
+ ## private attributes are explicitly declared.
84
+ #++
85
+
86
+ # entry data class
87
+ attr_reader :dbclass
88
+ private :dbclass
89
+
90
+ # input stream
91
+ attr_reader :stream
92
+ private :stream
93
+
94
+ # the last entry string read from the stream
95
+ attr_writer :entry
96
+ private :entry=
97
+
98
+ # the last entry as a parsed data object
99
+ attr_writer :parsed_entry
100
+ private :parsed_entry=
101
+
102
+ # start position of the entry
103
+ attr_writer :entry_start_pos
104
+ private :entry_start_pos=
105
+
106
+ # (end position of the entry) + 1
107
+ attr_writer :entry_ended_pos
108
+ private :entry_ended_pos=
109
+
110
+ # Does stream.pos if entry_pos_flag is not nil.
111
+ # Otherwise, returns nil.
112
+ def stream_pos
113
+ entry_pos_flag ? stream.pos : nil
114
+ end
115
+ private :stream_pos
116
+ end #class Template
117
+
118
+ # Default splitter.
119
+ # It sees following constants in the given class.
120
+ # DELIMITER:: (String) delimiter indicates the end of a entry.
121
+ # FLATFILE_HEADER:: (String) start of a entry, located on head of a line.
122
+ # DELIMITER_OVERRUN:: (Integer) excess read size included in DELIMITER.
123
+ #
124
+ class Default < Template
125
+ # Creates a new splitter.
126
+ # klass:: database class
127
+ # bstream:: input stream. It must be a BufferedInputStream object.
128
+ def initialize(klass, bstream)
129
+ super(klass, bstream)
130
+
131
+ @delimiter = klass::DELIMITER rescue nil
132
+ @header = klass::FLATFILE_HEADER rescue nil
133
+ # for specific classes' benefit
134
+ unless header
135
+ if (defined?(Bio::GenBank) and klass == Bio::GenBank) or
136
+ (defined?(Bio::GenPept) and klass == Bio::GenPept)
137
+ @header = 'LOCUS '
138
+ end
139
+ end
140
+ @delimiter_overrun = klass::DELIMITER_OVERRUN rescue nil
141
+ end
142
+
143
+ # (String) delimiter indicates the end of a entry.
144
+ attr_accessor :delimiter
145
+
146
+ # (String) start of a entry, located on head of a line.
147
+ attr_accessor :header
148
+
149
+ # (Integer) excess read data size included in delimiter.
150
+ attr_accessor :delimiter_overrun
151
+
152
+ # Skips leader of the entry.
153
+ #
154
+ # If @header is not nil, it reads till the contents of @header
155
+ # comes at the head of a line.
156
+ # If correct FLATFILE_HEADER is found, returns true.
157
+ # Otherwise, returns nil.
158
+ def skip_leader
159
+ if @header then
160
+ data = ''
161
+ while s = stream.gets(@header)
162
+ data << s
163
+ if data.split(/[\r\n]+/)[-1] == @header then
164
+ stream.ungets(@header)
165
+ return true
166
+ end
167
+ end
168
+ # @header was not found. For safety,
169
+ # pushes back data with removing white spaces in the head.
170
+ data.sub(/\A\s+/, '')
171
+ stream.ungets(data)
172
+ return nil
173
+ else
174
+ stream.skip_spaces
175
+ return nil
176
+ end
177
+ end
178
+
179
+ # gets a entry
180
+ def get_entry
181
+ p0 = stream_pos()
182
+ e = stream.gets(@delimiter)
183
+ if e and @delimiter_overrun then
184
+ if e[-@delimiter.size, @delimiter.size ] == @delimiter then
185
+ overrun = e[-@delimiter_overrun, @delimiter_overrun]
186
+ e[-@delimiter_overrun, @delimiter_overrun] = ''
187
+ stream.ungets(overrun)
188
+ end
189
+ end
190
+ p1 = stream_pos()
191
+ self.entry_start_pos = p0
192
+ self.entry = e
193
+ self.entry_ended_pos = p1
194
+ return entry
195
+ end
196
+ end #class Defalult
197
+
198
+
199
+ # A splitter for line oriented text data.
200
+ #
201
+ # The given class's object must have following methods.
202
+ # Klass#add_header_line(line)
203
+ # Klass#add_line(line)
204
+ # where 'line' is a string. They normally returns self.
205
+ # If the line is not suitable to add to the current entry,
206
+ # nil or false should be returned.
207
+ # Then, the line is treated as (for add_header_line) the entry data
208
+ # or (for add_line) the next entry's data.
209
+ #
210
+ class LineOriented < Template
211
+ # Creates a new splitter.
212
+ # klass:: database class
213
+ # bstream:: input stream. It must be a BufferedInputStream object.
214
+ def initialize(klass, bstream)
215
+ super(klass, bstream)
216
+ self.flag_to_fetch_header = true
217
+ end
218
+
219
+ # do nothing
220
+ def skip_leader
221
+ nil
222
+ end
223
+
224
+ # get an entry and return the entry as a string
225
+ def get_entry
226
+ if e = get_parsed_entry then
227
+ entry
228
+ else
229
+ e
230
+ end
231
+ end
232
+
233
+ # get an entry and return the entry as a data class object
234
+ def get_parsed_entry
235
+ p0 = stream_pos()
236
+ ent = @dbclass.new()
237
+
238
+ lines = []
239
+ line_overrun = nil
240
+
241
+ if flag_to_fetch_header then
242
+ while line = stream.gets("\n")
243
+ unless ent.add_header_line(line) then
244
+ line_overrun = line
245
+ break
246
+ end
247
+ lines.push line
248
+ end
249
+ stream.ungets(line_overrun) if line_overrun
250
+ line_overrun = nil
251
+ self.flag_to_fetch_header = false
252
+ end
253
+
254
+ while line = stream.gets("\n")
255
+ unless ent.add_line(line) then
256
+ line_overrun = line
257
+ break
258
+ end
259
+ lines.push line
260
+ end
261
+ stream.ungets(line_overrun) if line_overrun
262
+ p1 = stream_pos()
263
+
264
+ return nil if lines.empty?
265
+
266
+ self.entry_start_pos = p0
267
+ self.entry = lines.join('')
268
+ self.parsed_entry = ent
269
+ self.entry_ended_pos = p1
270
+
271
+ return ent
272
+ end
273
+
274
+ # rewinds the stream
275
+ def rewind
276
+ ret = super
277
+ self.flag_to_fetch_header = true
278
+ ret
279
+ end
280
+
281
+ #--
282
+ #private methods / attributes
283
+ #++
284
+
285
+ # flag to fetch header
286
+ attr_accessor :flag_to_fetch_header
287
+ private :flag_to_fetch_header
288
+ private :flag_to_fetch_header=
289
+
290
+ end #class LineOriented
291
+
292
+ end #module Splitter
293
+
294
+ end #class FlatFile
295
+ end #module Bio
296
+
297
+
@@ -0,0 +1,73 @@
1
+ #
2
+ # = bio/io/higet.rb - SOAP interface for HGC HiGet
3
+ #
4
+ # Copyright:: Copyright (C) 2005 Toshiaki Katayama <k@bioruby.org>
5
+ #
6
+ # $Id: higet.rb,v 1.3 2006/09/19 05:46:22 k Exp $
7
+ #
8
+
9
+ require 'bio/io/soapwsdl'
10
+
11
+ module Bio
12
+ class HGC
13
+
14
+ # == Description
15
+ #
16
+ # Interface for the HiGet service provided by Human Genome Center (HGC), Japan.
17
+ # HiGet performs full-text search against various biological databases.
18
+ #
19
+ # == References
20
+ #
21
+ # * http://higet.hgc.jp/
22
+ #
23
+ class HiGet < Bio::SOAPWSDL
24
+
25
+ SERVER_URI = "http://higet.hgc.jp/soap/higet.wsdl"
26
+
27
+ def initialize(wsdl = nil)
28
+ super(wsdl || SERVER_URI)
29
+ end
30
+
31
+ def higet_in_fasta(db, entries)
32
+ self.higet(db, entries, "-d fasta")
33
+ end
34
+
35
+ def higet_in_xml(db, entries)
36
+ self.higet(db, entries, "-d xml")
37
+ end
38
+
39
+ end
40
+
41
+ end # HGC
42
+ end # Bio
43
+
44
+
45
+ if __FILE__ == $0
46
+
47
+ begin
48
+ require 'pp'
49
+ alias p pp
50
+ rescue LoadError
51
+ end
52
+
53
+ puts ">>> Bio::HGC::HiGet"
54
+ serv = Bio::HGC::HiGet.new
55
+ serv.log = STDERR
56
+
57
+ puts "### HiFind"
58
+ puts serv.hifind("genbank", "human kinase", "-l 10")
59
+
60
+ puts "### HiGet"
61
+ puts serv.higet("genbank", "S40289", "")
62
+
63
+ puts "### HiGet (FASTA)"
64
+ puts serv.higet("genbank", "S40289", "-d fasta")
65
+
66
+ puts "### HiGet higet_in_fasta"
67
+ puts serv.higet_in_fasta("genbank", "S40289")
68
+
69
+ puts "### HiGet higet_in_xml"
70
+ puts serv.higet_in_xml("genbank", "S40289")
71
+
72
+ end
73
+
@@ -0,0 +1,442 @@
1
+ #
2
+ # = bio/io/hinv.rb - H-invDB web service (REST) client module
3
+ #
4
+ # Copyright:: Copyright (C) 2008 Toshiaki Katayama <k@bioruby.org>
5
+ # License:: The Ruby License
6
+ #
7
+ # $Id:$
8
+ #
9
+
10
+ require 'bio/command'
11
+ require 'rexml/document'
12
+
13
+ module Bio
14
+
15
+ # = Bio::Hinv
16
+ #
17
+ # Accessing the H-invDB web services.
18
+ #
19
+ # * http://www.h-invitational.jp/
20
+ # * http://www.jbirc.aist.go.jp/hinv/hws/doc/index.html
21
+ #
22
+ class Hinv
23
+
24
+ BASE_URI = "http://www.jbirc.aist.go.jp/hinv/hws/"
25
+
26
+ module Common
27
+ def query(options = nil)
28
+ response = Bio::Command.post_form(@url, options)
29
+ @result = response.body
30
+ @xml = REXML::Document.new(@result)
31
+ end
32
+ end
33
+
34
+
35
+ # Bio::Hinv.acc2hit("BC053657") # => "HIT000053961"
36
+ def self.acc2hit(acc)
37
+ serv = Acc2hit.new
38
+ serv.query("acc" => acc)
39
+ serv.result
40
+ end
41
+
42
+ # Bio::Hinv.hit2acc("HIT000022181") # => "AK097327"
43
+ def self.hit2acc(hit)
44
+ serv = Hit2acc.new
45
+ serv.query("hit" => hit)
46
+ serv.result
47
+ end
48
+
49
+ # Bio::Hinv.hit_cnt # => 187156
50
+ def self.hit_cnt
51
+ serv = HitCnt.new
52
+ serv.query
53
+ serv.result
54
+ end
55
+
56
+ # Bio::Hinv.hit_definition("HIT000000001") # => "Rho guanine ..."
57
+ def self.hit_definition(hit)
58
+ serv = HitDefinition.new
59
+ serv.query("hit" => hit)
60
+ serv.result
61
+ end
62
+
63
+ # Bio::Hinv.hit_pubmedid("HIT000053961") # => [7624364, 11279095, ... ]
64
+ def self.hit_pubmedid(hit)
65
+ serv = HitPubmedId.new
66
+ serv.query("hit" => hit)
67
+ serv.result
68
+ end
69
+
70
+ # Bio::Hinv.hit_xml("HIT000000001") # => "<?xml version="1.0" ..."
71
+ def self.hit_xml(hit)
72
+ serv = Bio::Hinv::HitXML.new
73
+ serv.query("hit" => hit)
74
+ puts serv.result
75
+ end
76
+
77
+ # Bio::Hinv.hix2hit("HIX0000004") # => ["HIT000012846", ... ]
78
+ def self.hix2hit(hix)
79
+ serv = Bio::Hinv::Hix2hit.new
80
+ serv.query("hix" => hix)
81
+ serv.result
82
+ end
83
+
84
+ # Bio::Hinv.hix_cnt # => 36073
85
+ def self.hix_cnt
86
+ serv = HixCnt.new
87
+ serv.query
88
+ serv.result
89
+ end
90
+
91
+ # Bio::Hinv.hix_represent("HIX0000001") # => "HIT000022181"
92
+ def self.hix_represent(hix)
93
+ serv = HixRepresent.new
94
+ serv.query("hix" => hix)
95
+ serv.result
96
+ end
97
+
98
+ # Bio::Hinv.id_search("HIT00002218*") # => ["HIT000022181", ... ]
99
+ def self.id_search(query)
100
+ serv = IdSearch.new
101
+ serv.query("query" => query)
102
+ serv.result
103
+ end
104
+
105
+ # Bio::Hinv.keyword_search("HIT00002218*") # => ["HIT000022181", ... ]
106
+ def self.keyword_search(query)
107
+ serv = KeywordSearch.new
108
+ serv.query("query" => query)
109
+ serv.result
110
+ end
111
+
112
+
113
+ # serv = Bio::Hinv::Acc2hit.new
114
+ # serv.query("acc" => "BC053657")
115
+ # puts serv.result
116
+ class Acc2hit
117
+ include Common
118
+
119
+ def initialize
120
+ @url = BASE_URI + "acc2hit.php"
121
+ end
122
+
123
+ # <?xml version='1.0' encoding='UTF-8'?>
124
+ # <H-Inv>
125
+ # <H-INVITATIONAL-ID>HIT000053961</H-INVITATIONAL-ID>
126
+ # </H-Inv>
127
+ def result
128
+ @xml.elements['//H-INVITATIONAL-ID'].text
129
+ end
130
+ end
131
+
132
+ # serv = Bio::Hinv::Hit2acc.new
133
+ # serv.query("hit" => "HIT000022181")
134
+ # puts serv.result
135
+ class Hit2acc
136
+ include Common
137
+
138
+ def initialize
139
+ @url = BASE_URI + "hit2acc.php"
140
+ end
141
+
142
+ # <?xml version='1.0' encoding='UTF-8'?>
143
+ # <H-Inv>
144
+ # <ACCESSION-NO>AK097327</ACCESSION-NO>
145
+ # </H-Inv>
146
+ def result
147
+ @xml.elements['//ACCESSION-NO'].text
148
+ end
149
+ end
150
+
151
+ # serv = Bio::Hinv::HitCnt.new
152
+ # serv.query
153
+ # puts serv.result
154
+ class HitCnt
155
+ include Common
156
+
157
+ def initialize
158
+ @url = BASE_URI + "hit_cnt.php"
159
+ end
160
+
161
+ # <?xml version='1.0' encoding='UTF-8'?>
162
+ # <H-Inv>
163
+ # <TRANSCRIPT_CNT>187156</TRANSCRIPT_CNT>
164
+ # </H-Inv>
165
+ def result
166
+ @xml.elements['//TRANSCRIPT_CNT'].text.to_i
167
+ end
168
+ end
169
+
170
+ # serv = Bio::Hinv::HitDefinition.new
171
+ # serv.query("hit" => "HIT000000001")
172
+ # puts serv.result
173
+ # puts serv.data_source_definition
174
+ # puts serv.cdna_rep_h_invitational
175
+ # puts serv.cdna_splicing_isoform_curation
176
+ # puts serv.data_source_db_reference_protein_motif_id
177
+ # puts serv.data_source_identity
178
+ # puts serv.data_source_coverage
179
+ # puts serv.data_source_homologous_species
180
+ # puts serv.data_source_similarity_category
181
+ class HitDefinition
182
+ include Common
183
+
184
+ def initialize
185
+ @url = BASE_URI + "hit_definition.php"
186
+ end
187
+
188
+ # <?xml version='1.0' encoding='UTF-8'?>
189
+ # <H-Inv>
190
+ # <HIT_FUNCTION>
191
+ # <H-INVITATIONAL-ID>HIT000000001</H-INVITATIONAL-ID>
192
+ # <DATA-SOURCE_DEFINITION>Rho guanine nucleotide exchange factor 10.</DATA-SOURCE_DEFINITION>
193
+ # <CDNA_REP-H-INVITATIONAL>Representative transcript</CDNA_REP-H-INVITATIONAL>
194
+ # <CDNA_SPLICING-ISOFORM_CURATION></CDNA_SPLICING-ISOFORM_CURATION>
195
+ # <DATA-SOURCE_DB-REFERENCE_PROTEIN-MOTIF-ID>NP_055444</DATA-SOURCE_DB-REFERENCE_PROTEIN-MOTIF-ID>
196
+ # <DATA-SOURCE_IDENTITY>100.0</DATA-SOURCE_IDENTITY>
197
+ # <DATA-SOURCE_COVERAGE>100.0</DATA-SOURCE_COVERAGE>
198
+ # <DATA-SOURCE_HOMOLOGOUS_SPECIES>Homo sapiens</DATA-SOURCE_HOMOLOGOUS_SPECIES>
199
+ # <DATA-SOURCE_SIMILARITY-CATEGORY>Identical to known human protein(Category I).</DATA-SOURCE_SIMILARITY-CATEGORY>
200
+ # </HIT_FUNCTION>
201
+ # </H-Inv>
202
+ def result
203
+ @xml.elements['//DATA-SOURCE_DEFINITION'].text
204
+ end
205
+ alias :data_source_definition :result
206
+
207
+ def cdna_rep_h_invitational
208
+ @xml.elements['//CDNA_REP-H-INVITATIONAL'].text
209
+ end
210
+ def cdna_splicing_isoform_curation
211
+ @xml.elements['//CDNA_SPLICING-ISOFORM_CURATION'].text
212
+ end
213
+ def data_source_db_reference_protein_motif_id
214
+ @xml.elements['//DATA-SOURCE_DB-REFERENCE_PROTEIN-MOTIF-ID'].text
215
+ end
216
+ def data_source_identity
217
+ @xml.elements['//DATA-SOURCE_IDENTITY'].text.to_f
218
+ end
219
+ def data_source_coverage
220
+ @xml.elements['//DATA-SOURCE_COVERAGE'].text.to_f
221
+ end
222
+ def data_source_homologous_species
223
+ @xml.elements['//DATA-SOURCE_HOMOLOGOUS_SPECIES'].text
224
+ end
225
+ def data_source_similarity_category
226
+ @xml.elements['//DATA-SOURCE_SIMILARITY-CATEGORY'].text
227
+ end
228
+ end
229
+
230
+ # serv = Bio::Hinv::HitPubmedId.new
231
+ # serv.query("hit" => "HIT000053961")
232
+ # puts serv.result
233
+ class HitPubmedId
234
+ include Common
235
+
236
+ def initialize
237
+ @url = BASE_URI + "hit_pubmedid.php"
238
+ end
239
+
240
+ # <?xml version='1.0' encoding='UTF-8'?>
241
+ # <H-Inv>
242
+ # <CDNA_DB-REFERENCE_PUBMED>7624364</CDNA_DB-REFERENCE_PUBMED>
243
+ # <CDNA_DB-REFERENCE_PUBMED>11279095</CDNA_DB-REFERENCE_PUBMED>
244
+ # <CDNA_DB-REFERENCE_PUBMED>15489334</CDNA_DB-REFERENCE_PUBMED>
245
+ # </H-Inv>
246
+ def result
247
+ list = []
248
+ @xml.elements.each('//CDNA_DB-REFERENCE_PUBMED') do |e|
249
+ list << e.text.to_i
250
+ end
251
+ return list
252
+ end
253
+ end
254
+
255
+ # serv = Bio::Hinv::HitXML.new
256
+ # serv.query("hit" => "HIT000000001")
257
+ # puts serv.result
258
+ class HitXML
259
+ include Common
260
+
261
+ def initialize
262
+ @url = BASE_URI + "hit_xml.php"
263
+ end
264
+
265
+ # <?xml version="1.0" standalone="yes" ?>
266
+ # <H-Inv>
267
+ # <cDNAXML>
268
+ # <CLUSTER-ID>HIX0021591</CLUSTER-ID>
269
+ # <CLUSTER-ID-VERSION>HIX0021591.11</CLUSTER-ID-VERSION>
270
+ # <H-INVITATIONAL-ID>HIT000000001</H-INVITATIONAL-ID>
271
+ # :
272
+ # </PROBE-MAPPING>
273
+ # </EXPRESSION>
274
+ # </cDNAXML>
275
+ # </H-Inv>
276
+ def result
277
+ @result
278
+ end
279
+ end
280
+
281
+ # serv = Bio::Hinv::Hix2hit.new
282
+ # serv.query("hix" => "HIX0000004")
283
+ # puts serv.result
284
+ class Hix2hit
285
+ include Common
286
+
287
+ def initialize
288
+ @url = BASE_URI + "hix2hit.php"
289
+ end
290
+
291
+ # <?xml version='1.0' encoding='UTF-8'?>
292
+ # <H-Inv>
293
+ # <H-INVITATIONAL-ID>HIT000012846</H-INVITATIONAL-ID>
294
+ # <H-INVITATIONAL-ID>HIT000022124</H-INVITATIONAL-ID>
295
+ # <H-INVITATIONAL-ID>HIT000007722</H-INVITATIONAL-ID>
296
+ # :
297
+ # <H-INVITATIONAL-ID>HIT000262478</H-INVITATIONAL-ID>
298
+ # </H-Inv>
299
+ def result
300
+ list = []
301
+ @xml.elements.each('//H-INVITATIONAL-ID') do |e|
302
+ list << e.text
303
+ end
304
+ return list
305
+ end
306
+ end
307
+
308
+ # serv = Bio::Hinv::HixCnt.new
309
+ # serv.query
310
+ # puts serv.result
311
+ class HixCnt
312
+ include Common
313
+
314
+ def initialize
315
+ @url = BASE_URI + "hix_cnt.php"
316
+ end
317
+
318
+ # <?xml version='1.0' encoding='UTF-8'?>
319
+ # <H-Inv>
320
+ # <LOCUS_CNT>36073</LOCUS_CNT>
321
+ # </H-Inv>
322
+ def result
323
+ @xml.elements['//LOCUS_CNT'].text.to_i
324
+ end
325
+ end
326
+
327
+ # serv = Bio::Hinv::HixRepresent.new
328
+ # serv.query("hix" => "HIX0000001")
329
+ # puts serv.result
330
+ # puts serv.rep_h_invitational_id
331
+ # puts serv.rep_accession_no
332
+ class HixRepresent
333
+ include Common
334
+
335
+ def initialize
336
+ @url = BASE_URI + "hix_represent.php"
337
+ end
338
+
339
+ # <?xml version='1.0' encoding='UTF-8'?>
340
+ # <H-Inv>
341
+ # <LOCUS>
342
+ # <CLUSTER-ID>HIX0000001</CLUSTER-ID>
343
+ # <REP-H-INVITATIONAL-ID>HIT000022181</REP-H-INVITATIONAL-ID>
344
+ # <REP-ACCESSION-NO>AK097327</REP-ACCESSION-NO>
345
+ # </LOCUS>
346
+ # </H-Inv>
347
+ def result
348
+ @xml.elements['//REP-H-INVITATIONAL-ID'].text
349
+ end
350
+ alias :rep_h_invitational_id :result
351
+
352
+ def rep_accession_no
353
+ @xml.elements['//REP-ACCESSION-NO'].text
354
+ end
355
+ end
356
+
357
+ # example at "http://www.jbirc.aist.go.jp/hinv/hws/doc/index_jp.html"
358
+ # is for hit_xml.php (not for hix_xml.php)
359
+ class HixXML
360
+ end
361
+
362
+ # serv = Bio::Hinv::KeywordSearch.new
363
+ # serv.query("query" => "HIT00002218*", "start" => 1, "end" => 100)
364
+ # puts serv.result
365
+ # puts serv.size
366
+ # puts serv.start
367
+ # puts serv.end
368
+ class KeywordSearch
369
+ include Common
370
+
371
+ def initialize
372
+ @url = BASE_URI + "keyword_search.php"
373
+ end
374
+
375
+ def query(hash = {})
376
+ default = {
377
+ "start" => 1,
378
+ "end" => 100
379
+ }
380
+ options = default.update(hash)
381
+ super(options)
382
+ end
383
+
384
+ # <?xml version='1.0' encoding='UTF-8'?>
385
+ # <HINVDB_SEARCH>
386
+ # <QUERY>HIT00002218*</QUERY>
387
+ # <SIZE>8</SIZE>
388
+ # <START>1</START>
389
+ # <END>8</END>
390
+ # <H-INVITATIONAL-ID>HIT000022180</H-INVITATIONAL-ID>
391
+ # <H-INVITATIONAL-ID>HIT000022181</H-INVITATIONAL-ID>
392
+ # <H-INVITATIONAL-ID>HIT000022183</H-INVITATIONAL-ID>
393
+ # <H-INVITATIONAL-ID>HIT000022184</H-INVITATIONAL-ID>
394
+ # <H-INVITATIONAL-ID>HIT000022185</H-INVITATIONAL-ID>
395
+ # <H-INVITATIONAL-ID>HIT000022186</H-INVITATIONAL-ID>
396
+ # <H-INVITATIONAL-ID>HIT000022188</H-INVITATIONAL-ID>
397
+ # <H-INVITATIONAL-ID>HIT000022189</H-INVITATIONAL-ID>
398
+ # </HINVDB_SEARCH>
399
+ def result
400
+ list = []
401
+ @xml.elements.each('//H-INVITATIONAL-ID') do |e|
402
+ list << e.text
403
+ end
404
+ return list
405
+ end
406
+
407
+ def size
408
+ @xml.elements['//SIZE'].text.to_i
409
+ end
410
+ def start
411
+ @xml.elements['//START'].text.to_i
412
+ end
413
+ def end
414
+ @xml.elements['//END'].text.to_i
415
+ end
416
+ end
417
+
418
+ # serv = Bio::Hinv::IdSearch.new
419
+ # serv.query("query" => "HIT00002218*", "id_type" => "H-INVITATIONAL-ID", "start" => 1, "end" => 100)
420
+ # puts serv.result
421
+ # puts serv.size
422
+ # puts serv.start
423
+ # puts serv.end
424
+ class IdSearch < KeywordSearch
425
+ def initialize
426
+ @url = BASE_URI + "id_search.php"
427
+ end
428
+
429
+ def query(hash = {})
430
+ default = {
431
+ "id_type" => "H-INVITATIONAL-ID",
432
+ "start" => 1,
433
+ "end" => 100
434
+ }
435
+ options = default.update(hash)
436
+ super(options)
437
+ end
438
+ end
439
+
440
+ end
441
+ end
442
+