> bay<-as.data.frame(scan("bay.dat",
+    what=list(wlc=0,wrc=0,wpy=0,blc=0,brc=0,bpy=0)))
> bay
   wlc  wrc     wpy  blc brc     bpy 
1 1429 1082 1625812 1006 999 1767995
2 1825 1394 1437511  809 913  996536
3 1048  667  565078  258 212  235442
4  484  273  229203  245 162  171292
5  176   87   84698   44  38   31789
> standard<-as.data.frame(list(
+    rate=(bay$wlc+bay$blc)/(bay$wpy+bay$bpy),
+    py=bay$wpy+bay$bpy))
> standard
          rate      py 
1 0.0007174833 3393807
2 0.0010821484 2434047
3 0.0016314396  800520
4 0.0018202474  400495
5 0.0018886228  116487
> expected<-standard$rate*bay$wpy
> smr<-sum(bay$wlc)/sum(expected)
> smr
[1] 1.175508
> crude.standard.rate<-sum(standard$py*standard$rate)/
+    sum(standard$py)
> crude.standard.rate
[1] 0.001025001
> indirect.standardized.rate<-crude.standard.rate*smr
> indirect.standardized.rate
[1] 0.001204898
> direct.standardized.rate<-sum((bay$wlc/bay$wpy)*standard$py
+    )/sum(standard$py)   
> direct.standardized.rate
[1] 0.001209954
> smr.ci<-smr+c(-1,1)*sqrt(smr)*1.96/sqrt(sum(expected))
> smr.ci
[1] 1.142800 1.208216
> smr.better.ci<-smr*exp(1.96*c(-1,1)/sqrt(sum(bay$wlc)))
> smr.better.ci
[1] 1.143251 1.208676