Stata DO script for cluster-level prevalence estimates