X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=getoturepcommand.h;h=ca3439d809f96f6eff93dcf3948abd1a8274a7f5;hp=a48cbf1d7fd195264dacf9c43ee12b27867a8ff4;hb=615301e57c25e241356a9c2380648d117709458d;hpb=f37a59b81c339f574d97042ff6e364146feb457a diff --git a/getoturepcommand.h b/getoturepcommand.h index a48cbf1..ca3439d 100644 --- a/getoturepcommand.h +++ b/getoturepcommand.h @@ -12,52 +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 "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 { + 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(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; - SparseMatrix* matrix; ListVector* list; - ReadOTUFile* read; - InputData* input; - FastaMap* fasta; GroupMap* groupMap; - string filename, fastafile, listfile, namesfile, groupfile, label; + 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; - bool groupError; - - bool abort, allLines; + 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 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*); - 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 - + 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