plot_inla_residuals: Plot residuals of observed vs predicted values for INLA model

In timcdlucas/INLAutils: Utility Functions for 'INLA'

Plot residuals of observed vs predicted values for INLA model

Description

Plot residuals of observed vs predicted values for INLA model

Usage

plot_inla_residuals(inla.model, observed)

Arguments

inla.model	An inla object
observed	The observed values

Examples

## 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),
               control.compute = list(return.marginals.predictor=TRUE))
 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),
                  control.compute = list(return.marginals.predictor=TRUE))
 plot_inla_residuals(mod.seeds, na.omit(Seeds$r / Seeds$n))
 
## End(Not run)

timcdlucas/INLAutils documentation built on Nov. 29, 2022, 5:41 a.m.