X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=bstick.cpp;h=67507a9967ef82effe923662298960ecf00618b8;hp=566099772e8cecc3d230dbb033aa1bebed1ab885;hb=050a3ff02473a3d4c0980964e1a9ebe52e55d6b8;hpb=510b1cfc25cd79391d6973ca20c5ec25fb1bb3b2 diff --git a/bstick.cpp b/bstick.cpp index 5660997..67507a9 100644 --- a/bstick.cpp +++ b/bstick.cpp @@ -18,7 +18,7 @@ double BStick::invSum(int index, double numSpec) sum += 1/(double)i; return sum; } - +/***********************************************************************/ RAbundVector BStick::getRAbundVector(SAbundVector* rank){ vector rData; int mr = 1; @@ -26,7 +26,7 @@ RAbundVector BStick::getRAbundVector(SAbundVector* rank){ int ns = 0; for(int i = rank->size()-1; i > 0; i--) { - int cur = rank->get(i); + double cur = rank->get(i); if(mr == 1 && cur > 0) mr = i; nb += cur; @@ -67,8 +67,8 @@ EstOutput BStick::getValues(SAbundVector* rank){ data[1] = 0.886/sqrt(rdata.size()); data[2] = 1.031/sqrt(rdata.size()); - /*mothurOut(critVal); mothurOutEndLine(); - mothurOut("If D-Statistic is less than the critical value then the data fits the Broken Stick model w/ 95% confidence.\n\n");*/ + /*m->mothurOut(critVal); m->mothurOutEndLine(); + m->mothurOut("If D-Statistic is less than the critical value then the data fits the Broken Stick model w/ 95% confidence.\n");*/ if (isnan(data[0]) || isinf(data[0])) { data[0] = 0; } @@ -78,7 +78,7 @@ EstOutput BStick::getValues(SAbundVector* rank){ return data; } catch(exception& e) { - errorOut(e, "BStick", "getValues"); + m->errorOut(e, "BStick", "getValues"); exit(1); } }