X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=mgclustercommand.h;h=3865bb20811a0ac29e29e0fec02392268ea0a657;hp=d5c784e27a28bb27e6b7cbdadbb8b94e25b199b5;hb=d1c97b8c04bb75faca1e76ffad60b37a4d789d3d;hpb=e150b0b0664caec517485ee6d69dcdade6dcae77

diff --git a/mgclustercommand.h b/mgclustercommand.h
index d5c784e..3865bb2 100644
--- 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,9 +31,12 @@ public:
 	vector<string> setParameters();
 	string getCommandName()			{ return "mgcluster";	}
 	string getCommandCategory()		{ return "Clustering";	}
-	string getHelpString();	
-	string getCitation() { return "http://www.mothur.org/wiki/Mgcluster"; }
 	
+	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()); }	
 	
@@ -44,21 +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&);
+    void createRabund(CountTable*&, ListVector*&, RAbundVector*&);
 
 };