# \donttest{
if (bru_safe_inla(multicore = FALSE)) {
# Generate some data
input.df <- data.frame(idx = 1:10, x = cos(1:10))
input.df <- within(
input.df,
y <- rpois(10, 5 + 2 * cos(1:10) + rnorm(10, mean = 0, sd = 0.1))
)
# Fit two models
fit1 <- bru(
y ~ x,
family = "poisson",
data = input.df,
options = list(control.compute = list(dic = TRUE))
)
fit2 <- bru(
y ~ x + rand(idx, model = "iid"),
family = "poisson",
data = input.df,
options = list(control.compute = list(dic = TRUE))
)
# Compare DIC
deltaIC(fit1, fit2)
}
# }
Run the code above in your browser using DataLab