- float maxbin = 0.0;
- for (int i = 0; i < rabund.size(); i++) {
- if (rabund.get(i) != 0) { //don't want log value of 0.
- if (scaler == "log10") {
- colorScale[(log10((rabund.get(i) / (float)rabund.getNumSeqs()) * 1000))] = "";
- if (maxbin < (log10((rabund.get(i) / (float)rabund.getNumSeqs()) * 1000))) { maxbin = (log10((rabund.get(i) / (float)rabund.getNumSeqs()) * 1000)); }
- }else if (scaler == "log2") {
- colorScale[(log2((rabund.get(i) / (float)rabund.getNumSeqs()) * 1000))] = "";
- if (maxbin < (log2((rabund.get(i) / (float)rabund.getNumSeqs()) * 1000))) { maxbin = (log2((rabund.get(i) / (float)rabund.getNumSeqs()) * 1000)); }
- }else if (scaler == "linear") {
- colorScale[rabund.get(i)] = "";
- if (maxbin < rabund.get(i)) { maxbin = rabund.get(i); }
- }else { //if user enters invalid scaler option.
- cout << scaler << " is not a valid scaler option. I will use log10." << endl;
- colorScale[(log10((rabund.get(i) / (float)rabund.getNumSeqs()) * 1000))] = "";
- if (maxbin < (log10((rabund.get(i) / (float)rabund.getNumSeqs()) * 1000))) { maxbin = (log10((rabund.get(i)) / (float)rabund.getNumSeqs()) * 1000); }
- }
- }else { colorScale[0] = "00"; }
+ float maxRelAbund = 0.0;
+
+ for(int i=0;i<rabund.size();i++){
+ float relAbund = rabund.get(i) / (float)rabund.getNumSeqs();
+ if(relAbund > maxRelAbund){ maxRelAbund = relAbund; }