spam (version 0.40-0)

makeprec: Create Precision Matrices

Description

Creates precision matrices for gridded GMRF.

Usage

precmat.GMRFreglat(n,m, par, model = "m1p1",  eps = .Spam$eps)

Arguments

n

first dimension of the grid.

m

second dimension of the grid.

par

parameters used to construct the matrix.

model

see details and examples.

eps

A tolerance parameter: elements of x such that abs(x) <= eps<="" code=""> set to zero. Defaults to eps = .Spam$eps

Value

• A spam matrix of dimension prod(dims)xprod(dims).

Details

The function should be illustrative on how to create precision matrices for gridded GMRF. Hence, no testing (positive definiteness is done). Please see the examples on the meaning of the different models.

See Also

precmat, toeplitz.spam, kronecker.spam

Examples

as.matrix(precmat.GMRFreglat(4, 3, c(.4),         'm1p1'))
as.matrix(precmat.GMRFreglat(4, 3, c(.4,.3),      'm1p2'))
as.matrix(precmat.GMRFreglat(4, 3, c(.4,.3,.2),   'm2p3'))
as.matrix(precmat.GMRFreglat(4, 3, c(.4,.3,.2,.1),'m2p4'))