5 * Created by westcott on 12/14/10.
6 * Copyright 2010 Schloss Lab. All rights reserved.
12 /***********************************************************************/
14 EstOutput Odum::getValues(vector<SharedRAbundVector*> shared) {
19 double sumDenom = 0.0;
21 for (int i = 0; i < shared[0]->getNumBins(); i++) {
23 int Aij = shared[0]->getAbundance(i);
24 int Bij = shared[1]->getAbundance(i);
26 sumNum += abs(Aij - Bij);
27 sumDenom += (Aij + Bij);
30 data[0] = sumNum / sumDenom;
32 if (isnan(data[0]) || isinf(data[0])) { data[0] = 0; }
37 m->errorOut(e, "Odum", "getValues");
42 /***********************************************************************/