lederhosen 0.1.6 → 0.1.7
Sign up to get free protection for your applications and to get access to all the features.
- 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
|