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 SharedMorHorn::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);
36 if ((tempA != 0) && (tempB != 0)) {
37 if (Atotal != 0) { sumSharedA = (tempA / (float)Atotal); }
38 if (Btotal != 0) { sumSharedB = (tempB / (float)Btotal); }
40 a += sumSharedA * sumSharedA;
41 b += sumSharedB * sumSharedB;
42 d += sumSharedA * sumSharedB;
46 morhorn = (2 * d) / (float) (a + b);
48 if (isnan(morhorn) || isinf(morhorn)) { morhorn = 0; }
55 cout << "Standard Error: " << e.what() << " has occurred in the SharedMorHorn class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
59 cout << "An unknown error has occurred in the SharedMorHorn class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
64 /***********************************************************************/