X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=classifyseqscommand.h;h=9c8dbcac9caf436fcf44b2c416a9e063509561aa;hb=ae9c5b48da9047548f86fc84897926d9f5102312;hp=ed2819bcc6c7a3adf335b8e7cd79ceb29e3ae01a;hpb=ef3f6d42fe720cd6d91419e5e32f8c04d8765010;p=mothur.git diff --git a/classifyseqscommand.h b/classifyseqscommand.h index ed2819b..9c8dbca 100644 --- a/classifyseqscommand.h +++ b/classifyseqscommand.h @@ -12,7 +12,6 @@ #include "mothur.h" #include "command.hpp" -#include "alignment.hpp" #include "classify.h" //KNN and Bayesian methods modeled from algorithms in @@ -28,30 +27,54 @@ class ClassifySeqsCommand : public Command { public: - ClassifySeqsCommand(string); + ClassifySeqsCommand(string); + ClassifySeqsCommand(); ~ClassifySeqsCommand(); + + vector setParameters(); + string getCommandName() { return "classify.seqs"; } + string getCommandCategory() { return "Phylotype Analysis"; } + string getHelpString(); + string getCitation() { return "Wang Q, Garrity GM, Tiedje JM, Cole JR (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol 73: 5261-7. [ for Bayesian classifier ] \nAltschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25: 3389-402. [ for BLAST ] \nDeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol 72: 5069-72. [ for kmer ] \nhttp://www.mothur.org/wiki/Classify.seqs"; } + string getDescription() { return "classify sequences"; } + int execute(); - void help(); + void help() { m->mothurOut(getHelpString()); } + + private: struct linePair { - int start; - int numSeqs; - linePair(int i, int j) : start(i), numSeqs(j) {} + unsigned long int start; + unsigned long int end; + linePair(unsigned long int i, unsigned long int j) : start(i), end(j) {} }; + vector processIDS; //processid vector lines; + vector fastaFileNames; + vector namefileNames; + vector groupfileNames; + vector outputNames; + map > nameMap; + map >::iterator itNames; Classify* classify; - string fastaFileName, templateFileName, distanceFileName, search, method, taxonomyFileName; - int processors, kmerSize, numWanted, cutoff; + string fastaFileName, templateFileName, distanceFileName, namefile, search, method, taxonomyFileName, outputDir, groupfile; + int processors, kmerSize, numWanted, cutoff, iters; float match, misMatch, gapOpen, gapExtend; - bool abort; + bool abort, probs; - int driver(linePair*, string, string); + int driver(linePair*, string, string, string); void appendTaxFiles(string, string); - void createProcesses(string, string); + int createProcesses(string, string, string); + string addUnclassifieds(string, int); + + int MPIReadNamesFile(string); + #ifdef USE_MPI + int driverMPI(int, int, MPI_File&, MPI_File&, MPI_File&, vector&); + #endif }; #endif