添加统计检验 P value

本期分享如何使用R为为分组数据差异添加统计检验P value。

使用R ggplot2自带的鸢尾花iris数据集，比较setosa、versicolor、virginica三类鸢尾花的Sepal.Width是否具有显著差异。

下文代码使用时，确保R中已安装了"ggplot2"，"ggsignif"，"ggpubr"，"ggbeeswarm"，"ggtext" 和"showtext"包。

---

数值 P value

详细注释代码，

library(ggplot2)      # 调用数据iris，使用geom_boxplot()绘制箱图
library(ggsignif)     # 使用geom_signif()添加p值
library(ggpubr)       # 使用theme_pubr()设置绘图主题
library(ggbeeswarm)   # 使用geom_quasirandom()添加蜂窝图
library(ggtext)       # 使用element_markdown()设置文字显示
library(showtext)     # 使用showtext_auto()解决R中中文字体显示「豆腐块」问题
showtext_auto()

ggplot(iris, aes(Species, Sepal.Width, fill = Species)) +
  geom_boxplot(fill = "transparent") +  # 绘制箱图
  geom_quasirandom(shape = 21, color = "black", 
                   alpha = 1, size = 5) +  # 绘制抖动图
  geom_signif(  # 添加p值