X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=getoturepcommand.h;h=ca3439d809f96f6eff93dcf3948abd1a8274a7f5;hp=7890d93db045be399a2854d5d17b314fb748358d;hb=615301e57c25e241356a9c2380648d117709458d;hpb=74844a60d80c6dd06e3fb02ee9b928424f9019b0 diff --git a/getoturepcommand.h b/getoturepcommand.h index 7890d93..ca3439d 100644 --- a/getoturepcommand.h +++ b/getoturepcommand.h @@ -12,16 +12,14 @@ /* The get.oturep command outputs a .fastarep file for each distance you specify, selecting one OTU representative for each bin. */ #include "command.hpp" -#include "globaldata.hpp" #include "listvector.hpp" #include "inputdata.h" -#include "readotu.h" #include "fastamap.h" #include "groupmap.h" #include "readmatrix.hpp" #include "formatmatrix.h" +#include "counttable.h" -typedef list::iterator MatData; typedef map SeqMap; struct repStruct { @@ -39,40 +37,54 @@ class GetOTURepCommand : public Command { public: GetOTURepCommand(string); - ~GetOTURepCommand(); - int execute(); - void help(); + GetOTURepCommand(); + ~GetOTURepCommand(){} + + vector setParameters(); + string getCommandName() { return "get.oturep"; } + string getCommandCategory() { return "OTU-Based Approaches"; } + + string getHelpString(); + string getOutputPattern(string); + string getCitation() { return "http://www.mothur.org/wiki/Get.oturep"; } + string getDescription() { return "gets a representative sequence for each OTU"; } + + int execute(); + void help() { m->mothurOut(getHelpString()); } + + private: - GlobalData* globaldata; ListVector* list; - ReadOTUFile* read; - InputData* input; - FastaMap* fasta; GroupMap* groupMap; ReadMatrix* readMatrix; FormatMatrix* formatMatrix; NameAssignment* nameMap; - string filename, fastafile, listfile, namefile, groupfile, label, sorted, phylipfile, columnfile, distFile, format, outputDir; + CountTable ct; + string filename, fastafile, listfile, namefile, groupfile, label, sorted, phylipfile, countfile, columnfile, distFile, format, outputDir, groups, method; ofstream out; ifstream in, inNames, inRow; - bool abort, allLines, groupError, large; + bool abort, allLines, groupError, large, weighted, hasGroups; set labels; //holds labels to be used map nameToIndex; //maps sequence name to index in sparsematrix - vector outputNames; + map nameFileMap; + vector outputNames, Groups; + map outputNameFiles; float cutoff; int precision; vector seqVec; // contains maps with sequence index and distance // for all distances related to a certain sequence vector rowPositions; - void readNamesFile(); + void readNamesFile(FastaMap*&); + void readNamesFile(bool); int process(ListVector*); SeqMap getMap(int); - string findRep(int, string&, ListVector*, int&); // returns the name of the "representative" sequence of given bin, - // fills a string containing the groups in that bin if a groupfile is given, - // and returns the number of sequences in the given bin - + string findRep(vector, string); // returns the name of the "representative" sequence of given bin or subset of a bin, for groups + string findRepAbund(vector, string); + int processNames(string, string); + int processFastaNames(string, string, FastaMap*&); + int readDist(); }; #endif