多组间比较散点图+误差棒（自备）

hx2024

已于 2024-05-07 21:39:49 修改

阅读量407

点赞数 3

文章标签： windows 服务器 linux 数据可视化

于 2024-05-07 00:48:34 首次发布

本文链接：https://blog.csdn.net/hx2024/article/details/138513760

版权

数据

计算四分位值

作图

数据

rm(list = ls())
library(tidyverse)
library(reshape2)
library(ggplot2)
library(dplyr)
library(ggpubr)
library(cowplot)
#library(Seurat)
data <- iris##鸢尾花数据集
dat <- data[,c(5,1)]#单个数据进行分析

计算四分位值

#根据分组计算四分位及中位数
dat1 <- dat %>% group_by(Species) %>% mutate(upper =  quantile(Sepal.Length, 0.75), 
                                                lower = quantile(Sepal.Length, 0.25),
                                                mean = mean(Sepal.Length),
                                                median = median(Sepal.Length))

> head(dat1)
# A tibble: 6 × 6
# Groups:   Species [1]
  Species Sepal.Length upper lower  mean median
  <fct>          <dbl> <dbl> <dbl> <dbl>  <dbl>
1 setosa           5.1   5.2   4.8  5.01      5
2 setosa           4.9   5.2   4.8  5.01      5
3 setosa           4.7   5.2   4.8  5.01      5

作图

#比较分组
my_comparisons =list( c("setosa","versicolor"),c("versicolor","virginica"),
                      c("setosa","virginica"))
P <- ggplot(dat1,aes(x=Species,y=Sepal.Length)) + #ggplot作图
  geom_jitter(shape = 21,aes(fill=Species),width = 0.25) + 
  stat_summary(fun=mean, geom="point", color="grey60") +
  theme_cowplot() +
  theme(axis.text = element_text(size = 10),axis.title = element_text(size = 10),legend.text = element_text(size = 10),
        legend.title = element_text(size = 10),plot.title = element_text(size = 10,face = 'plain'),legend.position = 'none') + 
  labs(title = "Species",y='Sepal.Lengthage') +
  geom_errorbar(aes(ymin = lower, ymax = upper),col = "grey60",width =  1)+#误差棒
  #差异检验
  stat_compare_means(comparisons=my_comparisons,
                   label.y = c(7.5, 8, 8.5),
                   method="t.test",#wilcox.test
                   label="p.signif")
P
dev.off()