hclustercommand.h

   1 #ifndef HCLUSTERCOMMAND_H
   2 #define HCLUSTERCOMMAND_H
   3
   4 /*
   5  *  hclustercommand.h
   6  *  Mothur
   7  *
   8  *  Created by westcott on 10/13/09.
   9  *  Copyright 2009 Schloss Lab. All rights reserved.
  10  *
  11  */
  12
  13 #include "command.hpp"
  14 #include "globaldata.hpp"
  15 #include "hcluster.h"
  16 #include "rabundvector.hpp"
  17 #include "sabundvector.hpp"
  18 #include "listvector.hpp"
  19 #include "readcluster.h"
  20
  21 /******************************************************************/
  22 //This command is an implementation of the HCluster algorythmn described in
  23 //ESPRIT: estimating species richness using large collections of 16S rRNA pyrosequences by
  24 //Yijun Sun1,2,*, Yunpeng Cai2, Li Liu1, Fahong Yu1, Michael L. Farrell3, William McKendree3
  25 //and William Farmerie1 1
  26
  27 //Interdisciplinary Center for Biotechnology Research, 2Department of Electrical and Computer Engineering,
  28 //University of Florida, Gainesville, FL 32610-3622 and 3Materials Technology Directorate, Air Force Technical
  29 //Applications Center, 1030 S. Highway A1A, Patrick AFB, FL 32925-3002, USA
  30 //Received January 28, 2009; Revised April 14, 2009; Accepted April 15, 2009
  31 /************************************************************/
  32 class HClusterCommand : public Command {
  33
  34 public:
  35         HClusterCommand(string);
  36         ~HClusterCommand();
  37         int execute();
  38         void help();
  39
  40 private:
  41         GlobalData* globaldata;
  42         HCluster* cluster;
  43         ListVector* list;
  44         RAbundVector* rabund;
  45         RAbundVector oldRAbund;
  46         ListVector oldList;
  47         ReadCluster* read;
  48
  49         bool abort, sorted, print_start;
  50         string method, fileroot, tag, distfile, format, phylipfile, columnfile, namefile, sort, showabund, timing, outputDir;
  51         double cutoff;
  52         int precision, length;
  53         ofstream sabundFile, rabundFile, listFile;
  54         time_t start;
  55         unsigned long loops;
  56
  57         void printData(string label);
  58 };
  59
  60 /************************************************************/
  61
  62 #endif