X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=bayesian.h;h=405fee3679df7b94ae8f392f62c27886e0b11e36;hp=0020962f4276c9deb2129a27b2d69974783aab82;hb=615301e57c25e241356a9c2380648d117709458d;hpb=a5afca18544555fba2d9c3670ad1f8574916b0a0 diff --git a/bayesian.h b/bayesian.h index 0020962..405fee3 100644 --- a/bayesian.h +++ b/bayesian.h @@ -18,30 +18,29 @@ class Bayesian : public Classify { public: - Bayesian(string, string, string, int, int, bool); - ~Bayesian() {}; + Bayesian(string, string, string, int, int, int, int, bool, bool); + ~Bayesian(); string getTaxonomy(Sequence*); - //map getConfidenceScores() { return taxConfidenceScore; } private: - //map > taxonomyProbability; //probability of a word being in a given taxonomy. - //taxonomyProbability[bacteria;][0] = probabtility that a sequence within bacteria; would contain kmer 0; - //taxonomyProbability[bacteria;][1] = probabtility that a sequence within bacteria; would contain kmer 1... - vector< vector > wordGenusProb; //vector of maps from genus to probability - //wordGenusProb[0][392] = probability that a sequence within genus that's index in the tree is 392 would contain kmer 0; + //wordGenusProb[0][392] = probability that a sequence within genus that's index in the tree is 392 would contain kmer 0; vector genusTotals; vector genusNodes; //indexes in phyloTree where genus' are located - int kmerSize, numKmers, confidenceThreshold; - bool probs; + vector WordPairDiffArr; + + int kmerSize, numKmers, confidenceThreshold, iters; string bootstrapResults(vector, int, int); int getMostProbableTaxonomy(vector); - void readProbFile(ifstream&, ifstream&); - //map parseTaxMap(string); + void readProbFile(ifstream&, ifstream&, string, string); + bool checkReleaseDate(ifstream&, ifstream&, ifstream&, ifstream&); + bool isReversed(vector&); + vector createWordIndexArr(Sequence*); + int generateWordPairDiffArr(); };