Tau在metafor包中为零

时间:2018-10-30 18:40:58

标签: metafor

考虑以下代码,以在软件包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赋零?对此有什么解释?

0 个答案:

没有答案