lederhosen 0.2.12 → 0.2.13
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- 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:
|