rSPDE (version 2.5.1)

precision.intrinsicCBrSPDEobj: Get the precision matrix of intrinsicCBrSPDEobj objects

Description

Function to get the precision matrix of a intrinsicCBrSPDEobj object

Usage

# S3 method for intrinsicCBrSPDEobj
precision(
  object,
  kappa = NULL,
  tau = NULL,
  alpha = NULL,
  beta = NULL,
  ld = FALSE,
  ...
)

Value

The precision matrix.

Arguments

object

The model object computed using intrinsic.matern.operators()

kappa

If non-null, update the range parameter.

tau

If non-null, update the precision parameter.

alpha

If non-null, update the alpha parameter.

beta

If non-null, update the beta parameter.

ld

If TRUE, return the log determinant of the precision matrix instead of the precision matrix. By default FALSE.

...

Currently not used.

See Also

simulate.intrinsicCBrSPDEobj(), intrinsic.matern.operators()

Examples

if (requireNamespace("RSpectra", quietly = TRUE)) {
  x <- seq(from = 0, to = 10, length.out = 201)
  beta <- 1
  alpha <- 1
  kappa <- 1
  op <- intrinsic.matern.operators(
    kappa = kappa, tau = 1, alpha = alpha,
    beta = beta, loc_mesh = x, d = 1
  )
Q <- precision(op) 
}