- int maxbin = 0;
- for (int i = 0; i < lookup.size(); i++) {
- for (int j = 0; j < lookup[i]->size(); j++) {
- //if (lookup[i]->getAbundance(j) != 0) { //don't want log value of 0.
- //if (scaler == "log10") {
- // colorScale[-log((log10(lookup[i]->getAbundance(j)) / (float)lookup[i]->getNumSeqs()))] = "";
- //cout << "abundance = " << lookup[i]->getAbundance(j) << '\t' << " relative adundance = " << (lookup[i]->getAbundance(j) / (float)lookup[i]->getNumSeqs()) << '\t';
- //cout << -log((log10(lookup[i]->getAbundance(j)) / lookup[i]->getNumSeqs())) << endl;
- //}else if (scaler == "log2") {
- //colorScale[-log((log2(lookup[i]->getAbundance(j)) / (float)lookup[i]->getNumSeqs()))] = ""; //cout << (int)log2(lookup[i]->getAbundance(j)) << endl;
- //cout << "abundance = " << lookup[i]->getAbundance(j) << '\t' << " relative adundance = " << lookup[i]->getAbundance(j) / (float)lookup[i]->getNumSeqs() << '\t';
- //cout << -log((log2(lookup[i]->getAbundance(j)) / lookup[i]->getNumSeqs())) << endl;
- // }else if (scaler == "linear") {
- colorScale[lookup[i]->getAbundance(j)] = "";
- if (maxbin < lookup[i]->getAbundance(j)) { maxbin = lookup[i]->getAbundance(j); }
- //cout << "abundance = " << lookup[i]->getAbundance(j) << '\t' << " relative adundance = " << lookup[i]->getAbundance(j) / (float)lookup[i]->getNumSeqs() << '\t';
- //cout << lookup[i]->getAbundance(j) /(float) lookup[i]->getNumSeqs() << endl;
- //}else { //if user enters invalid scaler option.
- // cout << scaler << " is not a valid scaler option. I will use log10." << endl;
- // colorScale[-log(log10(lookup[i]->getAbundance(j) / lookup[i]->getNumSeqs()))] = "";
- //}
- //}else { colorScale[0] = "00"; }