model
{
for(i in 1:N){
for(t in 1:S){
y[i,t]~dnorm(mu[i,t],sigmay)
mu[i,t]~dnorm(lamda[i,t],sigmamu)
lamda[i,t]<-alpha0+u+v+T
}
u~dnorm(0,sigma1)
RR<-exp(u+v)
T~dnorm(0,sigma2)
}
v[1:N]~car.normal(adj[],weights[],num[],tau)
for(k in 1:sumNumNeigh){
weights[k]<-1}
sigmay~dgamma(3.0E-3,1)
sigmamu~dgamma(2.0E-4,1)
sigma1~dunif(0,10)
sigma2~dunif(0,10)
alpha0~dflat()
tau~dgamma(0.5,5.0E-4)
}
#data
list(N=27,S=4,y=structure(
.Data=c(0.214285714,0.075,0.103448276,0.086206897,
0.308835027,0.328498912,0.264375877,0.24952381,
0.299898683,0.332673267,0.339412361,0.267364415,
0.068965517,0.066666667,0.135135135,0.2,
0.213754647,0.286516854,0.261595547,0.181019332,
0.293367347,0.355223881,0.325433526,0.268604139,
0.21,0.24137931,0.235955056,0.213836478,
0.256910569,0.32756324,0.290148448,0.254745255,
0.180851064,0.274428274,0.279461279,0.211911955,
0.345303867,0.380722892,0.333333333,0.242524917,
0.203937008,0.21141649,0.243523316,0.205426357,
0.263934426,0.269230769,0.294294294,0.269574944,
0.186138614,0.266787659,0.205769231,0.180159635,
0.254413292,0.285338016,0.252622378,0.214614878,
0.263265306,0.284313725,0.26827012,0.278653295,
0.263684913,0.274740484,0.251147541,0.208706786,
0.099518459,0.137535817,0.112275449,0.106309421,
0.248847926,0.310185185,0.264483627,0.227504244,
0.271505376,0.37514518,0.326708075,0.265886288,
0.241157556,0.30994152,0.334261838,0.203809524,
0.235457064,0.229828851,0.272727273,0.262435678,
0.257142857,0.234042553,0.256756757,0.221674877,
0.256345178,0.252427184,0.223844282,0.225742574,
0.278074866,0.327102804,0.216216216,0.256944444,
0.187155963,0.188821752,0.193006993,0.193548387,
0.349593496,0.292857143,0.253731343,0.286624204,
0.21799308,0.290909091,0.258302583,0.194366197),.Dim=c(27,4)),
adj=c(13,17,6,12,22,5,10,15,20,23,25,27,7,15,3,23,27,2,20,22,24,25,26,4,15,9,13,14,16,18,8,13,15,16,3,15,16,25,12,14,19,21,2,11,19,22,1,8,9,17,8,11,17,18,19,21,3,4,7,9,10,16,23,8,9,10,15,18,25,1,8,13,14,21,8,14,16,19,24,25,11,12,14,18,22,24,3,6,25,26,11,14,17,2,6,12,19,24,3,5,15,6,18,19,22,25,3,6,10,16,18,20,24,6,20,3,5),
num=c(2,3,7,2,3,6,2,5,4,4,4,4,4,6,7,6,5,6,6,4,3,5,3,5,7,2,2),sumNumNeigh = 117)
#Inits
list(alpha0=0,tau=1,sigmay=0,sigmamu=0,sigma1=0,sigma2=0,v=c(0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0))