Summary for posteriors of field parameters for an inla_rspde model from a rspde_result object
Source: R/inla_rspde.R
summary.rspde_result.RdSummary for posteriors of rSPDE field parameters in their original scales.
Usage
# S3 method for class 'rspde_result'
summary(object, digits = 6, ...)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)
}
#> Warning: the mean or mode of nu is very close to nu.upper.bound, please consider increasing nu.upper.bound, and refitting the model.
#> mean sd 0.025quant 0.5quant 0.975quant mode
#> tau 0.0123577 0.0311694 4.10394e-12 7.15796e-04 0.0986204 4.10394e-12
#> kappa 26.5330000 13.7946000 1.28897e+01 2.24129e+01 64.0107000 1.60940e+01
#> nu 1.5398200 0.4572140 5.65259e-01 1.71056e+00 1.9989900 1.99998e+00
# devel.tag
# }