Я следовал примеру с этого сайта о том, как создать лесной участок (https://designdatadecisions.wordpress.com/2016/07/02/forest-plot-with-horizontal-bands/). Однако, это не работает, и я не уверен, почему, поскольку сообщение об ошибке довольно неоднозначно, говоря, что есть неиспользуемые аргументы
top_results_poster <- read_excel("Z:/mr_top_results_poster.xlsx")
tabletext<-cbind(c("Trait","\n",top_results_poster$TRAIT),
c("No.SNPs","\n",top_results_poster$SNPs),
c("p for Egger","\n",top_results_poster$`p for Egger`),
c("p for Q statistic", "\n",top_results_poster$`P for Q statistic`))
library(forestplot)
png(file.path("Forestplot.png"),width=960, height=640)
forestplot(labeltext=tabletext, graph.pos=3,
mean=c(NA,NA,top_results_poster$mean),
lower=c(NA,NA,top_results_poster$LOWER),
upper=c(NA,NA,top_results_poster$UPPER),
title=" ",
xlab="Odds ratio (95% CI)",
txt_gp=fpTxtGp(label=gpar(cex=1.25),
ticks=gpar(cex=1.1),
xlab=gpar(cex=1.2)),
zero=1, cex=0.9, boxsize=0.5, colgap=unit(6,"mm"),
lwd.ci=2, ci.vertices=TRUE, ci.vertices.height = 0.4)
Error in forestplot(labeltext = tabletext, graph.pos = 3, mean = c(NA, :
unused arguments (labeltext = tabletext, graph.pos = 3, mean = c(NA, NA, top_results_poster$mean), lower = c(NA, NA, top_results_poster$LOWER), upper = c(NA, NA, top_results_poster$UPPER), title = " ", txt_gp = fpTxtGp(label = gpar(cex = 1.25), ticks = gpar(cex = 1.1), xlab = gpar(cex = 1.2)), zero = 1, cex = 0.9, boxsize = 0.5, colgap = unit(6, "mm"), lwd.ci = 2, ci.vertices = TRUE, ci.ve
rtices.height = 0,4)