探序基因肿瘤研究院 整理
构造输入文件:至少两列,第一列为group,值为1,2,3等;第二列为geneexp,值为某基因表达值。每一行代表一个样本
p <- ggboxplot(aaa, x = "group", y = "geneexp", add = "jitter",title = "the gene expression level in every group")
p <- p + stat_compare_means(method = "t.test",comparisons=list(c("1","2"),c("1","3"),c("2","3")))
print(p)
效果图: