base_data <- tibble::tibble(mean = c(0.578, 0.165, 0.246, 0.700, 0.348, 0.139, 1.017),
lower = c(0.372, 0.018, 0.072, 0.333, 0.083, 0.016, 0.365),
upper = c(0.898, 1.517, 0.833, 1.474, 1.455, 1.209, 2.831),
study = c("Auckland", "Block", "Doran", "Gamsu",
"Morrison", "Papageorgiou", "Tauesch"),
deaths_steroid = c("36", "1", "4", "14", "3", "1", "8"),
deaths_placebo = c("60", "5", "11", "20", "7", "7", "10"),
OR = c("0.58", "0.16", "0.25", "0.70", "0.35", "0.14", "1.02"))
base_data |>
forestplot(labeltext = c(study, deaths_steroid, deaths_placebo, OR),
clip = c(0.1, 2.5),
xlog = TRUE) |>
fp_add_header(study = c("", "Study"),
deaths_steroid = c("Deaths", "(steroid)"),
deaths_placebo = c("Deaths", "(placebo)"),
OR = c("", "OR")) |>
fp_set_style(box = "royalblue",
line = "darkblue") |>
fp_add_lines("steelblue")
base_data |>
forestplot(labeltext = c(study, deaths_steroid, deaths_placebo, OR),
clip = c(0.1, 2.5),
xlog = TRUE) |>
fp_add_header(study = c("", "Study"),
deaths_steroid = c("Deaths", "(steroid)"),
deaths_placebo = c("Deaths", "(placebo)"),
OR = c("", "OR")) |>
fp_set_style(box = "royalblue",
line = "darkblue") |>
# Add top line
fp_add_lines(h_3 = "darkred") |>
# Add surrounding box with fancy syntax
fp_add_lines(h_5 = gpar(col = "steelblue", columns = 1:4, lty = 2),
h_7 = gpar(col = "steelblue", columns = 1:4, lty = 2),
v_1 = gpar(col = "steelblue", rows = 5:6, lty = 3, lty = 2),
v_5 = gpar(col = "steelblue", rows = 5:6, lty = 3, lty = 2))
Run the code above in your browser using DataLab