modified some calculators and other stuff - pds

[mothur.git] / rarefactcommand.cpp

diff --git a/rarefactcommand.cpp b/rarefactcommand.cpp
index a83a750c5944aecfb3e740e5caca8bd231f11148..6e4197f7f8b4a9029709a9d4ab31f51c246d840b 100644 (file)
--- a/rarefactcommand.cpp
+++ b/rarefactcommand.cpp
@@ -17,6 +17,7 @@
 #include "npshannon.h"
 #include "shannon.h"
 #include "jackknife.h"
+#include "coverage.h"
 
 //**********************************************************************************************************************
 
@@ -50,6 +51,8 @@ RareFactCommand::RareFactCommand(){
                                        rDisplays.push_back(new RareDisplay(new Simpson(), new ThreeColumnFile(fileNameRoot+"r_simpson")));
                                }else if (globaldata->Estimators[i] == "bootstrap") { 
                                        rDisplays.push_back(new RareDisplay(new Bootstrap(), new ThreeColumnFile(fileNameRoot+"r_bootstrap")));
+                               }else if (globaldata->Estimators[i] == "coverage") { 
+                                       rDisplays.push_back(new RareDisplay(new Coverage(), new ThreeColumnFile(fileNameRoot+"r_coverage")));
                                }else if (globaldata->Estimators[i] == "nseqs") { 
                                        rDisplays.push_back(new RareDisplay(new NSeqs(), new ThreeColumnFile(fileNameRoot+"r_nseqs")));
                                }
@@ -89,7 +92,7 @@ int RareFactCommand::execute(){
                //if the users entered no valid calculators don't execute command
                if (rDisplays.size() == 0) { return 0; }
 
-               read = new ReadPhilFile(globaldata->inputFileName);     
+               read = new ReadOTUFile(globaldata->inputFileName);      
                read->read(&*globaldata); 
 
                order = globaldata->gorder;