Fit Besag Small Area Estimation model using TMB.

Random effects have an improper conditional autoregressive (ICAR) distribution with (generalised) precision matrix produced using the nb_to_precision function with input nb, the neighbourhood structure of sf.

Usage

besag_tmb(sf, its = 1000, ii = NULL)

Arguments

sf

A simple features object with some geometry.

its

Number of iterations in outer loop optimisation, passed to nlminb.

ii

The (zero-indexed) indices of the observations held-out.

Examples

besag_tmb(mw, its = 100)
#> Error in .class1(object): object 'mw' not found