5 * Created by Pat Schloss on 4/23/09.
6 * Copyright 2009 Patrick D. Schloss. All rights reserved.
10 #include "whittaker.h"
12 /***********************************************************************/
14 EstOutput Whittaker::getValues(vector<SharedRAbundVector*> shared){
20 int sTotal = shared[0]->getNumBins();
21 for(int i=0;i<sTotal;i++){
22 if(shared[0]->getAbundance(i) != 0){ countA++; }
23 if(shared[1]->getAbundance(i) != 0){ countB++; }
26 data[0] = 2*sTotal/(float)(countA+countB)-1;
30 cout << "Standard Error: " << e.what() << " has occurred in the Whittaker class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
34 cout << "An unknown error has occurred in the Whittaker class function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
39 /***********************************************************************/