def geneprobelistmatch(gene,probelist):
if 'Yeast6kArray' not in dir(): from TAMO.DataSources import Yeast6kArray
s_orf = SGD.gene2orf(gene)
orfs = []
for probe in probelist:
orfs.extend(Yeast6kArray.probe2orfs(probe))
#print "%4d --> %4d "%(len(probelist),len(orfs)),
return orforflistmatch(s_orf,orfs)
def probelist2categories(probelist,thresh=0.05):
"""
Which categories are overrepresented among probes bound in in a Yeast6k array experiment?
"""
if 'Yeast6kArray' not in dir(): from TAMO.DataSources import Yeast6kArray
orfs = []
preorfs = []
for probe in probelist:
preorfs.extend(Yeast6kArray.probe2orfs(probe))
for o in preorfs:
if (not o in orfs): orfs.append(o)
return orflist2categories(orfs,thresh)
