Bayesian inference under cluster sampling with probability proportional to size