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