import mir.math.common: approxEqual;

assert(0.5.betaProportionCCDF(0.5, 2).approxEqual(0.5));
assert(0.75.betaProportionCCDF((1.0 / 3), 3).approxEqual(0.0625));
assert(0.25.betaProportionCCDF((1.0 / 9), 4.5).approxEqual(0.1411133));