lladserconst

confirmedupstream bugPython · scikit-bioSystems & data structuresJuly 25, 2026

the question

Does scikit-bio's Lladser confidence interval use the correct table constant at rarefaction depth r=9?

what came out

scikit-bio's _CB_95 table stores 4.695227540 at r=9, a copy of _LOWER_CONFIDENCE_BOUND[9] instead of the true ~1.44. lladser_ci then returns a 95 percent interval 3.22x wider at r=9 than at r=8 and r=10, breaking three self-consistency invariants. Silent.

method & receipts

Result: confirmed
Reproducible — one script re-runs the whole thing from scratch.
Tested — a correctness/benchmark suite ships alongside the code.

→ read the code and re-run it

github.com/v-code01/lladserconst