Groups keyboard shortcuts have been updated
Dismiss
See shortcuts

INLA takes a long time even when the hyperparameters are fixed

17 views
Skip to first unread message

Alokesh Manna

unread,
Dec 7, 2024, 5:05:51 AM12/7/24
to R-inla discussion group
Hi INLA users,

In my computation when I fix the hyperparameters precision (log precision), INLA takes a long time to converge. Can you tell me what can be a possible reason? I write a snippet of code.

PS: I have a 15 random effect spatially varying coefficients. As INLA takes a long time to converge if I do not fix i am putting some initial value. Smaller precision sometime helps. But with an example of big data, without fixing the hyperparameters I got the initial guess.


f(as.numeric(spatial_ID_13),value_r1__value_u1,
                               constr=FALSE,
                               model = "besag",scale.model=TRUE,
                               graph = graph_adjacency_2,
                               hyper=list(theta = list(prior = "loggamma",
                                                       param=c(1,.0003),
                                                      initial=log(5.0),
                                                      fixed=TRUE,
                                                      scale.model=TRUE)
                               )
                              )


The other alignment is 

```
 family = "poisson" ,
                           
                            data = data.frame(x_matrix_2) ,
                           
                           
                            control.predictor = list(compute = TRUE),
                            control.compute = list(openmp.strategy = "huge",
                                                   dic = TRUE,config = TRUE,
                                                   save.memory = TRUE),
                            control.inla = list(
                              #diagonal = .01,
                              strategy="gaussian",
                              int.strategy = "eb"
                              ),
                            #verbose=TRUE,
                           
                            control.mode=list(restart=T)
```

Helpdesk (Haavard Rue)

unread,
Dec 7, 2024, 5:16:52 AM12/7/24
to Alokesh Manna, R-inla discussion group

can you upgrade to the most recent testing version and retry? INLA_24.12.06

> --
> 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/bef4fa35-9265-4712-9185-b15a91c0c9c9n%40googlegroups.com
> .

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

Reply all
Reply to author
Forward
0 new messages