如何在 ggplot2 中使用 stat_pvalue_manual() 将 t 和 p 值设为斜体

Question

如何使用 stat_pvalue_manual() 将 t 和

p

值从下面的箱线图中的

ggpubr()

斜体化？

• 剧情是：

• 代码是：

##
library(tidyverse) 
library(rstatix)   # Identify outliers and Stats Summary function
library(ggpubr)    # Plot correlation
library(RVAideMemoire) # Shapiro test on groups 
library(glue) # Round p value in the plots 
##
df %>%
      ggplot(., aes(x = TEST, y = VALUE)) +
      stat_boxplot(geom = "errorbar",
                   width = 0.15) +
      geom_boxplot(aes(fill = TEST), outlier.colour = "yellow", outlier.shape = 18,
                   outlier.size= 2, notch = F) +
      stat_pvalue_manual(df %>%
                           t_test(VALUE ~ TEST, paired = T) %>%
                           add_xy_position(),
                         label = "t ({df}) = {round(statistic, 2)}, p = {round(p, 3)}") 
         
               
               
               #stat_pvalue_manual(df %>%
                #                    wilcox_test(VALUE ~ TEST, paired = T) %>%
                 #                   add_xy_position(),
                  #                label = "W = {round(statistic, 2)}, p = {round(p, 3)}") +

问题

• 1：我见过

  annotate()

  here，但我找不到

  stat_pvalue_manual()

  的解决方案。 如何将 t 和 p 斜体化？ 非常感谢您的帮助。谢谢！
  

• 2：“奖金”：有谁知道如何制作

  stat_pvalue_manual()

  显示约定p< 0.05* , p < 0.01 ** and p<0.001***?

• 数据：

structure(list(ID = c(1L, 1L, 2L, 2L, 3L, 3L, 4L, 4L, 5L, 5L, 
6L, 6L, 7L, 7L, 8L, 8L, 9L, 9L, 10L, 10L, 11L, 11L, 12L, 12L, 
13L, 13L, 14L, 14L, 15L, 15L, 16L, 16L, 17L, 17L, 18L, 18L, 19L, 
19L, 20L, 20L, 21L, 21L, 22L, 22L, 23L, 23L, 24L, 24L, 25L, 25L, 
26L, 26L, 27L, 27L, 28L, 28L, 29L, 29L, 30L, 30L, 31L, 31L), 
    TEST = c("TEST1", "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", 
    "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", 
    "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", 
    "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", 
    "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", 
    "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", 
    "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", 
    "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", 
    "TEST2", "TEST1", "TEST2", "TEST1", "TEST2", "TEST1", "TEST2"
    ), VALUE = c(70, 73, 66, 140, 59, 59, 61, 70, 73, 107, 90, 
    91, 59, 95, 82, 133, 60, 80, 40, 45, 77, 96, 49, 50, 61, 
    56, 92, 120, 56, 110, 56, 110, 68, 83, 87, 113, 87, 114, 
    87, 134, 79, 57, 45, 65, 83, 71, 50, 52, 44, 56, 95, 138, 
    82, 126, 134, 118, 42, 72, 70, 87, 84, 112)), class = c("tbl_df", 
"tbl", "data.frame"), row.names = c(NA, -62L))