def main(): in_folder = sys.argv[1] out_dir = sys.argv[2] nickname = sys.argv[3] cutoffs = sys.argv[4] num_iter = int(sys.argv[5]) cutoffs = [float(x) for x in cutoffs.split(',')] cutoffs = sorted([1.01 if x == 1 else x for x in cutoffs]) poss_mat, col_headings = pt.get_pres_abs_mat(in_folder) results = simulate_reordering(poss_mat, col_headings, cutoffs, num_iter) cutoffs = [min(x, 1.0) for x in cutoffs] for result, cutoff in zip(results, cutoffs): file_name = out_dir + '/' + nickname + '_' + str(cutoff) + '.Rtab' make_rtab(result, file_name)
def main(): folder = sys.argv[1] poss_mat, cols = pt.get_pres_abs_mat(folder) gene_counts_dict = get_gene_counts_dict(poss_mat) print_c_for_fsgm(gene_counts_dict, len(cols))