# CAR, JAGS
model 
{
  for (j in 1:N.site){
    Y[j] ~ dpois(lambda[j])
    lambda[j] <- exp(r[j])
  }
  for (j in 2:N.site){
    r[j] ~ dnorm(r[j - 1], tau)
  }
  r[1] ~ dunif(0, 10)
  tau <- 1 / (s * s)
  s ~ dunif(0, 1.0E+4)
}