Program 2.16  Female Cancers

Model
{# Prior and likelihood: Multinomial response
 #      theta[] ~ ddirch(c[]);
 #       for (i in 1:I) { c[i] <- 1;}
 #       x[]  ~ dmulti( theta[] , X);

# Poisson analysis:priors
     for (i in 1:I) {  mu[i] ~ dgamma(0.001,0.001);
# likelihood
       x[i] ~ dpois(mu[i]);
# multinomial parameters
       theta[i] <- mu[i]/mu.sum}
       mu.sum <- sum(mu[])}  

National Model Data
list(x=c(14080,12990,6440,4350,3420,3190,2600,2420,1820,1760,23610), I   = 11,    # number of cancer types
   X   = 76680) # total deaths from cancer