基于空间滞后模型
y=ρWy+Xβ+ε
在WinBUGS中建一个模型如下：
model {
# Likelihood
for (i in 1 : N) {
  O ~ dpois(mu)
  log(mu) <-  alpha0 + alpha1 * X/10 + b
  RR <- exp(alpha0 + alpha1 * X/10 + b) # Area-specific relative risk (for maps)
}
# CAR prior distribution for random effects:
b[1:N] ~ car.normal(adj[], weights[], num[], tau)
for(k in 1:sumNumNeigh) {

  weights[k] <- 1
}
# Other priors:
alpha0 ~ dflat()
alpha1 ~ dnorm(0.0, 1.0E-5)
tau ~ dgamma(0.5, 0.0005) # prior on precision
sigma <- sqrt(1 / tau) # standard deviation
}

我的理解，b对应的是方程  y=ρWy+Xβ+ε  中的ρWy，对否？

如果是这样，那么如何得到ρ（也就是空间自相关系数）的值呢？