Hello,
I have the model:
form13 = y ~ -1 + intercept + temp +
f(idx, model="dmatern", locations = loc,
hyper = list(range = list(prior = "pc.range", param = c(30, 0.05)),
theta1 = list(prior="pc.prec", param=c(0.4, 0.05))))
fit_gpd_pred = inla(form13,
data = data_df_gpd_pred,
control.predictor = list(compute=TRUE, link=1),
family = "gp",
lincomb=lcs,
control.family = list(list(control.link = list(quantile = q_gpd))),
control.compute = list(dic = TRUE, waic = TRUE, cpo = TRUE, config = TRUE),
control.inla=list(strategy="laplace",int.strategy="ccd"),
verbose = FALSE,
num.threads = 2) which returns
Deviance Information Criterion (DIC) ...............: NA
Deviance Information Criterion (DIC, saturated) ....: NA
Effective number of parameters .....................: NA
WAIC is computed. If I remove the covariate and I have only an intercept in form13 then it is computed. WHy does this happen?