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()
