00001 #include "log.h" 00002 #include "format.h" 00003 #include "gfr.h" 00004 00005 00006 00007 int main (int argc, char *argv[]) 00008 { 00009 GfrEntry *currGE; 00010 double pvalueCutOff; 00011 int count; 00012 int countRemoved; 00013 00014 if (argc != 2) { 00015 usage ("%s <pvalueCutOff>",argv[0]); 00016 } 00017 count = 0; 00018 countRemoved = 0; 00019 pvalueCutOff = atof (argv[1]); 00020 gfr_init ("-"); 00021 puts (gfr_writeHeader ()); 00022 while (currGE = gfr_nextEntry ()){ 00023 if ( (MAX (currGE->pValueAB,currGE->pValueBA) < pvalueCutOff) && (currGE->pValueAB != -1) ) { 00024 countRemoved++; 00025 continue; 00026 } 00027 puts (gfr_writeGfrEntry (currGE)); 00028 count++; 00029 } 00030 gfr_deInit (); 00031 warn ("%s_pvalueCutOff: %f",argv[0],pvalueCutOff); 00032 warn ("%s_numRemoved: %d",argv[0],countRemoved); 00033 warn ("%s_numGfrEntries: %d",argv[0],count); 00034 return 0; 00035 } 00036