R中的漏斗图和regtest

Question

我正在 R 中学习元分析。我想从 metafor 库中计算 regtest。 请参阅我的数据示例：

library("metafor")

value <- c(0.0001,0.0003,0.0002,0.0009,0.01,0.002)
se <- c(1,3,2,4,0,3)
data <- data.frame(value, se)

regtest(x = data$value, sei = data$se, model="rma")

但是我Error in rma.uni(yi, vi, mods = X, intercept = FALSE, ...) : Division by zero when computing the inverse variance weights. In addition: Warning message: There are outcomes with non-positive sampling variances.得到一个错误Error in rma.uni(yi, vi, mods = X, intercept = FALSE, ...) : Division by zero when computing the inverse variance weights. In addition: Warning message: There are outcomes with non-positive sampling variances. Error in rma.uni(yi, vi, mods = X, intercept = FALSE, ...) : Division by zero when computing the inverse variance weights. In addition: Warning message: There are outcomes with non-positive sampling variances.

似乎值太小了，但我该如何解决？

Answer 1

不要再除以零了 ;)

value <- c(0.0001,0.0003,0.0002,0.0009,0.01,0.002)
se <- c(1,3,2,4,1,3)
data <- data.frame(value, se)

regtest(x = data$value, sei = data$se, model="rma")

---

Regression Test for Funnel Plot Asymmetry

Model:     mixed-effects meta-regression model
Predictor: standard error

Test for Funnel Plot Asymmetry: z = -0.0025, p = 0.9980
Limit Estimate (as sei -> 0):   b = 0.0069 (CI: -2.5466, 2.5605)