]> git.donarmstrong.com Git - mothur.git/blobdiff - unifracweightedcommand.h
changed unifrac.weighted()
[mothur.git] / unifracweightedcommand.h
index 047c92b2574515e20c439931808dce8dc5a34705..9deef732bfe100e27b99fb2e30551de0256e453e 100644 (file)
  *
  */
 
-#include <iostream>
-#include <fstream>
-#include <vector>
-#include <map>
-#include <cmath>
 #include "command.hpp"
 #include "weighted.h"
 #include "treemap.h"
@@ -39,7 +34,7 @@ class UnifracWeightedCommand : public Command {
                Tree* randT;  //random tree
                TreeMap* tmap;
                Weighted* weighted;
-               string weightedFile, sumFile, distFile;
+               string weightedFile, sumFile;
                int iters, numGroups, numComp;
                EstOutput userData;                     //weighted score info for user tree
                EstOutput randomData;           //weighted score info for random trees
@@ -47,16 +42,17 @@ class UnifracWeightedCommand : public Command {
                vector< map<float, float> > rscoreFreq;  //vector<weighted score, number of random trees with that score.> each group comb has an entry
                vector< map<float, float> > uscoreFreq;  //vector<weighted, number of user trees with that score.> each group comb has an entry
                vector< map<float, float> > totalrscoreFreq;  //vector<weighted score, number of random trees with that score.> each group comb has an entry
-               vector< map<float, float> > rCumul;  //vector<weighted score, number of random trees with that score.> each group comb has an entry
-               vector< map<float, float> > uCumul;  //vector<weighted, cumulative percentage of number of user trees with that score or lower.> each group comb has an entry
+               vector< map<float, float> > rCumul;  //vector<weighted score, number of random trees with that score or higher.> each group comb has an entry
+               vector< map<float, float> > uCumul;  //vector<weighted, cumulative percentage of number of user trees with that score or higher.> each group comb has an entry
                map<float, float>::iterator it;
                map<float, float>::iterator it2;
                
-               ofstream outSum, outDist, out;
+               ofstream outSum, out;
                
                void printWSummaryFile();
                void printWeightedFile();  
-               void saveRandomScores(); 
+               void saveRandomScores();
+               void setGroups(); 
 };