]> git.donarmstrong.com Git - mothur.git/commitdiff
adding misclassification output file
authorkdiverson <>
Wed, 11 Sep 2013 18:41:27 +0000 (14:41 -0400)
committerkd.iverson@gmail.com <>
Wed, 11 Sep 2013 18:41:27 +0000 (14:41 -0400)
classifysharedcommand.cpp

index c79e0e8f2aefcc39c7358fd7040681a8e0e7f924..c7eb6cd0daa18f5627824a94d297270a1ff32147 100755 (executable)
@@ -359,8 +359,11 @@ void ClassifySharedCommand::processSharedAndDesignData(vector<SharedRAbundVector
         
         vector< vector<int> > dataSet(numRows, vector<int>(numColumns, 0));
         
+        vector<string> names;
+        
         for (int i = 0; i < lookup.size(); i++) {
             string sharedGroupName = lookup[i]->getGroup();
+            names.push_back(sharedGroupName);
             string treatmentName = designMap.getGroup(sharedGroupName);
             
             int j = 0;
@@ -382,9 +385,17 @@ void ClassifySharedCommand::processSharedAndDesignData(vector<SharedRAbundVector
         variables["[distance]"] = lookup[0]->getLabel();
         string filename = getOutputFileName("summary", variables);
         outputNames.push_back(filename); outputTypes["summary"].push_back(filename);
-        
         randomForest.calcForrestVariableImportance(filename);
         
+        //
+        map<string, string> variable; 
+        variable["[filename]"] = outputDir + m->getRootName(m->getSimpleName(sharedfile)) + "misclassifications.";
+        variable["[distance]"] = lookup[0]->getLabel();
+        string mc_filename = getOutputFileName("summary", variable);
+        outputNames.push_back(mc_filename); outputTypes["summary"].push_back(mc_filename);
+        randomForest.getMissclassifications(mc_filename, intToTreatmentMap, names);
+        //
+        
         m->mothurOutEndLine();
     }
     catch(exception& e) {