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]]