X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=unifracweightedcommand.h;h=7317c8a43092a5712a06dda00589d1eb4528c649;hb=a98eb683e17d8e49583bf2d215ab7562a4cdca75;hp=ea2c41814ac1940a61bd8c48995ecebcff30f324;hpb=f18d75f9c996cb619a5051a7664620c0c5ae05fb;p=mothur.git diff --git a/unifracweightedcommand.h b/unifracweightedcommand.h index ea2c418..7317c8a 100644 --- a/unifracweightedcommand.h +++ b/unifracweightedcommand.h @@ -14,28 +14,38 @@ #include "weighted.h" #include "treemap.h" #include "progress.hpp" +#include "sharedutilities.h" +#include "fileoutput.h" -using namespace std; class GlobalData; class UnifracWeightedCommand : public Command { public: - UnifracWeightedCommand(); - ~UnifracWeightedCommand() { delete weighted; } + UnifracWeightedCommand(string); + ~UnifracWeightedCommand() { if (abort == false) { delete weighted; delete util; } } int execute(); + void help(); private: + struct linePair { + int start; + int num; + linePair(int i, int j) : start(i), num(j) {} + }; + vector lines; + GlobalData* globaldata; + SharedUtil* util; + FileOutput* output; vector T; //user trees vector utreeScores; //user tree unweighted scores vector WScoreSig; //tree weighted score signifigance when compared to random trees - percentage of random trees with that score or lower. vector groupComb; // AB. AC, BC... - Tree* randT; //random tree TreeMap* tmap; Weighted* weighted; - string weightedFile, weightedFileout, sumFile; + string sumFile, outputDir; int iters, numGroups, numComp, counter; EstOutput userData; //weighted score info for user tree EstOutput randomData; //weighted score info for random trees @@ -44,22 +54,24 @@ class UnifracWeightedCommand : public Command { vector< map > rScoreFreq; //map -vector entry for each combination. vector< map > rCumul; //map -vector entry for each c map validScores; //map contains scores from random - map::iterator it2; - map::iterator it; - ofstream outSum, out; - ifstream inFile; + bool abort, phylip, random; + string groups, itersString; + vector Groups, outputNames; //holds groups to be used + int processors; + + + ofstream outSum; void printWSummaryFile(); void printWeightedFile(); + void createPhylipFile(); //void removeValidScoresDuplicates(); int findIndex(float, int); - void setGroups(); void calculateFreqsCumuls(); - void initFile(string); - void output(vector); - void resetFile(); - + int createProcesses(Tree*, vector< vector >, vector&, vector< vector >&); + int driver(Tree*, vector< vector >, int, int, vector&, vector< vector >&); + };