lederhosen 0.2.12 → 0.2.13
Sign up to get free protection for your applications and to get access to all the features.
- data/.rspec +1 -1
- data/Gemfile +0 -1
- data/Rakefile +0 -8
- data/lederhosen.gemspec +1 -4
- data/lib/lederhosen/tasks/uc_filter.rb +12 -2
- data/lib/lederhosen/version.rb +1 -1
- data/spec/cli_spec.rb +7 -1
- data/spec/helpers_spec.rb +1 -1
- metadata +3 -17
data/.rspec
CHANGED
@@ -1 +1 @@
|
|
1
|
-
-c
|
1
|
+
-c --fail-fast
|
data/Gemfile
CHANGED
data/Rakefile
CHANGED
@@ -34,14 +34,6 @@ Rake::TestTask.new(:test) do |test|
|
|
34
34
|
test.verbose = true
|
35
35
|
end
|
36
36
|
|
37
|
-
require 'rcov/rcovtask'
|
38
|
-
Rcov::RcovTask.new do |test|
|
39
|
-
test.libs << 'test'
|
40
|
-
test.pattern = 'test/**/test_*.rb'
|
41
|
-
test.verbose = true
|
42
|
-
test.rcov_opts << '--exclude "gems/*"'
|
43
|
-
end
|
44
|
-
|
45
37
|
task :default => :test
|
46
38
|
|
47
39
|
require 'rdoc/task'
|
data/lederhosen.gemspec
CHANGED
@@ -5,7 +5,7 @@
|
|
5
5
|
|
6
6
|
Gem::Specification.new do |s|
|
7
7
|
s.name = "lederhosen"
|
8
|
-
s.version = "0.2.
|
8
|
+
s.version = "0.2.13"
|
9
9
|
|
10
10
|
s.required_rubygems_version = Gem::Requirement.new(">= 0") if s.respond_to? :required_rubygems_version=
|
11
11
|
s.authors = ["Austin G. Davis-Richardson"]
|
@@ -71,7 +71,6 @@ Gem::Specification.new do |s|
|
|
71
71
|
s.add_development_dependency(%q<rdoc>, ["~> 3.12"])
|
72
72
|
s.add_development_dependency(%q<bundler>, [">= 0"])
|
73
73
|
s.add_development_dependency(%q<jeweler>, [">= 0"])
|
74
|
-
s.add_development_dependency(%q<rcov>, [">= 0"])
|
75
74
|
else
|
76
75
|
s.add_dependency(%q<dna>, [">= 0"])
|
77
76
|
s.add_dependency(%q<progressbar>, [">= 0"])
|
@@ -79,7 +78,6 @@ Gem::Specification.new do |s|
|
|
79
78
|
s.add_dependency(%q<rdoc>, ["~> 3.12"])
|
80
79
|
s.add_dependency(%q<bundler>, [">= 0"])
|
81
80
|
s.add_dependency(%q<jeweler>, [">= 0"])
|
82
|
-
s.add_dependency(%q<rcov>, [">= 0"])
|
83
81
|
end
|
84
82
|
else
|
85
83
|
s.add_dependency(%q<dna>, [">= 0"])
|
@@ -88,7 +86,6 @@ Gem::Specification.new do |s|
|
|
88
86
|
s.add_dependency(%q<rdoc>, ["~> 3.12"])
|
89
87
|
s.add_dependency(%q<bundler>, [">= 0"])
|
90
88
|
s.add_dependency(%q<jeweler>, [">= 0"])
|
91
|
-
s.add_dependency(%q<rcov>, [">= 0"])
|
92
89
|
end
|
93
90
|
end
|
94
91
|
|
@@ -23,10 +23,12 @@ module Lederhosen
|
|
23
23
|
ohai "filtering #{input} to #{output}, reads = #{reads} & samples = #{samples}"
|
24
24
|
|
25
25
|
# load UC file
|
26
|
+
ohai "loading uc file"
|
26
27
|
clstr_info = Helpers.load_uc_file input
|
27
28
|
clstr_counts = clstr_info[:clstr_counts] # clstr_counts[:clstr][sample.to_i] = reads
|
28
29
|
|
29
30
|
# filter
|
31
|
+
ohai "filtering"
|
30
32
|
survivors = clstr_counts.reject do |a, b|
|
31
33
|
b.reject{ |i, j| j < reads }.length < samples
|
32
34
|
end
|
@@ -34,13 +36,18 @@ module Lederhosen
|
|
34
36
|
surviving_clusters = survivors.keys.to_set
|
35
37
|
|
36
38
|
# print filtered uc file
|
39
|
+
ohai "saving filtered table"
|
37
40
|
out = File.open(output, 'w')
|
38
|
-
|
41
|
+
|
39
42
|
lines = `wc -l #{input}`.split.first.to_i
|
43
|
+
|
40
44
|
pbar = ProgressBar.new 'saving', lines
|
45
|
+
kept, total = 1, 0
|
46
|
+
|
41
47
|
File.open(input) do |handle|
|
42
48
|
handle.each do |line|
|
43
49
|
pbar.inc
|
50
|
+
|
44
51
|
if line =~ /^#/
|
45
52
|
out.print line
|
46
53
|
next
|
@@ -48,15 +55,18 @@ module Lederhosen
|
|
48
55
|
|
49
56
|
total += 1
|
50
57
|
|
58
|
+
# check if cluster is in surviving clusters
|
51
59
|
if surviving_clusters.include? line.split[1].to_i
|
52
60
|
out.print line
|
53
61
|
kept += 1
|
54
62
|
end
|
63
|
+
|
55
64
|
end
|
56
65
|
end
|
66
|
+
|
57
67
|
pbar.finish
|
58
68
|
out.close
|
59
|
-
|
69
|
+
|
60
70
|
ohai "clusters: #{surviving_clusters.length}/#{clstr_counts.keys.length} = #{100*surviving_clusters.length/clstr_counts.keys.length.to_f}%"
|
61
71
|
ohai "reads: #{kept}/#{total} = #{100*kept/total.to_f}%"
|
62
72
|
end
|
data/lib/lederhosen/version.rb
CHANGED
data/spec/cli_spec.rb
CHANGED
@@ -54,7 +54,13 @@ describe Lederhosen::CLI do
|
|
54
54
|
it 'should identify clusters given a taxcollector database'
|
55
55
|
|
56
56
|
it 'should add names to otu abundance matrix given blat output' do
|
57
|
-
|
57
|
+
levels = %w{kingdom domain phylum class order genus speces}
|
58
|
+
level =
|
59
|
+
begin
|
60
|
+
levels.sample # 1.9 and up
|
61
|
+
rescue
|
62
|
+
levels.choice # 1.8 and lower
|
63
|
+
end
|
58
64
|
`./bin/lederhosen add_names --table=spec/data/otus.csv --blat=spec/data/blat.txt --level=#{level} --output=#{$test_dir}/named_otus.csv`
|
59
65
|
$?.success?.should be_true
|
60
66
|
end
|
data/spec/helpers_spec.rb
CHANGED
@@ -21,6 +21,6 @@ describe Lederhosen::Helpers do
|
|
21
21
|
reads = groups.values.first
|
22
22
|
Lederhosen::Helpers.trim_pairs reads[0], reads[1], "#{$test_dir}/munchen_trim_test.fasta"
|
23
23
|
# this test will break if trim parameters change
|
24
|
-
File.
|
24
|
+
File.readlines("#{$test_dir}/munchen_trim_test.fasta").grep(/^>/).length.should be_even
|
25
25
|
end
|
26
26
|
end
|
metadata
CHANGED
@@ -1,13 +1,13 @@
|
|
1
1
|
--- !ruby/object:Gem::Specification
|
2
2
|
name: lederhosen
|
3
3
|
version: !ruby/object:Gem::Version
|
4
|
-
hash:
|
4
|
+
hash: 13
|
5
5
|
prerelease:
|
6
6
|
segments:
|
7
7
|
- 0
|
8
8
|
- 2
|
9
|
-
-
|
10
|
-
version: 0.2.
|
9
|
+
- 13
|
10
|
+
version: 0.2.13
|
11
11
|
platform: ruby
|
12
12
|
authors:
|
13
13
|
- Austin G. Davis-Richardson
|
@@ -102,20 +102,6 @@ dependencies:
|
|
102
102
|
version_requirements: *id006
|
103
103
|
name: jeweler
|
104
104
|
prerelease: false
|
105
|
-
- !ruby/object:Gem::Dependency
|
106
|
-
type: :development
|
107
|
-
requirement: &id007 !ruby/object:Gem::Requirement
|
108
|
-
none: false
|
109
|
-
requirements:
|
110
|
-
- - ">="
|
111
|
-
- !ruby/object:Gem::Version
|
112
|
-
hash: 3
|
113
|
-
segments:
|
114
|
-
- 0
|
115
|
-
version: "0"
|
116
|
-
version_requirements: *id007
|
117
|
-
name: rcov
|
118
|
-
prerelease: false
|
119
105
|
description: Various tools for OTU clustering
|
120
106
|
email: harekrishna@gmail.com
|
121
107
|
executables:
|