1

我正在进行一项荟萃分析，检查三年级阅读理解干预的效果。我正在分析标准化的平均差异。我有 20 项研究，有 54 种影响（多种结果测量或治疗比较/研究）。鉴于样本量小且效果的相关性，我使用以下代码进行初步分析：

library(robumeta)
library(clubSandwich)
data (X3rdmetawithcolors)

CEmodel <- rma.mv(yi=adjg,
                  V = var_adjg, #variance
                  random = ~1 | study/eid, 
                  test= "t", 
                  data=X3rdmetaforR2, 
                  method="REML")
summary(CEmodel)
cecf <- coef_test(CEmodel,
                  cluster=X3rdmetaforR2$study_id, 
                  vcov = "CR2")
cecf

我也想做敏感性分析。我知道不建议将漏斗图和修剪和填充程序用于多级数据。但是，我想生成一个漏斗图，其中包含按研究着色的各个点。我有一个包含颜色、study_id 和内部引用的数据框，我已将其合并到我的主数据框中。我尝试使用以下代码来获取漏斗图：

res<-rma.mv(yi= adjg, 
      V = var_adjg, 
      random = ~1|study/eid, 
      data=X3rdmetacolor_intcite,
      method="REML")
funnel (res, col=X3rdmetacolor_intcite$color)

但是，这会产生一个漏斗图，其中所有点都以第一项研究的颜色（海蓝宝石）而不是不同的颜色着色。

关于我能做什么的任何想法？我还想创建一个显示颜色和内部引用的图例。

4

0 回答 0