Preparations before fitting with inla()
- create a stack with the data for estimation and prediction
# stack for estimation stk.e
stk.e <- inla.stack(
tag = "est",
data = list(y = d$value),
A = list(1, A),
effects = list(
data.frame(
b0 = rep(1, nrow(A)),
covtemp = d$covtemp,
covprec = d$covprec
),
s = indexs
)
)
# stack for prediction stk.p
stk.p <- inla.stack(
tag = "pred",
data = list(y = NA),
A = list(1, Ap),
effects = list(
data.frame(
b0 = rep(1, nrow(Ap)),
covtemp = dp$covtemp,
covprec = dp$covprec
),
s = indexs
)
)
# stk.full has stk.e and stk.p
stk.full <- inla.stack(stk.e, stk.p)