Description Usage Arguments Examples
View source: R/plot_inla_residuals.R
Plot residuals of observed vs predicted values for INLA model
1 | plot_inla_residuals(inla.model, observed)
|
inla.model |
An inla object |
observed |
The observed values |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | ## Not run:
library(INLA)
data(Epil)
observed <- Epil[1:30, 'y']
Epil <- rbind(Epil, Epil[1:30, ])
Epil[1:30, 'y'] <- NA
## make centered covariates
formula = y ~ Trt + Age + V4 +
f(Ind, model="iid") + f(rand,model="iid")
result = inla(formula, family="poisson", data = Epil,
control.predictor = list(compute = TRUE, link = 1))
plot_inla_residuals(result, observed)
data(Seeds)
l <- nrow(Seeds)
Seeds <- rbind(Seeds, Seeds)
Seeds$r[1:l] <- NA
formula = r ~ x1 * x2 + f(plate, model = "iid")
mod.seeds = inla(formula, data=Seeds, family = "binomial", Ntrials = n,
control.predictor = list(compute = TRUE, link = 1))
plot_inla_residuals(mod.seeds, na.omit(Seeds$r / Seeds$n))
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.