Question

我已经按照我想要的方式计算了效果大小并汇集了SE。唯一的事情是绘制森林图并让metafor计算摘要效果大小。我有超过30个.csv数据文件分别绘制。当我使用以下数据（如下）时，它可以平滑地绘制和计算摘要效果。

DeltaPI            Spooled
-75.35224985    7.618629848
-51.85221078    7.513461236
-37.77455275    7.164279414

我使用的是：

meta1<-rma(yi=mydata$DeltaPI, sei=mydata$Spooled)
forest(meta1,slab=paste(mydata$Study,mydata$Genotype..Experimental.),showweight=TRUE,alim=c(-100,25),at=c(-100,-50,0,25),xlab="Percentage Change of PI Score",cex=0.7,cex.lab=1,col="red")

然而，当我尝试用其他一些.csv文件做同样的事情时，rma会给出一个错误，并要求'measure'参数来绘制输出。由于该度量已经是手动计算的DeltaPI，我不想使用。

奇怪的是，即使我更改那些工作正常的文件（上面3个数据行）中的数据，它仍然会给出相同的错误。虽然，相同的数据在其他一些.csv文件中也能正常工作。

所以我不清楚为什么我会收到错误，解决方案是什么。

任何评论都将不胜感激！

Answer 1

我的猜测是这与绘图无关，但是在运行rma()命令时会发生。在我看来，在您正在阅读的数据中如何命名变量存在问题。现在您正在读取.csv文件中的数据，但这可能是正在发生的事情：

> library(metafor)
> dat <- data.frame(DeltaP1 = c(.2,.4), Spooled=c(.1,.1))
> rma(dat$DeltaPI, sei=dat$Spooled)

Error in rma(dat$DeltaPi, sei = dat$s) : 
  Specify the desired outcome measure via the 'measure' argument.

因此，实质上，您应该仔细检查变量名称。

Metafor测量参数错误

1 个答案: