1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
|
#!/usr/bin/env python
# File created on 09 Feb 2010
from __future__ import division
__author__ = "Justin Kuczynski"
__copyright__ = "Copyright 2011, The QIIME Project"
__credits__ = ["Justin Kuczynski"]
__license__ = "GPL"
__version__ = "1.4.0"
__maintainer__ = "Justin Kuczynski"
__email__ = "justinak@gmail.com"
__status__ = "Release"
from qiime.util import parse_command_line_parameters, create_dir
from qiime.util import make_option
import os.path
from qiime.rarefaction import RarefactionMaker
script_info={}
script_info['brief_description']="""Perform multiple subsamplings/rarefactions on an otu table"""
script_info['script_description']="""To perform bootstrap, jackknife, and rarefaction analyses, the otu table must be subsampled (rarefied). This script rarefies, or subsamples, OTU tables. This does not provide curves of diversity by number of sequences in a sample. Rather it creates a series of subsampled OTU tables by random sampling (without replacement) of the input OTU table. Samples that have fewer sequences then the requested rarefaction depth for a given output otu table are omitted from those ouput otu tables. The pseudo-random number generator used for rarefaction by subsampling is NumPy's default - an implementation of the Mersenne twister PRNG."""
script_info['script_usage']=[]
script_info['script_usage'].append(("""Examples:""","""An example of this script, where the user sets the minimum ("-m") and maximum ("-x") number of sequences per sample to 100 and 1200, respectively, while using steps ("-s") of 100, performing 2 iterations at each sampling depth ("-n"), and outputting the results to the directory "rarefaction_tables/" is shown by the following command:""","""multiple_rarefactions.py -i otu_table.txt -m 100 -x 1200 -s 100 -n 2 -o rarefaction_tables/"""))
script_info['script_usage'].append(("""""","""As a result, this command produces subsamples of the input otu_table.txt at 100 seqs per sample (twice), 200 seqs per sample (twice) ... 1200 seqs per sample (twice), which produces 24 rarefied otu talbes in the "rarefaction_tables" directory.""",""""""))
script_info['script_usage'].append(("""""","""Any sample containing fewer sequences in the input file than the requested number of sequences per sample is removed from the output rarefied otu table. To include samples with fewer than the requested number, you must manually add those samples to the resulting otu tables""",""""""))
script_info['output_description']="""The result of multiple_rarefactions.py consists of a number of files, which depend on the minimum/maximum number of sequences per samples, steps and iterations. The files have the same otu table format as the input otu_table.txt, and are named in the following way: rarefaction_100_0.txt, where "100" corresponds to the sequences per sample and "0" the iteration."""
script_info['required_options']=[
make_option('-i', '--input_path',
help='Input OTU table filepath.',
type='existing_filepath'),
make_option('-o', '--output_path',
help="Output directory.",
type='new_dirpath'),
make_option('-m', '--min', type=int,
help='Minimum number of seqs/sample for rarefaction.'),
make_option('-x', '--max', type=int,
help='Maximum number of seqs/sample (inclusive) for rarefaction. '),
make_option('-s', '--step', type=int,
help='Size of each steps between the min/max of' +\
' seqs/sample (e.g. min, min+step... for level <= max).')
]
script_info['optional_options']=[
### This option is screwed up: dest should equal the long form parameter name,
### but I'm not sure if we can do anything about it since '-' is not allowed
### in variable names... Hmmm... Changing the long-form parameter name
### would cause older parameter files not to work.
make_option('-n', '--num-reps', dest='num_reps', default=10, type=int,
help='The number of iterations at each step. [default: %default]'),
make_option('--lineages_included', default=False,
action="store_true",
help='Retain taxonomic (lineage) information for each OTU. Note:' +\
' this will only work if lineage information is in the input OTU' +\
' table. [default: %default]'),
make_option('-k', '--keep_empty_otus', default=False, action='store_true',
help='Retain OTUs of all zeros, which are usually omitted from' +\
' the output OTU tables. [default: %default]'),
]
script_info['option_label']={'input_path':'OTU table filepath',
'output_path': 'Output directory',
'min': 'Min # of seqs/sample',
'max': 'Max # of seqs/sample',
'step': 'Step size',
'num-reps':'# of iterations',
'lineages_included': 'Include lineages',
'keep_empty_otus':'Retain empty OTUs'}
script_info['version'] = __version__
def main():
option_parser, opts, args = parse_command_line_parameters(**script_info)
if opts.step <= 0:
option_parser.error("nonpositive step not allowed (%s was supplied)" % \
(opts.step,))
create_dir(opts.output_path, fail_on_exist=False)
maker = RarefactionMaker(opts.input_path, opts.min, opts.max,
opts.step, opts.num_reps)
maker.rarefy_to_files(opts.output_path, False,
include_lineages=opts.lineages_included,
empty_otus_removed=(not opts.keep_empty_otus))
if __name__ == "__main__":
main()
|
