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 CHANGED
@@ -1 +1 @@
1
- -c
1
+ -c --fail-fast
data/Gemfile CHANGED
@@ -12,5 +12,4 @@ group :development do
12
12
  gem 'rdoc', '~> 3.12'
13
13
  gem 'bundler'
14
14
  gem 'jeweler'
15
- gem 'rcov'
16
15
  end
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.12"
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
- kept, total = 1, 0
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
- ohai "Survivors"
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
@@ -2,7 +2,7 @@ module Lederhosen
2
2
  module Version
3
3
  MAJOR = 0
4
4
  MINOR = 2
5
- PATCH = 12
5
+ PATCH = 13
6
6
 
7
7
  STRING = [MAJOR, MINOR, PATCH].join('.')
8
8
  end
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
- level = %w{kingdom domain phylum class order genus speces}.choice
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.read("#{$test_dir}/munchen_trim_test.fasta").grep(/^>/).length.should be_even
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: 15
4
+ hash: 13
5
5
  prerelease:
6
6
  segments:
7
7
  - 0
8
8
  - 2
9
- - 12
10
- version: 0.2.12
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: