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(SharedRAbundVector* shared1, SharedRAbundVector* shared2) {
17 int Atotal, Btotal, tempA, tempB;
18 Atotal = 0; Btotal = 0;
19 float 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 < shared1->size(); i++) {
24 //store in temps to avoid multiple repetitive function calls
25 Atotal += shared1->getAbundance(i);
26 Btotal += shared2->getAbundance(i);
29 //calculate the theta denominator sums
30 for (int j = 0; j < shared1->size(); j++) {
31 //store in temps to avoid multiple repetitive function calls
32 tempA = shared1->getAbundance(j);
33 tempB = shared2->getAbundance(j);
40 a /= double(Atotal * Atotal);
41 b /= double(Btotal * Btotal);
42 d /= double(Atotal * Btotal);
44 morhorn = (2 * d) / (a + b);
46 if (isnan(morhorn) || isinf(morhorn)) { morhorn = 0; }
53 cout << "Standard Error: " << e.what() << " has occurred in the MorHorn class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
57 cout << "An unknown error has occurred in the MorHorn class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
62 /***********************************************************************/