我试图在metafor中创建一个funnel plot来测试是否存在发表偏倚。在漏斗图中，我试图从相同的研究中创建指示效应大小的相同颜色。我最近看到一篇论文，其中他们使用了以下代码应用了相同的方法：

mod1 <- rma(yi, vi, data = dat3)
pdf(NULL)
dev.control(displaylist="enable")
par(mar=c(4,4,0.1,0) + 0.1)
funnel(mod1, col =inferno(25)[pos],
       ylim = c(0, 0.26), xlim = c(-0.6, 0.8),
       ylab = "Standard error (SE)", xlab = "Effect size (lnRR)")

我只是想知道在col中，inferno(25)[pos]意味着什么。有人知道如何创建这样一个漏斗图，其中相同的颜色表示来自相同研究的效应量吗？