1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
|
//
// sharedrabundfloatvectors.hpp
// Mothur
//
// Created by Sarah Westcott on 5/15/17.
// Copyright © 2017 Schloss Lab. All rights reserved.
//
#ifndef sharedrabundfloatvectors_hpp
#define sharedrabundfloatvectors_hpp
#include "datavector.hpp"
#include "rabundvector.hpp"
#include "sharedrabundfloatvector.hpp"
#include "sharedordervector.h"
#include "ordervector.hpp"
/* DataStructure for a relabund file. */
//********************************************************************************************************************
inline bool compareRAbundFloats(SharedRAbundFloatVector* left, SharedRAbundFloatVector* right){ return (left->getGroup() < right->getGroup()); }
//********************************************************************************************************************
class SharedRAbundFloatVectors : public DataVector {
public:
SharedRAbundFloatVectors() : DataVector() { label = ""; numBins = 0; otuTag = "Otu"; }
SharedRAbundFloatVectors(string oTag) : DataVector() { label = ""; numBins = 0; otuTag = oTag; }
SharedRAbundFloatVectors(ifstream&, vector<string>&, string&, string&);
SharedRAbundFloatVectors(SharedRAbundFloatVectors& bv) : DataVector(bv), numBins(bv.numBins), otuTag(bv.otuTag) {
vector<SharedRAbundFloatVector*> data = bv.getSharedRAbundFloatVectors();
for (int i = 0; i < data.size(); i++) { push_back(data[i]); }
setLabels(bv.getLabel());
setOTUNames(bv.getOTUNames());
eliminateZeroOTUS();
}
~SharedRAbundFloatVectors() { clear(); }
void setLabels(string l);
float getOTUTotal(int bin);
vector<float> getOTU(int bin);
float removeOTU(int bin);
float removeOTUs(vector<int>, bool sorted=false); //bins to remove, sorted or not
float get(int bin, string group);
void set(int bin, float binSize, string group);
void setOTUNames(vector<string> names);
vector<string> getOTUNames();
string getOTUName(int);
void setOTUName(int, string);
int push_back(vector<float>, string binLabel=""); //add otu. mothur assumes abunds are in same order as groups.
int push_back(SharedRAbundFloatVector*);
void removeGroups(vector<string> g);
int removeGroups(int minSize, bool silent=false); // removes any groups with numSeqs < minSize
void resize(int n) { m->mothurOut("[ERROR]: can not use resize for SharedRAbundFloatVectors.\n"); }
void clear() { for (int i = 0; i < lookup.size(); i++) { if (lookup[i] != nullptr) { delete lookup[i]; lookup[i] = nullptr; } } lookup.clear(); groupNames.clear(); numBins = 0; currentLabels.clear(); }
int size() { return (int)lookup.size(); }
int getNumGroups() { return (int)lookup.size(); }
int getNumBins() { return numBins; }
float getNumSeqs(string); //group
float getNumSeqsSmallestGroup();
void print(ostream&, bool&);
vector<SharedRAbundVector*> getSharedRAbundVectors();
vector<SharedRAbundFloatVector*> getSharedRAbundFloatVectors();
RAbundVector getRAbundVector();
SAbundVector getSAbundVector();
OrderVector getOrderVector(map<string,int>*) { m->mothurOut("[ERROR]: can not convert SharedRAbundVectors to an ordervector, ordervectors assume no zero OTUS.\n"); m->setControl_pressed(true); OrderVector o; return o; }
vector<string> getNamesGroups();
void eliminateZeroOTUS(); //run after push_backs if groups are chosen
private:
void printHeaders(ostream&, bool&);
vector<SharedRAbundFloatVector*> lookup;
vector<string> currentLabels;
map<string, int> groupNames;
int numBins;
string otuTag;
};
#endif /* sharedrabundfloatvectors_hpp */
|
