-
Notifications
You must be signed in to change notification settings - Fork 109
/
demo_tmhmm_report.rb
68 lines (60 loc) · 1.36 KB
/
demo_tmhmm_report.rb
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
#
# = sample/demo_tmhmm_report.rb - demonstration of Bio::TMHMM::Report
#
# Copyright:: Copyright (C) 2003
# Mitsuteru C. Nakao <n@bioruby.org>
# License:: The Ruby License
#
#
# == Description
#
# Demonstration of Bio::TMHMM::Report, TMHMM output parser.
#
# == Usage
#
# Specify files containing SOSUI reports.
#
# $ ruby demo_tmhmm_report.rb files...
#
# Example usage using test data:
#
# $ ruby -Ilib sample/demo_tmhmm_report.rb test/data/TMHMM/sample.report
#
# == References
#
# * http://www.cbs.dtu.dk/services/TMHMM/
#
# == Development information
#
# The code was moved from lib/bio/appl/tmhmm/report.rb.
#
require 'bio'
#if __FILE__ == $0
begin
require 'pp'
alias p pp
rescue LoadError
end
Bio::TMHMM.reports(ARGF.read) do |ent|
puts '==>'
puts ent.to_s
pp ent
p [:entry_id, ent.entry_id]
p [:query_len, ent.query_len]
p [:predicted_tmhs, ent.predicted_tmhs]
p [:tmhs_size, ent.tmhs.size]
p [:exp_aas_in_tmhs, ent.exp_aas_in_tmhs]
p [:exp_first_60aa, ent.exp_first_60aa]
p [:total_prob_of_N_in, ent.total_prob_of_N_in]
ent.tmhs.each do |t|
p t
p [:entry_id, t.entry_id]
p [:version, t.version]
p [:status, t.status]
p [:range, t.range]
p [:pos, t.pos]
end
p [:helix, ent.helix]
p ent.tmhs.map {|t| t if t.status == 'TMhelix' }.compact
end
#end