# data from Table 3.14 of Breslow & Day Vol II, page 111

# 13 age-period strata

c0=c(2,0, 2,3,3, 2,7,10,1, 3,6,6,6)
c1=c(0,0, 4,3,1, 9,7, 3,1, 1,2,1,2)

pt0=c(3075.27, 936.75,
      2849.76,2195.59, 747.77,
      2085.43,1675.91,1501.73,440.21,
       833.61, 973.32,1027.12,674.44)

pt1=c(337.29,121.00,
      626.72,349.53,142.33,
      672.09,441.10,244.82,100.64,
      277.25,268.27,197.20, 92.75)
       
############

# conditional ... c1 | (c0+c1) ~ Binomial(c0+c1,p)
#  
#         p  = psi*pt1/(psi*pt1+pt0) 

#   Logit(p) = ... 

...

...

...


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

# unconditional 26 "cells" .. Poisson Regression 

p.t=c(pt0,pt1); cases=c(c0,c1); 
exposed = c(rep(0,13),rep(1,13)) ;
stratum=c(1:13,1:13)

...