X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=ace.cpp;h=54d0946eab9cf0b2230ef45ef8749b1feecb9da0;hb=e4827e0945cbda536064e5a345996b2a7dfcbb81;hp=d074310ab1678ca44abbb23372ee74b41321999b;hpb=addc715b6d7ea52440751cec10edad03e1a19b37;p=mothur.git diff --git a/ace.cpp b/ace.cpp index d074310..54d0946 100644 --- a/ace.cpp +++ b/ace.cpp @@ -13,9 +13,7 @@ EstOutput Ace::getValues(SAbundVector* rank) { try { - cout << "abund = " << abund << "\n"; data.resize(3,0); - // vector aceData(3,0); double ace, acelci, acehci; int nrare = 0; @@ -44,7 +42,7 @@ EstOutput Ace::getValues(SAbundVector* rank) { if(denom <= 0.0){ term1=0.0000; } else { term1 = (double)(srare * numsum)/(double)denom - 1.0; } if(term1 >= 0.0){ gamace = term1; } else { gamace = 0.0; } - + if(gamace >= 0.64){ gamace = gamace * (1 + (nrare * (1 - Cace) * numsum) / denom); if(gamace<0){ gamace = 0; }