本期分享如何使用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(f