森林地块聚合值的元数据

Question

第一次发帖

我在 metafor 中运行了一个多层次的 m-a，就像 Wolfgang 的想法一样，创建一个森林图来显示每项研究的总效应大小（例如，100 项研究的 m-a，其中一些具有单一效应，有些有 5 个效果，有些有 10 个，但森林图只会显示每个研究的一个效果）。教程（这里：https://www.metafor-project.org/doku.php/tips:forest_plot_with_aggregated_values）很清楚，我已经将它用于我自己的数据。

然而，我不能做的是让我第一个情节中的研究标题显示为作者姓名，而不是显示为“研究 1”、“研究 2”等

有人可以指导我如何将该信息（即作者姓名，大概在我运行聚合函数时未能包含的额外列中）添加到聚合数据或“escalc”对象吗？

所以，这是我使用的代码：

#通话数据 数据_RR<- read.csv("MA_RR_Data.csv") #create a subset data_RR<- subset(data_RR, pool==1) #create an escalc object dat <- escalc(measure = "RR", ai = bi_anx, bi = bi_not_anx, ci = not_bi_anx, di = not_bi_not_anx, data = data_RR, append = TRUE, add = 1/2, to = "only0", slab=paste("Study", First_Author, "Estimate", es.id, "obs", Obs))

假设研究中的效应量与 rho=0.8

相关

V<- vcalc(vi, cluster=Paper.ID, obs=es.id, data=dat, rho=0.8)

使用这个近似的 V 矩阵拟合多级模型

资源<- rma.mv(yi, V, random = ~ 1 | First_Author/es.id, data=dat, digits=3) res

#aggregate 来自研究内部的数据 聚合<- aggregate(dat, cluster=First_Author, V=vcov(res, type="obs"), addk=TRUE) #limit 'agg' to columns with headings: First_Author (2), yi (42), vi (43), ki (44) agg <- agg[c(2 42, 43, 44)] agg #meta-analyses pooled estimates res <- rma(yi, vi, method="EE", data=agg, digits=3) res #draw a forest plot forest(res, xlim=c(-4,5), mlab="Pooled Estimate", header=TRUE, ilab=ki, ilab.xpos=-2)

一切顺利， 皮特

PS 我知道情节特征的格式很草率，但是一旦我有了作者姓名，我就可以通过摆弄距离来解决这部分问题！

Answer 1

您将不得不再次使用

slab

参数，这次是将模型拟合到聚合数据（例如

rma()

）时作为

rma(yi, vi, method="EE", data=agg, digits=3, slab=First_Author

的一部分。