我有用于荟萃分析的数据集（df），看起来像这样。这是关于夜间与白天的认知测试。有些研究有多种结果，有些则只有一种。右侧，计算了每个研究结果的对冲 g (yi) 和方差 (vi)。
| 研究名称 | 时间_点 | 进入会话的时间 | 实验室领域 | 测试 | 测试结果 | 易 | vi |
|---|---|---|---|---|---|---|---|
| 研究1 | 23:00-03:00 | 2 | 1 | PVT | RT | -0.32 | 0.03 |
| 研究1 | 03:01-06:00 | 5 | 1 | PVT | RT | -0.27 | 0.03 |
| 研究1 | 06:01-06:59 | 7 | 1 | PVT | RT | 0.03 | 0.04 |
| 研究1 | 23:00-03:00 | 2 | 1 | PVT | 失效 | -0.88 | 0.16 |
| 研究1 | 03:01-06:00 | 5 | 1 | PVT | 失效 | -0.48 | 0.12 |
| 研究1 | 06:01-06:59 | 7 | 1 | PVT | 失效 | 0.10 | 0.05 |
| 研究2 | 23:00-03:00 | 3 | 1 | DSST | 正确 | -0.09 | 0.01 |
| 研究2 | 06:01-06:59 | 8 | 1 | DSST | 正确 | 0.22 | 0.10 |
| 研究3 | 06:01-06:59 | 4 | 2 | STM | 正确 | -1.53 | 0.49 |
| 研究 4 | 03:01-06:00 | 5 | 1 | PVT | RT | -0.71 | 0.15 |
由于一些研究有多个报告结果，我想将每项研究的数据汇总为一个效应大小。但是，当我这样做时，我收到以下错误（为了简单起见，我使用 struct="ID"）：
gg_studycluster <- aggregate(df, cluster=Study_Name, struct="ID", addk=TRUE) Error in aggregate.data.frame(df, cluster = Study_Name, struct = "ID", : argument "FUN" is missing, with no default
这应该与以下示例中的操作完全一样：https://wviechtb.github.io/metafor/reference/aggregate.escalc.html
大家知道是什么原因导致这个问题吗？
我已经使用 str(df) 检查了变量类型，其中 yi 和 vi 都是数字。
如果有人遇到同样的问题：
我想出了如何“解决”这个问题。我改用“MAd”包。
gg_studycluster <- agg(id = Study_name, es = yi, var = vi, method = "BHHR", cor = 0.50, data = df)
（显示适合我的数据的聚合方法）。