Re: [R-sig-Geo] BYM model in R-INLA - how to specify priors individually

2014-05-06 Thread James Rooney
model in R-INLA - how to specify priors individually Dear James, You should be able to set the prior of your random effects using hyper= in the call to f(). In your particular case, you can try: #Create areas IDs to match the values in nc.adj nc.sids$ID-1:100 nc.sids$ID2-1:100 hyper.besag -list

[R-sig-Geo] BYM model in R-INLA - how to specify priors individually

2014-05-05 Thread James Rooney
Hi all, I am attempting to implement the Besag-York-Mollie model in R-INLA. Specifically to recreate the WinBUGS model specified in ASDAR fig 11.10. However I cannot figure out how to specify the spatial and random priors to be different as they are in the WinBUGS model fig11.10. i.e. from

Re: [R-sig-Geo] BYM model in R-INLA - how to specify priors individually

2014-05-05 Thread VIRGILIO GOMEZ RUBIO
Dear James, You should be able to set the prior of your random effects using hyper= in the call to f(). In your particular case, you can try: #Create areas IDs to match the values in nc.adj nc.sids$ID-1:100 nc.sids$ID2-1:100 hyper.besag -list(prec=list(prior=loggamma, params=c(.1, .1)))