protk 1.4.4.beta1 → 1.4.4.beta2
Sign up to get free protection for your applications and to get access to all the features.
- checksums.yaml +4 -4
- data/bin/mzid_to_protxml.rb +2 -1
- data/lib/protk/mzidentml_doc.rb +2 -1
- data/lib/protk/protein_group.rb +2 -1
- metadata +1 -1
checksums.yaml
CHANGED
@@ -1,7 +1,7 @@
|
|
1
1
|
---
|
2
2
|
SHA1:
|
3
|
-
metadata.gz:
|
4
|
-
data.tar.gz:
|
3
|
+
metadata.gz: f278e7fe8e3a0955907a13952ced4fa0d772c204
|
4
|
+
data.tar.gz: 513fb214683486b9d3596a9b9f0956107d611170
|
5
5
|
SHA512:
|
6
|
-
metadata.gz:
|
7
|
-
data.tar.gz:
|
6
|
+
metadata.gz: e59fcb0724cbf42b8f63e65ca6d88fd91c7d45d3981964358c82e040f9235e2ce6768251d088966044827cadf5c9ce46f189345d3c27293d0b10b44daf018c10
|
7
|
+
data.tar.gz: 36329dca4cf416fc2b9bd6f0d395895ee2fd32fdba2f6ebd75c5d26ed21b497dc890cfc6bcf36b1f08bc7f3954bae76adea64450a431d6ac365a2930191997cc
|
data/bin/mzid_to_protxml.rb
CHANGED
@@ -71,7 +71,8 @@ protein_groups.each do |group_node|
|
|
71
71
|
if group_prob >= tool.minprob.to_f
|
72
72
|
$stdout.write "\n" if tool.debug
|
73
73
|
$protk.log "Writing group with probability #{group_prob}" , :info
|
74
|
-
group = ProteinGroup.from_mzid(group_node,mzid_doc,tool.minprob.to_f)
|
74
|
+
group = ProteinGroup.from_mzid(group_node,mzid_doc,tool.minprob.to_f)
|
75
|
+
# require 'byebug';byebug
|
75
76
|
prot_xml_writer.append_protein_group(group.as_protxml)
|
76
77
|
n_written+=1
|
77
78
|
end
|
data/lib/protk/mzidentml_doc.rb
CHANGED
@@ -127,8 +127,9 @@ class MzIdentMLDoc < Object
|
|
127
127
|
|
128
128
|
# Memoized because it gets called for every protein in a group
|
129
129
|
def get_proteins_for_group(group_node)
|
130
|
+
# puts group_node.attributes['group_number']
|
130
131
|
@proteins_for_group_cache ||= Hash.new do |h,key|
|
131
|
-
h[key] = self.find(
|
132
|
+
h[key] = self.find(key,"ProteinDetectionHypothesis")
|
132
133
|
end
|
133
134
|
@proteins_for_group_cache[group_node]
|
134
135
|
end
|
data/lib/protk/protein_group.rb
CHANGED
@@ -68,7 +68,8 @@ class ProteinGroup
|
|
68
68
|
def as_protxml()
|
69
69
|
node = XML::Node.new('protein_group')
|
70
70
|
node["group_number"] = self.group_number.to_s
|
71
|
-
|
71
|
+
# require 'byebug';byebug
|
72
|
+
node["group_probability"] = (self.group_probability * 0.01).to_s
|
72
73
|
self.proteins.each { |prot| node << prot.as_protxml }
|
73
74
|
node
|
74
75
|
end
|