1 #ifndef INDICATORCOMMAND_H
2 #define INDICATORCOMMAND_H
8 * Created by westcott on 11/12/10.
9 * Copyright 2010 Schloss Lab. All rights reserved.
13 #include "command.hpp"
15 #include "counttable.h"
16 #include "sharedrabundvector.h"
17 #include "sharedrabundfloatvector.h"
18 #include "inputdata.h"
20 class IndicatorCommand : public Command {
22 IndicatorCommand(string);
26 vector<string> setParameters();
27 string getCommandName() { return "indicator"; }
28 string getCommandCategory() { return "Hypothesis Testing"; }
29 string getOutputFileNameTag(string, string);
30 string getHelpString();
31 string getCitation() { return "Dufrene M, Legendre P (1997). Species assemblages and indicator species: The need for a flexible asymmetrical approach. Ecol Monogr 67: 345-66.\n McCune B, Grace JB, Urban DL (2002). Analysis of ecological communities. MjM Software Design: Gleneden Beach, OR. \nLegendre P, Legendre L (1998). Numerical Ecology. Elsevier: New York. \nhttp://www.mothur.org/wiki/Indicator"; }
32 string getDescription() { return "calculate the indicator value for each OTU"; }
35 void help() { m->mothurOut(getHelpString()); }
41 string treefile, sharedfile, relabundfile, groups, label, inputFileName, outputDir, designfile;
43 int iters, processors;
44 vector<string> outputNames, Groups;
45 vector<SharedRAbundVector*> lookup;
46 vector<SharedRAbundFloatVector*> lookupFloat;
50 int GetIndicatorSpecies(Tree*&);
51 int GetIndicatorSpecies();
52 set<string> getDescendantList(Tree*&, int, map<int, set<string> >, map<int, set<int> >&);
53 vector<float> getValues(vector< vector<SharedRAbundVector*> >&, map< vector<int>, vector<int> >);
54 vector<float> getValues(vector< vector<SharedRAbundFloatVector*> >&, map< vector<int>, vector<int> >);
55 map<int, float> getDistToRoot(Tree*&);
56 map< vector<int>, vector<int> > randomizeGroupings(vector< vector<SharedRAbundVector*> >&, int);
57 map< vector<int>, vector<int> > randomizeGroupings(vector< vector<SharedRAbundFloatVector*> >&, int);
58 vector<float> driver(vector< vector<SharedRAbundFloatVector*> >&, map< vector<int>, vector<int> >, int, vector<float>, int);
59 vector<float> driver(vector< vector<SharedRAbundVector*> >&, map< vector<int>, vector<int> >, int, vector<float>, int);
60 vector<float> getPValues(vector< vector<SharedRAbundFloatVector*> >&, map< vector<int>, vector<int> >, int, vector<float>);
61 vector<float> getPValues(vector< vector<SharedRAbundVector*> >&, map< vector<int>, vector<int> >, int, vector<float>);