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~~