def updateELevel(oFN, wigDir, rn=None, tn=None):
oNX = cgNexusFlat.Nexus(oFN, degPeak.degPeak)
oNX.load(['tcc', 'eLevel'], [rn, tn])
wigDict = cgWig.loadWigDict(wigDir)
for oID in oNX.eLevel:
expandTcc = bioLibCG.expandTcc(oNX.tcc[oID], 3)
coord_value = cgWig.getExpressionProfile(expandTcc, wigDict)
oNX.eLevel[oID] = max(coord_value.values())
oNX.save()
def updateELevel(oFN, wigDir, rn=None, tn=None):
'''Dont need to do it by chromosome because it is small enough'''
'''Also dont need to flip the strand because the wig is opposite as well'''
oNX = cgNexusFlat.Nexus(oFN, cgDegPeak.Peak)
oNX.load(['tcc', 'eLevel'], [rn, tn])
wigDict = cgWig.loadWigDict(wigDir)
for oID in oNX.eLevel:
expandTcc = bioLibCG.expandTcc(oNX.tcc[oID], 3)
coord_value = cgWig.getExpressionProfile(expandTcc, wigDict)
oNX.eLevel[oID] = max(coord_value.values())
oNX.save()
def updateELevel(oFN, wigDir, rn = None, tn = None):
oNX = cgNexusFlat.Nexus(oFN, degPeak.degPeak)
oNX.load(['tcc', 'eLevel'], [rn, tn])
wigDict = cgWig.loadWigDict(wigDir)
for oID in oNX.eLevel:
expandTcc = bioLibCG.expandTcc(oNX.tcc[oID], 3)
coord_value = cgWig.getExpressionProfile(expandTcc, wigDict)
oNX.eLevel[oID] = max(coord_value.values())
oNX.save()
def updateELevel(oFN, wigDir, rn = None, tn = None):
'''Dont need to do it by chromosome because it is small enough'''
'''Also dont need to flip the strand because the wig is opposite as well'''
oNX = cgNexusFlat.Nexus(oFN, cgDegPeak.Peak)
oNX.load(['tcc', 'eLevel'], [rn, tn])
wigDict = cgWig.loadWigDict(wigDir)
for oID in oNX.eLevel:
expandTcc = bioLibCG.expandTcc(oNX.tcc[oID], 3)
coord_value = cgWig.getExpressionProfile(expandTcc, wigDict)
oNX.eLevel[oID] = max(coord_value.values())
oNX.save()
