working on simulate(how=mean)