X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=unifracweightedcommand.h;h=9ea78baf17370d7c068e40d56014f2f4f6d49010;hb=cd37904452dc95b183ff313ff05720c562902487;hp=1b8b12bf789da9892a78017e7f888e0ad54ace01;hpb=bfbc55964f1977da72c2cea984288a427d370a59;p=mothur.git diff --git a/unifracweightedcommand.h b/unifracweightedcommand.h index 1b8b12b..9ea78ba 100644 --- a/unifracweightedcommand.h +++ b/unifracweightedcommand.h @@ -13,45 +13,55 @@ #include "command.hpp" #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() { delete weighted; delete util; } int execute(); + void help(); private: 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 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, sumFile, distFile; - int iters, numGroups, numComp; + string sumFile; + int iters, numGroups, numComp, counter; EstOutput userData; //weighted score info for user tree EstOutput randomData; //weighted score info for random trees - vector< map > validScores; //vector each group comb has an entry - vector< map > rscoreFreq; //vector each group comb has an entry - vector< map > uscoreFreq; //vector each group comb has an entry - vector< map > totalrscoreFreq; //vector each group comb has an entry - vector< map > rCumul; //vector each group comb has an entry - vector< map > uCumul; //vector each group comb has an entry - map::iterator it; - map::iterator it2; + vector< vector > rScores; //vector each group comb has an entry + vector< vector > uScores; //vector each group comb has an entry + 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 - ofstream outSum, outDist, out; + bool abort; + string groups, itersString; + vector Groups; //holds groups to be used + + + ofstream outSum; void printWSummaryFile(); void printWeightedFile(); - void saveRandomScores(); + //void removeValidScoresDuplicates(); + int findIndex(float, int); + void calculateFreqsCumuls(); + };