Лесной участок в R с линиями эффекта и без эффекта

Я сделал лесной участок в r, используя пакет forestplot. Код выглядит следующим образом:

#forestplot 


labeltext2 <- c("A", "B", "C", "D", "E", "F", "G", "H", "I", "J", "K", "L", "M", "N", "O", "P", "Q", "R", "S", "Summary Effect")

effect2 <- c(0.12, 0.61, 0.11, 0.25, 0.24, 0.63, 0.33, 0.41, 0.38, 0.52, 0.23, 0.47,
         0.32, 0.36, 0.03, 0.15, 0.25, 0.67, 0.003, 0.32)


lower_2 <- c(0.08, 0.51, .03, 0.11, 0.06, 0.62, 0.11, 0.29, 0.18, 
         0.4, 0.19, 0.39, 0.24, 0.39, -0.13, 0.01, 0.09, 0.43,
         -0.08, 0.19)

higher_2 <- c(0.16, 0.71, .19, 0.39, 0.42, 0.64, 0.55, 0.53, 0.58, 
          0.64, 0.27, 0.55, 0.4, 0.4, 0.19, 0.29, 0.41, 0.91,
          0.08, 0.44)

forestplot(labeltext2, effect2, lower_2, higher_2, zero = .32,
       cex = 2, 
       lineheight = "auto",
       xlab = "effect size",
       xticks = c(-.5, 0, .5, 1, 1.5),
       title = "ForestPlot",
       new_page = TRUE)

Что позволяет мне получить это изображение: Forestplot

У этого изображения есть строка эффекта, идущая к 0,32, суммарный эффект, использующий нулевой аргумент. Я хотел бы добавить дополнительную толстую черную линию на 0, чтобы показать линию «без эффекта». Кто-нибудь знает, как это сделать? Я открыт для использования другого пакета/функции.

Благодарю вас!


r statistics forestplot
person smd2179    schedule 01.08.2016    source источник
comment
Вы можете попробовать abline(v=0). Если пакет использует базовую графику, это должно работать.   -  person lmo    schedule 01.08.2016

Ответы (1)


arrow_upward
0
arrow_downward

Я не могу понять, как лесной участок использует сетку, но вот ручное решение:

forestplot(labeltext2, effect2, lower_2, higher_2, zero = .32,
           cex = 2, 
           lineheight = "auto",
           xlab = "effect size",
           xticks = c(-.5, 0, .5, 1, 1.5),
           title = "ForestPlot",
           new_page = TRUE)

# Add line manually
x = .397
y0 = .12
y1 = .875
grid.lines(c(x, x), c(y0, y1), default.units = "npc",
           gp = gpar(lwd = 2))
person Vandenman    schedule 01.08.2016