囤优质画图链接,有需要自取哈!!!也会随时补充 哈哈哈
链接1:https://mp.weixin.qq.com/s?src=11×tamp=1710320556&ver=5135&signature=RBCFiyhK6o1hnTE1GB8N4GAAInGM7qSuJKkFhi7P7AsU4c4JhV6Mq22SsinSPlIfbnfIQZGyFcmASFD4E4fpkypcISNrc-giz3WKbYW034A41QjyJePWe3Gg8ldhIq&new=1
library(ggplot2)#绘图
library(ggsignif)
library(ggpubr)#计算差异
rm(list=ls())#清除环境
library(readxl)
data1 <- read_excel("CRISPR1A.xlsx")
table(data1$Genotype) #查看品种
# 像这种多组之间两两比较,需要提前设置一个list,里面包含的是两两比较的对象
list = list(c("setosa","versicolor"),c("virginica","versicolor"),c("setosa","virginica"))
ggplot(data1, # 用来画图的数据集名称
aes(Genotype, # 数据集的Genotype列作为横坐标
Leaf1, # 数据集的value列作为纵坐标
color = Genotype, # 根据数据集的Genotype列设置箱线图的边框颜色
fill = Genotype))+ # 根据数据集的Genotype列设置散点图和箱线图的填充颜色
geom_jitter(width=.15)+ # 将散点图的宽度设置为0.15,缩窄一些
geom_boxplot(width=.4, alpha=.2)+ # 将箱线图的宽度设置为0.4,缩窄了一些,不透明度设置为0.2,填充颜色变得很浅
scale_color_manual(values = c("#BF5960","#6F99AD","#F9A363"))+ # 自定义边框颜色
scale_fill_manual(values = c("#BF5960","#6F99AD","#F9A363"))+ # 自定义填充颜色,可以和边框颜色不一样,我觉得配套的比较好看,就用了一样的颜色
theme_test()+ # 把背景颜色去掉并添加边框,可以修改为其他背景,比如theme_classic
labs(y="Leaf1", x="")+ # 修改横纵坐标名称
stat_compare_means(comparisons = list, method = "wilcox.test", #像这种非配对的非正态分布的数据,可以用非参数检验方法,这里选的是非参的其中一种,叫秩和检验
size=8, # 将p值字体改小一些
label = "p.signif")+ # 把p值改为星号
scale_y_continuous(expand = c(0.1,0.1))+ # 将y轴上下拉宽一点,这样被遮住的p值就可以呈现出来啦~
theme(axis.text = element_text(color="black", size=10),
axis.line = element_blank()) # 把横纵坐标的字体颜色改为黑色,看起来更沉稳
上述代码搬运自链接1