Рассмотрим код ниже для запуска мета-анализа в пакете 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")
Таким образом, значение tau2 равно нулю res$tau2
. Я пробовал другие методы, такие как "ML"
, "DL"
, "HE"
, "EB"
и "HS"
, чтобы проверить результаты, и все они дают одинаковое нулевое значение для тау. Единственными методами с ненулевым значением для тау являются "REML"
и "SJ"
, которые очень близки к нулю. Почему упомянутые методы дают ноль для тау? Чем это объясняется?