Add na.rm=T to region functions
Add simulation results, remove n=10 from simulation (bad behavior for small n)