def get_normalized_probeset(locus, inbredsetid): normalized_probesets = [] probesetxrefs = get_probesetxref_inbredsetid(locus, inbredsetid) for probesetxref in probesetxrefs: normalized_probeset = [] # probesetfreezeid = probesetxref[4] probesetfreeze = datastructure.get_probesetfreeze(probesetfreezeid) normalized_probeset.append(probesetfreeze[0]) normalized_probeset.append(probesetfreeze[1]) normalized_probeset.append(probesetfreeze[2]) # probesetid = probesetxref[0] probeset = get_probeset(probesetid) normalized_probeset.append(probeset[1]) normalized_probeset.append(probeset[2]) normalized_probeset.append(probeset[3]) normalized_probeset.append(probeset[4]) normalized_probeset.append(probeset[5]) normalized_probeset.append(probeset[6]) # normalized_probeset.append(probesetxref[1]) normalized_probeset.append(probesetxref[2]) # locus = probesetxref[3] geno = genotypes.get_geno(inbredsetid=inbredsetid, name=locus) normalized_probeset.append(geno[2]) normalized_probeset.append(geno[3]) # normalized_probesets.append(normalized_probeset) return normalized_probesets
def generate_probesets_2(probesetfreezesfile, outputdir): file = open(probesetfreezesfile, 'r') for line in file: line = line.strip() cells = line.split() probesetfreezeid = cells[0] probesetfreeze = datastructure.get_probesetfreeze(probesetfreezeid) probesetfreezeid = probesetfreeze[0] probesetfreezename = probesetfreeze[1] inbredset = datastructure.get_inbredset(probesetfreezeid) inbredsetid = inbredset[0] # outputfile = open("%s/%d_%s.txt" % (outputdir, probesetfreezeid, probesetfreezename), "w+") outputfile.write("%s\t" % "ProbeSetId") outputfile.write("%s\t" % "Symbol") outputfile.write("%s\t" % "Description") outputfile.write("%s\t" % "Chr") outputfile.write("%s\t" % "MB") outputfile.write("%s\t" % "Marker_Chr") outputfile.write("%s\t" % "Marker_MB") outputfile.write("%s\t" % "Mean_Expression") outputfile.write("%s\t" % "SE") outputfile.write("%s\t" % "LRS") outputfile.write("%s\t" % "pValue") outputfile.write("\n") outputfile.flush() # probesetxrefs = probesets.get_probesetxref(probesetfreezeid) print("%s:\n\t%d probesetxrefs" % (probesetfreeze, len(probesetxrefs))) for probesetxref in probesetxrefs: # probesetid = probesetxref[0] locus = probesetxref[2] lrs = probesetxref[3] pvalue = probesetxref[4] mean = probesetxref[5] se = probesetxref[6] # probeset = probesets.get_probeset(probesetid) probesetname = probeset[1] probesetsymbol = probeset[2] probesetdescription = probeset[3] probesetchr = probeset[5] probesetmb = probeset[6] # if locus is None or not locus: genochr = "" genomb = "" else: geno = genotypes.get_geno(inbredsetid=inbredsetid, name=locus) genochr = geno[2] genomb = geno[3] # outputfile.write("%s\t" % probesetname) outputfile.write("%s\t" % probesetsymbol) outputfile.write("%s\t" % probesetdescription) outputfile.write("%s\t" % probesetchr) outputfile.write("%s\t" % probesetmb) outputfile.write("%s\t" % genochr) outputfile.write("%s\t" % genomb) outputfile.write("%s\t" % mean) outputfile.write("%s\t" % se) outputfile.write("%s\t" % lrs) outputfile.write("%s\t" % pvalue) outputfile.write("\n") outputfile.flush() # outputfile.close() file.close()