#  Confidence interval for a marginal odds
#  ratio in a 1-1 matched pair study

x<-matrix(c(9,37,16,82),ncol=2, byrow=T)
xs1<-sum(x[1,])
xs2<-sum(x[2,])
x1s<-sum(x[,1])
x2s<-sum(x[,2])
xss<-sum(x)
or<-xs1*x2s/(xs2*x1s)
logor<-log(or)
selogor<-sqrt(1/(xs1*xs2)+1/(x1s*x2s)
              -2*(x[1,1]*x[2,2]-x[1,2]*x[2,1])/
                 (xs1*xs2*x1s*x2s))
logcil<-logor-1.96*selogor
logciu<-logor+1.96*selogor
cil<-exp(logcil)
ciu<-exp(logciu)

or
cil
ciu