except Exception, inst: sys.stderr.write("Unexpected error opening %s: %s\n" % (infile, inst)) sys.exit() if (verbose): print "Reading " + str(argv[0]) elems = tree.getroot() thresholds = [ 0.0015, 0.0025, 0.005, 0.0075, 0.01, 0.015, 0.020, 0.030, 0.050, 0.060, 0.075, 0.090, 0.1, 0.15, 0.2 ] proteins_percolator = per.getProteins(elems) peptides = per.getPeptides(elems) psms = per.getPSMs(elems) percolator_proteins_target = list() percolator_proteins_decoy = list() twohit_proteins_target = list() twohit_proteins_decoy = list() twohit_proteins_target_unique = list() twohit_proteins_decoy_unique = list() parsimony_proteins_target = list() parsimony_proteins_decoy = list() mayu_proteins_target = list() mayu_proteins_decoy = list() mayu_proteins_target_1psm = list()
parser = etree.XMLParser(ns_clean=False, huge_tree=True) try: tree = etree.parse(infile,parser) except Exception, inst: sys.stderr.write("Unexpected error opening %s: %s\n" % (infile, inst)) sys.exit() if(verbose): print "Reading " + str(argv[0]) elems = tree.getroot() percolatorPSMs = per.getPSMs(elems) percolatorPeptides = per.getPeptides(elems) percolatorProteins = per.getProteins(elems) if(verbose): print "Read " + str(len(percolatorPSMs)) + " PSMs" print "Read " + str(len(percolatorPeptides)) + " Peptides" print "Read " + str(len(percolatorProteins)) + " Proteins" if(fdr < 1.0 and fdr > 0.0): percolatorPSMs = [psm for psm in percolatorPSMs if psm.qvalue <= fdr] percolatorPeptides = [pep for pep in percolatorPeptides if pep.qvalue <= fdr] percolatorProteins = [prot for prot in percolatorProteins if prot.qvalue <= fdr] if(len(percolatorPSMs) > 0): if(verbose): print "writing in " + str(OutputFile) utils.writePsms(percolatorPSMs, "psms_" + OutputFile)
parser = etree.XMLParser(ns_clean=False, huge_tree=True) try: tree = etree.parse(infile, parser) except Exception, inst: sys.stderr.write("Unexpected error opening %s: %s\n" % (infile, inst)) sys.exit() if (verbose): print "Reading " + str(argv[0]) elems = tree.getroot() percolatorPSMs = per.getPSMs(elems) percolatorPeptides = per.getPeptides(elems) percolatorProteins = per.getProteins(elems) if (verbose): print "Read " + str(len(percolatorPSMs)) + " PSMs" print "Read " + str(len(percolatorPeptides)) + " Peptides" print "Read " + str(len(percolatorProteins)) + " Proteins" if (fdr < 1.0 and fdr > 0.0): percolatorPSMs = [ psm for psm in percolatorPSMs if psm.qvalue <= fdr ] percolatorPeptides = [ pep for pep in percolatorPeptides if pep.qvalue <= fdr ] percolatorProteins = [ prot for prot in percolatorProteins if prot.qvalue <= fdr
parser = etree.XMLParser(ns_clean=False, huge_tree=True) try: tree = etree.parse(infile,parser) except Exception, inst: sys.stderr.write("Unexpected error opening %s: %s\n" % (infile, inst)) sys.exit() if(verbose): print "Reading " + str(argv[0]) elems = tree.getroot() thresholds = [0.0015,0.0025,0.005,0.0075,0.01,0.015,0.020,0.030,0.050,0.060,0.075,0.090,0.1,0.15,0.2] proteins_percolator= per.getProteins(elems) peptides = per.getPeptides(elems) psms = per.getPSMs(elems) percolator_proteins_target = list() percolator_proteins_decoy = list() twohit_proteins_target = list() twohit_proteins_decoy = list() twohit_proteins_target_unique = list() twohit_proteins_decoy_unique = list() parsimony_proteins_target = list() parsimony_proteins_decoy = list() mayu_proteins_target = list() mayu_proteins_decoy = list() mayu_proteins_target_1psm = list()