# NOT RUN {
### copy data into 'dat' and examine data
dat <- dat.nielweise2008
dat
### standard (inverse-variance) random-effects model
res <- rma(measure="IRR", x1i=x1i, t1i=t1i, x2i=x2i, t2i=t2i, data=dat)
print(res, digits=3)
predict(res, transf=exp, digits=2)
### random-effects conditional Poisson model
# }
# NOT RUN {
res <- rma.glmm(measure="IRR", x1i=x1i, t1i=t1i, x2i=x2i, t2i=t2i, data=dat, model="CM.EL")
print(res, digits=3)
predict(res, transf=exp, digits=2)
# }
Run the code above in your browser using DataLab