write script to compare mean/median regions
run small simulation to compare estimators