Return the Q matrix from INLA

[Matt Suttton]

Hi all,

I am currently working on a new method and was hoping to gain access to the Q(theta) matrix for specific theta and x. For this run theta would be held fixed and I am looking at returning a function:

f(x) = log pi(theta) + 0.5*log(|Q(theta)|) - 0.5*x'Q(theta)x

I've had a toy with the rgeneric function and can see some useful functions there but I would like to ideally return this from a fitted INLA object. Is this possible?

Best,

Matt

[Helpdesk]

would that be Q(theta) for the whole model?

r = inla(..., control.compute=list(config=TRUE))

will return the 'configurations' in

r$misc$configs

and

r$misc$config[[1]]

etc.

is this what you want ?

> --
> 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 on the web, visit
> https://groups.google.com/d/msgid/r-inla-discussion-group/5c03f111-0b2e-472a-a88d-bc75bff67797n%40googlegroups.com
> .

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