Question

我正在尝试对7项随机研究进行广义线性混合效应（二项式-正态）荟萃分析，其中每项研究都记录了治疗和安慰剂人群（暴露和对照）中不良事件的存在。 / p>

为此，我希望使用metabin function（元软件包）。但是，我遇到了一个错误，并且不确定为什么。例如。运行此代码：

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进行此操作），但是会得到相同的错误消息。我也尝试阅读source code for 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")

这是一个收敛问题，有人知道是否有解决办法？关于此错误消息，我唯一能找到的是a problem (and thus solution) which does not apply to me。

任何帮助将不胜感激：）

GLMM进行荟萃分析-使用Metabin出错

0 个答案: