powered by
jsmc(mod1, mod2, nsamp = 5000)
temp <- jsmc(gmmpen[[1]], gmmpen[[2]])
Run the code above in your browser using DataLab