]> git.donarmstrong.com Git - mothur.git/blobdiff - whittaker.cpp
changing command name classify.shared to classifyrf.shared
[mothur.git] / whittaker.cpp
index 95f1c9c5fcb789a19a7496a83177ab639d0067e6..5ab1e5dd88eedb8e0fb0dbe52d0b42e0c61d3c44 100644 (file)
 
 /***********************************************************************/
 
-EstOutput Whittaker::getValues(SharedRAbundVector* shared1, SharedRAbundVector* shared2){
+EstOutput Whittaker::getValues(vector<SharedRAbundVector*> shared){
        try{
                data.resize(1);
 
                int countA = 0;
                int countB = 0;
-               int sTotal = shared1->getNumBins();
+               int sTotal = shared[0]->getNumBins();
                for(int i=0;i<sTotal;i++){
-                       if(shared1->getAbundance(i) != 0){      countA++;       }
-                       if(shared2->getAbundance(i) != 0){      countB++;       }               
+                       if(shared[0]->getAbundance(i) != 0){    countA++;       }
+                       if(shared[1]->getAbundance(i) != 0){    countB++;       }               
                }
                
-               data[0] = 2*sTotal/(float)(countA+countB)-1;
+               data[0] = 2-2*sTotal/(float)(countA+countB);
                return data;
        }
        catch(exception& e) {
-               cout << "Standard Error: " << e.what() << " has occurred in the Whittaker class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
+               m->errorOut(e, "Whittaker", "getValues");
                exit(1);
        }
-       catch(...) {
-               cout << "An unknown error has occurred in the Whittaker class function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
-               exit(1);
-       }       
 }
 
-/***********************************************************************/
\ No newline at end of file
+/***********************************************************************/