X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=qstat.cpp;h=e07cdb3fb860b96a17ba2972474f0abb8bf05d7e;hp=ed4471a8d18b983294bd3824832d093316fe40fa;hb=cf9987b67aa49777a4c91c2d21f96e58bf17aa82;hpb=c196b6b4768ccb84955d773ff0f22e4994d1ba7b diff --git a/qstat.cpp b/qstat.cpp index ed4471a..e07cdb3 100644 --- a/qstat.cpp +++ b/qstat.cpp @@ -31,8 +31,8 @@ EstOutput QStat::getValues(SAbundVector* rank){ int r3 = -1; int r1Ind = 0; int r3Ind = 0; - int sumSpec = 0; - int iqSum = 0; + double sumSpec = 0; + double iqSum = 0; for(int i = 1; i < rank->size(); i++) { if(r1 != -1 && r3 != -1) i = rank->size(); @@ -52,7 +52,6 @@ EstOutput QStat::getValues(SAbundVector* rank){ } double qstat = (.5*r1 + iqSum + .5*r3)/log((double)r3Ind/r1Ind); - //cout << "QStat:\nQStatistic = " << qstat << "\n\n"; data[0] = qstat; if (isnan(data[0]) || isinf(data[0])) { data[0] = 0; } @@ -60,13 +59,9 @@ EstOutput QStat::getValues(SAbundVector* rank){ return data; } catch(exception& e) { - cout << "Standard Error: " << e.what() << " has occurred in the QStat class Function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; + m->errorOut(e, "QStat", "getValues"); exit(1); } - catch(...) { - cout << "An unknown error has occurred in the QStat class function getValues. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n"; - exit(1); - } } /***********************************************************************/