我是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,但我真的不知道我在做什么。
1条答案
按热度按时间dldeef671#
问题是
at=transf.rtoz(c(-4,-2,0,.2,.4,.6))
产生-Inf
和Inf
值。我想你的意思是使用at=transf.rtoz(c(-.4,-.2,0,.2,.4,.6))
,然后它应该工作。