use scipy inverse beta function for exact confidence intervals, if available