INLAutils

INLAutils copied to clipboard

x <- rnorm(100)
y <- x + rnorm(100, sd = 0.2)
d <- data.frame(x, y)

result1 = inla(y ~ x, family="normal", data = d, 
               control.predictor = list(compute = TRUE, link = 1),
               control.compute = list(return.marginals.predictor=TRUE))

plot_inla_residuals(result1, y)

looks horrible even with correct model

Maybe I should just switch to PIT values.