ms-sequest 0.1.1 → 0.1.2
Sign up to get free protection for your applications and to get access to all the features.
- data/History +4 -0
- data/Rakefile +4 -4
- data/VERSION +1 -1
- data/lib/ms/sequest/sqt.rb +1 -2
- data/lib/ms/sequest/srf.rb +10 -7
- data/lib/ms/sequest/srf/pepxml.rb +8 -8
- data/spec/ms/sequest/srf/sqt_spec.rb +1 -0
- metadata +8 -19
data/History
CHANGED
data/Rakefile
CHANGED
@@ -11,13 +11,13 @@ Jeweler::Tasks.new do |gem|
|
|
11
11
|
gem.email = "jtprince@gmail.com"
|
12
12
|
gem.authors = ["John T. Prince"]
|
13
13
|
gem.rubyforge_project = 'mspire'
|
14
|
-
gem.add_runtime_dependency "ms-ident", ">= 0.
|
15
|
-
gem.add_runtime_dependency "ms-core", ">= 0.0.
|
14
|
+
gem.add_runtime_dependency "ms-ident", ">= 0.1.1"
|
15
|
+
gem.add_runtime_dependency "ms-core", ">= 0.0.19"
|
16
16
|
#gem.add_runtime_dependency "ms-msrun", ">= 0.3.4"
|
17
|
-
gem.add_runtime_dependency "trollop", "~> 1.16"
|
17
|
+
gem.add_runtime_dependency "trollop", "~> 1.16.2"
|
18
18
|
gem.add_development_dependency "jeweler", "~> 1.5.2"
|
19
19
|
gem.add_development_dependency "ms-testdata", ">= 0.1.1"
|
20
|
-
gem.add_development_dependency "spec-more", ">= 0"
|
20
|
+
gem.add_development_dependency "spec-more", ">= 0.0.4"
|
21
21
|
end
|
22
22
|
Jeweler::RubygemsDotOrgTasks.new
|
23
23
|
|
data/VERSION
CHANGED
@@ -1 +1 @@
|
|
1
|
-
0.1.
|
1
|
+
0.1.2
|
data/lib/ms/sequest/sqt.rb
CHANGED
data/lib/ms/sequest/srf.rb
CHANGED
@@ -14,8 +14,7 @@ require 'ms/sequest/params'
|
|
14
14
|
module Ms ; end
|
15
15
|
module Ms::Sequest ; end
|
16
16
|
|
17
|
-
class Ms::Sequest::Srf
|
18
|
-
include Ms::Ident::Search
|
17
|
+
class Ms::Sequest::Srf < Ms::Ident::Search
|
19
18
|
|
20
19
|
class NoSequestParamsError < ArgumentError
|
21
20
|
end
|
@@ -102,7 +101,7 @@ class Ms::Sequest::Srf
|
|
102
101
|
# # will attempt to read peptide hit information (equivalent to .out
|
103
102
|
# # files), otherwise, just reads the dta information.
|
104
103
|
def initialize(filename=nil, opts={})
|
105
|
-
@
|
104
|
+
@peptide_hits = []
|
106
105
|
|
107
106
|
@dta_files = []
|
108
107
|
@out_files = []
|
@@ -270,7 +269,7 @@ class Ms::Sequest::Srf
|
|
270
269
|
outfile.charge = ind[2]
|
271
270
|
|
272
271
|
pep_hits = @out_files[i].hits
|
273
|
-
@
|
272
|
+
@peptide_hits.push( *pep_hits )
|
274
273
|
pep_hits.each do |pep_hit|
|
275
274
|
pep_hit[15] = @base_name
|
276
275
|
pep_hit[16] = ind[0]
|
@@ -286,7 +285,7 @@ class Ms::Sequest::Srf
|
|
286
285
|
filter_by_precursor_mass_tolerance! if params
|
287
286
|
|
288
287
|
if opts[:link_protein_hits]
|
289
|
-
(@
|
288
|
+
(@peptide_hits, @proteins) = merge!([self.peptide_hits]) do |_protein, _peptides|
|
290
289
|
Ms::Sequest::Srf::Out::Protein.new(_protein.reference, _peptides)
|
291
290
|
end
|
292
291
|
end
|
@@ -615,7 +614,7 @@ end
|
|
615
614
|
# deltacn should also be changed to reflect it.
|
616
615
|
# mh - the theoretical mass + h
|
617
616
|
# proteins are created as SRF prot objects with a reference and linked to their
|
618
|
-
#
|
617
|
+
# peptide_hits (from global hash by reference)
|
619
618
|
# ppm = 10^6 * ∆m_accuracy / mass_measured [ where ∆m_accuracy = mass_real – mass_measured ]
|
620
619
|
# This is calculated for the M+H mass!
|
621
620
|
# num_other_loci is the number of other loci that the peptide matches beyond
|
@@ -728,7 +727,7 @@ end
|
|
728
727
|
Ms::Sequest::Srf::Out::Protein = Struct.new( *%w(reference peptides).map(&:to_sym) )
|
729
728
|
|
730
729
|
class Ms::Sequest::Srf::Out::Protein
|
731
|
-
include Ms::Ident::Protein
|
730
|
+
#include Ms::Ident::Protein
|
732
731
|
|
733
732
|
## we shouldn't have to do this because this is inlcuded in SpecID::Protein, but
|
734
733
|
## under some circumstances it won't work without explicitly calling it.
|
@@ -741,6 +740,10 @@ class Ms::Sequest::Srf::Out::Protein
|
|
741
740
|
end
|
742
741
|
$VERBOSE = tmp
|
743
742
|
|
743
|
+
def first_entry
|
744
|
+
reference.split(/[\s\r]/)[0]
|
745
|
+
end
|
746
|
+
|
744
747
|
# "<Ms::Sequest::Srf::Out::Protein reference=\"#{@reference}\">"
|
745
748
|
|
746
749
|
undef_method :inspect
|
@@ -155,24 +155,24 @@ class Ms::Sequest::Srf
|
|
155
155
|
msms_pipeline_analysis.merge!(:summary_xml => summary_xml_filename, :pepxml_version => opt[:pepxml_version]) do |msms_run_summary|
|
156
156
|
# prep the sample enzyme and search_summary
|
157
157
|
msms_run_summary.merge!(
|
158
|
-
:base_name => File.join(opt[:mz_dir], srf.base_name_noext),
|
158
|
+
{ :base_name => File.join(opt[:mz_dir], srf.base_name_noext),
|
159
159
|
:ms_manufacturer => opt[:ms_manufacturer],
|
160
160
|
:ms_model => opt[:ms_model],
|
161
161
|
:ms_ionization => opt[:ms_ionization],
|
162
162
|
:ms_mass_analyzer => opt[:ms_mass_analyzer],
|
163
163
|
:ms_detector => opt[:ms_detector],
|
164
164
|
:raw_data => opt[:raw_data].first,
|
165
|
-
:raw_data_type => opt[:raw_data].first,
|
165
|
+
:raw_data_type => opt[:raw_data].first, }
|
166
166
|
) do |sample_enzyme, search_summary, spectrum_queries|
|
167
167
|
sample_enzyme.merge!(params.sample_enzyme_hash)
|
168
168
|
sample_enzyme.name = opt[:enzyme] if opt[:enzyme]
|
169
169
|
search_summary.merge!(
|
170
|
-
:base_name=> srf.resident_dir + '/' + srf.base_name_noext,
|
170
|
+
{ :base_name=> srf.resident_dir + '/' + srf.base_name_noext,
|
171
171
|
:search_engine => 'SEQUEST',
|
172
172
|
:precursor_mass_type => params.precursor_mass_type,
|
173
173
|
:fragment_mass_type => params.fragment_mass_type,
|
174
174
|
:out_data_type => opt[:out_data_type],
|
175
|
-
:out_data => opt[:out_data],
|
175
|
+
:out_data => opt[:out_data], }
|
176
176
|
) do |search_database, enzymatic_search_constraint, modifications_ar, parameters_hash|
|
177
177
|
search_database.merge!(:local_path => db_filename, :seq_type => opt[:db_seq_type], :database_name => opt[:db_name], :orig_database_url => opt[:db_orig_url], :database_release_date => opt[:db_release_date], :database_release_identifier => opt[:db_release_id])
|
178
178
|
|
@@ -184,9 +184,9 @@ class Ms::Sequest::Srf
|
|
184
184
|
end
|
185
185
|
|
186
186
|
enzymatic_search_constraint.merge!(
|
187
|
-
:enzyme => opt[:enzyme] ? opt[:enzyme] : params.enzyme,
|
187
|
+
{ :enzyme => opt[:enzyme] ? opt[:enzyme] : params.enzyme,
|
188
188
|
:max_num_internal_cleavages => params.max_num_internal_cleavages,
|
189
|
-
:min_number_termini => params.min_number_termini,
|
189
|
+
:min_number_termini => params.min_number_termini, }
|
190
190
|
)
|
191
191
|
modifications_ar.replace(modifications_obj.modifications)
|
192
192
|
parameters_hash.merge!(params.opts)
|
@@ -240,10 +240,10 @@ class Ms::Sequest::Srf
|
|
240
240
|
end
|
241
241
|
end
|
242
242
|
Ms::Ident::Pepxml::SpectrumQuery.new(
|
243
|
-
:spectrum => [srf.base_name_noext, *i_ar].join('.'), :start_scan => i_ar[0], :end_scan => i_ar[1],
|
243
|
+
{ :spectrum => [srf.base_name_noext, *i_ar].join('.'), :start_scan => i_ar[0], :end_scan => i_ar[1],
|
244
244
|
:precursor_neutral_mass => dta_file.mh - h_plus, :assumed_charge => i_ar[2],
|
245
245
|
:retention_time_sec => ret_time,
|
246
|
-
:search_results => [sr],
|
246
|
+
:search_results => [sr], }
|
247
247
|
)
|
248
248
|
end
|
249
249
|
spectrum_queries.replace(spec_queries)
|
metadata
CHANGED
@@ -2,7 +2,7 @@
|
|
2
2
|
name: ms-sequest
|
3
3
|
version: !ruby/object:Gem::Version
|
4
4
|
prerelease:
|
5
|
-
version: 0.1.
|
5
|
+
version: 0.1.2
|
6
6
|
platform: ruby
|
7
7
|
authors:
|
8
8
|
- John T. Prince
|
@@ -10,7 +10,7 @@ autorequire:
|
|
10
10
|
bindir: bin
|
11
11
|
cert_chain: []
|
12
12
|
|
13
|
-
date: 2011-
|
13
|
+
date: 2011-07-05 00:00:00 -06:00
|
14
14
|
default_executable:
|
15
15
|
dependencies:
|
16
16
|
- !ruby/object:Gem::Dependency
|
@@ -21,7 +21,7 @@ dependencies:
|
|
21
21
|
requirements:
|
22
22
|
- - ">="
|
23
23
|
- !ruby/object:Gem::Version
|
24
|
-
version: 0.
|
24
|
+
version: 0.1.1
|
25
25
|
type: :runtime
|
26
26
|
version_requirements: *id001
|
27
27
|
- !ruby/object:Gem::Dependency
|
@@ -32,7 +32,7 @@ dependencies:
|
|
32
32
|
requirements:
|
33
33
|
- - ">="
|
34
34
|
- !ruby/object:Gem::Version
|
35
|
-
version: 0.0.
|
35
|
+
version: 0.0.19
|
36
36
|
type: :runtime
|
37
37
|
version_requirements: *id002
|
38
38
|
- !ruby/object:Gem::Dependency
|
@@ -43,7 +43,7 @@ dependencies:
|
|
43
43
|
requirements:
|
44
44
|
- - ~>
|
45
45
|
- !ruby/object:Gem::Version
|
46
|
-
version:
|
46
|
+
version: 1.16.2
|
47
47
|
type: :runtime
|
48
48
|
version_requirements: *id003
|
49
49
|
- !ruby/object:Gem::Dependency
|
@@ -76,7 +76,7 @@ dependencies:
|
|
76
76
|
requirements:
|
77
77
|
- - ">="
|
78
78
|
- !ruby/object:Gem::Version
|
79
|
-
version:
|
79
|
+
version: 0.0.4
|
80
80
|
type: :development
|
81
81
|
version_requirements: *id006
|
82
82
|
description: reads .SRF, .SQT and supports conversions
|
@@ -161,16 +161,5 @@ rubygems_version: 1.6.2
|
|
161
161
|
signing_key:
|
162
162
|
specification_version: 3
|
163
163
|
summary: An mspire library supporting SEQUEST, Bioworks, SQT, etc
|
164
|
-
test_files:
|
165
|
-
|
166
|
-
- spec/ms/sequest/params_spec.rb
|
167
|
-
- spec/ms/sequest/pepxml/modifications_spec.rb
|
168
|
-
- spec/ms/sequest/pepxml_spec.rb
|
169
|
-
- spec/ms/sequest/sqt_spec.rb
|
170
|
-
- spec/ms/sequest/sqt_spec_helper.rb
|
171
|
-
- spec/ms/sequest/srf/pepxml_spec.rb
|
172
|
-
- spec/ms/sequest/srf/search_spec.rb
|
173
|
-
- spec/ms/sequest/srf/sqt_spec.rb
|
174
|
-
- spec/ms/sequest/srf_spec.rb
|
175
|
-
- spec/ms/sequest/srf_spec_helper.rb
|
176
|
-
- spec/spec_helper.rb
|
164
|
+
test_files: []
|
165
|
+
|