Neatly annotate forest plot using ggplot2

Question

How can I neatly annotate the four points in this forest plot using the "estimate" column? I would like them vertically aligned if possible.

library(ggplot2)
library(wesanderson)

## make data

df1 <- data.frame(matrix(, nrow=4, ncol=5))
colnames(df1) <- c("Metabolite", "OR", "CI_low", "CI_high", "Method")
df1$Metabolite <- c("Isoleucine", "Leucine", "Isoleucine", "Leucine")
df1$OR <- c(0.75, 0.64, 0.88, 0.95)
df1$CI_low <- c(0.61, 0.48, 0.8, 0.85)
df1$CI_high <- c(0.89, 0.8, 0.96, 1.05)
df1$Method <- c("Metabolomic", "Metabolomic", "GRS", "GRS")
df1$estimate <- c("0.75 (0.61, 0.89)", "0.64 (0.48, 0.8)",  "0.88 (0.8, 0.96)",  "0.95 (0.85, 1.05)")

# Plot using ggplot ------------------------------

ggplot(df1, aes(x=Metabolite, y=OR, ymin=CI_low, ymax=CI_high,col=Method,fill=Method)) + 
  #specify position here
  geom_linerange(size=1,position=position_dodge(width = 0.5), show.legend = FALSE) +
  geom_hline(yintercept=1, lty=2) +
  #specify position here too
  geom_point(size=3, shape=21, colour="white", stroke = 0.5,position=position_dodge(width = 0.5)) +
  # scale_fill_manual(values=barCOLS)+
  # scale_color_manual(values=dotCOLS)+
  scale_fill_manual(values = wes_palette(name = "GrandBudapest2")) +
  scale_color_manual(values = wes_palette(name = "GrandBudapest2")) +
  scale_x_discrete(name="Metabolite") +
  scale_y_continuous(name="Odds ratio (95% CI)", 
                     limits = c(0.4, 1.5), 
                     breaks = c(0.4, 1, 1.5),
                     labels=c("0.4", "1", "1.5")) +
  coord_flip() +
  guides(fill = guide_legend(reverse = TRUE)) +
  theme_classic()

Who's W. E. Sanderson? ;-)

qdread
– qdread

2021-01-08 14:21:49 +00:00
Commented Jan 8, 2021 at 14:21 — qdread
– qdread, Commented Jan 8, 2021 at 14:21
google.com/…

krtbris
– krtbris

2021-01-08 14:23:02 +00:00
Commented Jan 8, 2021 at 14:23 — krtbris
– krtbris, Commented Jan 8, 2021 at 14:23

teunbrand · Accepted Answer · 2021-01-08 14:10:58Z

How about a typical geom_text() with a dodge position?

library(ggplot2)
library(wesanderson)

## make data

df1 <- data.frame(matrix(nrow=4, ncol=5))
colnames(df1) <- c("Metabolite", "OR", "CI_low", "CI_high", "Method")
df1$Metabolite <- c("Isoleucine", "Leucine", "Isoleucine", "Leucine")
df1$OR <- c(0.75, 0.64, 0.88, 0.95)
df1$CI_low <- c(0.61, 0.48, 0.8, 0.85)
df1$CI_high <- c(0.89, 0.8, 0.96, 1.05)
df1$Method <- c("Metabolomic", "Metabolomic", "GRS", "GRS")
df1$estimate <- c("0.75 (0.61, 0.89)", "0.64 (0.48, 0.8)",  "0.88 (0.8, 0.96)",  "0.95 (0.85, 1.05)")

# Plot using ggplot ------------------------------

ggplot(df1, aes(x=Metabolite, y=OR, ymin=CI_low, ymax=CI_high,col=Method,fill=Method)) + 
  geom_linerange(size=1,position=position_dodge(width = 0.5), show.legend = FALSE) +
  geom_hline(yintercept=1, lty=2) +
  geom_point(size=3, shape=21, colour="white", stroke = 0.5,position=position_dodge(width = 0.5)) +
  geom_text(aes(y = 0.8,label = estimate), 
            position = position_dodge(width = 1),
            show.legend = FALSE, hjust = 0.5) +
  scale_fill_manual(values = wes_palette(name = "GrandBudapest2")) +
  scale_color_manual(values = wes_palette(name = "GrandBudapest2")) +
  scale_x_discrete(name="Metabolite") +
  scale_y_continuous(name="Odds ratio (95% CI)", 
                     limits = c(0.4, 1.5), 
                     breaks = c(0.4, 1, 1.5),
                     labels=c("0.4", "1", "1.5")) +
  coord_flip() +
  guides(fill = guide_legend(reverse = TRUE)) +
  theme_classic()

^{Created on 2021-01-08 by the reprex package (v0.3.0)}

Collectives™ on Stack Overflow

Neatly annotate forest plot using ggplot2

1 Answer 1

Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

Comments

Your Answer

Sign up or log in

Post as a guest

Related