Skip to contents

Summary for posteriors of field parameters for an inla_rspde model from a rspde_result object

Source: R/inla_rspde.R
summary.rspde_result.Rd

Summary for posteriors of rSPDE field parameters in their original scales.

Usage

# S3 method for rspde_result
summary(object, digits = 6, ...)

Arguments

object

A rspde_result object.

digits

integer, used for number formatting with signif()

...

Currently not used.

Value

Returns a data.frame

containing the summary.

Examples

# \donttest{
# devel version
if (requireNamespace("INLA", quietly = TRUE)) {
  library(INLA)

  set.seed(123)

  m <- 100
  loc_2d_mesh <- matrix(runif(m * 2), m, 2)
  mesh_2d <- inla.mesh.2d(
    loc = loc_2d_mesh,
    cutoff = 0.05,
    max.edge = c(0.1, 0.5)
  )
  sigma <- 1
  range <- 0.2
  nu <- 0.8
  kappa <- sqrt(8 * nu) / range
  op <- matern.operators(
    mesh = mesh_2d, nu = nu,
    range = range, sigma = sigma, m = 2,
    parameterization = "matern"
  )
  u <- simulate(op)
  A <- inla.spde.make.A(
    mesh = mesh_2d,
    loc = loc_2d_mesh
  )
  sigma.e <- 0.1
  y <- A %*% u + rnorm(m) * sigma.e
  Abar <- rspde.make.A(mesh = mesh_2d, loc = loc_2d_mesh)
  mesh.index <- rspde.make.index(name = "field", mesh = mesh_2d)
  st.dat <- inla.stack(
    data = list(y = as.vector(y)),
    A = Abar,
    effects = mesh.index
  )
  rspde_model <- rspde.matern(
    mesh = mesh_2d,
    nu.upper.bound = 2
  )
  f <- y ~ -1 + f(field, model = rspde_model)
  rspde_fit <- inla(f,
    data = inla.stack.data(st.dat),
    family = "gaussian",
    control.predictor =
      list(A = inla.stack.A(st.dat))
  )
  result <- rspde.result(rspde_fit, "field", rspde_model)
  summary(result)
}
#>             mean        sd 0.025quant   0.5quant 0.975quant       mode
#> tau    0.0251086 0.0109143  0.0080407  0.0239145  0.0494529  0.0199181
#> kappa 16.8878000 3.4370600 11.3481000 16.4601000 24.7903000 15.5779000
#> nu     0.9383100 0.1241800  0.7309400  0.9249410  1.2100900  0.8777150
# devel.tag
# }