2015-02-23 58 views
0

我以我想要的方式计算了效果大小和汇总SE。唯一的事情是绘制一个森林情节,让metafor计算总结效果的大小。我有超过30个.csv数据文件分开绘制。当我用下面的数据做到这一点时(如下),它平滑地绘制和计算汇总效果。Metafor测量参数错误

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给出了一个错误,并要求“措施”参数绘制输出。由于该手段已经是我手动计算的DeltaPI,我不想使用。

奇怪的是,即使我改变那些不工作的文件中的数据与正常工作的文件(上面的3个数据行),它仍然会给出相同的错误。虽然,相同的数据在其他.csv文件中正常工作。

所以我不清楚为什么我得到的错误和解决方案是什么。

如有任何意见,将不胜感激!

回答

0

我的猜测是,这与绘图无关,但在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. 

所以,本质上,你应该仔细检查变量名称。