boxplot_pair_显著性

setwd('D:\\11_NM_DNAmeth\\06_all_meth/')
library(ggplot2)
library(ggsignif)
library(ggpubr)
data<-read.table('./4_twobed_sin.txt',header = 0)
head(data)
colnames(data)<-c('sample2','level','type','value')#'sample',
table(data$sample2)
table(data$type)
data$sample2<- factor(data$sample2,
                     levels= c('BCD',"CK"),ordered = TRUE)

# compare_means(log2TPM~sample, data=data, group.by = "type",method = "wilcox.test" )
head(data)
> head(data)
  sample2 level type     value
1     BCD   TEs    C 0.6131419
2     BCD   TEs    C 0.3783133
3     BCD   TEs    C 0.1697793
4     BCD   TEs    C