added heatmap command

[mothur.git] / summarysharedcommand.cpp

diff --git a/summarysharedcommand.cpp b/summarysharedcommand.cpp
index 44bd581929a28ae16890530d4963fba1e730d328..8b3a63f8688e557f290aa9db1acda63ae4de5426 100644 (file)
--- a/summarysharedcommand.cpp
+++ b/summarysharedcommand.cpp
@@ -20,8 +20,10 @@
 #include "sharedsorest.h"
 #include "sharedthetayc.h"
 #include "sharedthetan.h"
+#include "sharedkstest.h"
+#include "sharedbdiversity.h"
 #include "sharedochiai.h"
-#include "sharedanderberg.h"
+//#include "sharedanderberg.h"
 #include "sharedkulczynski.h"
 #include "sharedkulczynskicody.h"
 #include "sharedlennon.h"
@@ -64,12 +66,14 @@ SummarySharedCommand::SummarySharedCommand(){
                                        sumCalculators.push_back(new SharedThetaYC());
                                }else if (globaldata->Estimators[i] == "sharedthetan") { 
                                        sumCalculators.push_back(new SharedThetaN());
+                               }else if (globaldata->Estimators[i] == "sharedkstest") { 
+                                       sumCalculators.push_back(new SharedKSTest());
                                }else if (globaldata->Estimators[i] == "sharednseqs") { 
                                        sumCalculators.push_back(new SharedNSeqs());
                                }else if (globaldata->Estimators[i] == "sharedochiai") { 
                                        sumCalculators.push_back(new SharedOchiai());
-                               }else if (globaldata->Estimators[i] == "sharedanderberg") { 
-                                       sumCalculators.push_back(new SharedAnderberg());
+                               //}else if (globaldata->Estimators[i] == "sharedanderberg") { 
+                                       //sumCalculators.push_back(new SharedAnderberg());
                                }else if (globaldata->Estimators[i] == "sharedkulczynski") { 
                                        sumCalculators.push_back(new SharedKulczynski());
                                }else if (globaldata->Estimators[i] == "sharedkulczynskicody") { 
@@ -81,9 +85,12 @@ SummarySharedCommand::SummarySharedCommand(){
                                }else if (globaldata->Estimators[i] == "sharedbraycurtis") { 
                                        sumCalculators.push_back(new SharedBrayCurtis());
                                }
+                               else if (globaldata->Estimators[i] == "sharedbdiversity") { 
+                                       sumCalculators.push_back(new SharedBDiversity());
+                               }
+                               
                        }
                }
-               
                //reset calc for next command
                globaldata->setCalc("");
 
@@ -152,6 +159,7 @@ int SummarySharedCommand::execute(){
                                int n = 1; 
                                for (int k = 0; k < (lookup.size() - 1); k++) { // pass cdd each set of groups to commpare
                                        for (int l = n; l < lookup.size(); l++) {
+                                               outputFileHandle << order->getLabel() << '\t' << (lookup[k]->getGroup() + lookup[l]->getGroup()) << '\t' << '\t'; //print out label and group
                                                outputFileHandle << order->getLabel() << '\t';
                                                
                                                //sort groups to be alphanumeric