X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=getoturepcommand.cpp;h=c64549334ded9f4d32140ad64cbe898c5c41a023;hp=61092372b18ae196a652d0b0f66c6b04d22d0266;hb=a2cde58c1e72199498a2142983ef040dce36da10;hpb=ac03f1f6c27b5bfdf2cfb6d45c3667c3e0281f51 diff --git a/getoturepcommand.cpp b/getoturepcommand.cpp index 6109237..c645493 100644 --- a/getoturepcommand.cpp +++ b/getoturepcommand.cpp @@ -760,6 +760,8 @@ string GetOTURepCommand::findRepAbund(vector names, string group) { try{ vector reps; string rep = "notFound"; + + if (m->debug) { m->mothurOut("[DEBUG]: group=" + group + " names.size() = " + toString(names.size()) + " " + names[0] + "\n"); } if ((names.size() == 1)) { return names[0]; @@ -773,7 +775,7 @@ string GetOTURepCommand::findRepAbund(vector names, string group) { if (countfile != "") { //if countfile is not blank then we can assume the list file contains only uniques, otherwise we assume list file contains everyone. int numRep = 0; if (group != "") { numRep = ct.getGroupCount(names[i], group); } - else { numRep = ct.getGroupCount(names[i]); } + else { numRep = ct.getNumSeqs(names[i]); } if (numRep > maxAbund) { reps.clear(); reps.push_back(names[i]); @@ -834,7 +836,7 @@ string GetOTURepCommand::findRep(vector names, string group) { if (countfile != "") { //if countfile is not blank then we can assume the list file contains only uniques, otherwise we assume list file contains everyone. int numRep = 0; if (group != "") { numRep = ct.getGroupCount(names[i], group); } - else { numRep = ct.getGroupCount(names[i]); } + else { numRep = ct.getNumSeqs(names[i]); } for (int j = 1; j < numRep; j++) { //don't add yourself again seqIndex.push_back(nameToIndex[names[i]]); }