working on pam

[mothur.git] / mgclustercommand.h

diff --git a/mgclustercommand.h b/mgclustercommand.h
index 889586ae305ac47bdb0f2d15880507aa68655dff..3865bb20811a0ac29e29e0fec02392268ea0a657 100644 (file)
--- a/mgclustercommand.h
+++ b/mgclustercommand.h
@@ -12,12 +12,12 @@
 
 #include "command.hpp"
 #include "readblast.h"
-#include "sparsematrix.hpp"
 #include "nameassignment.hpp"
 #include "cluster.hpp"
 #include "hcluster.h"
 #include "rabundvector.hpp"
 #include "sabundvector.hpp"
+#include "counttable.h"
 
 /**********************************************************************/
 
@@ -31,8 +31,9 @@ public:
        vector<string> setParameters();
        string getCommandName()                 { return "mgcluster";   }
        string getCommandCategory()             { return "Clustering";  }
-       string getOutputFileNameTag(string, string);
+       
        string getHelpString(); 
+    string getOutputPattern(string);   
        string getCitation() { return "Schloss PD, Handelsman J (2008). A statistical toolbox for metagenomics. BMC Bioinformatics 9: 34. \nhttp://www.mothur.org/wiki/Mgcluster"; }
        string getDescription()         { return "cluster your sequences into OTUs using a blast file"; }
 
@@ -46,22 +47,24 @@ private:
        Cluster* cluster;
        HCluster* hcluster;
        ListVector* list;
+    CountTable* ct;
        ListVector oldList;
+    RAbundVector rav;
        vector<seqDist> overlapMatrix;
        vector<string> outputNames;
        
-       string blastfile, method, namefile, overlapFile, distFile, outputDir;
+       string blastfile, method, namefile, countfile, overlapFile, distFile, outputDir;
        ofstream sabundFile, rabundFile, listFile;
        double cutoff;
-       float penalty;
+       float penalty, adjust;
        int precision, length, precisionLength;
-       bool abort, minWanted, hclusterWanted, merge, hard;
+       bool abort, minWanted, hclusterWanted, merge, hard, cutoffSet;
        
        void printData(ListVector*);
        ListVector* mergeOPFs(map<string, int>, float);
        void sortHclusterFiles(string, string);
        vector<seqDist> getSeqs(ifstream&);
-    RAbundVector createRabund(ListVector, map<string, int>);
+    void createRabund(CountTable*&, ListVector*&, RAbundVector*&);
 
 };