Início / coding / Perguntas / 78357543
Accepted
zhang
Asked: 2024-04-20 16:35:43 +0800 CST

Adicionando marcadores de significância para comparações de grupos em um boxplot usando ggplot2

  • 772

Eu tenho dados como o seguinte:

set.seed(123)
expr_data <- data.frame(
    cell_line_name = rep(c("CL1", "CL2"), each = 50),
    GSE_id = rep(c("GSE1", "GSE2"), each = 50),
    model = rep(c("M1", "M2"), each = 50),
    expr = rnorm(100),
    age_group = sample(c('A', 'B'), 100, replace = TRUE)
)

    cell_line_name GSE_id model        expr age_group
1              CL1   GSE1    M1 -0.34391723         B
2              CL1   GSE1    M1  0.09049665         A
3              CL1   GSE1    M1  1.59850877         A
4              CL1   GSE1    M1 -0.08856511         A
5              CL1   GSE1    M1  1.08079950         B
6              CL1   GSE1    M1  0.63075412         B
7              CL1   GSE1    M1 -0.11363990         A
8              CL1   GSE1    M1 -1.53290200         B
9              CL1   GSE1    M1 -0.52111732         A
10             CL1   GSE1    M1 -0.48987045         B
11             CL1   GSE1    M1  0.04715443         A

Eu o que polt um boxplot com cell_line_name, GSE_id, e modelar como eixo x, expr como eixo y, e mostrar a comparação entre o grupo A e B

ggplot(expr_data.gene, aes(x = paste0(cell_line_name, '_', GSE_id, '_', model), y = expr, color = age_group)) + 
    geom_boxplot() + 
    geom_point(position = position_jitterdodge(), alpha = 0.5)

insira a descrição da imagem aqui

Também o que adicionar rótulo significativo para cada comparação no valor x, tentei ggsignife ggpubr, mas não funcionou.

ggplot(expr_data.gene, aes(x = paste0(cell_line_name, '_', GSE_id, '_', model), y = expr, color = age_group)) + 
    geom_boxplot() + 
    geom_point(position = position_jitterdodge(), alpha = 0.5) + 
    ggpubr::stat_compare_means(comparisons = list(c("A", "B")))

resultado esperado: insira a descrição da imagem aqui

1 respostas

  1. Best Answer

    A comparação precisa ser os nomes ou o índice de 2 valores no eixo x, então você pode:

    • Use facet_wrapem vez de se esquivar e ocultar os painéis, ou
    • Incluir age_groupno eixo x:
    library(tidyverse)
library(ggsignif)

set.seed(123)

expr_data <- data.frame(
  cell_line_name = rep(c("CL1", "CL2"), each = 50),
  GSE_id = rep(c("GSE1", "GSE2"), each = 50),
  model = rep(c("M1", "M2"), each = 50),
  expr = rnorm(100),
  age_group = sample(c('A', 'B'), 100, replace = TRUE)
)

# Using facet_wrap with hidden panels
expr_data |> 
  mutate(x = paste0(cell_line_name, '_', GSE_id, '_', model)) |> 
  ggplot(aes(x = age_group, y = expr, color = age_group)) + 
  geom_boxplot() + 
  geom_point(position = position_jitter(), alpha = 0.5) + 
  geom_signif(
    comparisons = list(c("A", "B")),
    map_signif_level = TRUE, textsize = 6
  ) +
  facet_wrap(~x, scales = "free_x", strip.position = "bottom") +
  ylim(c(NA, 3)) +
  labs(x = "Axis Title") +
  theme_bw() +
  theme(
    axis.text.x = element_blank(),
    axis.ticks.x = element_blank(),
    panel.border = element_blank(),
    strip.background = element_blank()
  )


    
# Adding age-group into the x-axis
expr_data |> 
  mutate(x = paste0(cell_line_name, '_', GSE_id, '_', model, "_", age_group)) |> 
  ggplot(aes(x = x, y = expr, color = age_group)) + 
  geom_boxplot() + 
  geom_point(position = position_jitter(), alpha = 0.5) + 
  geom_signif(
    comparisons = list(c("CL1_GSE1_M1_A", "CL1_GSE1_M1_B"), c("CL2_GSE2_M2_A", "CL2_GSE2_M2_B")),
    map_signif_level = TRUE, textsize = 6
  ) +
  ylim(c(NA, 3)) +
  labs(x = "Axis Title") +
  theme_bw() +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

    Criado em 20/04/2024 com reprex v2.1.0

    • 1

relate perguntas