# \donttest{
# Define two models and sample them
# options(mc.cores = 3) # Speed up sampling
ex = mcp_example("intercepts") # Get some simulated data.
model1 = list(y ~ 1 + x, ~ 1)
model2 = list(y ~ 1 + x) # Without a change point
fit1 = mcp(model1, ex$data)
fit2 = mcp(model2, ex$data)
# Compute LOO for each and compare (works for waic(fit) too)
fit1$loo = loo(fit1)
fit2$loo = loo(fit2)
loo::loo_compare(fit1$loo, fit2$loo)
# }
Run the code above in your browser using DataLab