X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=treemap.h;h=7ffd1e7003154d672d544aa8a993f371a7f9a28a;hp=d9617f5455e03f0f5fc702fffaf3516335351b15;hb=615301e57c25e241356a9c2380648d117709458d;hpb=58cf1d08fee8c64334979075fa57bcafb035a2ed

diff --git a/treemap.h b/treemap.h
index d9617f5..7ffd1e7 100644
--- a/treemap.h
+++ b/treemap.h
@@ -9,11 +9,8 @@
  *
  */
 
-#include <map>
-#include <string>
-#include <iostream>
-#include <fstream>
-#include "utilities.hpp"
+#include "mothur.h"
+#include "listvector.hpp"
 
 /* This class is used by the read.tree command to build the tree container. */
 
@@ -22,29 +19,51 @@ struct GroupIndex {
 	int		vectorIndex;
 };
 
-
-
 class TreeMap {
 public:
+	TreeMap() { m = MothurOut::getInstance(); }
 	TreeMap(string);
 	~TreeMap();
-	void readMap();
+    
+	int readMap();
+    int readMap(string);
 	int getNumGroups();
 	int getNumSeqs();
-	void setIndex(string, int);  //sequencename, index
-	int getIndex(string);		//returns vector index of sequence
+	//void setIndex(string, int);  //sequencename, index
+	//int getIndex(string);		//returns vector index of sequence
+	bool isValidGroup(string);  //return true if string is a valid group
+	void removeSeq(string);  //removes a sequence, this is to accomadate trees that do not contain all the seqs in your groupfile
 	string getGroup(string);
-	vector<string> namesOfGroups;
-	vector<string> namesOfSeqs;
-	void print(ostream&);
-	
+	void addSeq(string, string);
+	void addGroup(string s) { setNamesOfGroups(s); }
+	vector<string> getNamesOfGroups() {
+		sort(namesOfGroups.begin(), namesOfGroups.end());
+		return namesOfGroups;
+	}
+    
+    void print(ostream&);
+	void makeSim(vector<string>);  //takes groupmap info and fills treemap for use by tree.shared command.
+	void makeSim(ListVector*);  //takes listvector info and fills treemap for use by tree.shared command.	
+    vector<string> getNamesSeqs();
+	vector<string> getNamesSeqs(vector<string>); //get names of seqs belonging to a group or set of groups
+    int getCopy(TreeMap&);
+    
+    vector<string> namesOfSeqs;
+    map<string,int> seqsPerGroup;	//groupname, number of seqs in that group.
+	map<string, GroupIndex> treemap; //sequence name and <groupname, vector index>
+
+    
 private:
+	vector<string> namesOfGroups;
 	ifstream fileHandle;
 	string groupFileName;
 	int numGroups;
 	map<string, GroupIndex>::iterator it;
+	map<string, int>::iterator it2;
 	void setNamesOfGroups(string); 
-	map<string, GroupIndex> treemap; //sequence name and groupname
+	MothurOut* m;
+	
+	
 };
 
 #endif