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
87
88
|
//
// sharedrabundfloatvector.hpp
// Mothur
//
// Created by Sarah Westcott on 7/25/17.
// Copyright © 2017 Schloss Lab. All rights reserved.
//
#ifndef sharedrabundfloatvector_hpp
#define sharedrabundfloatvector_hpp
#include "datavector.hpp"
#include "rabundvector.hpp"
#include "sabundvector.hpp"
#include "ordervector.hpp"
/* Data Structure for a rabund file.
This class is a child to datavector. It represents OTU information at a certain distance.
A rabundvector can be converted into and ordervector, listvector or sabundvector.
Each member of the internal container "data" represents an individual OTU.
So data[0] = 6, because there are six member in that OTU.
example: listvector = a,b,c,d,e,f g,h,i j,k l m
rabundvector = 6 3 2 1 1
sabundvector = 2 1 1 0 0 1
ordervector = 1 1 1 1 1 1 2 2 2 3 3 4 5
a,b,g,j,l = sample A
c,d,h,i = sample B
e,f,k,m = sample C
The sharedRabund class is very similar to rabund. SharedRabund allows for 0 otus in a sample, rabunds do not. SharedRabund also know their group. SharedRabunds are stored as floats, but printed as integers or floats depending on whether it represents a shared or relabund file.
sharedrabund = A 5 2 1 1 1 0
B 5 2 2 0 0 0
C 5 2 0 1 0 1
*/
//class RAbundFloatVector;
//class OrderVector;
class SharedRAbundFloatVector : public DataVector {
public:
SharedRAbundFloatVector();
SharedRAbundFloatVector(int);
SharedRAbundFloatVector(vector<float>, float, int, float); //maxRank, numbins, numSeqs
SharedRAbundFloatVector(vector<float>);
SharedRAbundFloatVector(const SharedRAbundFloatVector& bv) : DataVector(bv), data(bv.data), maxRank(bv.maxRank), numBins(bv.numBins), numSeqs(bv.numSeqs), group(bv.group) {};
SharedRAbundFloatVector(ifstream&);
SharedRAbundFloatVector(ifstream& f, string l, string g, int); //filehandle, label
~SharedRAbundFloatVector();
int getNumBins();
float getNumSeqs();
float getMaxRank();
float remove(int);
float remove(vector<int>);
void set(int, float);
float get(int);
vector<float> get() { return data; }
void push_back(float);
void resize(int);
int size();
void clear();
void print(ostream&); //nonsorted
RAbundVector getRAbundVector();
SharedRAbundVector getSharedRAbundVector();
RAbundFloatVector getRAbundFloatVector();
SAbundVector getSAbundVector();
OrderVector getOrderVector(map<string,int>*);
string getGroup() { return group; } //group = "" for rabunds without groupInfo
void setGroup(string g) { group = g; }
private:
vector<float> data;
float maxRank;
int numBins;
float numSeqs;
string group;
};
#endif /* sharedrabundfloatvector_hpp */
|
