lederhosen 0.1.6 → 0.1.7
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/Gemfile +1 -1
- data/bin/lederhosen +1 -1
- data/examples/hierarchical_clustering.sh +51 -0
- data/examples/pipeline.sh +1 -1
- data/lederhosen.gemspec +1 -1
- data/lib/lederhosen/buffer.rb +2 -2
- data/lib/lederhosen/helpers.rb +10 -5
- data/lib/lederhosen/tasks/add_names.rb +80 -80
- data/lib/lederhosen/tasks/cluster.rb +1 -1
- data/lib/lederhosen/tasks/name.rb +2 -2
- data/lib/lederhosen/tasks/otu_filter.rb +43 -43
- data/lib/lederhosen/tasks/otu_table.rb +13 -13
- data/lib/lederhosen/tasks/rep_reads.rb +2 -2
- data/lib/lederhosen/tasks/squish.rb +42 -42
- data/lib/lederhosen.rb +2 -1
- data/lib/version.rb +1 -1
- data/readme.md +7 -3
- data/spec/data/ILT_L_9_B_001_1.txt.gz +0 -0
- data/spec/data/ILT_L_9_B_001_3.txt.gz +0 -0
- data/spec/data/ILT_L_9_B_002_1.txt.gz +0 -0
- data/spec/data/ILT_L_9_B_002_3.txt.gz +0 -0
- data/spec/helpers_spec.rb +7 -7
- data/spec/misc_spec.rb +1 -1
- data/spec/pipeline_spec.rb +18 -22
- data/spec/spec_helper.rb +5 -1
- metadata +13 -16
- data/spec/data/ILT_L_9_B_001_1.txt +0 -400
- data/spec/data/ILT_L_9_B_001_3.txt +0 -400
- data/spec/data/ILT_L_9_B_002_1.txt +0 -400
- data/spec/data/ILT_L_9_B_002_3.txt +0 -400
- data/spec/data/blast_out.txt +0 -10
- data/spec/data/blat.txt +0 -86
data/spec/pipeline_spec.rb
CHANGED
@@ -1,8 +1,4 @@
|
|
1
|
-
|
2
|
-
|
3
|
-
`mkdir -p #{test_dir}`
|
4
|
-
|
5
|
-
$stderr.puts "test dir: #{test_dir}"; sleep 1
|
1
|
+
require 'spec_helper'
|
6
2
|
|
7
3
|
describe 'the pipeline' do
|
8
4
|
|
@@ -12,44 +8,44 @@ describe 'the pipeline' do
|
|
12
8
|
end
|
13
9
|
|
14
10
|
it 'should trim reads' do
|
15
|
-
`./bin/lederhosen trim --reads-dir=spec/data/IL*.txt --out-dir=#{test_dir}/trimmed`
|
11
|
+
`./bin/lederhosen trim --reads-dir=spec/data/IL*.txt.gz --out-dir=#{$test_dir}/trimmed`
|
16
12
|
$?.success?.should be_true
|
17
13
|
end
|
18
14
|
|
19
15
|
it 'should join reads' do
|
20
|
-
`./bin/lederhosen join --trimmed=#{test_dir}/trimmed/*.fasta --output=#{test_dir}/joined.fasta`
|
16
|
+
`./bin/lederhosen join --trimmed=#{$test_dir}/trimmed/*.fasta --output=#{$test_dir}/joined.fasta`
|
21
17
|
$?.success?.should be_true
|
22
18
|
end
|
23
19
|
|
24
20
|
it 'should sort reads' do
|
25
|
-
`./bin/lederhosen sort --input=#{test_dir}/joined.fasta --output=#{test_dir}/sorted.fasta`
|
21
|
+
`./bin/lederhosen sort --input=#{$test_dir}/joined.fasta --output=#{$test_dir}/sorted.fasta`
|
26
22
|
$?.success?.should be_true
|
27
23
|
end
|
28
24
|
|
29
25
|
it 'should cluster reads' do
|
30
|
-
`./bin/lederhosen cluster --identity=0.80 --input=#{test_dir}/sorted.fasta --output=#{test_dir}/clusters.uc`
|
26
|
+
`./bin/lederhosen cluster --identity=0.80 --input=#{$test_dir}/sorted.fasta --output=#{$test_dir}/clusters.uc`
|
31
27
|
$?.success?.should be_true
|
32
28
|
end
|
33
29
|
|
34
30
|
it 'should build OTU abundance matrices' do
|
35
|
-
`./bin/lederhosen otu_table --clusters=#{test_dir}/clusters.uc --output=#{test_dir}/test_tables --joined=#{test_dir}/joined.fasta`
|
31
|
+
`./bin/lederhosen otu_table --clusters=#{$test_dir}/clusters.uc --output=#{$test_dir}/test_tables --joined=#{$test_dir}/joined.fasta`
|
36
32
|
$?.success?.should be_true
|
37
33
|
end
|
38
34
|
|
39
35
|
it 'should split joined.fasta into reads for each cluster' do
|
40
|
-
`./bin/lederhosen split --reads=#{test_dir}/joined.fasta --clusters=#{test_dir}/clusters.uc --out-dir=#{test_dir}/split --min-clst-size=1`
|
36
|
+
`./bin/lederhosen split --reads=#{$test_dir}/joined.fasta --clusters=#{$test_dir}/clusters.uc --out-dir=#{$test_dir}/split --min-clst-size=1`
|
41
37
|
end
|
42
38
|
|
43
|
-
|
44
|
-
|
45
|
-
it 'should add names to otu abundance matrix given blat output' do
|
46
|
-
level = %w{kingdom domain phylum class order genus speces}.choice
|
47
|
-
`./bin/lederhosen add_names --table=spec/data/otus.csv --blat=spec/data/blat.txt --level=#{level} --output=#{test_dir}/named_otus.csv`
|
48
|
-
$?.success?.should be_true
|
49
|
-
end
|
39
|
+
it 'should name clusters given a taxcollector database'
|
50
40
|
|
51
|
-
|
52
|
-
|
53
|
-
|
54
|
-
|
41
|
+
it 'should add names to otu abundance matrix given blat output' do
|
42
|
+
level = %w{kingdom domain phylum class order genus speces}.choice
|
43
|
+
`./bin/lederhosen add_names --table=spec/data/otus.csv --blat=spec/data/blat.txt --level=#{level} --output=#{$test_dir}/named_otus.csv`
|
44
|
+
$?.success?.should be_true
|
45
|
+
end
|
46
|
+
|
47
|
+
it 'should squish otu abundance matrix by same name' do
|
48
|
+
`./bin/lederhosen squish --csv-file=#{$test_dir}/named_otus.csv --output=#{$test_dir}/squished.csv`
|
49
|
+
$?.success?.should be_true
|
50
|
+
end
|
55
51
|
end
|
data/spec/spec_helper.rb
CHANGED
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: 21
|
5
5
|
prerelease:
|
6
6
|
segments:
|
7
7
|
- 0
|
8
8
|
- 1
|
9
|
-
-
|
10
|
-
version: 0.1.
|
9
|
+
- 7
|
10
|
+
version: 0.1.7
|
11
11
|
platform: ruby
|
12
12
|
authors:
|
13
13
|
- Austin G. Davis-Richardson
|
@@ -15,7 +15,7 @@ autorequire:
|
|
15
15
|
bindir: bin
|
16
16
|
cert_chain: []
|
17
17
|
|
18
|
-
date: 2012-07-
|
18
|
+
date: 2012-07-17 00:00:00 Z
|
19
19
|
dependencies:
|
20
20
|
- !ruby/object:Gem::Dependency
|
21
21
|
name: dna
|
@@ -129,6 +129,7 @@ files:
|
|
129
129
|
- .rvmrc
|
130
130
|
- Gemfile
|
131
131
|
- bin/lederhosen
|
132
|
+
- examples/hierarchical_clustering.sh
|
132
133
|
- examples/pipeline.sh
|
133
134
|
- lederhosen.gemspec
|
134
135
|
- lib/lederhosen.rb
|
@@ -150,12 +151,10 @@ files:
|
|
150
151
|
- lib/lederhosen/tasks/uc_filter.rb
|
151
152
|
- lib/version.rb
|
152
153
|
- readme.md
|
153
|
-
- spec/data/ILT_L_9_B_001_1.txt
|
154
|
-
- spec/data/ILT_L_9_B_001_3.txt
|
155
|
-
- spec/data/ILT_L_9_B_002_1.txt
|
156
|
-
- spec/data/ILT_L_9_B_002_3.txt
|
157
|
-
- spec/data/blast_out.txt
|
158
|
-
- spec/data/blat.txt
|
154
|
+
- spec/data/ILT_L_9_B_001_1.txt.gz
|
155
|
+
- spec/data/ILT_L_9_B_001_3.txt.gz
|
156
|
+
- spec/data/ILT_L_9_B_002_1.txt.gz
|
157
|
+
- spec/data/ILT_L_9_B_002_3.txt.gz
|
159
158
|
- spec/data/otus.csv
|
160
159
|
- spec/helpers_spec.rb
|
161
160
|
- spec/misc_spec.rb
|
@@ -195,12 +194,10 @@ signing_key:
|
|
195
194
|
specification_version: 3
|
196
195
|
summary: 16S rRNA clustering for paired-end Illumina
|
197
196
|
test_files:
|
198
|
-
- spec/data/ILT_L_9_B_001_1.txt
|
199
|
-
- spec/data/ILT_L_9_B_001_3.txt
|
200
|
-
- spec/data/ILT_L_9_B_002_1.txt
|
201
|
-
- spec/data/ILT_L_9_B_002_3.txt
|
202
|
-
- spec/data/blast_out.txt
|
203
|
-
- spec/data/blat.txt
|
197
|
+
- spec/data/ILT_L_9_B_001_1.txt.gz
|
198
|
+
- spec/data/ILT_L_9_B_001_3.txt.gz
|
199
|
+
- spec/data/ILT_L_9_B_002_1.txt.gz
|
200
|
+
- spec/data/ILT_L_9_B_002_3.txt.gz
|
204
201
|
- spec/data/otus.csv
|
205
202
|
- spec/helpers_spec.rb
|
206
203
|
- spec/misc_spec.rb
|