def getSRTSByAllele(alleles,pw,samples):
"""
Returns a dict containing sample rank templates for each alleles given the pathway(pw)
"""
srts = {}
for allele in alleles:
srts[allele] = dirac.getSRT(sd.getExpression(pw,[s for a,s in samples[allele]]))
return srts
def genSRTs(self, cstrain, pw):
#self.p.start("genSRTs")
if (cstrain, pw) not in self._srt_cache:
srts = None
sd = self._sd
mi = self._mi
samples = []
for allele in self.getAlleles(cstrain):
samples += [s for a,s in self._sample_x_allele[cstrain][allele]]
#self.p.start("getExpression")
expFrame = sd.getExpression( pw, samples)
srts = dirac.getSRT( expFrame )
self._srt_cache[(cstrain,pw)] = srts
return self._srt_cache[(cstrain, pw)]
def genSRTs(self, cstrain, pw):
#self.p.start("genSRTs")
if (cstrain, pw) not in self._srt_cache:
srts = None
sd = self._sd
mi = self._mi
samples = []
for allele in self.getAlleles(cstrain):
samples += [
s for a, s in self._sample_x_allele[cstrain][allele]
]
#self.p.start("getExpression")
expFrame = sd.getExpression(pw, samples)
srts = dirac.getSRT(expFrame)
self._srt_cache[(cstrain, pw)] = srts
return self._srt_cache[(cstrain, pw)]
