model {
for(i in 1:N) {
   h[i] ~ dpois(lambda[i])
   mu[i] <- b0 + b1*SOI[i] + b2*SST[i] + b3*RI[i] + eta[i]
   lambda[i] <- exp(mu[i])
   tt[i] <- tau[RI[i] + 1]
   eta[i] ~ dnorm(0, tt[i])
   }
  b0 ~ dnorm(0, .0001)
  b1 ~ dnorm(0, .0001)
  b2 ~ dnorm(0, .0001)
  b3 <- log(pm)
  pm ~ dunif(lo, hi)
  tau[1] ~ dgamma(.001, .001)
  tau[2] ~ dgamma(.001, .001)
}