Added option to control n.sims, default n.sims=10,000