Skip to contents

Define a spacetime model object for the f() call.

Usage

barrierModel.define(
  mesh,
  barrier.triangles,
  prior.range,
  prior.sigma,
  range.fraction = 0.1,
  constr = FALSE,
  debug = FALSE,
  verbose = FALSE,
  useINLAprecomp = TRUE,
  libpath = NULL
)

Arguments

mesh

a spatial mesh

barrier.triangles

a integer vector to specify which triangles centers are in the barrier domain, or a list with integer vector if more than one.

prior.range

numeric vector containing U and a to define the probability statements P(range < U) = a used to setup the PC-prior for range. If a = 0 then U is taken to be the fixed value for the range.

prior.sigma

numeric vector containing U and a to define the probability statements P(range > U) = a used to setup the PC-prior for sigma. If a = 0 then U is taken to be the fixed value for sigma.

range.fraction

numeric to specify the fraction of the range for the barrier domain. Default value is 0.1. This has to be specified with care in order to have it small enough to make it act as barrier but not too small in order to prevent numerical issues.

constr

logical to indicate if the integral of the field over the domain is to be constrained to zero. Default value is FALSE.

debug

logical indicating if to run in debug mode.

verbose

logical indicating if to print parameter values.

useINLAprecomp

logical indicating if is to be used shared object pre-compiled by INLA. This will not be considered if the argument libpath is provided.

libpath

string to the shared object. Default is NULL.

Value

objects to be used in the f() formula term in INLA.

Details

See the paper.