add p to sim.factor do.call