bay<-as.data.frame(scan("bay.dat",
   what=list(wlc=0,wrc=0,wpy=0,blc=0,brc=0,bpy=0)))
standard<-as.data.frame(list(rate=(bay$wlc+bay$blc)/(bay$wpy+bay$bpy),
   py=bay$wpy+bay$bpy))
expected<-standard$rate*bay$wpy
smr<-sum(bay$wlc)/sum(expected)
smr
crude.standard.rate<-sum(standard$py*standard$rate)/sum(standard$py)
crude.standard.rate
indirect.standardized.rate<-crude.standard.rate*smr
indirect.standardized.rate
direct.standardized.rate<-sum((bay$wlc/bay$wpy)*standard$py
   )/sum(standard$py)
direct.standardized.rate
smr.ci<-smr+c(-1,1)*sqrt(smr)*1.96/sqrt(sum(expected))
smr.ci
smr.better.ci<-smr*exp(1.96*c(-1,1)/sqrt(sum(bay$wlc)))
smr.better.ci