SMR

library(epitools)
o=2500; e=1250; scale=1; 
ageadjust.indirect(count=o,pop=1,stdcount=scale*e,stdpop=scale,stdrate=e)
c((o-1.96*sqrt(o))/e,(o+1.96*sqrt(o))/e)

Mantel-Haenszel

library(rmeta)
c1=c( 1, 0, 1, 0, 1, 0, 0, 1); 
c0=c( 0, 1, 0, 1, 0, 1, 1, 0); 
n1=c(24,10,17, 7,10,10, 7, 8) ; 
n0=c(48,14,32,16,15,16,18,16);
meta.MH(n1,n0,c1,c0)