data <- data.frame(matrix(rnorm(600), 100))
result <- mjmcmc.parallel(runs = 2,
cores = 1,
y = matrix(rnorm(100), 100),
x = data,
loglik.pi = gaussian.loglik)
compute_effects(result,labels = names(data))
Run the code above in your browser using DataLab