dpsdPosSim(I = 30, J = 200, K = 6, muN = -0.7, s2aN = 0.2, s2bN = 0.2, muD = 0, s2aD = 0.2, s2bD = 0.2, muR = qnorm(0.25), s2aR = 0.2, s2bR = 0.2, crit = matrix(rep(c(-1.6, -0.5, 0, 0.5, 1.6), each = I), ncol = (K - 1)))I here~~J here~~K here~~muN here~~s2aN here~~s2bN here~~muD here~~s2aD here~~s2bD here~~muR here~~s2aR here~~s2bR here~~crit here~~