zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Let's go the other way. 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@5422 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
chartl 2011-03-11 17:36:41 +00:00
parent 572b2707f2
commit be1f6af815
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
python/recalAssociation.py
View File

@ -21,11 +21,11 @@ def recalibrate(fname,hist,out):
denom = sum(hist.values()) denom = sum(hist.values())
cumQuals = dict() cumQuals = dict()
for key1 in hist: for key1 in hist:
sumLess = 0 sumMore = 1
for key2 in hist: for key2 in hist:
if ( key2 <= key1 ): if ( key2 > key1 ):
sumLess += hist[key2] sumMore += hist[key2]
cumQuals[key1] = int(-10*log10((0.0+sumLess)/denom)) cumQuals[key1] = min(int(-10*log10((0.0+sumMore)/denom)),150)
use = parseTDF use = parseTDF
if ( fname.endswith(".wig") ): if ( fname.endswith(".wig") ):
use = parseWiggle use = parseWiggle