meta包metacont()函数忽略自定义权重问题求助

在使用meta包的metacont()做连续数据的meta分析时，尝试通过w.random参数指定自定义权重，但发现权重并未生效，默认权重和指定后的权重完全一致：

> anx_ma$w.random
[1] 0.2942658 0.2919072 0.2951889 0.2941619 0.2855660 0.2957122
> anx_ma_weighted$w.random
[1] 0.2942658 0.2919072 0.2951889 0.2941619 0.2855660 0.2957122

可复现代码如下：

label <- c(
  "AAA",
  "BBB",
  "CCC",
  "DDD",
  "EEE",
  "FFF"
)

num.tre <- c(39, 40, 45, 35, 59, 50)
anxiety.treat.mean <- c(1.46, 1.35, 5.32, 1.74, 1.94, 5.10)
anxiety.treat.sd   <- c(0.51, 0.60, 1.70, 2.62, 1.41, 1.50)
num.cont <- c(39, 40, 45, 35, 52, 49)
anxiety.cont.mean <- c(2.07, 2.82, 6.38, 2.74, 9.19, 5.50)
anxiety.cont.sd   <- c(0.73, 0.66, 1.56, 1.66, 1.37, 1.20)
new.weights <- c(0.5, 0.1, 0.1, 0.3, 0.2, 0.2)

# Combine into a data frame
DB0 <- data.frame(
  label = label,
  num.tre = num.tre,
  anxiety.treat.mean = anxiety.treat.mean,
  anxiety.treat.sd = anxiety.treat.sd,
  num.cont = num.cont,
  anxiety.cont.mean = anxiety.cont.mean,
  anxiety.cont.sd = anxiety.cont.sd,
  new.weights = new.weights,
  stringsAsFactors = FALSE
)

### 默认权重模型
anx_ma <- meta::metacont(num.tre, anxiety.treat.mean, anxiety.treat.sd,
                                num.cont, anxiety.cont.mean, anxiety.cont.sd,
                                data = DB0, studlab = label,
                                common = FALSE,
                                sm = "SMD", method.tau = "REML")

问题原因 #

metacont()函数中的w.random参数是内部用于存储计算出的随机效应权重，并非用来接收用户自定义权重的参数。当你传入自定义值时，函数会忽略该输入，依然按照默认逻辑计算随机效应权重，因此权重不会改变。

如果需要使用自定义权重，需要使用meta包的底层函数metagen()，它支持直接指定研究的权重值。

解决方法 #

1. 先通过metacont()计算出每个研究的效应量（SMD）及其标准误；
2. 调用metagen()，传入效应量、标准误以及自定义权重；
3. 自定义权重需满足非负，推荐标准化为总和为1（函数也支持原始权重值，会自动处理）。

修改后的代码示例：

# 第一步：用metacont计算效应量和标准误
anx_ma <- meta::metacont(num.tre, anxiety.treat.mean, anxiety.treat.sd,
                                num.cont, anxiety.cont.mean, anxiety.cont.sd,
                                data = DB0, studlab = label,
                                common = FALSE,
                                sm = "SMD", method.tau = "REML")

# 提取效应量和标准误
smd <- anx_ma$TE
se_smd <- anx_ma$seTE

# 第二步：标准化自定义权重（可选，确保权重和为1）
std_weights <- DB0$new.weights / sum(DB0$new.weights)

# 使用metagen指定自定义权重
anx_ma_weighted <- meta::metagen(
  TE = smd,
  seTE = se_smd,
  studlab = DB0$label,
  data = DB0,
  common = FALSE,
  random = TRUE,
  method.tau = "REML",
  w.random = std_weights
)

# 查看自定义权重是否生效
anx_ma_weighted$w.random

运行后即可看到w.random变为自定义的标准化权重：

> anx_ma_weighted$w.random
[1] 0.3571429 0.0714286 0.0714286 0.2142857 0.1428571 0.1428571

内容的提问来源于stack exchange，提问作者cccnrc