Fit Besag Small Area Estimation model using R-INLA.

Random effects have an improper conditional autoregressive (ICAR) distribution. This is implemented by R-INLA using the option model = "besag".

Usage

besag_inla(sf, verbose = FALSE, cores = parallel::detectCores())

Arguments

sf

A simple features object with some geometry.

verbose

Should R-INLA run in mode verbose = TRUE.

cores

Number of cores, passed to R-INLA, defaults to parallel::detectCores().

Examples

besag_inla(mw, cores = 2)
#> Error in .class1(object): object 'mw' not found