Я пытаюсь провести обобщенный метаанализ линейных смешанных эффектов (биномиальный-нормальный) для 7 рандомизированных исследований, где каждое исследование регистрирует наличие неблагоприятного события в группах лечения и плацебо (воздействие и контроль).
Для этого я надеюсь использовать метабиновую функцию (метапакет). Тем не менее, я получаю сообщение об ошибке, и я не уверен, почему. Например. работает этот код:
install.packages('meta')
# Data
data<-data.frame(exposure.events=c(11,34,152,4,60,3,25), exposure.population=c(184,152,9500,77,2012,15,60), control.events=c(3,33,4729,133,1441,1,25), control.population=c(184,375,613978,15865,480485,105,238), Study=c("1","2","3","4","5","6","7"))
# Calling metabin
metabin(event.e=exposure.events, n.e=exposure.population, event.c=control.events, n.c=control.population, studlab=Study, data=data, method="GLMM",model.glmm = "CM.AL",method.tau = "ML")
Я получаю этот вывод:
Ошибка в metafor :: rma.glmm (ai = event.e [! Exclude], n1i = n.e [! Exclude],:
Не подходит для модели ML.
Я также пытался вызвать функцию rma.glmm напрямую (вместо того, чтобы делать это через metabin), но получил то же сообщение об ошибке. Я также попытался прочитать исходный код для rma.glmm , но я не уверен, что понимаю, что происходит. Тем не менее, я думаю, что проблема связана с третьим (самым крупным) исследованием, и, в частности, с размером контрольной группы населения, так как оба следующих шага проходят гладко:
# Modifying 3rd row's control population
data<-data.frame(exposure.events=c(11,34,152,4,60,3,25), exposure.population=c(184,152,9500,77,2012,15,60), control.events=c(3,33,4729,133,1441,1,25), control.population=c(184,375,61378,15865,480485,105,238), Study=c("1","2","3","4","5","6","7"))
metabin(event.e=exposure.events, n.e=exposure.population, event.c=control.events, n.c=control.population, studlab=Study, data=data, method="GLMM",model.glmm = "CM.AL",method.tau = "ML")
# Deleting 3rd row
data<-data.frame(exposure.events=c(11,34,4,60,3,25), exposure.population=c(184,152,77,2012,15,60), control.events=c(3,33,133,1441,1,25), control.population=c(184,375,15865,480485,105,238), Study=c("1","2","3","4","5","6"))
metabin(event.e=exposure.events, n.e=exposure.population, event.c=control.events, n.c=control.population, studlab=Study, data=data, method="GLMM",model.glmm = "CM.AL",method.tau = "ML")
Это проблема сходимости, и кто-нибудь знает, есть ли способ обойти это? Единственное, что я могу найти в этом сообщении об ошибке, это проблема (и, следовательно, решение), которая не относится ко мне .
Любая помощь будет очень признательна :) 1023 *