1 #ifndef MGCLUSTERCOMMAND_H
2 #define MGCLUSTERCOMMAND_H
8 * Created by westcott on 12/11/09.
9 * Copyright 2009 Schloss Lab. All rights reserved.
13 #include "command.hpp"
14 #include "readblast.h"
15 #include "sparsematrix.hpp"
16 #include "nameassignment.hpp"
17 #include "globaldata.hpp"
18 #include "cluster.hpp"
21 /**********************************************************************/
23 class MGClusterCommand : public Command {
26 MGClusterCommand(string);
32 GlobalData* globaldata;
34 NameAssignment* nameMap;
39 vector<DistNode> overlapMatrix;
43 string blastfile, method, namefile, overlapFile, distFile;
44 ofstream sabundFile, rabundFile, listFile;
45 float cutoff, penalty;
46 int precision, length, precisionLength;
47 bool abort, minWanted, hclusterWanted, exitedBreak, merge;
49 void printData(ListVector*);
50 ListVector* mergeOPFs(map<string, int>, float);
51 void sortHclusterFiles(string, string);
52 vector<DistNode> getSeqs(ifstream&);
56 /**********************************************************************/