def main(self, args, command): ''' Parameters ---------- Output ------ ''' self._check_general(args) self._logging_setup(args) logging.info("Running command: %s" % ' '.join(command)) if args.subparser_name == self.DATA: self._check_data(args) d = Data() d.do() if args.subparser_name == self.ANNOTATE: self._check_annotate(args) a = Annotate(# Define inputs and outputs args.output, # Define type of annotation to be carried out args.ko, args.pfam, args.tigrfam, args.cog, args.hypothetical, # Cutoffs args.evalue, args.bit, args.id, args.aln_query, args.aln_reference, args.cascaded, args.c, # Parameters args.threads, args.parallel, args.suffix) a.do(args.genome_directory, args.protein_directory, args.genome_files, args.protein_files) elif args.subparser_name == self.CLASSIFY: self._check_classify(args) c = Classify() c.do(args.custom_modules, args.cutoff, args.genome_and_annotation_file, args.genome_and_annotation_matrix, args.output) elif args.subparser_name == self.ENRICHMENT: self._check_enrichment(args) e = Enrichment() e.do(# Inputs args.annotation_matrix, args.annotation_file, args.metadata, args.modules, args.abundances, args.do_all, args.do_ivi, args.do_gvg, args.do_ivg, args.pval_cutoff, args.proportions_cutoff, args.threshold, args.multi_test_correction, args.output) elif args.subparser_name == self.COMPARE: self._check_compare(args) c = Compare() c.do(args.enrichm_annotate_output) elif(args.subparser_name == NetworkAnalyser.PATHWAY or args.subparser_name == NetworkAnalyser.EXPLORE or args.subparser_name == NetworkAnalyser.TRAVERSE): self._check_network(args) na=NetworkAnalyser(args.metadata) na.do(args.matrix, args.transcriptome, args.metabolome, args.depth, args.filter, args.limit, args.queries, args.subparser_name, args.starting_compounds, args.steps, args.number_of_queries, args.output) logging.info('Done!')