Computing Marginal Probability of Fixed Effect Over Time

[Alexander Jack]

Hello,

I am working on a spatiotemporal model where I have interactions between fixed effects and time. These interactions are significant and represent something I would like to discuss in a paper. So I would like to be able to show readers how the marginal probability of a given fixed effect changes over time. I am able to generate a marginal for the entirety of the time series using inla.rmarginal or something similar but how would I do that for each specific timepoint in my series so I can get confidence intervals and estimates over time as well?

Alex 

[Helpdesk (Haavard Rue)]

vignette("svc", package="INLA") ???


On Wed, 2024-07-17 at 11:21 -0700, 'Alexander Jack' via R-inla discussion group
wrote:

> --
> 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/e56f323c-d932-4c85-ae3f-c679c431e05an%40googlegroups.com
> .

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

[Alexander Jack]

Would this work for lattice data?

[Helpdesk (Haavard Rue)]

the concept will work for essentially anything.

if you have a covariate x and want to check its dependence in time, then

y ~ x + ...

is replaced with

y ~ f(time, x, model="...") + ..

where you have to add the model for the regression coefficient, like

model="rw1",constr=FALSE

or

model="rw2",constr=FALSE




On Wed, 2024-07-17 at 12:30 -0700, 'Alexander Jack' via R-inla discussion group