Truncation errors (falling back to rejection sampling despite logcdf present)

Logcdf isn’t enough, icdf is also needed to avoid rejection sampling. According to this issue, that method is still missing for BetaBinomial: Add icdf functions for distributions · Issue #6612 · pymc-devs/pymc · GitHub