X-Git-Url: https://git.donarmstrong.com/?p=mothur.git;a=blobdiff_plain;f=lefsecommand.h;h=cb1bbb6da80f7da1567979dd158fb5861e74a685;hp=c9c4092e99f6bb402145f2e3432753ce642c7c6e;hb=d1c97b8c04bb75faca1e76ffad60b37a4d789d3d;hpb=b25ede2ad307ae76f8a610443e0ec3ec69621ce7 diff --git a/lefsecommand.h b/lefsecommand.h index c9c4092..cb1bbb6 100644 --- a/lefsecommand.h +++ b/lefsecommand.h @@ -41,23 +41,34 @@ public: string getOutputPattern(string); string getHelpString(); - string getCitation() { return "http://www.mothur.org/wiki/Lefse"; } + string getCitation() { return "Segata, N., J. Izard, L. Waldron, D. Gevers, L. Miropolsky, W. S. Garrett, and C. Huttenhower. 2011. Metagenomic biomarker discovery and explanation. Genome Biol 12:R60, http://www.mothur.org/wiki/Lefse"; } string getDescription() { return "brief description"; } int execute(); void help() { m->mothurOut(getHelpString()); } private: - bool abort, allLines; - string outputDir, sharedfile, designfile, mclass, subclass, classes; + bool abort, allLines, wilc, wilcsamename, curv, subject, normMillion; + string outputDir, sharedfile, designfile, mclass, subclass, rankTec, multiClassStrat; vector outputNames; set labels; - float anovaAlpha, wilcoxonAlpha; + double anovaAlpha, wilcoxonAlpha, fBoots, ldaThreshold; + int nlogs, iters, strict, minC; - int process(vector&, DesignMap&); - vector runKruskalWallis(vector&, DesignMap&); - vector runWilcoxon(vector&, DesignMap&, vector); - + int process(vector&, DesignMap&); + int normalize(vector&); + map runKruskalWallis(vector&, DesignMap&); + map runWilcoxon(vector&, DesignMap&, map, map >& class2SubClasses, map >& subClass2GroupIndex, map); + bool testOTUWilcoxon(map >& class2SubClasses, vector abunds, map >& subClass2GroupIndex, map); + map testLDA(vector&, map, map >& class2GroupIndex, map >&); + bool contastWithinClassesOrFewPerClass(vector< vector >&, vector rands, int minCl, map > class2GroupIndex, map indexToClass); + vector< vector > lda(vector< vector >& adjustedLookup, vector rand_s, map& indexToClass, vector); + vector< vector > getMeans(vector& lookup, map >& class2GroupIndex); + int printResults(vector< vector >, map, map, string, vector); + + //for testing + bool printToCoutForRTesting(vector< vector >& adjustedLookup, vector rand_s, map >& class2GroupIndex, map bins, map >&, vector); + int makeShared(int); }; /**************************************************************************************************/