/***********************************************************************/
-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;
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
+/***********************************************************************/