X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=mgclustercommand.cpp;h=e76dd4860b03826db116323af7e42528a591c526;hb=70491a12902e89b85cfa6b44a7b7fbe066ee2ac1;hp=b9c9d1eb845543bc38893068081b9367b913980b;hpb=1a5c2356c1b955c6ec024b2baf9f46377ee7c72e;p=mothur.git diff --git a/mgclustercommand.cpp b/mgclustercommand.cpp index b9c9d1e..e76dd48 100644 --- a/mgclustercommand.cpp +++ b/mgclustercommand.cpp @@ -247,7 +247,7 @@ int MGClusterCommand::execute(){ if(countfile != "") { //map nameMapCounts = m->readNames(namefile); ct = new CountTable(); - ct->readTable(countfile); + ct->readTable(countfile, false); rabund = new RAbundVector(); createRabund(ct, list, rabund); }else { @@ -271,10 +271,10 @@ int MGClusterCommand::execute(){ map variables; variables["[filename]"] = fileroot; - if (countfile != "") { variables["[tag2]"] = "unique_list"; } variables["[clustertag]"] = tag; string sabundFileName = getOutputFileName("sabund", variables); string rabundFileName = getOutputFileName("rabund", variables); + if (countfile != "") { variables["[tag2]"] = "unique_list"; } string listFileName = getOutputFileName("list", variables); if (countfile == "") { @@ -312,10 +312,13 @@ int MGClusterCommand::execute(){ outputTypes.clear(); return 0; } - + + //cluster using cluster classes while (distMatrix->getSmallDist() < cutoff && distMatrix->getNNodes() > 0){ + if (m->debug) { cout << "numNodes=" << distMatrix->getNNodes() << " smallDist = " << distMatrix->getSmallDist() << endl; } + cluster->update(cutoff); if (m->control_pressed) {