# Table 4.3 Daily Alcohol Consumption / Esophageal Ca.

# Breslow and Day Volume I page 137

# jh 2008.02.20

# c1 and c0: numbers of exposed/unexposed cases

# d1 and d0: numbers of exposed/unexposed in denominator series

#     Alcohol 1: > 80 g/day 0: < 80

# 6 age strata, mid ages 30:80

c1=c(1,4,25,42,19,5)
c0=c(0,5,21,34,36,8)

d1=c(  9, 26, 29, 27,18, 0)
d0=c(106,164,138,139,88,31)

n=c1+c0+d1+d0

# m-h

sum(c1*d0/n)/sum(c0*d1/n)

# woolf

# unconditional logistic regression

# ...


#####  ML fitting

## unconditional, via GLM

exposure=c(rep(1,6),rep(0,6))
cases=c(c1,c0);
controls=c(d1,d0);
strata=c(1:6,1:6)

...

  
## conditional, via GLM