做贝叶斯时空交互模型时,有一事不明,特来请教各位大神!
for(i in 1:N) {
Onset[i] ~ dpois(mu[i])
log(mu[i]) <- offs[i]
+ beta[1] * cons[i]
+ beta[2] *x[i]
+ u1[county[i]]*equals(year[i],1) + u2[county[i]]*equals(year[i],2)+ u3[county[i]]*equals(year[i],3)
}
# CAR prior distribution for random effects
u1[1:181] ~ car.normal(adj[],weights[],num[],tau.u[1])
u2[1:181] ~ car.normal(adj[],weights[],num[],tau.u[2])
u3[1:181] ~ car.normal(adj[],weights[],num[],tau.u[3])
这个模型中,使用winbugs只能得到beta、tau.u,sigma,但相关文章中报告了不同年代的空间相关系数,不知道是哪个参数?特此请教!