X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=tree.h;h=0b316c3042f01cbb000a3aed975a778f8af5f9fc;hb=fca3f55d5ded10c3dc77856f3cc4a1c53b02bb6f;hp=6404db784eec791783e609dfc8ff8d8a971c61aa;hpb=9099173ea17b13be96563c7579ecc8dc2ef65d20;p=mothur.git diff --git a/tree.h b/tree.h index 6404db7..0b316c3 100644 --- a/tree.h +++ b/tree.h @@ -10,34 +10,42 @@ * */ -using namespace std; - #include "treenode.h" #include "globaldata.hpp" +class GlobalData; /* This class represents the treefile. */ class Tree { public: + Tree(string); + Tree(int); Tree(); //to generate a tree from a file - ~Tree() {}; - + ~Tree(); void getCopy(Tree*); //makes tree a copy of the one passed in. + void getSubTree(Tree*, vector); //makes tree a that contains only the names passed in. void assembleRandomTree(); - void assembleRandomUnifracTree(); - void createNewickFile(); + void assembleRandomUnifracTree(vector); + void assembleRandomUnifracTree(string, string); + void createNewickFile(string); int getIndex(string); void setIndex(string, int); int getNumNodes() { return numNodes; } int getNumLeaves(){ return numLeaves; } - map mergeUserGroups(int); //returns a map with a groupname and the number of times that group was seen in the children + map mergeUserGroups(int, vector); //returns a map with a groupname and the number of times that group was seen in the children void printTree(); + void print(ostream&); + void print(ostream&, string); + int findRoot(); //return index of root node //this function takes the leaf info and populates the non leaf nodes - void assembleTree(); + int assembleTree(); + int assembleTree(string); vector tree; //the first n nodes are the leaves, where n is the number of sequences. + map< string, vector > groupNodeInfo; //maps group to indexes of leaf nodes with that group, different groups may contain same node because of names file. + private: GlobalData* globaldata; int numNodes, numLeaves; @@ -46,14 +54,24 @@ private: map::iterator it, it2; map mergeGroups(int); //returns a map with a groupname and the number of times that group was seen in the children + map mergeGcounts(int); - map Tree::mergeGcounts(int); + void addNamesToCounts(); void randomTopology(); void randomBlengths(); - void randomLabels(); - int findRoot(); //return index of root node - void printBranch(int); //recursively print out tree - void setGroups(); + void randomLabels(vector); + //void randomLabels(string, string); + void printBranch(int, ostream&, string); //recursively print out tree + void parseTreeFile(); //parses through tree file to find names of nodes and number of them + //this is required in case user has sequences in the names file that are + //not included in the tree. + //only takes names from the first tree in the tree file and assumes that all trees use the same names. + int readTreeString(ifstream&); + int populateNewTree(vector&, int, int&); + void printBranch(int, ostream&, string, vector&); + + MothurOut* m; + }; #endif