void bigWigCorrelateList(char *listFile)
/* Correlate all files in list to each other */
{
char **fileNames = NULL;
int fileCount = 0;
char *buf = NULL;
readAllWords(listFile, &fileNames, &fileCount, &buf);
int i;
for (i=0; i<fileCount; ++i)
{
char *aPath = fileNames[i];
char aName[FILENAME_LEN];
if (rootNames)
splitPath(aPath, NULL, aName, NULL);
else
safef(aName, sizeof(aName), "%s", aPath);
int j;
for (j=i+1; j<fileCount; ++j)
{
char *bPath = fileNames[j];
char bName[FILENAME_LEN];
if (rootNames)
splitPath(bPath, NULL, bName, NULL);
else
safef(bName, sizeof(bName), "%s", bPath);
struct correlate *c = bigWigCorrelate(aPath, bPath);
printf("%s\t%s\t%g\n", aName, bName, correlateResult(c));
correlateFree(&c);
}
}
}
int main(int argc, char *argv[])
/* Process command line. */
{
optionInit(&argc, argv, options);
if (argc != 3)
usage();
restrictFile = optionVal("restrict", restrictFile);
threshold = optionDouble("threshold", threshold);
bigWigCorrelate(argv[1], argv[2]);
return 0;
}
void bigWigCorrelatePair(char *aFileName, char *bFileName)
/* Correlate a pair of bigWigs and print result */
{
struct correlate *c = bigWigCorrelate(aFileName, bFileName);
printf("%g\n", correlateResult(c));
}