Fit Centroid MVN Small Area Estimation model using aghq.

Random effects have a multivariate Gaussian distribution with covariance matrix calculated using centroid_covariance. Kernel hyper-parameters are fixed.

Usage

fck_aghq(sf, k = 3, kernel = matern, ii = NULL, ...)

Arguments

sf

A simple features object with some geometry.

k

Number of quadrature points per hyperparameter dimension.

kernel

A kernel function, defaults to matern.

ii

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

...

Additional arguments to kernel.

Examples

fck_aghq(mw)
#> Error in sf::st_centroid(sf): object 'mw' not found