fc724a23c9a0ad34a6ff6e530eb93945fd94a5c5
kate
Sun Mar 11 18:24:53 2018 -0700
Add score filtering. refs #17512
diff --git src/hg/hgc/interactClick.c src/hg/hgc/interactClick.c
index a65a7ca0..6bab38f 100644
--- src/hg/hgc/interactClick.c
+++ src/hg/hgc/interactClick.c
@@ -1,174 +1,178 @@
/* Details page for interact type tracks */
/* Copyright (C) 2018 The Regents of the University of California
* See README in this or parent directory for licensing information. */
#include "common.h"
#include "obscure.h"
#include "hdb.h"
#include "hgc.h"
#include "interact.h"
#include "interactUi.h"
static struct interact *getInteractFromTable(struct trackDb *tdb, char *item,
char *chrom, int start, int end)
/* Retrieve this item from track table */
{
struct sqlConnection *conn = NULL;
struct customTrack *ct = lookupCt(tdb->track);
char *table;
if (ct != NULL)
{
conn = hAllocConn(CUSTOM_TRASH);
table = ct->dbTableName;
}
else
{
conn = hAllocConnTrack(database, tdb);
table = tdb->table;
}
if (conn == NULL)
return NULL;
struct interact *inters = NULL, *inter;
char **row;
int offset;
char where[512];
if (isNotEmpty(item))
sqlSafefFrag(where, sizeof(where), "name='%s'", item);
struct sqlResult *sr = hRangeQuery(conn, table, chrom, start, end, where, &offset);
while ((row = sqlNextRow(sr)) != NULL)
{
inter = interactLoad(row+offset);
slAddHead(&inters, inter);
}
slReverse(&inters);
sqlFreeResult(&sr);
hFreeConn(&conn);
return inters;
}
static struct interact *getInteractFromFile(char *file, char *item, char *chrom, int start, int end)
/* Retrieve interact BED item from big file */
{
struct bbiFile *bbi = bigBedFileOpen(file);
struct lm *lm = lmInit(0);
struct bigBedInterval *bb, *bbList = bigBedIntervalQuery(bbi, chrom, start, end, 0, lm);
for (bb = bbList; bb != NULL; bb = bb->next)
{
char startBuf[16], endBuf[16];
char *row[32];
bigBedIntervalToRow(bb, chrom, startBuf, endBuf, row, ArraySize(row));
struct interact *interact = interactLoad(row);
if (interact == NULL)
continue;
if (sameString(interact->name, item))
return interact;
}
return NULL;
}
static struct interact *getInteract(struct trackDb *tdb, char *item,
char *chrom, int start, int end)
/* Retrieve interact BED item from track */
{
struct interact *interact = NULL;
char *file = trackDbSetting(tdb, "bigDataUrl");
if (file != NULL)
interact = getInteractFromFile(file, item, chrom, start, end);
else
interact = getInteractFromTable(tdb, item, chrom, start, end);
return interact;
}
void doInteractDetails(struct trackDb *tdb, char *item)
/* Details of interaction item */
{
char *chrom = cartString(cart, "c");
int start = cartInt(cart, "o");
int end = cartInt(cart, "t");
struct interact *inter = getInteract(tdb, item, chrom, start, end);
if (inter == NULL)
errAbort("Can't find interaction '%s'", item);
genericHeader(tdb, item);
char startBuf[1024], endBuf[1024], sizeBuf[1024];
-if (isNotEmpty(inter->name))
+if (isNotEmpty(inter->name) && differentString(".", inter->name))
printf("Interaction name: %s
\n", inter->name);
printf("Interaction region: ");
if (interactOtherChrom(inter))
printf("across chromosomes
");
else
{
sprintLongWithCommas(startBuf, inter->chromStart+1);
sprintLongWithCommas(endBuf, inter->chromEnd);
sprintLongWithCommas(sizeBuf, inter->chromEnd - inter->chromStart);
printf("%s:%s-%s",
inter->chrom, inter->chromStart, inter->chromEnd,
inter->chrom, startBuf, endBuf);
printf(" %s bp
\n", sizeBuf);
}
printf("Score: %d
\n", inter->score);
printf("Value: %0.3f
\n", inter->value);
-if (isNotEmpty(inter->exp))
+if (isNotEmpty(inter->exp) && differentString(".", inter->name))
printf("Experiment: %s
\n", inter->exp);
puts("
");
/* print info for both regions */
/* Use different labels:
1) directional (source/target)
2) non-directional same chrom (lower/upper)
3) non-directional other chrom (this/other)
*/
char *region1Label = "Source";
char *region1Chrom = inter->sourceChrom;
int region1Start = inter->sourceStart;
int region1End = inter->sourceEnd;
char *region1Name = inter->sourceName;
+if (sameString(".", inter->sourceName))
+ region1Name = "";
char *region2Label = "Target";
char *region2Chrom = inter->targetChrom;
int region2Start = inter->targetStart;
int region2End = inter->targetEnd;
char *region2Name = inter->targetName;
+if (sameString(".", inter->targetName))
+ region2Name = "";
if (!interactUiDirectional(tdb))
{
if (interactOtherChrom(inter))
{
region1Label = "This";
region2Label = "Other";
}
else
{
region1Label = "Lower";
region2Label = "Upper";
}
}
sprintLongWithCommas(startBuf, region1Start + 1);
sprintLongWithCommas(endBuf, region1End);
sprintLongWithCommas(sizeBuf, region1End - region1Start);
printf("%s region: %s "
"%s:%s-%s",
region1Label, region1Name, region1Chrom, region1Start+1, region1End,
region1Chrom, startBuf, endBuf);
printf(" %s bp
\n", sizeBuf);
sprintLongWithCommas(startBuf, region2Start+1);
sprintLongWithCommas(endBuf, region2End);
sprintLongWithCommas(sizeBuf, region2End - region2Start);
printf("%s region: %s "
"%s:%s-%s",
region2Label, region2Name, region2Chrom, region2Start+1, region2End,
region2Chrom, startBuf, endBuf);
printf(" %s bp
\n", sizeBuf);
#ifdef TODO /* TODO: get count and score stats of all interactions in window ?*/
double *scores;
AllocArray(scores, count);
#endif
}