如何更改ggerrorplot（ggpubr）中特定组的颜色？

Answer 1

我通过使用group 和color 美学设置和geom_errorbar 来尝试使用ggplot2。可以使用geom_text 添加文本注释。

library(dplyr)
library(ggplot2)

# dummy data
SUSC_TRIG_subset = data.frame(
  Habitat_Year_Cohort = c("Seagrass (hard-bottom) 2005","Seagrass (soft-bottom) 2005", "Mussel (mixed substrate) 2005",
                          "After Hatch Year Male 2021", "After Hatch Year Female 2021", "Hatch Year Male 2021", "Hatch Year Female 2021", 
                          "Seagrass (hard-bottom) 2006","Seagrass (soft-bottom) 2006", "Mussel (mixed substrate) 2006"),
  label = c("Birch Bay", "Padilla Bay", "Penn Cove", "Adult Male", "Adult Female", "Juvenile Male", "Juvenile Female", "Birch Bay", "Padilla Bay", "Penn Cove"),
  Month = c("Dec", "Dec", "Dec", "Dec", "Dec", "Dec", "Dec", "Mar", "Mar", "Mar"),
  True.TRIG = rnorm(30, sd = 0.5) + 1
) %>%
  mutate(
    Habitat_Year_Cohort = factor(Habitat_Year_Cohort, levels = unique(Habitat_Year_Cohort))
  )

# ggplot2
g <- SUSC_TRIG_subset %>%
  group_by(
    Habitat_Year_Cohort
  ) %>%
  summarise( # for Mean-SD error bar
    ymin = mean(True.TRIG) - sd(True.TRIG),
    ymax = mean(True.TRIG) + sd(True.TRIG),
    TRIG = mean(True.TRIG),
    Month = first(Month),
    label = first(label)
  ) %>%
  mutate( # create group for color
    group = case_when(
      grepl("Seagrass \\(hard-bottom\\)", Habitat_Year_Cohort) ~ "Seagrass (hard-bottom)",
      grepl("Seagrass \\(soft-bottom\\)", Habitat_Year_Cohort) ~ "Seagrass (soft-bottom)",
      TRUE ~ "Mussel"
    )) %>%
  ggplot(aes(Month, TRIG, 
             group = Habitat_Year_Cohort, # group for dodge position
             color = group)) +
  geom_errorbar(aes(ymin = ymin, ymax = ymax), 
                position = position_dodge(0.8), 
                width = 0.2) +
  geom_point(position = position_dodge(0.8), 
             size = 4) + # size of point
  geom_text(aes(label = label), color = "black", angle = 90, 
            position = position_dodge(0.8), vjust = 1.2, size = 3, hjust = -0.2) +
  scale_color_manual(values = c("Seagrass (hard-bottom)" = "#FFFFFF", 
                                "Seagrass (soft-bottom)" = "#999999", 
                                "Mussel" = "#000000")) +
  theme(legend.title=element_blank(), # remove legend title
        legend.position="top") +
  labs(x = NULL) # remove x label

g

手动覆盖一个独立的水平误差条geom_errorbarh 用于年份注释。

df_year = data.frame(
  y = c(0.1, 0.1, 0.1),
  x = c(1 - 0.22, 1 + 0.16, 2),
  Year = c("2005", "2021", "2006")
) %>%
  cbind(
    xmin = .$x - c(0.45, 0.6, 1) * 0.3,
    xmax = .$x + c(0.45, 0.65, 1) * 0.3
  )

g +
  geom_errorbarh(data = df_year, aes(x = x, y = y, xmin = xmin, xmax = xmax), 
                 inherit.aes = F, height = 0.1) +
  geom_text(data = df_year, aes(x = x, y = y, label = Year), 
            vjust = 1.3, size = 3, inherit.aes = F)
#> Warning: Ignoring unknown aesthetics: x

^{创建于 2022-09-29，reprex v2.0.2}