n.i.1 <- c( 20, 208, 24, 190, 58, 36, 51)
a.i <- c( 2,79,0,98,15,34,9)
n.i.2 <- c( 20, 119, 22, 185, 29, 51, 47)
c.i <- c(9,106,14,98,12,49,9)
m1 <- metabin( event.e = a.i,n.e = n.i.1,
event.c = c.i,n.c = n.i.2,
studlab = paste('Study',1:7), sm = 'OR',
common = FALSE, random = TRUE )
find_umax(m1 , common.effect = FALSE, alternative = 'two-sided',
t = 0.05 , confidence = 0.95 )
Run the code above in your browser using DataLab