get_interpolation_plot: Plot the interpolated surface

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

Description

Plot the interpolated surface

Usage

get_interpolation_plot(obs_coord, dpc_grid, fit, fineness = 16,
  contour_binwidth = NULL, do_grid = TRUE)

Arguments

obs_coord	Coordinates of original observations
dpc_grid	Discrete Process Convolution grid
fit	Fitted model
fineness	Resolution of interpolation (1 = coarse, 10 = fine)
contour_binwidth	Width of bins for contour plot
do_grid	Plot DPC gridpoints?

Value

Plot of gridded posterior predictive means

Examples

s = expand.grid(lat = -1:1, lon=-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_mcmc(s, dpc_grid, y, 10, 1000, priors, 100, 1)
get_interpolation_plot(s, dpc_grid, fit, contour_binwidth = 0.1)

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