Fit Integrated MVN Small Area Estimation model using aghq.

Random effects have a multivariate Gaussian distribution with covariance matrix calculated using integrated_covariance.

Usage

fik_aghq(
  sf,
  k = 3,
  L = 10,
  type = "hexagonal",
  kernel = matern,
  ii = NULL,
  ...
)

Arguments

sf

A simple features object with some geometry.

k

Number of quadrature points per hyperparameter dimension.

L

The number of Monte Carlo samples to draw from each area.

type

The type argument of sf::st_sample, defaults to "hexagonal"

kernel

A kernel function, defaults to matern.

ii

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

...

Additional arguments to kernel.

Examples

fik_tmb(mw)
#> Error in nrow(sf): object 'mw' not found