EstOutput Ace::getValues(SAbundVector* rank) {
try {
data.resize(3,0);
- // vector<double> aceData(3,0);
double ace, acelci, acehci;
int nrare = 0;
if(denom <= 0.0){ term1=0.0000; } else { term1 = (double)(srare * numsum)/(double)denom - 1.0; }
if(term1 >= 0.0){ gamace = term1; } else { gamace = 0.0; }
-
+
if(gamace >= 0.64){
gamace = gamace * (1 + (nrare * (1 - Cace) * numsum) / denom);
if(gamace<0){ gamace = 0; }
return data;
}
catch(exception& e) {
- cout << "Standard Error: " << e.what() << " has occurred in the Ace class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
+ errorOut(e, "Ace", "getValues");
exit(1);
}
- catch(...) {
- cout << "An unknown error has occurred in the Ace class function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
- exit(1);
- }
}
/***********************************************************************/