Bayesian vs Frequentist Statistical Model for Randomized Benchmarking

008_demo.py

with h_model:
    BoundedGamma = pm.Bound(pm.Gamma,
                            lower= 0.0005,
                            upper= 0.004 + 0.001 * (N - 1) )
    sigma = BoundedGamma("σ_Beta",
                           alpha = 10 - 5 * (N - 1), 
                           beta = 10000 - 8000 * (N - 1),
                           testval = 0.001 + 0.0015 * (N - 1), 
                           shape = counts.shape[1])