added classify.shared command and random forest files. added count file to pcr.seqs...

[mothur.git] / mothurout.cpp

diff --git a/mothurout.cpp b/mothurout.cpp
index c5eb0dca828a47023973f1684986ba05d45bb0fb..7d40e80ead727205f4bcc00269f78bd7f81057ca 100644 (file)
--- a/mothurout.cpp
+++ b/mothurout.cpp
@@ -2912,8 +2912,31 @@ string MothurOut::removeQuotes(string tax) {
        }
 }
 /**************************************************************************************************/
-
-
+// function for calculating standard deviation
+double MothurOut::getStandardDeviation(vector<int>& featureVector){
+    try {
+        //finds sum
+        double average = 0; 
+        for (int i = 0; i < featureVector.size(); i++) { average += featureVector[i]; }
+        average /= (double) featureVector.size();
+        
+        //find standard deviation
+        double stdDev = 0;
+        for (int i = 0; i < featureVector.size(); i++) { //compute the difference of each dist from the mean, and square the result of each
+            stdDev += ((featureVector[i] - average) * (featureVector[i] - average));
+        }
+          
+        stdDev /= (double) featureVector.size(); 
+        stdDev = sqrt(stdDev);
+        
+        return stdDev;
+    }
+       catch(exception& e) {
+               errorOut(e, "MothurOut", "getStandardDeviation");
+               exit(1);
+       }
+}
+/**************************************************************************************************/