gammaLikeSim(NN = 1, NS = 1, I = 30, J = 200, K = 6, muS = log(0.5), s2aS = 0.2, s2bS = 0.2, lagEffect = 0, shape = 2, crit = matrix(rep(c(0.75, 0.8, 1, 1.35, 1.6), each = I), ncol = (K - 1)))NN here~~NS here~~I here~~J here~~K here~~muS here~~s2aS here~~s2bS here~~lagEffect here~~shape here~~crit here~~