src/hg/makeDb/hgBedsToBedExps/hgBedsToBedExps.c 1.6
1.6 2010/05/06 18:02:38 kent
Making it so scores go from 0-1000
Index: src/hg/makeDb/hgBedsToBedExps/hgBedsToBedExps.c
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/makeDb/hgBedsToBedExps/hgBedsToBedExps.c,v
retrieving revision 1.5
retrieving revision 1.6
diff -b -B -U 4 -r1.5 -r1.6
--- src/hg/makeDb/hgBedsToBedExps/hgBedsToBedExps.c 6 May 2010 17:53:41 -0000 1.5
+++ src/hg/makeDb/hgBedsToBedExps/hgBedsToBedExps.c 6 May 2010 18:02:38 -0000 1.6
@@ -185,17 +185,23 @@
for (i=0; i<sourceCount; ++i)
levels[i] = 0;
/* Set levels according to sourceSite list. */
+ double maxLevel = 0.0;
for (site = range->val; site != NULL; site = site->next)
{
- levels[site->source->sourceIx] = site->site->dataValue;
+ double level = site->site->dataValue;
+ if (level > 1000)
+ level = 1000;
+ if (level > maxLevel)
+ maxLevel = level;
+ levels[site->source->sourceIx] = level;
}
/* Output. */
fprintf(f, "%s\t%d\t%d\t", chrom, range->start, range->end);
fprintf(f, "%s\t", factor->factor);
- fprintf(f, "0\t");
+ fprintf(f, "%d\t", round(maxLevel)); /* score */
fprintf(f, "+\t"); /* strand.... */
fprintf(f, "%d\t%d\t", range->start, range->end);
fprintf(f, "0\t"); /* itemRgb */
fprintf(f, "1\t"); /* block count */