Susan Elias
2018-Aug-23 22:09 UTC
[R] How does Markov random field (bs=mrf) in mgvc gam handle repeated measures on the spatial units?
Hello R,
I am attempting a spatio-temporal model in mgcv gam.
I am using a factor smooth to define each of 27 areal units in a shapefile
("id") as subjects (essentially) which have undergone 23 repeated
measurements in time ("year").
I am using the Markov random field to define the neighbor relationships of
the area units.
I have combed through Wood (2017) and other resources and have not found an
answer to my question, which is: how is the mrf handling the repeated
measures? My model is:
model<-gam(response) ~ s(year) +
ti(year,id, bs="fs", m=1) +
s(id, bs = 'mrf', xt = list(nb = neighbors),
k=5),
family=tw(), data=data)
Susan
[[alternative HTML version deleted]]