n1 <- c(204, 201, 932, 130, 77)
n2 <- c(106, 103, 415, 132, 83)
f1 <- c(24, 40, 93, 14, 5)
f2 <- c(12, 9, 28, 3, 1)
x1 <- c(4, 4, 5, 3, 26)
x2 <- c(1, 1, 1, 0, 0)
X <- matrix(cbind(x1, x2), 5, 2)
meta.lm.oddsratio(.05, f1, f2, n1, n2, X)
# Should return:
# Estimate SE z p LL UL
# b0 1.541895013 0.69815801 2.20851868 0.027 0.1735305 2.91025958
# b1 -0.004417932 0.04840623 -0.09126784 0.927 -0.0992924 0.09045653
# b2 -1.071122269 0.60582695 -1.76803337 0.077 -2.2585213 0.11627674
# exp(Estimate) exp(LL) exp(UL)
# b0 4.6734381 1.1894969 18.361564
# b1 0.9955918 0.9054779 1.094674
# b2 0.3426238 0.1045049 1.123307
Run the code above in your browser using DataLab