src/utils/bedToBigBed/bedToBigBed.c 1.13
1.13 2009/11/05 20:00:23 kent
Fixing bug where sumSquares was actually just sum.
Index: src/utils/bedToBigBed/bedToBigBed.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/utils/bedToBigBed/bedToBigBed.c,v
retrieving revision 1.12
retrieving revision 1.13
diff -b -B -U 4 -r1.12 -r1.13
--- src/utils/bedToBigBed/bedToBigBed.c 5 Nov 2009 19:35:01 -0000 1.12
+++ src/utils/bedToBigBed/bedToBigBed.c 5 Nov 2009 20:00:23 -0000 1.13
@@ -312,9 +312,9 @@
sum->validCount += overlap;
if (sum->minVal > val) sum->minVal = val;
if (sum->maxVal < val) sum->maxVal = val;
sum->sumData += val * overlap;
- sum->sumSquares += val * overlap;
+ sum->sumSquares += val*val * overlap;
bbiOutputOneSummaryFurtherReduce(sum, &twiceReducedList, doubleReductionSize,
&boundsPt, boundsEnd, usage->size, lm, f);
/* Move summary to next part. */
@@ -331,9 +331,9 @@
sum->validCount += size;
if (sum->minVal > val) sum->minVal = val;
if (sum->maxVal < val) sum->maxVal = val;
sum->sumData += val * size;
- sum->sumSquares += val * size;
+ sum->sumSquares += val*val * size;
}
if (sum != NULL)
{
bbiOutputOneSummaryFurtherReduce(sum, &twiceReducedList, doubleReductionSize,