+ double morhorn, sumSharedA, sumSharedB, a, b, d;
morhorn = 0.0; sumSharedA = 0.0; sumSharedB = 0.0; a = 0.0; b = 0.0; d = 0.0;
//get the total values we need to calculate the theta denominator sums
morhorn = 0.0; sumSharedA = 0.0; sumSharedB = 0.0; a = 0.0; b = 0.0; d = 0.0;
//get the total values we need to calculate the theta denominator sums
- for (int i = 0; i < shared1->size(); i++) {
+ for (int i = 0; i < shared[0]->size(); i++) {
//store in temps to avoid multiple repetitive function calls
//store in temps to avoid multiple repetitive function calls
- Atotal += shared1->getAbundance(i);
- Btotal += shared2->getAbundance(i);
+ Atotal += shared[0]->getAbundance(i);
+ Btotal += shared[1]->getAbundance(i);
}
}
- //calculate the theta denominator sums
- for (int j = 0; j < shared1->size(); j++) {
+ //calculate the denominator sums
+ for (int j = 0; j < shared[0]->size(); j++) {
//store in temps to avoid multiple repetitive function calls
//store in temps to avoid multiple repetitive function calls
- tempA = shared1->getAbundance(j);
- tempB = shared2->getAbundance(j);
+ tempA = shared[0]->getAbundance(j);
+ tempB = shared[1]->getAbundance(j);
+ float relA = tempA / Atotal;
+ float relB = tempB / Btotal;
- a += tempA * tempA;
- b += tempB * tempB;
- d += tempA * tempB;
+ a += relA * relA;
+ b += relB * relB;
+ d += relA * relB;
}
}
- a /= double(Atotal * Atotal);
- b /= double(Btotal * Btotal);
- d /= double(Atotal * Btotal);
-
morhorn = (2 * d) / (a + b);
morhorn = (2 * d) / (a + b);
-
+
if (isnan(morhorn) || isinf(morhorn)) { morhorn = 0; }
data[0] = morhorn;
if (isnan(morhorn) || isinf(morhorn)) { morhorn = 0; }
data[0] = morhorn;
@@ -50,13+48,9 @@ EstOutput MorHorn::getValues(SharedRAbundVector* shared1, SharedRAbundVector* sh
return data;
}
catch(exception& e) {
return data;
}
catch(exception& e) {
- cout << "Standard Error: " << e.what() << " has occurred in the MorHorn class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
+ m->errorOut(e, "MorHorn", "getValues");
exit(1);
}
exit(1);
}
- catch(...) {
- cout << "An unknown error has occurred in the MorHorn class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";