Я провожу метаанализ в метафоре.После преобразования лесного графика обратно из логарифмического масштаба в нормальный масштаб, я надеюсь добавить исходные доверительные интервалы (указанные в ci.lb и ci.ub в данных ниже).Это связано с тем, что доверительные интервалы, которые были преобразованы метафорой обратно из логарифмической шкалы, отличаются от оригиналов, и я хочу, чтобы на моем лесном участке отображались исходные доверительные интервалы (как показано в моем анализе пропорциональных рисков Кокса)
Код, который я использую:
library(metafor)
library(tidyverse)
#create data [this is just an example- I have real data]
dat <- tibble(study = c(1, 2, 3, 4, 5, 6),
yi = runif(6, 0.001, 0.5),
sei = runif(6, 0.2, 0.4),
ni = runif(6, 70, 400),
outcome = "DFS",
author = "author",
year = c("2001","2002", "2003","2004", "2005", "2006"),
ci.lb = runif(6, 0.6, 1.5),
ci.ub = runif(6, 1.3, 4))
# random-effects model analysis for outcome DFS
res <- rma(yi, sei, ni, data=dat)
res
predict(res, transf=exp, digits=2)
#forest plot
par(mar=c(4,4,1,2))
res <- rma(yi, sei=sei, data=dat, slab=paste(author, year, sep=", "))
tmp <- cumul(res, order=order(dat$year))
forest(tmp, xlim=c(-4,2), at=log(c(0.125, 0.25, 0.5, 1, 2)),
atransf=exp, digits=c(2,3), cex=0.75, annotate=FALSE)
par(cex=0.75, font=2)
text(-4, 15, "Author(s) and Year", pos=4)
text( 2, 15, "Risk Ratio [95% CI]", pos=2)
forest(res, atransf= exp, showweights = T)
Кто-нибудь знает, как я могу это сделать?Боролись с этим в течение некоторого времени.Буду очень признателен за любые ваши идеи :) спасибо
участок, выведенный этим кодом