shannoneven.cpp

   1 /*
   2  *  shannoneven.cpp
   3  *  Mothur
   4  *
   5  *  Created by Pat Schloss on 8/21/10.
   6  *  Copyright 2010 Schloss Lab. All rights reserved.
   7  *
   8  */
   9
  10 #include "shannoneven.h"
  11 #include "shannon.h"
  12
  13 /***********************************************************************/
  14
  15 EstOutput ShannonEven::getValues(SAbundVector* rank){
  16         try {
  17                 //vector<double> simpsonData(3,0);
  18                 data.resize(1,0);
  19                 vector<double> shanData(3,0);
  20                 Shannon* shannon = new Shannon();
  21                 shanData = shannon->getValues(rank);
  22
  23                 long int sobs = rank->getNumBins();
  24                 if(sobs > 1){
  25                         data[0] = shanData[0] / log(sobs);
  26                 }
  27                 else{
  28                         data[0] = 1;
  29                 }
  30
  31                 delete shannon;
  32                 return data;
  33         }
  34         catch(exception& e) {
  35                 m->errorOut(e, "ShannonEven", "getValues");
  36                 exit(1);
  37         }
  38 }
  39
  40 /***********************************************************************/