Try dplyr do() for the simulations