model{
for(i in 1:n)
{
y[i] ~ dpois(theta)
}
theta ~ dgamma(1,0.1)
}
list(y=c(6,5,4,6,6,3,12,7,4,2,6,7,4),n=13)
list(theta=1)

# Table 2.1 example model with two theta's
model{
for(i in 1:k1)
{
y1[i] ~ dpois(theta1)
}

for(i in 1:k2)
{
y2[i] ~ dpois(theta2)
}
theta1~dgamma(1,0.1)
theta2~dgamma(1,0.1)
diff <- theta2-theta1
}
list(y1=c(0,1,1,0,2,3,0,1,1,1,1,2,0,1,3,0,1,2,1,3,3,4,1,3,2,0), 
y2=c(2,0,3,0,0,1,1,1,1,0,0,2,2,0,1,2,0,0,1,1,1,0,2), k1=26,k2=23)
list(theta1=0.5,theta2=0.5)

# Table 2.1 model with one theta
model{
for(i in 1:k1)
{
y1[i] ~ dpois(theta)
}
for(i in 1:k2)
{
y2[i] ~ dpois(theta)
}
theta~dgamma(1,0.1)
}
list(y1=c(0,1,1,0,2,3,0,1,1,1,1,2,0,1,3,0,1,2,1,3,3,4,1,3,2,0), 
y2=c(2,0,3,0,0,1,1,1,1,0,0,2,2,0,1,2,0,0,1,1,1,0,2), k1=26,k2=23)
list(theta=0.5)