if (square == true) { readSquareMatrix(filehandle); }
else { readLTMatrix(filehandle); }
-
-
- printMatrix(cout);
//sort sequences so they are gathered in groups for processing
sortGroups(0, numSeqs-1);
- cout << "after sort" << endl;
- printMatrix(cout);
-
+
}
catch(exception& e) {
cout << "Standard Error: " << e.what() << " has occurred in the FullMatrix class Function FullMatrix. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";
i++;
j--;
-cout << "swapping rows " << i << " " << j << endl;
-printMatrix(cout); cout << endl;
}
} while(i <= j);
order = SharedList->getSharedOrderVector();
}
+ //set users groups
+ setGroups();
+
//output estimator names as column headers
outputFileHandle << "label" <<'\t' << "comparison" << '\t';
for(int i=0;i<sumCalculators.size();i++){
}
outputFileHandle << endl;
- //set users groups
- setGroups();
-
while(order != NULL){
if(globaldata->allLines == 1 || globaldata->lines.count(count) == 1 || globaldata->labels.count(order->getLabel()) == 1){
try {
lookup.clear();
//create and initialize vector of sharedvectors, one for each group
- for (int i = 0; i < globaldata->gGroupmap->getNumGroups(); i++) {
+ for (int i = 0; i < globaldata->Groups.size(); i++) {
SharedRAbundVector* temp = new SharedRAbundVector(order->getNumBins());
temp->setLabel(order->getLabel());
- temp->setGroup(globaldata->gGroupmap->namesOfGroups[i]);
+ temp->setGroup(globaldata->Groups[i]);
lookup.push_back(temp);
}
}
}
}
-
- //get rid of vectors from groups you don't want to analyze
- for (int r = 0; r < lookup.size(); r++) {
- if (inUsersGroups(lookup[r]->getGroup(), globaldata->Groups) != true) {
- lookup.erase(lookup.begin()+r);
- }
- }
-
}
catch(exception& e) {
cout << "Standard Error: " << e.what() << " has occurred in the SummarySharedCommand class Function getSharedVectors. Please contact Pat Schloss at pschloss@microbio.umass.edu." << "\n";