be4311c07e14feb728abc6425ee606ffaa611a58 markd Fri Jan 22 06:46:58 2021 -0800 merge with master diff --git src/hg/utils/otto/isca/validateISCA.sh src/hg/utils/otto/isca/validateISCA.sh index c1def64..2c0f717 100755 --- src/hg/utils/otto/isca/validateISCA.sh +++ src/hg/utils/otto/isca/validateISCA.sh @@ -18,28 +18,28 @@ -o $i == "iscaBenignGainCum" -o $i == "iscaBenignLossCum" then fields='*' else fields="chrom,chromStart,chromEnd,name,score,strand,thickStart,thickEnd" fi echo "select $fields from $i" | hgsql $db | tail -n +2 | sort > $i.out f=$i"New" echo "select $fields from $f" |hgsql $db | tail -n +2 | sort > $f.out oldCount=`cat $i.out | wc -l` newCount=`cat $f.out | wc -l` common=`join -t '\001' $i.out $f.out | wc -l` onlyOld=`join -t '\001' -v 1 $i.out $f.out | wc -l` onlyNew=`join -t '\001' -v 2 $i.out $f.out | wc -l` echo $i $newCount "-" $onlyNew "=" $common "=" $oldCount "-" $onlyOld -# rm $i.out $f.out + rm $i.out $f.out done > newISCA.stats cat newISCA.stats | awk -v db=$db -v tooMuch=$tooMuch ' { if (($4/$6 > tooMuch) || ($10/$6 > tooMuch)) { print "validate on " db "." $1 " failed:" $4,$6,$4/$6,$10,$6,$10/$6; exit 1 } }' exit 0