老师你好,我写了个循环画图的代码,可是现在图上面需要加上text函数可是只有最后一个图有text函数我想问下有什么解决的办法么

# 读取数据

data_1051312 <- read.csv("RS1051312.csv", header = T)

data_362549 <- read.csv("RS362549.csv", header = T)

data_362998 <- read.csv("RS362998.csv", header =T)

data_363006 <- read.csv("RS363006.csv", header = T)

data_6077690 <- read.csv("RS6077690.csv", header = T)

data_8636 <- read.csv("RS8636.csv", header =T)

data_RS3746544 <- read.csv("RS3746544.csv", header = T)



# 导入metafor包

library(metafor)



data <- list(data_1051312, data_362549, data_362998,

          data_363006, data_6077690,  data_8636, data_RS3746544)

library(tidyverse)

ouput <- list()

q <- list()

p <- list()

regtest1 <- list()

ranktest1 <- list()

leve <-list()

par(op)

for(i in 1:7){

  # 计算出escalc函数结果

  ouput[[i]] <- escalc(measure="RR", 

                  ai=tpos, bi=tneg, ci=cpos, di=cneg, data=data[i])

  #计算rma函数结果

  q[[i]]<-rma(measure="OR", ai=tpos, bi=tneg, ci=cpos, di=cneg, data=data[i])

  #得到森林图

  forest(q[[i]],atransf = exp,slab=paste(data[[i]]$author,data[[i]]$year,sep=","),

         annotate = F)

  就是这里我需要加上

text(-16,                15, "Author(s) and Year",   pos=4)

text(6,                  15, "Risk Ratio [95% CI]", pos=2)

可是只有最后一个图加上了 我想老师帮我完善一下这个循环

  # 得到漏斗图

  p[[i]] <- funnel(q[[i]])

  # 得到Egger检验

  regtest1[[i]] <- regtest(q[[i]])

  # 得到Begg检验

  ranktest1[[i]] <- ranktest(q[[i]])

}

请先 登录 后评论

1 个回答

王莹 - 基因课管家

forest函数前面加上pdf(file=paste( i, "pdf", sep = ".") )函数创建pdf文件, 再绘图代码结束时,加上dev.off()关闭pdf画布。

请先 登录 后评论
  • 1 关注
  • 0 收藏,31 浏览
  • 王平 提出于 2019-11-29 20:02

相似问题