如何将p值添加到ggplot中的分组图表中？

Question

当我想在我的图中添加 p 值时：

library(tidyverse)
library(ggpubr)
library(rstatix)

 stat.test3 <- MP %>%
    group_by(TBI) %>%
    wilcox_test(age ~ mp_1) %>%
    adjust_pvalue(method = "bonferroni") %>%
    add_significance("p.adj")%>%
  mutate(y.position = 35)

 C2<-  ggplot(data=MP, aes(x=TBI, y=age, fill=mp_1))+
    geom_violin()+
    geom_boxplot(width=.2, fatten=NULL, position = position_dodge(0.9))+
    stat_summary(fun="median", geom="point", position = position_dodge(0.9))+
    stat_summary(fun.data = "mean_se", geom = "errorbar", width=.1, position = position_dodge(0.9))+
      scale_fill_brewer(name="Mind-pop", palette = "Accent")

  C2+ stat_pvalue_manual(stat.test3, xmin = "TBI",xmax = NULL)

它给了我这个错误：

FUN(X[[i]], ...) 中的错误：找不到对象“mp_1” 将 stat_pvalue 添加到对象后会显示此错误。 我该如何解决？

Answer 1

我不熟悉ggpubr，所以不能说我了解根本问题，但似乎color=mp_1 而不是fill=mp_1 可能会解决您的问题。这是在以下行中：

C2 <- ggplot(data=MP, aes(x=TBI, y=age, color=mp_1)).

完整的代码如下。我还更改了y.position，以便重要性位于情节的顶部。

MP <- structure(list(age = c(55L, 54L, 56L, 60L, 55L, 53L, 61L, 56L, 
                             58L, 58L, 56L, 58L, 58L, 58L, 59L, 57L, 56L, 60L, 57L, 58L, 61L, 
                             60L),
                     mp_1 = structure(c(2L, 2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 
                            2L, 2L, 2L, 2L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("No", "Yes"), class = "factor"), 
                     TBI = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("HC", 
                "TBI"), class = "factor")), class = "data.frame", row.names = c(NA, -22L))

library(tidyverse)
library(ggpubr)
library(rstatix)
    
stat.test3 <- MP %>%
  group_by(TBI) %>%
  wilcox_test(age ~ mp_1) %>%
  adjust_pvalue(method = "bonferroni") %>%
  add_significance("p.adj") %>%
  mutate(y.position = 61.5)

C2 <-  ggplot(data=MP, aes(x=TBI, y=age, color=mp_1))+
  geom_violin() +
  geom_boxplot(width=.2, fatten=NULL, position = position_dodge(0.9))+
 stat_summary(fun="median", geom="point", position = position_dodge(0.9))+
 stat_summary(fun.data = "mean_se", geom = "errorbar", width=.1, position = position_dodge(0.9))+
 scale_fill_brewer(name="Mind-pop", palette = "Accent")

C2 + stat_pvalue_manual(stat.test3, xmin = "TBI",xmax = NULL)