5 * Created by westcott on 12/15/10.
6 * Copyright 2010 Schloss Lab. All rights reserved.
10 #include "structkulczynski.h"
12 /***********************************************************************/
13 EstOutput StructKulczynski::getValues(vector<SharedRAbundVector*> shared) {
21 for (int i = 0; i < shared[0]->getNumBins(); i++) {
23 int A = shared[0]->getAbundance(i);
24 int B = shared[1]->getAbundance(i);
31 data[0] = 1.0 - (0.5 * ((sumMin / sumA) + (sumMin / sumB)));
33 if (isnan(data[0]) || isinf(data[0])) { data[0] = 0; }
38 m->errorOut(e, "StructKulczynski", "getValues");
42 /***********************************************************************/