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["default"]                   = "default";
 168         }
 169         catch(exception& e) {
 170                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialShared. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 171                 exit(1);
 172         }
 173         catch(...) {
 174                 cout << "An unknown error has occurred in the ValidCalculator class function initialShared. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 175                 exit(1);
 176         }
 177 }
 178
 179 /********************************************************************/
 180 void ValidCalculators::initialRarefaction() {
 181         try {
 182                 rarefaction["sobs"]                     = "sobs";
 183                 rarefaction["chao"]                     = "chao";
 184                 rarefaction["ace"]                      = "ace";
 185                 rarefaction["jack"]                     = "jack";
 186                 rarefaction["shannon"]          = "shannon";
 187                 rarefaction["npshannon"]        = "npshannon";
 188                 rarefaction["simpson"]          = "simpson";
 189                 rarefaction["bootstrap"]        = "bootstrap";
 190                 rarefaction["nseqs"]            = "nseqs";
 191                 rarefaction["default"]      = "default";
 192         }
 193         catch(exception& e) {
 194                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialRarefaction. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 195                 exit(1);
 196         }
 197         catch(...) {
 198                 cout << "An unknown error has occurred in the ValidCalculator class function initialRarefaction. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 199                 exit(1);
 200         }
 201 }
 202
 203 /********************************************************************/
 204
 205 void ValidCalculators::initialSummary() {
 206         try {
 207                 summary["sobs"]                 = "sobs";
 208                 summary["chao"]                 = "chao";
 209                 summary["ace"]                  = "ace";
 210                 summary["jack"]                 = "jack";
 211                 summary["shannon"]              = "shannon";
 212                 summary["npshannon"]    = "npshannon";
 213                 summary["simpson"]              = "simpson";
 214                 summary["bootstrap"]    = "bootstrap";
 215                 summary["nseqs"]                = "nseqs";
 216                 summary["default"]          = "default";
 217         }
 218         catch(exception& e) {
 219                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialSummary. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 220                 exit(1);
 221         }
 222         catch(...) {
 223                 cout << "An unknown error has occurred in the ValidCalculator class function initialSummary. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 224                 exit(1);
 225         }
 226 }
 227
 228 /********************************************************************/
 229 void ValidCalculators::initialSharedSummary() {
 230         try {
 231                 sharedsummary["sharedsobs"]                             = "sharedsobs";
 232                 sharedsummary["sharedchao"]                             = "sharedchao";
 233                 sharedsummary["sharedace"]                              = "sharedace";
 234                 sharedsummary["sharedjabund"]                   = "sharedjabund";
 235                 sharedsummary["sharedsorensonabund"]    = "sharedsorensonabund";
 236                 sharedsummary["sharedjclass"]                   = "sharedjclass";
 237                 sharedsummary["sharedsorclass"]                 = "sharedsorclass";
 238                 sharedsummary["sharedjest"]                             = "sharedjest";
 239                 sharedsummary["sharedsorest"]                   = "sharedsorest";
 240                 sharedsummary["sharedthetayc"]                  = "sharedthetayc";
 241                 sharedsummary["sharedthetan"]                   = "sharedthetan";
 242                 sharedsummary["sharednseqs"]                    = "sharednseqs";
 243                 sharedsummary["sharedochiai"]                   = "sharedochiai";
 244                 sharedsummary["sharedanderberg"]                = "sharedanderberg";
 245                 sharedsummary["default"]                                = "default";
 246         }
 247         catch(exception& e) {
 248                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialSharedSummary. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 249                 exit(1);
 250         }
 251         catch(...) {
 252                 cout << "An unknown error has occurred in the ValidCalculator class function initialSharedSummary. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 253                 exit(1);
 254         }
 255 }
 256
 257
 258 /********************************************************************/
 259
 260 void ValidCalculators::initialSharedRarefact() {
 261         try {
 262                 sharedrarefaction["sharedobserved"]     = "sharedobserved";
 263                 sharedrarefaction["sharednseqs"]        = "sharednseqs";
 264                 sharedrarefaction["default"]        = "default";
 265         }
 266         catch(exception& e) {
 267                 cout << "Standard Error: " << e.what() << " has occurred in the ValidCalculator class Function initialSharedRarefact. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 268                 exit(1);
 269         }
 270         catch(...) {
 271                 cout << "An unknown error has occurred in the ValidCalculator class function initialSharedRarefact. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
 272                 exit(1);
 273         }
 274 }
 275
 276 /********************************************************************/
 277
 278
 279