Вычисление оценок с помощью rma.glmm: ошибка, когда 2 оценки идентичны

Question

Я использую rma.glmm для мета пропорций из 2 разных исследований. Например, доля лиц с (xi) и без (xii) нежелательным явлением через 6 месяцев:

library(metafor)

#6 months
study=c("1", "2")
ni = c(41, 19)
xi = c(26, 14)
xii = c(15, 5)
NP6monthAT <- data.frame(study, xi, xii, ni)
res2 <- rma.glmm(measure="PLO", xi = xi, ni = ni, data = NP6monthAT)
predict (res2, transf = transf.ilogit, digits=2)
res2 <- rma.glmm(measure = "PLO", xi = xii, ni = ni, data = NP6monthAT)
predict (res2, transf = transf.ilogit, digits=2)

ОДНАКО, случайно в определенный момент времени обе пропорции в 2 разных популяциях идентичны (11,1%), и я получаю ошибку:

#12 months
study=c("1", "2")
ni=c(81, 45)
xi=c(9, 5)
xii=c(72, 40)
NNPNNP12monthAT<-data.frame(study, xi, xii, ni)
res4<-rma.glmm(measure="PLO", xi=xi, ni=ni, data=NNPNNP12monthAT)
predict(res4, transf=transf.ilogit, digits=2)

Ошибка в файле rma.glmm(measure = "PLO", xi = xi, ni = ni, data = NNPNNP12monthAT): не может соответствовать модели ML.

Я понимаю, что оценка будет равна 11,1 (поскольку это то, что происходит в обеих группах)... но я хочу получить результат с доверительным интервалом, какие-либо советы о том, что я могу сделать, чтобы получить эту информацию?

Answer 1

Эта проблема проистекает из lme4::glmer() который используется, чтобы соответствовать этой модели:

> glmer(cbind(xi,ni-xi) ~ 1 + (1 | study), family=binomial, data=NNPNNP12monthAT)
Error: Response is constant

Очевидно, что не может быть какой-либо неоднородности между исследованиями, когда шансы во всех исследованиях одинаковы. Простое решение состоит в том, чтобы соответствовать модели FE (так как с tau^2 равно 0, это то, что вы получите в любом случае). Итак, используйте:

res4 <- rma.glmm(measure="PLO", xi=xi, ni=ni, data=NNPNNP12monthAT, method="FE")