### copy data into 'dat' and examine data
dat <- dat.anand1999
dat
if (FALSE) {
### load metafor package
library(metafor)
### High-Intensity OA vs Control
rma.mh(measure="OR", ai=ai, n1i=n1i, ci=ci, n2i=n2i, data=dat,
subset=(intensity=="high" & asp.t==0 & asp.c==0), digits=2)
### High- or Moderate-Intensity OA vs Aspirin
rma.mh(measure="OR", ai=ai, n1i=n1i, ci=ci, n2i=n2i, data=dat,
subset=(intensity %in% c("high","moderate") & asp.t==0 & asp.c==1), digits=2)
### Moderate-Intensity OA vs Control
rma.mh(measure="OR", ai=ai, n1i=n1i, ci=ci, n2i=n2i, data=dat,
subset=(intensity=="moderate" & asp.t==0 & asp.c==0), digits=2)
### High- or Moderate-Intensity OA and Aspirin vs Aspirin
rma.mh(measure="OR", ai=ai, n1i=n1i, ci=ci, n2i=n2i, data=dat,
subset=(intensity %in% c("high","moderate") & asp.t==1 & asp.c==1), digits=2)
### Low-Intensity OA and Aspirin vs Aspirin
rma.mh(measure="OR", ai=ai, n1i=n1i, ci=ci, n2i=n2i, data=dat,
subset=(intensity=="low" & asp.t==1 & asp.c==1), digits=2)
}
Run the code above in your browser using DataLab