 Solution

Use a hierarchical Bayesian model with binomial observation errors and a common distribution for the underlying risks:Give $(\alpha ,\beta )$ a shared hyperprior, or parameterize them by mean $m=\alpha /(\alpha +\beta )$ and concentration $\kappa =\alpha +\beta$ and put priors on $m\in (0,1)$ and $\kappa >0$. Conditional on these hyperparameters, the risks are independent and identically distributed; after integrating them out they are exchangeable random variables with shared uncertainty. The data estimate the overall risk level and between-hospital variation. Conditional Beta-binomial conjugacy givesand full Bayes point estimates under squared-error loss areAn Empirical Bayes method instead estimates the hyperparameters from the marginal beta-binomial distribution likelihoodand inserts their fitted values in the conditional posterior means. Both approaches provide partial pooling, stronger for smaller hospitals. Full Bayes also propagates hyperparameter uncertainty into the credible intervals. If systematic covariates are needed, one can place exchangeable residual effects in a logistic regression model rather than assuming a common unconditional risk distribution.