simMPRM function

simulate data according to MPRM