########################################################################

# Poisson Distribution.

# Hospital A

########################################################################

model

{

y ~ dpois(lambda.y) # likelihood

lambda.y <- eo*lambda

lambda ~ dgamma(alpha, beta) # prior

}

# Data

list( y=1, eo=66, alpha=16, beta=15174)

# Initial values

list( lambda = 0.001)