X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=ace.cpp;h=92c023f663fa11e8e2ac9bc0041a15162a1abdb1;hp=54d0946eab9cf0b2230ef45ef8749b1feecb9da0;hb=d1c97b8c04bb75faca1e76ffad60b37a4d789d3d;hpb=e4827e0945cbda536064e5a345996b2a7dfcbb81

diff --git a/ace.cpp b/ace.cpp
index 54d0946..92c023f 100644
--- a/ace.cpp
+++ b/ace.cpp
@@ -16,12 +16,12 @@ EstOutput Ace::getValues(SAbundVector* rank) {
 		data.resize(3,0);
 		double ace, acelci, acehci;
 	
-		int nrare = 0;
-		int srare = 0;
-		int sabund = 0;
+		double nrare = 0;
+		double srare = 0;
+		double sabund = 0;
 	
 		double Cace, term1, gamace;
-		int numsum = 0;
+		double numsum = 0;
 	
 		double maxRank = (double)rank->getMaxRank();
 	
@@ -33,7 +33,7 @@ EstOutput Ace::getValues(SAbundVector* rank) {
 			}
 			else if(i>abund)	{sabund += rank->get(i);}
 		}
-		int sobs = srare + sabund;
+		double sobs = srare + sabund;
 	
 		if (nrare == 0){ Cace = 0.0000; }
 		else { Cace = 1.0000 -(double)rank->get(1)/(double)nrare; }
@@ -62,7 +62,8 @@ EstOutput Ace::getValues(SAbundVector* rank) {
 		I have also added the forumlae to calculate the 95% confidence intervals.
 		*/
 	
-		int j,D_s=0,nn=0,ww=0,Max_Index=rank->getMaxRank()+1;
+		double j,D_s=0,nn=0,ww=0;
+		int Max_Index=rank->getMaxRank()+1;
 		double pp, temp1, temp2;
 		vector<double> Part_N_Part_F(Max_Index+1,0.0);
     
@@ -140,13 +141,9 @@ EstOutput Ace::getValues(SAbundVector* rank) {
 		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";
+		m->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);
-	}	
 }
 
 /***********************************************************************/