Как я могу добавить исходные доверительные интервалы отношения рисков обратно в участок метафора леса после преобразования из преобразования журнала? - PullRequest
0 голосов
/ 02 декабря 2018

Я провожу метаанализ в метафоре.После преобразования лесного графика обратно из логарифмического масштаба в нормальный масштаб, я надеюсь добавить исходные доверительные интервалы (указанные в 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)

Кто-нибудь знает, как я могу это сделать?Боролись с этим в течение некоторого времени.Буду очень признателен за любые ваши идеи :) спасибо

участок, выведенный этим кодом

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...