考虑以下代码,以在软件包metafor
中运行元分析。首先是数据:
TP1 = scan()
105 54 152 22 145 13
FP1 = scan()
65 11 49 40 51 14
FN1 = scan()
62 59 91 10 50 9
TN1 = scan()
131 89 134 47 54 37
data = data.frame("TP"=TP1, "FP"=FP1, "FN"=FN1, "TN"=TN1)
然后运行元分析:
dat <- escalc(measure = "RR", ai = TP, bi = FP, ci = FN, di = TN,
data = data, append =TRUE)
res <- rma(measure = "RR", ai = TP, bi = FP, ci = FN, di = TN,
data = data, method="PM")
因此res$tau2
,tau2的值为零。我尝试了其他方法,例如"ML"
,"DL"
,"HE"
,"EB"
和"HS"
来测试结果,并且都给出相同的tau零值。 tau唯一具有非零值的方法是"REML"
和"SJ"
,它们非常接近零。为什么上述方法给tau赋零?对此有什么解释?