Hi CGS,
From a cursory inspection of the batch file, it looks like the matrix variable
doneSites got renamed to
fullSites in the fprintf statement without having ever been declared. Try reverting these to the original name:
Code:/*------------------------------------------------------------------------*/
function ReportSite2 (siteI, siteM)
{
if (siteI != siteM)
{
for (k=0; k < 13; k=k+1)
{
doneSites[siteI][k] = doneSites[siteM][k];
}
}
/*fprintf (stdout, "Site ", Format(siteI+1,4,0)," ", Format(doneSites[siteI][0],5,3) ,",", Format(doneSites[siteI][1],5,3), "/", Format(doneSites[siteI][2],5,3),
" :: ", Format(doneSites[siteI][3],5,3) ,",", Format(doneSites[siteI][4],5,3), "/", Format(doneSites[siteI][5],5,3),
" :: ", Format(doneSites[siteI][6],5,3) ,",", Format(doneSites[siteI][7],5,3), ",", Format(doneSites[siteI][8],5,3), ",", Format(doneSites[siteI][9],5,3), "/", Format(doneSites[siteI][10],5,3),
" LRT = ",Format(doneSites[siteM][11],5,3),
" p = ",Format(doneSites[siteM][12],5,3), "\n"); */
fprintf (stdout, "| Codon: ", Format(siteI+1,4,0),
"| Sample 1 (dS,dN Int, dN Leaf): ",
Format(doneSites[siteI][0],6,2),
Format(doneSites[siteI][1],6,2),
Format(doneSites[siteI][2],8,2),
"| Sample 2 (dS,dN Int, dN Leaf): ",
Format(doneSites[siteI][3],6,2),
Format(doneSites[siteI][4],6,2),
Format(doneSites[siteI][5],8,2),
"| LRT: ", Format(doneSites[siteI][11],7,2),
"| p: ", Format(doneSites[siteI][12],5,2), "\n");
return 0;
}
Cheers,
- Art.
Home
Pages: 1
CompareSelectivePressureIVL bug? (Read 2889 times)
