X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=unifracunweightedcommand.h;h=fd39ae43852bee7ad970d77035fa839c50f62422;hb=ffc44592ff7ae94f14f9e21f87198e33d323cd1d;hp=0a925385e3d6ef3f8621029038ebaafb1bfbec49;hpb=74c78f9abd9e733f0c2f812efec97a76632fcbf8;p=mothur.git diff --git a/unifracunweightedcommand.h b/unifracunweightedcommand.h index 0a92538..fd39ae4 100644 --- a/unifracunweightedcommand.h +++ b/unifracunweightedcommand.h @@ -15,46 +15,54 @@ #include "treemap.h" #include "sharedutilities.h" #include "fileoutput.h" - - -class GlobalData; +#include "readtree.h" class UnifracUnweightedCommand : public Command { public: - UnifracUnweightedCommand(); - ~UnifracUnweightedCommand() { delete unweighted; delete util; } - int execute(); + UnifracUnweightedCommand(string); + UnifracUnweightedCommand(); + ~UnifracUnweightedCommand() {} + + vector setParameters(); + string getCommandName() { return "unifrac.unweighted"; } + string getCommandCategory() { return "Hypothesis Testing"; } + string getHelpString(); + string getCitation() { return "Lozupone C, Knight R (2005). UniFrac: a new phylogenetic method for comparing microbial communities. Appl Environ Microbiol 71: 8228-35. \nhttp://www.mothur.org/wiki/Unifrac.unweighted"; } + string getDescription() { return "generic tests that describes whether two or more communities have the same structure"; } + + int execute(); + void help() { m->mothurOut(getHelpString()); } + private: - GlobalData* globaldata; - SharedUtil* util; FileOutput* output; vector T; //user trees TreeMap* tmap; - Unweighted* unweighted; string sumFile, allGroups; vector groupComb; // AB. AC, BC... - int iters, numGroups, numComp, counter; - EstOutput userData; //unweighted score info for user tree - EstOutput randomData; //unweighted score info for random trees + int iters, numGroups, numComp, counter, processors, subsampleSize, subsampleIters; vector< vector > utreeScores; //scores for users trees for each comb. vector< vector > UWScoreSig; //tree score signifigance when compared to random trees - percentage of random trees with that score or higher. map validScores; //map contains scores from random vector< map > rscoreFreq; //map -vector entry for each combination. vector< map > rCumul; //map -vector entry for each combination. - map::iterator it2; - map::iterator it; + bool abort, phylip, random, includeRoot, consensus, subsample; + string groups, itersString, outputDir, outputForm, treefile, groupfile, namefile; + vector Groups, outputNames; //holds groups to be used + ofstream outSum, out; ifstream inFile; + int runRandomCalcs(Tree*, vector); void printUWSummaryFile(int); void printUnweightedFile(); - + void createPhylipFile(int); + vector buildTrees(vector< vector >&, int, TreeMap&); + int getConsensusTrees(vector< vector >&, int); + int getAverageSTDMatrices(vector< vector >&, int); }; - - -#endif \ No newline at end of file +#endif