X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=indicatorcommand.cpp;h=8d1d7f7a6864b261a0a31611a1e02270730cfd32;hb=e0ce7cbc93d7d2fbb753ca694182db092a0ea0e7;hp=1864f7473b0e21ff06d84ce76968f81374189954;hpb=0ca63a8165baa0afa459e644ebe140ba496d5ba0;p=mothur.git

diff --git a/indicatorcommand.cpp b/indicatorcommand.cpp
index 1864f74..8d1d7f7 100644
--- a/indicatorcommand.cpp
+++ b/indicatorcommand.cpp
@@ -99,7 +99,6 @@ IndicatorCommand::IndicatorCommand(string option)  {
 			m->clearGroups();
 			m->clearAllGroups();
 			m->Treenames.clear();
-			m->names.clear();
 			
 			vector<string> tempOutNames;
 			outputTypes["tree"] = tempOutNames;
@@ -236,13 +235,11 @@ int IndicatorCommand::execute(){
 			designMap->readDesignMap();
 			
 			//fill Groups - checks for "all" and for any typo groups
-			SharedUtil* util = new SharedUtil();
+			SharedUtil util;
 			vector<string> nameGroups = designMap->getNamesOfGroups();
-			util->setGroups(Groups, nameGroups);
+			util.setGroups(Groups, nameGroups);
 			designMap->setNamesOfGroups(nameGroups);
-			delete util;
 			
-			//loop through the Groups and fill Globaldata's Groups with the design file info
 			vector<string> namesSeqs = designMap->getNamesSeqs(Groups);
 			m->setGroups(namesSeqs);
 		}
@@ -320,8 +317,9 @@ int IndicatorCommand::execute(){
 				else { for (int i = 0; i < lookupFloat.size(); i++) {  delete lookupFloat[i];  } }
 				for (int i = 0; i < T.size(); i++) {  delete T[i];  }  delete treeMap; return 0; 
 			}
-				
-			T[0]->assembleTree();
+            
+			map<string, string> nameMap;	
+			T[0]->assembleTree(nameMap);
 					
 			/***************************************************/
 			//    create ouptut tree - respecting pickedGroups //
@@ -329,12 +327,11 @@ int IndicatorCommand::execute(){
 			Tree* outputTree = new Tree(m->getNumGroups(), treeMap); 
 			
 			outputTree->getSubTree(T[0], m->getGroups());
-			outputTree->assembleTree();
+			outputTree->assembleTree(nameMap);
 				
 			//no longer need original tree, we have output tree to use and label
 			for (int i = 0; i < T.size(); i++) {  delete T[i];  } 
 			
-					
 			if (m->control_pressed) { 
 				if (designfile != "") { delete designMap; }
 				if (sharedfile != "") {  for (int i = 0; i < lookup.size(); i++) {  delete lookup[i];  } }
@@ -470,17 +467,17 @@ int IndicatorCommand::GetIndicatorSpecies(){
 				
 			if (m->control_pressed) { out.close(); return 0; }
 			
-			out << (j+1) << '\t' << indicatorValues[j] << '\t'; 
+			out << m->currentBinLabels[j] << '\t' << indicatorValues[j] << '\t'; 
 			
 			if (pValues[j] > (1/(float)iters)) { out << pValues[j] << endl; } 
 			else { out << "<" << (1/(float)iters) << endl; }
 			
 			if (pValues[j] <= 0.05) {
-				cout << "OTU" << j+1 << '\t' << indicatorValues[j]  << '\t';
+				cout << m->currentBinLabels[j] << '\t' << indicatorValues[j]  << '\t';
 				string pValueString = "<" + toString((1/(float)iters)); 
 				if (pValues[j] > (1/(float)iters)) { pValueString = toString(pValues[j]); cout << pValues[j];} 
 				else { cout << "<" << (1/(float)iters); }
-				m->mothurOutJustToLog("OTU" + toString(j+1) + "\t" + toString(indicatorValues[j]) + "\t" + pValueString); 
+				m->mothurOutJustToLog(m->currentBinLabels[j] + "\t" + toString(indicatorValues[j]) + "\t" + pValueString); 
 				m->mothurOutEndLine(); 
 			}
 		}
@@ -516,7 +513,7 @@ int IndicatorCommand::GetIndicatorSpecies(Tree*& T){
 		
 		//print headings
 		out << "TreeNode\t";
-		for (int i = 0; i < numBins; i++) { out << "OTU" << (i+1) << "_IndValue" << '\t' << "pValue" << '\t'; }
+		for (int i = 0; i < numBins; i++) { out << m->currentBinLabels[i] << "_IndValue" << '\t' << "pValue" << '\t'; }
 		out << endl;
 		
 		m->mothurOutEndLine(); m->mothurOut("Node\tSpecies\tIndicatorValue\tpValue\n");
@@ -673,11 +670,11 @@ int IndicatorCommand::GetIndicatorSpecies(Tree*& T){
 				}
 				
 				if (pValues[j] <= 0.05) {
-					cout << i+1 << "\tOTU" << j+1 << '\t' << indicatorValues[j]  << '\t';
+					cout << i+1 << '\t' << m->currentBinLabels[j] << '\t' << indicatorValues[j]  << '\t';
 					string pValueString = "<" + toString((1/(float)iters)); 
 					if (pValues[j] > (1/(float)iters)) { pValueString = toString(pValues[j]); cout << pValues[j];} 
 					else { cout << "<" << (1/(float)iters); }
-					m->mothurOutJustToLog(toString(i) + "\tOTU" + toString(j+1) + "\t" + toString(indicatorValues[j]) + "\t" + pValueString); 
+					m->mothurOutJustToLog(toString(i) + "\t" + m->currentBinLabels[j] + "\t" + toString(indicatorValues[j]) + "\t" + pValueString); 
 					m->mothurOutEndLine(); 
 				}
 			}
@@ -1120,7 +1117,7 @@ vector<float> IndicatorCommand::getPValues(vector< vector<SharedRAbundFloatVecto
 	try {
 		vector<float> pvalues;
 		
-#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
+#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux) || (__linux__) || (__unix__) || (__unix)
 		if(processors == 1){
 			pvalues = driver(groupings, groupingsMap, num, indicatorValues, iters);
 			for (int i = 0; i < pvalues.size(); i++) { pvalues[i] /= (double)iters; }
@@ -1234,7 +1231,7 @@ vector<float> IndicatorCommand::getPValues(vector< vector<SharedRAbundVector*> >
 	try {
 		vector<float> pvalues;
 		
-#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux)
+#if defined (__APPLE__) || (__MACH__) || (linux) || (__linux) || (__linux__) || (__unix__) || (__unix)
 		if(processors == 1){
 			pvalues = driver(groupings, groupingsMap, num, indicatorValues, iters);
 			for (int i = 0; i < pvalues.size(); i++) { pvalues[i] /= (double)iters; }