2015-02-23 56 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. 

所以,本質上,你應該仔細檢查變量名稱。