X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=mgclustercommand.h;h=3865bb20811a0ac29e29e0fec02392268ea0a657;hp=a1dac221b891a8551a36b7493d5bb7a603204cf5;hb=cf9987b67aa49777a4c91c2d21f96e58bf17aa82;hpb=9c23307c583d4e8595f75278c13e708788f2f058

diff --git a/mgclustercommand.h b/mgclustercommand.h
index a1dac22..3865bb2 100644
--- a/mgclustercommand.h
+++ b/mgclustercommand.h
@@ -12,11 +12,12 @@
 
 #include "command.hpp"
 #include "readblast.h"
-#include "sparsematrix.hpp"
 #include "nameassignment.hpp"
-#include "globaldata.hpp"
 #include "cluster.hpp"
 #include "hcluster.h"
+#include "rabundvector.hpp"
+#include "sabundvector.hpp"
+#include "counttable.h"
 
 /**********************************************************************/
 
@@ -24,31 +25,46 @@ class MGClusterCommand : public Command {
 
 public:
 	MGClusterCommand(string);
-	~MGClusterCommand();
-	int execute();
-	void help();
+	MGClusterCommand();
+	~MGClusterCommand(){}
+	
+	vector<string> setParameters();
+	string getCommandName()			{ return "mgcluster";	}
+	string getCommandCategory()		{ return "Clustering";	}
+	
+	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"; }
+
+	int execute(); 
+	void help() { m->mothurOut(getHelpString()); }	
+	
 	
 private:
-	GlobalData* globaldata;
 	ReadBlast* read;
 	NameAssignment* nameMap;
 	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;
+	bool abort, minWanted, hclusterWanted, merge, hard, cutoffSet;
 	
 	void printData(ListVector*);
 	ListVector* mergeOPFs(map<string, int>, float);
 	void sortHclusterFiles(string, string);
 	vector<seqDist> getSeqs(ifstream&);
+    void createRabund(CountTable*&, ListVector*&, RAbundVector*&);
 
 };