validcalculator.cpp

   1 /*
   2  *  validcalculator.cpp
   3  *  Dotur
   4  *
   5  *  Created by Sarah Westcott on 1/5/09.
   6  *  Copyright 2009 Schloss Lab UMASS Amherst. All rights reserved.
   7  *
   8  */
   9
  10 #include "validcalculator.h"
  11
  12 /********************************************************************/
  13 ValidCalculators::ValidCalculators() {
  14         try {
  15                  initialSingle();
  16                  initialShared();
  17                  initialRarefaction();
  18                  initialSharedRarefact();
  19                  initialSummary();
  20                  initialSharedSummary();
  21         }
  22         catch(exception& e) {
  23                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function ValidCalculator. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
  24                 exit(1);
  25         }
  26         catch(...) {
  27                 cout << "An unknown error has occurred in the ValidCalculator class function ValidCalculator. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
  28                 exit(1);
  29         }
  30 }
  31
  32 /********************************************************************/
  33
  34 ValidCalculators::~ValidCalculators() {}
  35
  36 /********************************************************************/
  37
  38 bool ValidCalculators::isValidCalculator(string parameter, string calculator) {
  39         try {
  40                 //are you looking for a calculator for a single parameter
  41                 if (parameter == "single") {
  42                         //is it valid
  43                         if ((single.find(calculator)) != (single.end())) {
  44                                 return true;
  45                         }else {
  46                                 cout << calculator << " is not a valid estimator for the collect.single command and will be disregarded. Valid estimators are ";
  47                                 for (it = single.begin(); it != single.end(); it++) {
  48                                         cout << it->first << ", ";
  49                                 }
  50                                 cout << endl;
  51                                 return false; }
  52                 //are you looking for a calculator for a shared parameter
  53                 }else if (parameter == "shared") {
  54                         //is it valid
  55                         if ((shared.find(calculator)) != (shared.end())) {
  56                                 return true;
  57                         }else {
  58                                 cout << calculator << " is not a valid estimator for the collect.shared command and will be disregarded.  Valid estimators are ";
  59                                 for (it = shared.begin(); it != shared.end(); it++) {
  60                                         cout << it->first << ", ";
  61                                 }
  62                                 cout << endl;
  63                                 return false; }
  64                 //are you looking for a calculator for a rarefaction parameter
  65                 }else if (parameter == "rarefaction") {
  66                         //is it valid
  67                         if ((rarefaction.find(calculator)) != (rarefaction.end())) {
  68                                 return true;
  69                         }else {
  70                                 cout << calculator << " is not a valid estimator for the rarefaction.single command and will be disregarded. Valid estimators are ";
  71                                 for (it = rarefaction.begin(); it != rarefaction.end(); it++) {
  72                                         cout << it->first << ", ";
  73                                 }
  74                                 cout << endl;
  75                                 return false; }
  76                 //are you looking for a calculator for a summary parameter
  77                 }else if (parameter == "summary") {
  78                         //is it valid
  79                         if ((summary.find(calculator)) != (summary.end())) {
  80                                 return true;
  81                         }else {
  82                                 cout << calculator << " is not a valid estimator for the summary.shared command and will be disregarded. Valid estimators are ";
  83                                 for (it = summary.begin(); it != summary.end(); it++) {
  84                                         cout << it->first << ", ";
  85                                 }
  86                                 cout << endl;
  87                                 return false; }
  88                 //are you looking for a calculator for a sharedsummary parameter
  89                 }else if (parameter == "sharedsummary") {
  90                         //is it valid
  91                         if ((sharedsummary.find(calculator)) != (sharedsummary.end())) {
  92                                 return true;
  93                         }else {
  94                                 cout << calculator << " is not a valid estimator for the summary.shared command and will be disregarded. Valid estimators are ";
  95                                 for (it = sharedsummary.begin(); it != sharedsummary.end(); it++) {
  96                                         cout << it->first << ", ";
  97                                 }
  98                                 cout << endl;
  99                                 return false; }
 100                 }else if (parameter == "sharedrarefaction") {
 101                         //is it valid
 102                         if ((sharedrarefaction.find(calculator)) != (sharedrarefaction.end())) {
 103                                 return true;
 104                         }else {
 105                                 cout << calculator << " is not a valid estimator for the rarefaction.shared command and will be disregarded. Valid estimator is ";
 106                                 for (it = sharedrarefaction.begin(); it != sharedrarefaction.end(); it++) {
 107                                         cout << it->first << ", ";
 108                                 }
 109                                 cout << endl;
 110                                 return false; }
 111                 //not a valid paramter
 112                 }else { return false; }
 113
 114         }
 115         catch(exception& e) {
 116                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function isValidCalculator. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 117                 exit(1);
 118         }
 119         catch(...) {
 120                 cout << "An unknown error has occurred in the ValidCalculator class function isValidCalculator. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 121                 exit(1);
 122         }
 123 }
 124
 125 /********************************************************************/
 126 void ValidCalculators::initialSingle() {
 127         try {
 128
 129                 single["sobs"]          = "sobs";
 130                 single["chao"]          = "chao";
 131                 single["ace"]           = "ace";
 132                 single["jack"]          = "jack";
 133                 single["shannon"]       = "shannon";
 134                 single["npshannon"]     = "npshannon";
 135                 single["simpson"]       = "simpson";
 136                 single["bootstrap"]     = "bootstrap";
 137                 single["default"]       = "default";
 138                 single["nseqs"]         = "nseqs";
 139         }
 140         catch(exception& e) {
 141                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialSingle. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 142                 exit(1);
 143         }
 144         catch(...) {
 145                 cout << "An unknown error has occurred in the ValidCalculator class function initialSingle. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 146                 exit(1);
 147         }
 148 }
 149
 150 /********************************************************************/
 151 void ValidCalculators::initialShared() {
 152         try {
 153                 shared["sharedsobs"]                    = "sharedsobs";
 154                 shared["sharedchao"]                    = "sharedchao";
 155                 shared["sharedace"]                             = "sharedace";
 156                 shared["sharedjabund"]                  = "sharedjabund";
 157                 shared["sharedsorensonabund"]   = "sharedsorensonabund";
 158                 shared["sharedjclass"]                  = "sharedjclass";
 159                 shared["sharedsorclass"]                = "sharedsorclass";
 160                 shared["sharedjest"]                    = "sharedjest";
 161                 shared["sharedsorest"]                  = "sharedsorest";
 162                 shared["sharedthetayc"]                 = "sharedthetayc";
 163                 shared["sharedthetan"]                  = "sharedthetan";
 164                 shared["sharednseqs"]                   = "sharednseqs";
 165                 shared["sharedochiai"]                  = "sharedochiai";
 166                 shared["sharedanderberg"]               = "sharedanderberg";
 167                 shared["sharedkulczynski"]              = "sharedkulczynski";
 168                 shared["sharedkulczynskicody"]  = "sharedkulczynskicody";
 169                 shared["sharedlennon"]                  = "sharedlennon";
 170                 shared["sharedmorisitahorn"]    = "sharedmorisitahorn";
 171                 shared["sharedbraycurtis"]              = "sharedbraycurtis";
 172                 shared["default"]                   = "default";
 173         }
 174         catch(exception& e) {
 175                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialShared. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 176                 exit(1);
 177         }
 178         catch(...) {
 179                 cout << "An unknown error has occurred in the ValidCalculator class function initialShared. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 180                 exit(1);
 181         }
 182 }
 183
 184 /********************************************************************/
 185 void ValidCalculators::initialRarefaction() {
 186         try {
 187                 rarefaction["sobs"]                     = "sobs";
 188                 rarefaction["chao"]                     = "chao";
 189                 rarefaction["ace"]                      = "ace";
 190                 rarefaction["jack"]                     = "jack";
 191                 rarefaction["shannon"]          = "shannon";
 192                 rarefaction["npshannon"]        = "npshannon";
 193                 rarefaction["simpson"]          = "simpson";
 194                 rarefaction["bootstrap"]        = "bootstrap";
 195                 rarefaction["nseqs"]            = "nseqs";
 196                 rarefaction["default"]      = "default";
 197         }
 198         catch(exception& e) {
 199                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialRarefaction. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 200                 exit(1);
 201         }
 202         catch(...) {
 203                 cout << "An unknown error has occurred in the ValidCalculator class function initialRarefaction. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 204                 exit(1);
 205         }
 206 }
 207
 208 /********************************************************************/
 209
 210 void ValidCalculators::initialSummary() {
 211         try {
 212                 summary["sobs"]                 = "sobs";
 213                 summary["chao"]                 = "chao";
 214                 summary["ace"]                  = "ace";
 215                 summary["jack"]                 = "jack";
 216                 summary["shannon"]              = "shannon";
 217                 summary["npshannon"]    = "npshannon";
 218                 summary["simpson"]              = "simpson";
 219                 summary["bootstrap"]    = "bootstrap";
 220                 summary["nseqs"]                = "nseqs";
 221                 summary["default"]          = "default";
 222         }
 223         catch(exception& e) {
 224                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialSummary. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 225                 exit(1);
 226         }
 227         catch(...) {
 228                 cout << "An unknown error has occurred in the ValidCalculator class function initialSummary. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 229                 exit(1);
 230         }
 231 }
 232
 233 /********************************************************************/
 234 void ValidCalculators::initialSharedSummary() {
 235         try {
 236                 sharedsummary["sharedsobs"]                             = "sharedsobs";
 237                 sharedsummary["sharedchao"]                             = "sharedchao";
 238                 sharedsummary["sharedace"]                              = "sharedace";
 239                 sharedsummary["sharedjabund"]                   = "sharedjabund";
 240                 sharedsummary["sharedsorensonabund"]    = "sharedsorensonabund";
 241                 sharedsummary["sharedjclass"]                   = "sharedjclass";
 242                 sharedsummary["sharedsorclass"]                 = "sharedsorclass";
 243                 sharedsummary["sharedjest"]                             = "sharedjest";
 244                 sharedsummary["sharedsorest"]                   = "sharedsorest";
 245                 sharedsummary["sharedthetayc"]                  = "sharedthetayc";
 246                 sharedsummary["sharedthetan"]                   = "sharedthetan";
 247                 sharedsummary["sharednseqs"]                    = "sharednseqs";
 248                 sharedsummary["sharedochiai"]                   = "sharedochiai";
 249                 sharedsummary["sharedanderberg"]                = "sharedanderberg";
 250                 sharedsummary["sharedkulczynski"]               = "sharedkulczynski";
 251                 sharedsummary["sharedkulczynskicody"]   = "sharedkulczynskicody";
 252                 sharedsummary["sharedlennon"]                   = "sharedlennon";
 253                 sharedsummary["sharedmorisitahorn"]             = "sharedmorisitahorn";
 254                 sharedsummary["sharedbraycurtis"]               = "sharedbraycurtis";
 255                 sharedsummary["default"]                                = "default";
 256         }
 257         catch(exception& e) {
 258                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialSharedSummary. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 259                 exit(1);
 260         }
 261         catch(...) {
 262                 cout << "An unknown error has occurred in the ValidCalculator class function initialSharedSummary. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 263                 exit(1);
 264         }
 265 }
 266
 267
 268 /********************************************************************/
 269
 270 void ValidCalculators::initialSharedRarefact() {
 271         try {
 272                 sharedrarefaction["sharedobserved"]     = "sharedobserved";
 273                 sharedrarefaction["sharednseqs"]        = "sharednseqs";
 274                 sharedrarefaction["default"]        = "default";
 275         }
 276         catch(exception& e) {
 277                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialSharedRarefact. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 278                 exit(1);
 279         }
 280         catch(...) {
 281                 cout << "An unknown error has occurred in the ValidCalculator class function initialSharedRarefact. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 282                 exit(1);
 283         }
 284 }
 285
 286 /********************************************************************/
 287
 288
 289