[R] Simulation of 1PL and 2PL IRT data