本文整理汇总了Python中qiime.parse.parse_qiime_parameters函数的典型用法代码示例。如果您正苦于以下问题:Python parse_qiime_parameters函数的具体用法?Python parse_qiime_parameters怎么用?Python parse_qiime_parameters使用的例子?那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了parse_qiime_parameters函数的20个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的Python代码示例。
示例1: test_run_core_diversity_analyses_error_on_invalid_category
def test_run_core_diversity_analyses_error_on_invalid_category(self):
"""run_core_diversity_analyses: error raised on invalid categories
"""
# too few values in 'month' category
self.assertRaises(ValueError,
run_core_diversity_analyses,
self.test_data['biom'][0],
self.test_data['map'][0],
20,
output_dir=self.test_out,
params=parse_qiime_parameters({}),
qiime_config=self.qiime_config,
categories=['SampleType', 'month'],
tree_fp=self.test_data['tree'][0],
parallel=False,
status_update_callback=no_status_updates)
# invalid category name
self.assertRaises(ValueError,
run_core_diversity_analyses,
self.test_data['biom'][0],
self.test_data['map'][0],
20,
output_dir=self.test_out,
params=parse_qiime_parameters({}),
qiime_config=self.qiime_config,
categories=['not-a-real-category'],
tree_fp=self.test_data['tree'][0],
parallel=False,
status_update_callback=no_status_updates)
开发者ID:Springbudder,项目名称:qiime,代码行数:30,代码来源:test_core_diversity_analyses.py
示例2: run_command
def run_command(self, options, arguments):
verbose = options["verbose"]
input_fp = options["input_fp"]
output_dir = options["output_dir"]
verbose = options["verbose"]
print_only = options["print_only"]
parallel = options["parallel"]
# No longer checking that jobs_to_start > 2, but
# commenting as we may change our minds about this.
# if parallel: raise_error_on_parallel_unavailable()
if options["parameter_fp"]:
try:
parameter_f = open(options["parameter_fp"])
except IOError:
raise QiimeCommandError, "Can't open parameters file (%s). Does it exist? Do you have read access?" % options[
"parameter_fp"
]
params = parse_qiime_parameters(parameter_f)
else:
params = parse_qiime_parameters([])
# empty list returns empty defaultdict for now
params["parallel"]["jobs_to_start"] = self._validate_jobs_to_start(
options["jobs_to_start"], qiime_config["jobs_to_start"], parallel
)
try:
makedirs(output_dir)
except OSError:
if options["force"]:
pass
else:
# Since the analysis can take quite a while, I put this check
# in to help users avoid overwriting previous output.
print "Output directory already exists. Please choose " + "a different directory, or force overwrite with -f."
exit(1)
if print_only:
command_handler = print_commands
else:
command_handler = call_commands_serially
if verbose:
status_update_callback = print_to_stdout
else:
status_update_callback = no_status_updates
run_qiime_data_preparation(
input_fp,
output_dir,
command_handler=command_handler,
params=params,
qiime_config=qiime_config,
parallel=parallel,
status_update_callback=status_update_callback,
)
开发者ID:gregcaporaso,项目名称:cmd-abstraction,代码行数:60,代码来源:interfaces.py
示例3: main
def main():
option_parser, opts, args = parse_command_line_parameters(**script_info)
verbose = opts.verbose
input_fp = opts.input_fp
output_dir = opts.output_dir
verbose = opts.verbose
print_only = opts.print_only
parallel = opts.parallel
# No longer checking that jobs_to_start > 2, but
# commenting as we may change our minds about this.
# if parallel: raise_error_on_parallel_unavailable()
if opts.parameter_fp:
try:
parameter_f = open(opts.parameter_fp)
except IOError:
raise IOError, "Can't open parameters file (%s). Does it exist? Do you have read access?" % opts.parameter_fp
params = parse_qiime_parameters(parameter_f)
else:
params = parse_qiime_parameters([])
# empty list returns empty defaultdict for now
jobs_to_start = opts.jobs_to_start
default_jobs_to_start = qiime_config["jobs_to_start"]
validate_and_set_jobs_to_start(params, jobs_to_start, default_jobs_to_start, parallel, option_parser)
try:
makedirs(output_dir)
except OSError:
if opts.force:
pass
else:
# Since the analysis can take quite a while, I put this check
# in to help users avoid overwriting previous output.
print "Output directory already exists. Please choose " + "a different directory, or force overwrite with -f."
exit(1)
if print_only:
command_handler = print_commands
else:
command_handler = call_commands_serially
if verbose:
status_update_callback = print_to_stdout
else:
status_update_callback = no_status_updates
run_qiime_data_preparation(
input_fp,
output_dir,
command_handler=command_handler,
params=params,
qiime_config=qiime_config,
parallel=parallel,
status_update_callback=status_update_callback,
)
开发者ID:ranjit58,项目名称:qiime,代码行数:59,代码来源:pick_otus_through_otu_table.py
示例4: main
def main():
option_parser, opts, args = parse_command_line_parameters(**script_info)
verbose = opts.verbose
otu_table_fp = opts.otu_table_fp
output_dir = opts.output_dir
mapping_fp = opts.mapping_fp
verbose = opts.verbose
print_only = opts.print_only
mapping_cat = opts.mapping_category
sort=opts.sort
if mapping_cat != None and mapping_fp == None:
option_parser.error("If passing -c must also pass -m.")
if opts.parameter_fp:
try:
parameter_f = open(opts.parameter_fp, 'U')
except IOError:
raise IOError,\
"Can't open parameters file (%s). Does it exist? Do you have read access?" \
% opts.parameter_fp
params = parse_qiime_parameters(parameter_f)
parameter_f.close()
else:
params = parse_qiime_parameters([])
# empty list returns empty defaultdict for now
try:
makedirs(output_dir)
except OSError:
if opts.force:
pass
else:
# Since the analysis can take quite a while, I put this check
# in to help users avoid overwriting previous output.
option_parser.error("Output directory already exists. Please choose"
" a different directory, or force overwrite with -f.")
if print_only:
command_handler = print_commands
else:
command_handler = call_commands_serially
if verbose:
status_update_callback = print_to_stdout
else:
status_update_callback = no_status_updates
run_summarize_taxa_through_plots(
otu_table_fp=otu_table_fp,
mapping_fp=mapping_fp,
output_dir=output_dir,
mapping_cat=mapping_cat,
sort=sort,
command_handler=command_handler,
params=params,
qiime_config=qiime_config,
status_update_callback=status_update_callback)
开发者ID:gxenomics,项目名称:qiime,代码行数:59,代码来源:summarize_taxa_through_plots.py
示例5: main
def main():
option_parser, opts, args = \
parse_command_line_parameters(suppress_verbose=True, **script_info)
input_dir = opts.input_dir
parameter_fp = opts.parameter_fp
read1_indicator = opts.read1_indicator
read2_indicator = opts.read2_indicator
match_barcodes = opts.match_barcodes
barcode_indicator = opts.barcode_indicator
leading_text = opts.leading_text
trailing_text = opts.trailing_text
include_input_dir_path = opts.include_input_dir_path
output_dir = abspath(opts.output_dir)
remove_filepath_in_name = opts.remove_filepath_in_name
print_only = opts.print_only
if remove_filepath_in_name and not include_input_dir_path:
option_parser.error("If --remove_filepath_in_name is enabled, "
"--include_input_dir_path must also be enabled.")
if opts.parameter_fp:
with open(opts.parameter_fp, 'U') as parameter_f:
params_dict = parse_qiime_parameters(parameter_f)
params_str = get_params_str(params_dict['join_paired_ends'])
else:
params_dict = {}
params_str = ""
create_dir(output_dir)
all_files = []
extensions = ['.fastq.gz', '.fastq', '.fq.gz', '.fq']
for root, dir, fps in walk(input_dir):
for fp in fps:
for extension in extensions:
if fp.endswith(extension):
all_files += [abspath(join(root, fp))]
pairs, bc_pairs = get_pairs(all_files, read1_indicator, read2_indicator,
match_barcodes, barcode_indicator)
commands = create_commands_jpe(pairs, output_dir,
params_str, leading_text, trailing_text, include_input_dir_path,
remove_filepath_in_name, match_barcodes, bc_pairs)
qiime_config = load_qiime_config()
if print_only:
command_handler = print_commands
else:
command_handler = call_commands_serially
logger = WorkflowLogger(generate_log_fp(output_dir),
params=params_dict,
qiime_config=qiime_config)
# Call the command handler on the list of commands
command_handler(commands,
status_update_callback=no_status_updates,
logger=logger,
close_logger_on_success=True)
开发者ID:Springbudder,项目名称:qiime,代码行数:60,代码来源:multiple_join_paired_ends.py
示例6: test_run_core_diversity_analyses_parallel
def test_run_core_diversity_analyses_parallel(self):
"""run_core_diversity_analyses functions with categories in parallel
"""
run_core_diversity_analyses(
self.test_data['biom'][0],
self.test_data['map'][0],
20,
output_dir=self.test_out,
params=parse_qiime_parameters({}),
arare_num_steps=3,
qiime_config=self.qiime_config,
categories=['SampleType', 'days_since_epoch'],
tree_fp=self.test_data['tree'][0],
parallel=True,
status_update_callback=no_status_updates)
# Basic sanity test that output directories and files exist
fps = [
'%s/bdiv_even20' % self.test_out,
'%s/arare_max20' % self.test_out,
'%s/taxa_plots' % self.test_out,
'%s/bdiv_even20/unweighted_unifrac_dm.txt' % self.test_out,
'%s/bdiv_even20/weighted_unifrac_pc.txt' % self.test_out,
'%s/arare_max20/compare_chao1/days_since_epoch_stats.txt' % self.test_out,
'%s/arare_max20/compare_PD_whole_tree/SampleType_boxplots.pdf' % self.test_out,
'%s/index.html' % self.test_out,
'%s/table_mc%d.biom.gz' % (self.test_out, 20)
]
for fp in fps:
self.assertTrue(exists(fp))
开发者ID:Springbudder,项目名称:qiime,代码行数:30,代码来源:test_core_diversity_analyses.py
示例7: setUp
def setUp(self):
"""setup the test values"""
self.qiime_config = load_qiime_config()
self.dirs_to_remove = []
self.files_to_remove = []
#this is specific to the web-apps only
test_dir = abspath(dirname(__file__))
self.fna_original_fp = os.path.join(test_dir, 'support_files', \
'test.fna')
tmp_dir = self.qiime_config['temp_dir'] or '/tmp/'
if not exists(tmp_dir):
makedirs(tmp_dir)
# if test creates the temp dir, also remove it
self.dirs_to_remove.append(tmp_dir)
self.wf_out = get_tmp_filename(tmp_dir=tmp_dir,
prefix='qiime_wf_out',suffix='',result_constructor=str)
if not exists(self.wf_out):
makedirs(self.wf_out)
self.dirs_to_remove.append(self.wf_out)
#print self.wf_out
working_dir = self.qiime_config['working_dir'] or './'
jobs_dir = join(working_dir,'jobs')
if not exists(jobs_dir):
# only clean up the jobs dir if it doesn't already exist
self.dirs_to_remove.append(jobs_dir)
self.params = parse_qiime_parameters(qiime_parameters_f.split('\n'))
signal.signal(signal.SIGALRM, timeout)
# set the 'alarm' to go off in allowed_seconds seconds
signal.alarm(allowed_seconds_per_test)
开发者ID:ElDeveloper,项目名称:qiime_web_app,代码行数:35,代码来源:test_chain_pick_otus.py
示例8: setUp
def setUp(self):
""" """
self.test_data = get_test_data_fps()
self.files_to_remove = []
self.dirs_to_remove = []
# Create example output directory
tmp_dir = get_qiime_temp_dir()
self.test_out = mkdtemp(dir=tmp_dir,
prefix='core_qiime_analyses_test_',
suffix='')
self.dirs_to_remove.append(self.test_out)
self.qiime_config = load_qiime_config()
self.params = parse_qiime_parameters(params_f1)
# suppress stderr during tests (one of the systems calls in the
# workflow prints a warning, and we can't suppress that warning with
# warnings.filterwarnings) here because it comes from within the code
# executed through the system call. Found this trick here:
# http://stackoverflow.com/questions/9949633/suppressing-print-as-stdout-python
self.saved_stderr = sys.stderr
sys.stderr = StringIO()
initiate_timeout(180)
开发者ID:Kleptobismol,项目名称:qiime,代码行数:25,代码来源:test_downstream.py
示例9: test_run_core_diversity_analyses
def test_run_core_diversity_analyses(self):
"""run_core_diversity_analyses functions with categories
"""
run_core_diversity_analyses(
self.test_data["biom"][0],
self.test_data["map"][0],
20,
output_dir=self.test_out,
params=parse_qiime_parameters({}),
qiime_config=self.qiime_config,
categories=["SampleType", "days_since_epoch"],
tree_fp=self.test_data["tree"][0],
parallel=False,
status_update_callback=no_status_updates,
)
# Basic sanity test that output directories and files exist
fps = [
"%s/bdiv_even20" % self.test_out,
"%s/arare_max20" % self.test_out,
"%s/taxa_plots" % self.test_out,
"%s/bdiv_even20/unweighted_unifrac_dm.txt" % self.test_out,
"%s/bdiv_even20/weighted_unifrac_pc.txt" % self.test_out,
"%s/arare_max20/compare_chao1/days_since_epoch_stats.txt" % self.test_out,
"%s/arare_max20/compare_PD_whole_tree/SampleType_boxplots.pdf" % self.test_out,
"%s/index.html" % self.test_out,
"%s/table_mc%d.biom.gz" % (self.test_out, 20),
]
for fp in fps:
self.assertTrue(exists(fp))
开发者ID:kylebittinger,项目名称:qiime,代码行数:30,代码来源:test_core_diversity_analyses.py
示例10: test_run_core_diversity_analyses_no_categories
def test_run_core_diversity_analyses_no_categories(self):
"""run_core_diversity_analyses functions without categories
"""
# this takes a long time, so use a longer sigalrm
run_core_diversity_analyses(
self.test_data['biom'][0],
self.test_data['map'][0],
20,
output_dir=self.test_out,
params=parse_qiime_parameters({}),
qiime_config=self.qiime_config,
categories=None,
tree_fp=self.test_data['tree'][0],
parallel=False,
status_update_callback=no_status_updates)
# Basic sanity test that output directories and files exist
fps = [
'%s/bdiv_even20' % self.test_out,
'%s/arare_max20' % self.test_out,
'%s/taxa_plots' % self.test_out,
'%s/bdiv_even20/unweighted_unifrac_dm.txt' % self.test_out,
'%s/bdiv_even20/weighted_unifrac_pc.txt' % self.test_out,
'%s/index.html' % self.test_out,
'%s/table_mc%d.biom.gz' % (self.test_out, 20)
]
for fp in fps:
self.assertTrue(exists(fp))
# categorical output files don't exist
self.assertFalse(exists(
'%s/arare_max20/compare_chao1/days_since_epoch_stats.txt' % self.test_out))
self.assertFalse(exists(
'%s/arare_max20/compare_PD_whole_tree/SampleType_boxplots.pdf' % self.test_out))
开发者ID:Springbudder,项目名称:qiime,代码行数:34,代码来源:test_core_diversity_analyses.py
示例11: test_run_core_diversity_analyses_no_tree
def test_run_core_diversity_analyses_no_tree(self):
"""run_core_diversity_analyses functions without tree
"""
# this takes a long time, so use a longer sigalrm
run_core_diversity_analyses(
self.test_data['biom'][0],
self.test_data['map'][0],
20,
output_dir=self.test_out,
params=parse_qiime_parameters(
['beta_diversity:metrics bray_curtis',
'alpha_diversity:metrics observed_species,chao1']),
qiime_config=self.qiime_config,
categories=['SampleType'],
tree_fp=None,
parallel=False,
status_update_callback=no_status_updates)
# Basic sanity test that output directories and files exist
fps = [
'%s/bdiv_even20' % self.test_out,
'%s/arare_max20' % self.test_out,
'%s/taxa_plots' % self.test_out,
'%s/bdiv_even20/bray_curtis_dm.txt' % self.test_out,
'%s/arare_max20/compare_observed_species/SampleType_boxplots.pdf' % self.test_out,
'%s/index.html' % self.test_out,
'%s/table_mc%d.biom.gz' % (self.test_out, 20)
]
for fp in fps:
self.assertTrue(exists(fp))
# phylogenetic diversity output files do not exist
self.assertFalse(exists(
'%s/bdiv_even20/unweighted_unifrac_dm.txt' % self.test_out))
开发者ID:Springbudder,项目名称:qiime,代码行数:34,代码来源:test_core_diversity_analyses.py
示例12: setUp
def setUp(self):
""" """
self.qiime_config = load_qiime_config()
self.dirs_to_remove = []
self.files_to_remove = []
#this is specific to the web-apps only
test_dir = abspath(dirname(__file__))
sff_original_fp = os.path.join(test_dir, 'support_files', \
'Fasting_subset.sff')
self.sff_fp = os.path.join('/%s/' % environ['HOME'],
'Fasting_subset.sff')
self.files_to_remove.append(self.sff_fp)
copy(sff_original_fp, self.sff_fp)
self.illumina_fps = [os.path.join(test_dir, 'support_files', \
's_8_1_sequence_100_records.txt'),
os.path.join(test_dir, 'support_files', \
's_8_2_sequence_100_records.txt')]
self.illumina_map_fp = os.path.join(test_dir, 'support_files', \
's8_map_incomplete.txt')
self.fasta_fps=[os.path.join(test_dir,'support_files',
'test_split_lib_seqs.fasta')]
self.fasta_map_fp = os.path.join(test_dir, 'support_files', \
'fasta_mapping_file.txt')
tmp_dir = "/%s/test_wf" % environ['HOME']
self.dirs_to_remove.append(tmp_dir)
#self.qiime_config['temp_dir'] or '/tmp/'
if not exists(tmp_dir):
makedirs(tmp_dir)
# if test creates the temp dir, also remove it
#self.dirs_to_remove.append(tmp_dir)
self.wf_out="/%s/test_processed_data" % environ['HOME']
#print self.wf_out
self.dirs_to_remove.append(self.wf_out)
self.gg_out=os.path.join(self.wf_out,'gg_97_otus')
if not exists(self.gg_out):
makedirs(self.gg_out)
#self.dirs_to_remove.append(self.gg_out)
self.fasting_mapping_fp = get_tmp_filename(tmp_dir=tmp_dir,
prefix='qiime_wf_mapping',suffix='.txt')
fasting_mapping_f = open(self.fasting_mapping_fp,'w')
fasting_mapping_f.write(fasting_map)
fasting_mapping_f.close()
self.files_to_remove.append(self.fasting_mapping_fp)
self.params = parse_qiime_parameters(qiime_parameters_f)
signal.signal(signal.SIGALRM, timeout)
# set the 'alarm' to go off in allowed_seconds seconds
signal.alarm(allowed_seconds_per_test)
开发者ID:ElDeveloper,项目名称:qiime_web_app,代码行数:58,代码来源:test_submit_sff_and_metadata_to_db.py
示例13: setUp
def setUp(self):
""" """
self.test_data = get_test_data_fps()
self.files_to_remove = []
self.dirs_to_remove = []
# Create example output directory
tmp_dir = get_qiime_temp_dir()
self.test_out = mkdtemp(dir=tmp_dir,
prefix='core_qiime_analyses_test_',
suffix='')
self.dirs_to_remove.append(self.test_out)
self.qiime_config = load_qiime_config()
self.params = parse_qiime_parameters([])
self.params_sortmerna = parse_qiime_parameters(
['pick_otus:otu_picking_method\tsortmerna'])
self.params_sumaclust = parse_qiime_parameters(
['pick_otus:otu_picking_method\tsumaclust'])
initiate_timeout(60)
开发者ID:nbresnick,项目名称:qiime,代码行数:21,代码来源:test_upstream.py
示例14: setUp
def setUp(self):
""" """
self.qiime_config = load_qiime_config()
self.dirs_to_remove = []
self.files_to_remove = []
#this is specific to the web-apps only
test_dir = abspath(dirname(__file__))
sff_original_fp = os.path.join(test_dir, 'support_files', \
'Fasting_subset.sff')
# copy sff file to working directory
self.sff_dir = tempfile.mkdtemp()
self.dirs_to_remove.append(self.sff_dir)
self.sff_fp = os.path.join(self.sff_dir, 'Fasting_subset.sff')
copy(sff_original_fp, self.sff_fp)
self.files_to_remove.append(self.sff_fp)
tmp_dir = self.qiime_config['temp_dir'] or '/tmp/'
if not exists(tmp_dir):
makedirs(tmp_dir)
# if test creates the temp dir, also remove it
self.dirs_to_remove.append(tmp_dir)
self.wf_out = get_tmp_filename(tmp_dir=tmp_dir,
prefix='qiime_wf_out',suffix='',result_constructor=str)
self.dirs_to_remove.append(self.wf_out)
self.fasting_mapping_fp = get_tmp_filename(tmp_dir=tmp_dir,
prefix='qiime_wf_mapping',suffix='.txt')
fasting_mapping_f = open(self.fasting_mapping_fp,'w')
fasting_mapping_f.write(fasting_map)
fasting_mapping_f.close()
self.files_to_remove.append(self.fasting_mapping_fp)
working_dir = self.qiime_config['working_dir'] or './'
jobs_dir = join(working_dir,'jobs')
if not exists(jobs_dir):
# only clean up the jobs dir if it doesn't already exist
self.dirs_to_remove.append(jobs_dir)
self.params = parse_qiime_parameters(qiime_parameters_f.split('\n'))
signal.signal(signal.SIGALRM, timeout)
# set the 'alarm' to go off in allowed_seconds seconds
signal.alarm(allowed_seconds_per_test)
开发者ID:ElDeveloper,项目名称:qiime_web_app,代码行数:48,代码来源:test_per_study_sff.py
示例15: setUp
def setUp(self):
""" """
self.test_data = get_test_data_fps()
self.files_to_remove = []
self.dirs_to_remove = []
# Create example output directory
tmp_dir = get_qiime_temp_dir()
self.test_out = mkdtemp(dir=tmp_dir,
prefix='core_qiime_analyses_test_',
suffix='')
self.dirs_to_remove.append(self.test_out)
self.qiime_config = load_qiime_config()
self.params = parse_qiime_parameters([])
initiate_timeout(60)
开发者ID:Bonder-MJ,项目名称:qiime,代码行数:17,代码来源:test_upstream.py
示例16: test_parse_qiime_parameters
def test_parse_qiime_parameters(self):
"""parse_qiime_parameters: functions with valid input """
lines = ["#Don't edit this file!",\
"pick_otus:similarity 0.94",\
"pick_otus:otu_picking_method\tcdhit",\
"align_seqs:verbose",\
"assign_taxonomy:use_rdp\ttRuE",\
"assign_taxonomy:something\tNone",\
"",\
"#some_script:fake_parameter\t99.0"]
actual = parse_qiime_parameters(lines)
expected = {'pick_otus':\
{'similarity':'0.94', 'otu_picking_method':'cdhit'},\
'assign_taxonomy':\
{'use_rdp':None}}
self.assertEqual(actual,expected)
# default dict functions as expected -- looking up non-existant key
# returns empty dict
self.assertEqual(actual['some_other_script'],{})
开发者ID:Ecogenomics,项目名称:FrankenQIIME,代码行数:20,代码来源:test_parse.py
示例17: main
def main():
option_parser, opts, args =\
parse_command_line_parameters(**script_info)
parameters = {}
# get the tree insertion method to use
module = opts.insertion_method
# create output directory
output_dir = opts.output_dir
create_dir(output_dir)
# list of tree insertion methods
tree_insertion_module_names = \
{'raxml_v730': brokit.raxml_v730,
'parsinsert': brokit.parsinsert,
'pplacer': brokit.pplacer}
# load input sequences and convert to phylip since the tools require
# the query sequences to phylip-compliant names
load_aln = parse_fasta(open(opts.input_fasta_fp, 'U'))
aln = DenseAlignment(load_aln)
seqs, align_map = aln.toPhylip()
if opts.method_params_fp:
param_dict = parse_qiime_parameters(open(opts.method_params_fp, 'U'))
if module == 'raxml_v730':
# load the reference sequences
load_ref_aln = \
DenseAlignment(parse_fasta(open(opts.refseq_fp, 'U')))
# combine and load the reference plus query
combined_aln = parse_fasta(StringIO(load_ref_aln.toFasta() +
'\n' + aln.toFasta()))
# overwrite the alignment map
aln = DenseAlignment(combined_aln)
seqs, align_map = aln.toPhylip()
try:
parameters = param_dict['raxml']
except:
parameters = {}
tree = convert_tree_tips(align_map, opts.starting_tree_fp)
# write out the tree with phylip labels
updated_tree_fp = join(output_dir,
'%s_phylip_named_tree.tre' % (module))
write_updated_tree_file(updated_tree_fp, tree)
# set the primary parameters for raxml
parameters['-w'] = abspath(output_dir) + '/'
parameters["-n"] = split(splitext(get_tmp_filename())[0])[-1]
parameters["-t"] = updated_tree_fp
if "-f" not in parameters:
parameters["-f"] = 'v'
if "-m" not in parameters:
parameters["-m"] = 'GTRGAMMA'
elif module == 'pplacer':
try:
parameters = param_dict['pplacer']
except:
parameters = {}
# make sure stats file is passed
if not opts.stats_fp:
raise IOError(
'When using pplacer, the RAxML produced info file is required.')
# set the primary parameters for pplacer - allow for user-defined
parameters['--out-dir'] = abspath(output_dir) + '/'
parameters["-t"] = opts.starting_tree_fp
parameters['-r'] = opts.refseq_fp
parameters['-s'] = opts.stats_fp
elif module == 'parsinsert':
try:
parameters = param_dict['parsinsert']
except:
parameters = {}
# define log fp
log_fp = join(output_dir, 'parsinsert.log')
# define tax assignment values fp
tax_assign_fp = join(output_dir, 'parsinsert_assignments.log')
parameters["-l"] = log_fp
parameters["-o"] = tax_assign_fp
parameters["-s"] = opts.refseq_fp
parameters["-t"] = opts.starting_tree_fp
# call the module and return a tree object
result = \
tree_insertion_module_names[module].insert_sequences_into_tree(seqs,
moltype=DNA, params=parameters)
#.........这里部分代码省略.........
开发者ID:MicrobiomeResearch,项目名称:qiime,代码行数:101,代码来源:insert_seqs_into_tree.py
示例18: main
def main():
option_parser, opts, args =\
parse_command_line_parameters(**script_info)
#get all the options
cd_dir=path.join(opts.fs_fp,'arare')
tmp_prefix=get_tmp_filename('',suffix='').strip()
output_dir=path.join(opts.fs_fp,'arare','arare_'+tmp_prefix)
web_fp=path.join(opts.web_fp,'arare','arare_'+tmp_prefix)
otu_table_fp=opts.otu_table_fp
mapping_file_fp=opts.mapping_file_fp
file_name_prefix=opts.fname_prefix
user_id=int(opts.user_id)
meta_id=int(opts.meta_id)
bdiv_rarefied_at=int(opts.bdiv_rarefied_at)
jobs_to_start=opts.jobs_to_start
tree_fp=opts.tree_fp
command_handler=call_commands_serially
status_update_callback=no_status_updates
zip_fpath=opts.zip_fpath
zip_fpath_db=opts.zip_fpath_db
run_date=opts.run_date
force=True
try:
from data_access_connections import data_access_factory
from enums import ServerConfig
import cx_Oracle
data_access = data_access_factory(ServerConfig.data_access_type)
except ImportError:
print "NOT IMPORTING QIIMEDATAACCESS"
pass
try:
parameter_f = open(opts.params_path)
except IOError:
raise IOError,\
"Can't open parameters file (%s). Does it exist? Do you have read access?"\
% opts.params_path
params=parse_qiime_parameters(parameter_f)
try:
makedirs(output_dir)
except OSError:
if force:
pass
else:
# Since the analysis can take quite a while, I put this check
# in to help users avoid overwriting previous output.
print "Output directory already exists. Please choose "+\
"a different directory, or force overwrite with -f."
exit(1)
commands=[]
python_exe_fp = qiime_config['python_exe_fp']
script_dir = get_qiime_scripts_dir()
logger = WorkflowLogger(generate_log_fp(output_dir),
params=params,
qiime_config=qiime_config)
# determine whether to run alpha-diversity in serial or parallel
serial_or_parallel = params['serial_or_parallel']['method']
if serial_or_parallel=='Serial':
arare_cmd='%s %s/alpha_rarefaction.py -i %s -m %s -o %s -t %s -p %s -f' %\
(python_exe_fp, script_dir, otu_table_fp, mapping_file_fp, \
output_dir,tree_fp,opts.params_path)
else:
arare_cmd='%s %s/alpha_rarefaction.py -i %s -m %s -o %s -t %s -a -O 50 -p %s -f' %\
(python_exe_fp, script_dir, otu_table_fp, mapping_file_fp, \
output_dir,tree_fp,opts.params_path)
commands.append([('Alpha-Rarefaction',arare_cmd)])
command_handler(commands, status_update_callback, logger)
#zip the distance matrices
cmd_call='cd %s; zip -r %s %s' % (cd_dir,zip_fpath,'arare_'+tmp_prefix)
system(cmd_call)
#convert link into web-link
web_link=path.join(web_fp, 'alpha_rarefaction_plots',
'rarefaction_plots.html')
#add the distance matrices
valid=data_access.addMetaAnalysisFiles(True, int(meta_id), web_link,
'ARARE', run_date, 'ARARE')
if not valid:
raise ValueError, 'There was an issue uploading the filepaths to the DB!'
开发者ID:ElDeveloper,项目名称:qiime_web_app,代码行数:89,代码来源:alpha_rarefaction.py
示例19: main
def main():
option_parser, opts, args =\
parse_command_line_parameters(**script_info)
study_id = opts.study_id
for i in range(1,25):
output_dir='/home/wwwdevuser/user_data/studies/study_%s/processed_data_%s/' % (study_id,i)
if not exists(output_dir):
break
else:
continue
print output_dir
sff_fname=opts.sff_fname
map_fname = opts.map_fname
verbose = opts.verbose
print_only = opts.print_only
parallel = opts.parallel
denoise=opts.denoise
if parallel:
raise_error_on_parallel_unavailable()
try:
parameter_f = open(opts.parameter_fp)
except IOError:
raise IOError,\
"Can't open parameters file (%s). Does it exist? Do you have read access?"\
% opts.parameter_fp
try:
makedirs(output_dir)
except OSError:
if opts.force:
pass
else:
# Since the analysis can take quite a while, I put this check
# in to help users avoid overwriting previous output.
print "Output directory already exists. Please choose "+\
"a different directory, or force overwrite with -f."
exit(1)
if print_only:
command_handler = print_commands
else:
command_handler = web_app_call_commands_serially
if verbose:
status_update_callback = print_to_stdout
else:
status_update_callback = no_status_updates
run_process_sff_through_pick_otus(sff_input_fp=sff_fname,\
mapping_fp=map_fname,\
output_dir=output_dir,\
denoise=denoise,\
command_handler=command_handler,\
params=parse_qiime_parameters(parameter_f),\
qiime_config=qiime_config,\
parallel=parallel,\
status_update_callback=status_update_callback)
开发者ID:ElDeveloper,项目名称:qiime_web_app,代码行数:62,代码来源:process_sff_and_metadata_for_db_orig.py
示例20: main
def main():
option_parser, opts, args = parse_command_line_parameters(**script_info)
if opts.output_dir is None:
opts.output_dir = opts.output_filepath + '_dir'
if opts.parameter_fp:
try:
parameter_f = open(opts.parameter_fp, 'U')
except IOError:
raise IOError("Can't open parameters file (%s). Does it exist? Do you have read access?"
% opts.parameter_fp)
params = parse_qiime_parameters(parameter_f)
parameter_f.close()
else:
params = parse_qiime_parameters([])
# empty list returns empty defaultdict for now
try:
makedirs(opts.output_dir)
except OSError:
if opts.force:
pass
else:
# Since the analysis can take quite a while, I put this check
# in to help users avoid overwriting previous output.
option_parser.error("Output directory already exists. Please choose"
" a different directory, or force overwrite with -f.")
if opts.print_only:
command_handler = print_commands
else:
command_handler = call_commands_serially
if opts.verbose:
status_update_callback = print_to_stdout
else:
status_update_callback = no_status_updates
# set env variable
if opts.platform == 'flx':
existing_pyro_fp = os.environ['PYRO_LOOKUP_FILE']
new_pyro_fp = os.path.join(
os.path.split(existing_pyro_fp)[0],
'LookUp_E123.dat')
os.environ['PYRO_LOOKUP_FILE'] = new_pyro_fp
elif opts.platform == 'titanium':
existing_pyro_fp = os.environ['PYRO_LOOKUP_FILE']
new_pyro_fp = os.path.join(
os.path.split(existing_pyro_fp)[0],
'LookUp_Titanium.dat')
os.environ['PYRO_LOOKUP_FILE'] = new_pyro_fp
else:
raise RuntimeError(
'could not find PYRO_LOOKUP_FILE for platform ' +
platform)
if opts.truncate_len:
try:
truncate_len_int_check = int(opts.truncate_len)
truncate_len = str(truncate_len_int
|
请发表评论