如何在 metaanalysis of correlation中创建一个forestplot(关于xlim的错误信息)?

vzgqcmou  于 2023-05-04  发布在  其他
关注(0)|答案(1)|浏览(139)

我是R语言的初学者,我试图运行相关性的元分析,但当我试图将数据绘制成图时,我得到了一个错误。我正在学习代码:我真的很感激你能给予我任何帮助或建议。这是我的代码:

install.packages("metafor")
install.packages("robumeta")
install.packages("dplyr")

library(metafor)
library(robumeta)
library(dplyr)

Main<- read.csv("Meta-main.csv")
Main <- mutate(Main, study_id = 1:23)
Main <- Main %>% select(study_id, Authur:Quality)
View(Main)

Main <-escalc(measure = "ZCOR", ri=EffSize, ni=Ssize, data = Main, slab = paste(Authur, Year, sep = ", "))
View(Main)

res <- rma(yi, vi, data = Main)
res
predict(res, digits = 3, transf = transf.ztor)
confint(res)

b_res <- rma(yi, vi, data = Main, slab = study_id)
baujat(b_res)

inf<-influence(res)
print(inf)
plot(inf)

forest(res, xlim=c(-1.6,1.6), atransf = transf.ztor,
       at = transf.rtoz(c(-4,-2,0,.2,.4,.6)), digits = c(2,1),cex = .8)

这就是错误:
plot.window(...)中出错:需要有限的“xlim”值
我试着自己解x〈-1.6:1.6,但我真的不知道我在做什么。

dldeef67

dldeef671#

问题是at=transf.rtoz(c(-4,-2,0,.2,.4,.6))产生-InfInf值。我想你的意思是使用at=transf.rtoz(c(-.4,-.2,0,.2,.4,.6)),然后它应该工作。

相关问题