Extract results

• res$summary.fixed
• res$summary.fitted.values
• res$marginals.fitted.values

Extract needed subsets with: inla.stack.index(stack = , tag =)

Compute the probability of exceeding a threshold (PM2.5 >= 10), by using inla.pmarginal():

1 - inla.pmarginal(q = 10, marginal = marg)