X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=getoturepcommand.h;h=ca3439d809f96f6eff93dcf3948abd1a8274a7f5;hp=4e2b3a0e685aec843624613df7bda5c7e184b4c1;hb=615301e57c25e241356a9c2380648d117709458d;hpb=c5c7502f435e1413c19e373dab1dfebcaa67588d diff --git a/getoturepcommand.h b/getoturepcommand.h index 4e2b3a0..ca3439d 100644 --- a/getoturepcommand.h +++ b/getoturepcommand.h @@ -12,46 +12,79 @@ /* 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 "sparsematrix.hpp" #include "listvector.hpp" #include "inputdata.h" -#include "readmatrix.hpp" #include "fastamap.h" +#include "groupmap.h" +#include "readmatrix.hpp" +#include "formatmatrix.h" +#include "counttable.h" +typedef map SeqMap; -class GlobalData; - -typedef list::iterator MatData; +struct repStruct { + string name; + int bin; + int size; + string group; + + repStruct(){} + repStruct(string n, int b, int s, string g) : name(n), bin(b), size(s), group(g) {} + ~repStruct() {} +}; class GetOTURepCommand : public Command { - + public: - GetOTURepCommand(); - ~GetOTURepCommand(); - int execute(); + GetOTURepCommand(string); + 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; - SparseMatrix* matrix; ListVector* list; - ListVector* listOfNames; - ReadMatrix* read; - InputData* input; - FastaMap* fasta; - string filename, fastafile, namesfile; + GroupMap* groupMap; + ReadMatrix* readMatrix; + FormatMatrix* formatMatrix; + NameAssignment* nameMap; + CountTable ct; + string filename, fastafile, listfile, namefile, groupfile, label, sorted, phylipfile, countfile, columnfile, distFile, format, outputDir, groups, method; ofstream out; - ifstream in, inNames; - - + ifstream in, inNames, inRow; + bool abort, allLines, groupError, large, weighted, hasGroups; + set labels; //holds labels to be used map nameToIndex; //maps sequence name to index in sparsematrix - map::iterator it; - map::iterator it2; - map::iterator it3; - - void readNamesFile(); - string FindRep(int); // returns name of "representative" sequence of given bin. + 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(FastaMap*&); + void readNamesFile(bool); + int process(ListVector*); + SeqMap getMap(int); + 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