HyPhy message board
http://www.hyphy.org/cgi-bin/hyphy_forums/YaBB.pl
HYPHY Package >> Contributions >> CodonModelCompare - save results to file
http://www.hyphy.org/cgi-bin/hyphy_forums/YaBB.pl?num=1125942257

Message started by avilella on Sep 5th, 2005 at 10:44am

Title: CodonModelCompare - save results to file
Post by avilella on Sep 5th, 2005 at 10:44am
Dears,

I found out that CodonModelCompare.bf differs from NucModelCompare.bf in that the former won't ask for a results file, while the latter does.

As I find more convenient for batched jobs to save the results in a separate file, I took to liberty of adding the results file prompt in CodonModelCompare.bf.

In hopes of being useful to others, here is the diff:
-----
--- CodonModelCompare.bf.ori    2005-09-05 19:24:51.967872600 +0200
+++ CodonModelCompare.bf        2005-09-05 19:28:44.062588840 +0200
@@ -631,6 +631,10 @@
       fscanf  (stdin,"Number", rejectAt);
}

+SetDialogPrompt ("Save results to:");
+
+fprintf (PROMPT_FOR_FILE, CLEAR_FILE);
+BASE_PATH = LAST_FILE_PATH;


KEEP_OPTIMAL_ORDER = 1;
@@ -942,6 +946,9 @@
               }
       }

+       fprintf (BASE_PATH,CLEAR_FILE,"\n\nRemaining models:\n\n#   |  Model   | # prm |    lnL    |      LRT       |    AIC     |   P-Value        |");
+       fprintf (BASE_PATH,"\n----|----------|-------|-----------|----------------|------------|------------------|");
+
       fprintf (stdout,"\n\nRemaining models:\n\n#   |  Model   | # prm |    lnL    |      LRT       |    AIC     |   P-Value        |");
       fprintf (stdout,"\n----|----------|-------|-----------|----------------|------------|------------------|");

@@ -970,11 +977,15 @@
                       modelNum = modelNum + 1;
                       PRINT_DIGITS = 3;
                       fprintf (stdout,"\n",v2);
+                       fprintf (BASE_PATH,"\n",v2);
                       PRINT_DIGITS = 1;
                       fprintf (stdout," | (",0,resultCache[v2][0],resultCache[v2][1],resultCache[v2][2],resultCache[v2][3],resultCache[v2][4],") | ");
                       fprintf (stdout,Format (np,5,0));
+                       fprintf (BASE_PATH," | (",0,resultCache[v2][0],resultCache[v2][1],resultCache[v2][2],resultCache[v2][3],resultCache[v2][4],") | ");
+                       fprintf (BASE_PATH,Format (np,5,0));
                       PRINT_DIGITS = 8;
                       fprintf (stdout, " |  ",lnL," | ",Format(LRT,14,3), " |  ", AIC, "  |  ", );
+                       fprintf (BASE_PATH, " |  ",lnL," | ",Format(LRT,14,3), " |  ", AIC, "  |  ", );
                       PRINT_DIGITS = 15;
                       if (LRT==0)
                       {
@@ -990,6 +1001,7 @@
                               v4 = v2;
                       }
                       fprintf (stdout,pValue," |");
+                       fprintf (BASE_PATH,pValue," |");

               }
       }
@@ -1002,6 +1014,7 @@
       }

       fprintf (stdout, "\n\nAIC based winner: (", modelString, ") with AIC = ", v5, "\n\n");
+       fprintf (BASE_PATH, "\n\nAIC based winner: (", modelString, ") with AIC = ", v5, "\n\n");

       dummy = printModelMatrix (modelString);

@@ -1021,12 +1034,14 @@
       if (Abs(modelString2))
       {
               fprintf (stdout, "\nThis model is better known as:", modelString2, "\n");
+               fprintf (BASE_PATH, "\nThis model is better known as:", modelString2, "\n");
       }

}
else
{
       fprintf (stdout, "\nGeneral Reversible Model is the winner!\n");
+       fprintf (BASE_PATH, "\nGeneral Reversible Model is the winner!\n");
}


Title: Re: CodonModelCompare - save results to file
Post by Sergei on Sep 6th, 2005 at 12:01pm
Dear Albert,

Thanks for the modified version of the analysis. I will include the additions in the next build of HyPhy.

Cheers,
Sergei

HyPhy message board » Powered by YaBB 2.5.2!
YaBB Forum Software © 2000-2020. All Rights Reserved.