箱线图加看两者之间的显著性

library(ggplot2)
library(ggpmisc)
library(ggpubr)
library(readxl)
setwd(“F:/论文/博士文章1/数据”)
plot_data <- read_excel(“10-数据.xlsx”,sheet = “Sheet10”)

p<- ggplot(data=plot_data)+geom_boxplot(mapping=aes(x=type,y=SR,colour = type ),alpha = 0.2,size=1.5,width = 0.6)+
scale_color_manual(limits=c(“BB”,“BQ”,“CH”,“DC”,“HTT”,“JDX”,“JSYC”,“LM”,“QY”,“SH”,“TM”,“WLFY”,“WLJM”,“WLJS”),
values=c(“#85B22E”,“#5F80B4”,“#E29827”,“#922927”,“#9FA0EB”,“#5F87B4”,“#E29895”,“#34DB97”,“#85B36E”,“#6DEAE9”,“#E25267”,“#928567”,“#E6BFE2”,“#C0E6A2”))+

theme_classic( # 主题设置，这个是无线条主题
base_line_size = 1 # 坐标轴的粗细
)+
labs(x=“type”,y=“SR”) +# 添加标题，x轴，y轴内容
theme(plot.title = element_text(size = 15,
colour = “black”,
hjust = 0.5),
axis.title.y = element_text(size = 15, #Y轴坐标名字
# family = “myFont”,
color = “black”,
face = “bold”,
vjust = 1.9,
hjust = 0.5,
angle = 90),
axis.title.x = element_text(size = 15, #Y轴坐标名字
# family = “myFont”,
color = “black”,
face = “bold”,
vjust = 1.9,
hjust = 0.5),
legend.title = element_text(color=“black”, # 修改图例的标题
size=15,
face=“bold”),
legend.text = element_text(color=“black”, # 设置图例标签文字
size = 15,
face = “bold”),
axis.text.x = element_text(size = 15, # 修改X轴上字体大小，
color = “black”, # 颜色
face = “bold”, # face取值：plain普通，bold加粗，italic斜体，bold.italic斜体加粗
vjust = 0.5, # 位置
hjust = 0.5,
angle = 0), #角度
axis.text.y = element_text(size = 13,
color = “black”, # 颜色
face = “bold”,
vjust = 0.5,
hjust = 0.5,
angle = 0) )
p

ggsave(“SR关系.tiff”,path = “F:/论文/博士文章1/0tu数据/模型”,width = 15,height =13.5)
tiff(“SR关系.tiff”,width = 200,height = 180,units = “px”)
dev.off()

    #查看两者之间的关系                             
    library("PMCMRplus") 
    plot_data$type <- as.factor(plot_data$type)
    b <- bwsAllPairsTest(SR~type,data = plot_data)
    
    summary(b)