inla.mesh.project: Methods for projecting to/from an inla.mesh

In andrewzm/INLA: Functions which allow to perform full Bayesian analysis of latent Gaussian models using Integrated Nested Laplace Approximaxion

Description

Calculate a lattice projection to/from an inla.mesh

Usage

inla.mesh.project(...)
inla.mesh.projector(...)

## S3 method for class 'inla.mesh'
inla.mesh.projector(mesh,
             loc = NULL,
             lattice = NULL,
             xlim = range(mesh$loc[,1]),
             ylim = range(mesh$loc[,2]),
             dims = c(100,100),
             projection = NULL,
             ...)

## S3 method for class 'inla.mesh.1d'
inla.mesh.projector(mesh,
             loc = NULL,
             xlim = mesh$interval,
             dims = 100, ...)

## S3 method for class 'inla.mesh.projector'
inla.mesh.project(projector, field, ...)

## S3 method for class 'inla.mesh'
inla.mesh.project(mesh, loc, field = NULL, ...)
## S3 method for class 'inla.mesh.1d'
inla.mesh.project(mesh, loc, field = NULL, ...)

Arguments

mesh	An inla.mesh or inla.mesh.1d object.
loc	Projection locations.
lattice	An inla.mesh.lattice object.
xlim	X-axis limits for a lattice.
ylim	Y-axis limits for a lattice.
dims	Lattice dimensions.
projector	An inla.mesh.projector object.
field	Basis function weights, one per mesh basis function, describing the function to be avaluated at the projection locationssFunction values for on the mesh
projection	One of c("default", "longlat", "longsinlat", "mollweide").
...	Additional arguments passed on to methods.

Details

The call inla.mesh.project(mesh, loc, field=..., ...), is a shortcut to inla.mesh.project(inla.mesh.projector(mesh, loc), field).

Value

For inla.mesh.project(mesh, ...), a list with projection information. For inla.mesh.projector(mesh, ...), an inla.mesh.projector object. For inla.mesh.project(projector, field, ...), a field projected from the mesh onto the locations given by the projector object.

Author(s)

Finn Lindgren [e-mail]

See Also

inla.mesh, inla.mesh.1d, inla.mesh.lattice

Examples

n = 20
loc = matrix(runif(n*2), n, 2)
mesh = inla.mesh.create(loc, refine=list(max.edge=0.05))
proj = inla.mesh.projector(mesh)
field = cos(mesh$loc[,1]*2*pi*3)*sin(mesh$loc[,2]*2*pi*7)
image(proj$x, proj$y, inla.mesh.project(proj, field))

if (require(rgl)) {
  plot(mesh, rgl=TRUE, col=field, draw.edges=FALSE, draw.vertices=FALSE)
}

andrewzm/INLA documentation built on May 10, 2019, 11:12 a.m.