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-2*sTotal/(float)(countA+countB);
30 m->errorOut(e, "Whittaker", "getValues");
35 /***********************************************************************/