program define bssim
  version 6.0
  if "`1'" == "?" {
    global S_1 "lb_stan ub_stan lb_n ub_n lb_p ub_p lb_bc ub_bc"
    exit
  }
  use welfsub, clear
  gen x = uniform()
  sort x
  keep if _n <= `2'
  ci agrc
  local ub_stan = `r(ub)'
  local lb_stan = `r(lb)'
  bs "summarize agrc" "r(mean)", reps(1000)
  post `1' `lb_stan' `ub_stan' r(lb_n) r(ub_n) r(lb_p) r(ub_p) r(lb_bc) r(ub_bc)
end