1d:T6da,# Summary for posteriors of field parameters for an `inla_rspde` model from a `rspde_result` object Summary for posteriors of rSPDE field parameters in their original scales. ```r ## S3 method for class 'rspde_result' summary(object, digits = 6, ...) ``` ## Arguments - `object`: A `rspde_result` object. - `digits`: integer, used for number formatting with signif() - `...`: Currently not used. ## Returns Returns a `data.frame` containing the summary. ## Examples ```r #tryCatch version tryCatch({ 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) } #stable.tryCatch }, error = function(e){print("Could not run the example")}) ```1e:T49a,M7.49933 0.25C3.49635 0.25 0.25 3.49593 0.25 7.50024C0.25 10.703 2.32715 13.4206 5.2081 14.3797C5.57084 14.446 5.70302 14.2222 5.70302 14.0299C5.70302 13.8576 5.69679 13.4019 5.69323 12.797C3.67661 13.235 3.25112 11.825 3.25112 11.825C2.92132 10.9874 2.44599 10.7644 2.44599 10.

summary.rspde_result function