sffinfo bug with flow grams right index when clipQualRight=0

[mothur.git] / nseqs.h

diff --git a/nseqs.h b/nseqs.h
index 42a1f57b57f59737d91d99d9546659c5730a5a2b..c0f9549c87584bbf99e503388ff737c9f5343e8e 100644 (file)
--- a/nseqs.h
+++ b/nseqs.h
@@ -25,7 +25,26 @@ public:
                data[0] = (double)rank->getNumSeqs();
                return data;
        }
-       EstOutput getValues(vector<SharedRAbundVector*>) {return data;};
+       
+       EstOutput getValues(vector<SharedRAbundVector*> shared) { //return number of sequences in the sharedotus
+               
+               int numGroups = shared.size();
+               data.clear(); data.resize(numGroups,0);
+
+               for (int i = 0; i < shared[0]->getNumBins(); i++) {
+                       //get bin values and set sharedByAll 
+                       bool sharedByAll = true;
+                       for (int j = 0; j < numGroups; j++) {
+                               if (shared[j]->getAbundance(i) == 0) { sharedByAll = false; }
+                       }
+                       
+                       //they are shared
+                       if (sharedByAll == true) {  for (int j = 0; j < numGroups; j++) {  data[j] += shared[j]->getAbundance(i);  } }
+               }
+
+               return data;
+       }
+       string getCitation() { return "http://www.mothur.org/wiki/Nseqs"; }
 };
 
 /***********************************************************************/