Effect size
(lnRR=escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=ex1,append=TRUE,slab=paste(author,'(',year,')')) )
## trial author year tpos tneg cpos cneg yi vi
## 1 1 Aronson 1948 4 119 11 128 -0.8893 0.3256
## 2 2 Ferguson & Simes 1949 6 300 29 274 -1.5854 0.1946
## 3 3 Rosenthal et al 1960 3 228 11 209 -1.3481 0.4154
## 4 4 Hart & Sutherland 1977 62 13536 248 12619 -1.4416 0.0200
## 5 5 Frimodt-Moller et al 1973 33 5036 47 5761 -0.2175 0.0512
## 6 6 Stein & Aronson 1953 180 1361 372 1079 -0.7861 0.0069
## 7 7 Vandiviere et al 1973 8 2537 10 619 -1.6209 0.2230
## 8 8 TPT Madras 1980 505 87886 499 87892 0.0120 0.0040
## 9 9 Coetzee & Berjak 1968 29 7470 45 7232 -0.4694 0.0564
## 10 10 Rosenthal et al 1961 17 1699 65 1600 -1.3713 0.0730
## 11 11 Comstock et al 1974 186 50448 141 27197 -0.3394 0.0124
## 12 12 Comstock & Webster 1969 5 2493 3 2338 0.4459 0.5325
## 13 13 Comstock et al 1976 27 16886 29 17825 -0.0173 0.0714
Random effect model
lnRR=escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=ex1,append=TRUE,slab=paste(author,'(',year,')'))
result=rma(yi, vi, data=lnRR, method="HE",slab=paste(author,'(',year,')'))
forest(result,atransf=exp,cex=0.8,showweight=T)
text(8,14.5, "wight", pos=2, cex=0.75,font=4)
text(10.3,14.6, "\t\t Relative Risk\n[95% CI]", pos=2, cex=0.75,font=4)
text(-11,14.5, "Author(s) and Year", pos=4, cex=0.75,font=4)
Fix effect model
lnRR=escalc(measure="RR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=ex1,append=TRUE,slab=paste(author,'(',year,')'))
result=rma.mh(ai=tpos, bi=tneg, ci=cpos, di=cneg, data=ex1, measure="RR",slab=paste(author,'(',year,')')) #MH#
forest(result,atransf=exp,cex=0.8,showweight=T)
text(8,14.5, "wight", pos=2, cex=0.75,font=4)
text(10.3,14.6, "\t\t Relative Risk\n[95% CI]", pos=2, cex=0.75,font=4)
text(-11,14.5, "Author(s) and Year", pos=4, cex=0.75,font=4)
