def main(): global options # setup Tanimoto class tanimoto = Tanimoto(options.check, options.verbose) tanimoto.import_data(options.drugs_file, options.targets_file) # setup PValue class p_value = PValue(options.check, options.verbose) p_value.set_tanimoto_class(tanimoto) p_value.set_n_bootstrap(options.n_bootstrap) p_value.set_tc_cutoff(options.tc_cutoff) p_value.set_use_random_seed(options.use_random_seed) p_value.import_targets_file(options.targets_file) # execute p_value class value = p_value.compute(options.protA_ID, options.protB_ID) print 'p-value: %.3f' % value