result <- gmjmcmc.parallel(
runs = 1,
cores = 1,
list(populations = "best", complex.measure = 2, tol = 0.0000001),
matrix(rnorm(600), 100),
P = 2,
gaussian.loglik,
loglik.alpha = gaussian.loglik.alpha,
c("p0", "exp_dbl")
)
plot(result)
Run the code above in your browser using DataLab