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

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

#inits;
list(theta=0)