1. ## 代码来源：生信技能树
    3. cg="ENSG00000152583" #以ENSG00000152583基因为例
    4. deg_anno[match(cg, deg_anno$ENSEMBL),]
    5. exp <- c(t(express_cpm[match(cg,rownames(express_cpm)),]))
    6. test <- data.frame(value=exp,group=group_list)
    7. library(ggplot2) 
    8. library(ggpubr) 
    9. df <- data.frame(value=exp,group=group_list)
    10. head(df)
    11. p <- ggboxplot(df, x = "group", y = "value",
    12.                color = "group", palette = "jco",
    13.                add = "jitter")
    14. p
    15. #  Add p-value
    16. p + stat_compare_means()

    代码来源：https://mp.weixin.qq.com/s/QskjnJBwtQhA9DRa82qY9Q