Estimating prior sd for the parameter p of a beta binomial regression

You can also use

pm.sample_prior_predictive()