自由度为30的泊松分布
program a,reclassified
drop _all
set obs 30
gen x=rpoisson
sum x
return scalar mean_a=r(mean)
end

simulate xbar=r(mean_a),seed(1234)reps(10000)nodots:a

hist xbar,normal xtitle("xbar from a")