model {
    theta ~ dnorm(0.0,1.0E-6)
  for (i in 1:n) {
    zeros[i] <- 0
    zeros[i] ~ dpois(phi[i])
    phi[i] <- log(1+pow(x[i]-theta,2))
  }
}

data;
list(n=9,x=c(-0.774,0.597,7.575,0.397,-0.865,-0.318,-0.125,0.961,1.039))

inits;
list(theta=0)