*
*/
-using namespace std;
-
-
#include "sharedordervector.h"
-#include "utilities.hpp"
-
+#include "sharedutilities.h"
/***********************************************************************/
//This function is used to read a .shared file for the collect.shared, rarefaction.shared and summary.shared commands
//if you don't use a list and groupfile.
-SharedOrderVector::SharedOrderVector(ifstream& f) : DataVector() {
+SharedOrderVector::SharedOrderVector(ifstream& f) : DataVector() { //reads in a shared file
try {
globaldata = GlobalData::getInstance();
+ maxRank = 0; numBins = 0; numSeqs = 0;
if (globaldata->gGroupmap == NULL) { groupmap = new GroupMap(); }
for(int i=0;i<num;i++){
f >> inputData;
+
for (int j = 0; j < inputData; j++) {
- push_back(i+1, i+1, groupN);
+ push_back(i, i, groupN);
numSeqs++;
}
}
for(int i=0;i<num;i++){
f >> inputData;
+
for (int j = 0; j < inputData; j++) {
- push_back(i+1, i+1, groupN);
+ push_back(i, i, groupN);
numSeqs++;
}
}
f.seekg(pos, ios::beg);
if (globaldata->gGroupmap == NULL) { globaldata->gGroupmap = groupmap; }
-
+
updateStats();
+
}
catch(exception& e) {
cout << "Standard Error: " << e.what() << " has occurred in the SharedOrderVector class Function SharedOrderVector. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
/***********************************************************************/
int SharedOrderVector::getNumBins(){
- if(needToUpdate == 1){ updateStats(); }
return numBins;
}
/***********************************************************************/
int SharedOrderVector::getNumSeqs(){
- if(needToUpdate == 1){ updateStats(); }
return numSeqs;
}
/***********************************************************************/
int SharedOrderVector::getMaxRank(){
- if(needToUpdate == 1){ updateStats(); }
return maxRank;
}
data[index].group = groupName;
data[index].bin = binNumber;
data[index].abundance = abund;
- needToUpdate = 1;
-
+ //if (abund > maxRank) { maxRank = abund; }
+ updateStats();
}
/***********************************************************************/
/***********************************************************************/
-
+//commented updateStats out to improve speed, but whoever calls this must remember to update when they are done with all the pushbacks they are doing
void SharedOrderVector::push_back(int binNumber, int abund, string groupName){
individual newGuy;
newGuy.group = groupName;
newGuy.abundance = abund;
newGuy.bin = binNumber;
data.push_back(newGuy);
- needToUpdate = 1;
+ //numSeqs++;
+ //numBins++;
+ //if (abund > maxRank) { maxRank = abund; }
+ //updateStats();
}
/***********************************************************************/
exit(1);
}
+}
+/***********************************************************************/
+vector<SharedRAbundVector*> SharedOrderVector::getSharedRAbundVector() {
+ try {
+ SharedUtil* util;
+ util = new SharedUtil();
+ vector<SharedRAbundVector*> lookup;
+
+ util->setGroups(globaldata->Groups, globaldata->gGroupmap->namesOfGroups);
+ util->getSharedVectors(globaldata->Groups, lookup, this);
+
+ return lookup;
+ }
+ catch(exception& e) {
+ cout << "Standard Error: " << e.what() << " has occurred in the SharedOrderVector class Function getSharedRAbundVector. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
+ exit(1);
+ }
+ catch(...) {
+ cout << "An unknown error has occurred in the SharedOrderVector class function getSharedRAbundVector. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
+ exit(1);
+ }
+
}
/***********************************************************************/
SharedSAbundVector SharedOrderVector::getSharedSAbundVector(string group) {
/***********************************************************************/
SharedOrderVector SharedOrderVector::getSharedOrderVector(){
+ random_shuffle(data.begin(), data.end());
return *this;
}
numBins = 0;
maxRank = 0;
- for(int i=0;i<data.size();i++){
- if(data[i].bin != -1){
- numSeqs++;
- }
- }
-
- vector<individual> hold(numSeqs);
-
+ numSeqs = data.size();
+
+ vector<int> hold(numSeqs, 0);
for(int i=0;i<numSeqs;i++){
- if(data[i].bin != -1){
- hold[data[i].bin].bin = hold[data[i].bin].bin+1;
- }
+ hold[data[i].bin] = hold[data[i].bin]+1;
}
-
+
for(int i=0;i<numSeqs;i++){
- if(data[i].bin > numBins) { numBins = data[i].bin; }
- if(data[i].abundance > maxRank) { maxRank = data[i].abundance; }
+ if(hold[i] > 0) { numBins++; }
+ if(hold[i] > maxRank) { maxRank = hold[i]; }
}
- numBins++; //if you have 10 bins largest .bin is 9 since we start at 0.
+
}
catch(exception& e) {
cout << "Standard Error: " << e.what() << " has occurred in the SharedOrderVector class Function updateStats. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";