#include "qstat.h"
#include "bergerparker.h"
#include "bstick.h"
+#include "goodscoverage.h"
+#include "coverage.h"
//**********************************************************************************************************************
for (i=0; i<globaldata->Estimators.size(); i++) {
if (validCalculator->isValidCalculator("summary", globaldata->Estimators[i]) == true) {
-
if(globaldata->Estimators[i] == "sobs"){
sumCalculators.push_back(new Sobs());
}else if(globaldata->Estimators[i] == "chao"){
sumCalculators.push_back(new Chao1());
- }else if(globaldata->Estimators[i] == "geom"){
+ }else if(globaldata->Estimators[i] == "coverage"){
+ sumCalculators.push_back(new Coverage());
+ }else if(globaldata->Estimators[i] == "geometric"){
sumCalculators.push_back(new Geom());
- }else if(globaldata->Estimators[i] == "logsd"){
+ }else if(globaldata->Estimators[i] == "logseries"){
sumCalculators.push_back(new LogSD());
}else if(globaldata->Estimators[i] == "qstat"){
sumCalculators.push_back(new QStat());
sumCalculators.push_back(new Bootstrap());
}else if (globaldata->Estimators[i] == "nseqs") {
sumCalculators.push_back(new NSeqs());
+ }else if (globaldata->Estimators[i] == "goodscoverage") {
+ sumCalculators.push_back(new GoodsCoverage());
}
}
}
openOutputFile(outputFileName, outputFileHandle);
outputFileHandle << "label";
- read = new ReadPhilFile(globaldata->inputFileName);
+ read = new ReadOTUFile(globaldata->inputFileName);
read->read(&*globaldata);
for(int i=0;i<sumCalculators.size();i++){
for(int i=0;i<sumCalculators.size();i++){
vector<double> data = sumCalculators[i]->getValues(sabund);
outputFileHandle << '\t';
- //sumCalculators[i]->print(outputFileHandle);
+ sumCalculators[i]->print(outputFileHandle);
}
outputFileHandle << endl;