getwd()
library(readr)
datct<-read_delim("boxplot_clock_tissue.txt",delim = "\t")
datct
df = melt(datct)
#预览数据
head(df)
install.packages("ggplot")
install.packages("RColorBrewer")
install.packages("path/to/ggplot2.tar.gz", repos = NULL)
install.packages("path/to/RColorSrewer.tar.gz", repos = NULL)
library(ggplot2)
library(RColorBrewer)
display.brewer.all()
mycol<-brewer.pal(n=3,name = "Set3")
mycol
mycol<-c("#8DD3C7","#FB8072","#80B1D3")#要查看代号运行mycol
library(ggpubr)
p<-ggboxplot(df,x="variable",y="value",color="variable",palette = mycol,legend="none",
add ="jitter",add.params = list(size=3,jitter=0.1,alpha=0.5),outlier.shape = NA)+
stat_compare_means(method = "kruskal.test",label.y = 4,label.x = 0.8,label = "p.format")+
stat_compare_means(comparisons = list(c("G","F"),c("G","H"),c("F","H")),method="wilcox",label="p.signif")+
theme(panel.grid = element_line(color = NA),panel.grid.minor = element_line(color = NA),
panel.border = element_rect(fill = NA,color = NA),
axis.text.x= element_text(size = 10,face ="italic",hjust = 0.5),
axis.text.y = element_text(size = 10,face ="plain",hjust = 0.5),
axis.title.y =element_text(vjust=0.2,size=12,face="bold"))
p
library(gridExtra)
library(cowplot)
ggsave("123.pdf",width=9,height=6)