st_coordinates()
sf
inla()
res$summary.fixed
res$summary.fitted.values
res$marginals.fitted.values
Extract needed subsets with: inla.stack.index(stack = , tag =)
inla.stack.index(stack = , tag =)
Compute the probability of exceeding a threshold (PM2.5 >= 10), by using inla.pmarginal():
inla.pmarginal()
1 - inla.pmarginal(q = 10, marginal = marg)