# Analyze data in Table 5.2 B&D Vol II (BrCa estrogen) pages 187-190

ds=read.table("BD-II-5-2-BrCa-estrogen.txt",header=T)
ds
attach(ds)

nc.H0=nc.H0+nc.H2

# M-H summary OR

  sum(c.H1*nc.H0/n) / sum(c.H0*nc.H1/n)
  
# Conditional MLE

 ... without the age covariate
 
 ... with the age covariate ???