Smoothness parameter in dmatern

[Nafsika Antoniadou]

Hello,

I have the model:

form12 = y ~ -1 + intercept +
  f(idx, model="dmatern", locations = loc,
    hyper = list(range = list(prior = "pc.range", param = c(30, 0.5)),
                 theta1 = list(prior="pc.prec", param=c(0.4, 0.05)),
                 theta3 = list(prior="loggamma", param=c(1, 0.05))))

# Fit the model
fit <- inla(form12,
            data = data_df_gpd,
            control.predictor = list(compute = TRUE, link = 1),
            family = "gp",
            control.family = list(list(control.link = list(quantile = q_gpd))),
            control.compute = list(dic = TRUE, waic = TRUE, cpo = TRUE),
            control.inla = list(strategy = "laplace", int.strategy = "ccd"),
            verbose = FALSE,
            num.threads = 2)
But in: 
fit$summary.hyperpar I cannot see the posterior of the smoothness parameter (theta3 ).

Is there another way? Thank you!

[Helpdesk (Haavard Rue)]

the parameter is by default fixed, but you can set fixed=FALSE, and it will be
estimated.

Alternative, is to use the rSPDE package

https://davidbolin.github.io/rSPDE/

> --
> You received this message because you are subscribed to the Google Groups "R-
> inla discussion group" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to r-inla-discussion...@googlegroups.com.
> To view this discussion, visit
> https://groups.google.com/d/msgid/r-inla-discussion-group/d1ba6eb6-d76b-4477-ad88-174ff220f27fn%40googlegroups.com
> .

--
Håvard Rue
he...@r-inla.org