Hi All,
I have successfully run the CompareSelectivePressure bf with my data, but when I try the CompareSelectivePressureIVL bf I always get the following message:
An error occurred:
Operation MAccess is not defined for 0
Current BL Command: Print to the File stout:| Codon:
Format(site+1,4,0)
| Sample 1 (dS,dN Int, dN Leaf):
Format(fullSites[site|][0],6,2)
Format(fullSites[site|][1],6,2)
Format(fullSites[site|][2],6,2)
| Sample 2 (dS, d
Current operation/job has been terminated

Any guess whats goin' on?

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.

Dear CGS,

It is indeed a bug in the batch file; Art's fix should take care of it and I will update the code to work properly in the next build.

Thanks for pointing this out.

Cheers,
Sergei

Thanks for the quick reply Sergei!

Dear CGS,

No problem. The credit goes to Art though, as he valiantly manned the trenches while I was out of the country:)

Cheers,
Sergei