get_em: Optimize kernel parameters

In rtlemos/scallops: Lemos and Sanso' (2012) model, edited for speed

Description

Optimize kernel parameters

Usage

get_em(nburn, priors, s, y, iso_kernel_matrix)

Arguments

nburn	Number of iterations
priors	Prior distributions for gamma (Gaussian) and tau (Gamma)
s	Coordinates of observations
y	Values of observations
iso_kernel_matrix	Isotropic kernel matrix

Value

List with optimized kernel parameters, as well as corresponding mean values of gamma and tau

Examples

s = expand.grid(lat = c(-1:1), lon = c(-1:1))
y = c(1,0,0,0,1,0,0,0,1)
dpc_grid = get_grid(c(-1,1), c(-1,1), spacing = 2)
priors = get_priors(dpc_grid)
iso_kernel_matrix = get_kernel_matrix(s, dpc_grid)
fit = get_em(10, priors, s, y, iso_kernel_matrix)
cbind(get_estimates(s, dpc_grid, fit), obs = y)

rtlemos/scallops documentation built on May 4, 2019, 7:43 p.m.