2 * sharedmorisitahorn.cpp
5 * Created by Sarah Westcott on 3/24/09.
6 * Copyright 2009 Schloss Lab UMASS Amherst. All rights reserved.
10 #include "sharedmorisitahorn.h"
12 /***********************************************************************/
13 EstOutput MorHorn::getValues(vector<SharedRAbundVector*> shared) {
17 double Atotal, Btotal, tempA, tempB;
18 Atotal = 0; Btotal = 0;
19 double morhorn, sumSharedA, sumSharedB, a, b, d;
20 morhorn = 0.0; sumSharedA = 0.0; sumSharedB = 0.0; a = 0.0; b = 0.0; d = 0.0;
22 //get the total values we need to calculate the theta denominator sums
23 for (int i = 0; i < shared[0]->size(); i++) {
24 //store in temps to avoid multiple repetitive function calls
25 Atotal += shared[0]->getAbundance(i);
26 Btotal += shared[1]->getAbundance(i);
29 //calculate the denominator sums
30 for (int j = 0; j < shared[0]->size(); j++) {
31 //store in temps to avoid multiple repetitive function calls
32 tempA = shared[0]->getAbundance(j);
33 tempB = shared[1]->getAbundance(j);
34 float relA = tempA / Atotal;
35 float relB = tempB / Btotal;
42 morhorn = (2 * d) / (a + b);
44 if (isnan(morhorn) || isinf(morhorn)) { morhorn = 0; }
51 m->errorOut(e, "MorHorn", "getValues");
56 /***********************************************************************/