get_point_estimate: Get a point estimate for a single location

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

Description

Get a point estimate for a single location

Usage

get_point_estimate(lon, lat, dpc_grid, fit)

Arguments

lon	Longitude of point
lat	Latitude of point
dpc_grid	Discrete Process Convolution grid
fit	Fitted model

Value

Posterior predictive mean at given location

Examples

s = data.frame(lat = rep(0, 3), lon = c(-1:1))
y = c(1,1,0)
dpc_grid = get_grid(c(-1,1), c(0,0), spacing = 2)
priors = get_priors(dpc_grid)
iso_kernel_matrix = get_kernel_matrix(s, dpc_grid)
fit = get_mcmc(s, dpc_grid, y, 10, 1000, priors, 100, 1)
get_point_estimate(lon = s$lon[1], lat = s$lat[1], dpc_grid = dpc_grid, fit = fit)

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