Can we manually create the stopping criteria of INLA?

18 views
Skip to first unread message

Alokesh Manna

unread,
Dec 9, 2024, 11:47:11 AM12/9/24
to R-inla discussion group
Hi INLA group,

Can we manually create the stopping criteria like "it will run 10000 steps"? In MCMC we can do that. Is there any default value for INLA?

Alokesh

Helpdesk (Haavard Rue)

unread,
Dec 13, 2024, 2:04:11 AM12/13/24
to Alokesh Manna, R-inla discussion group
not from R, no.
> --
> 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/b928f0ae-e669-43eb-a1d2-9ab1c1288a4en%40googlegroups.com
> .

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

Alokesh Manna

unread,
Dec 13, 2024, 8:28:15 AM12/13/24
to Helpdesk, R-inla discussion group
Thanks for the feedback!

Any way without R?

In future it maybe great to see some traceplot etc.

Finn Lindgren

unread,
Dec 13, 2024, 8:34:24 AM12/13/24
to Alokesh Manna, Helpdesk, R-inla discussion group
While seeing the optimiser trace might be useful for internal debugging, I don’t think it would be of general use (optimisers don’t necessarily take an humanly logical path towards the optimum) and the amount of data storage would be huge, as would the added time of actually writing all those values to disk; storing hundreds of copies of the model would be really bad. That said, nonlinear inlabru models do store the trace (usually the number of iterations is small, so the memory requirement is less than it would be to store the full inla optimisation) of its iterated inla method for nonlinear predictor models; bru_convergence_plot() displays a summary (in light of my comments above I’ll likely at least make it an option to only store the convergence summaries).
Finn

On 13 Dec 2024, at 13:28, Alokesh Manna <alokesh...@gmail.com> wrote:



Helpdesk (Haavard Rue)

unread,
Dec 13, 2024, 8:57:39 AM12/13/24
to Finn Lindgren, Alokesh Manna, R-inla discussion group

If you want to see the optimization path, you can do

results <- inla(...)

plot(result, plot.opt.trace=TRUE)

which will give two plots; as attached
Screenshot from 2024-12-13 16-55-49.png

Alokesh Manna

unread,
Dec 13, 2024, 9:03:57 AM12/13/24
to Helpdesk, Finn Lindgren, R-inla discussion group
Very helpful!

Finn Lindgren

unread,
Dec 13, 2024, 9:04:53 AM12/13/24
to Helpdesk, Alokesh Manna, R-inla discussion group
Just in case someone thinks my comment doesn't match Haavard's, I was
thinking of the full latent variable trace, not just the theta; but
good to know that the theta trace is there!
Finn
--
Finn Lindgren
email: finn.l...@gmail.com

Alokesh Manna

unread,
Dec 13, 2024, 12:21:59 PM12/13/24
to Finn Lindgren, Helpdesk, R-inla discussion group
Thank you!
Reply all
Reply to author
Forward
0 new messages