- SAbundVector sabund1 = shared[0]->getSAbundVector();
- SAbundVector sabund2 = shared[1]->getSAbundVector();
+ double aux = ceil(pow(sampleSize+1, (1/(double)3)));
+ double est0 = rank->get(1)+1;
+ if (aux > est0) { est0 = aux; } //est0 = max(rank->get(1)+1, aux)
+
+ vector<double> ests;
+ double numr = 0.0;
+ for (int i = 1; i < rank->getNumBins()-1; i++) {
+
+ if (m->control_pressed) { break; }
+
+ int abund = rank->get(i);
+
+ if (abund != 0) {
+
+ int abundNext = rank->get(i+1);
+ if (abundNext == 0) { numr = aux; }
+ else {
+ if (abundNext+1 > aux) { numr = abundNext+1; } //numr = max(abundNext+1, aux)
+ else { numr = aux; }
+ }
+ double denr = aux;
+ if (abund > aux) { denr = abund; } //denr = max(abund, aux)
+ ests.push_back((abund+1)*numr/denr);
+ }
+ }
+ numr = aux;