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")
)
preds <- predict(result$results, matrix(rnorm(600), 100))
Run the code above in your browser using DataLab