R语言 boxplot作图图内展示校正后的P值(padj)

HenryPeng21

已于 2022-07-29 23:43:05 修改

阅读量4.1k

点赞数 2

分类专栏： R boxplot padj 文章标签： r语言

于 2022-07-25 20:58:40 首次发布

本文链接：https://blog.csdn.net/HenryPeng21/article/details/125983036

版权

Boxplot绘制步骤

1. 导入数据--> 2. 数据预处理 --> 3. 统计检验+pvalue校正 --> 4. pvalue加到boxplot中

导入数据

library(openxlsx)
setwd("/Users/mac/Desktop/Imperial\ College\ of\ London/1--data")

2. 数据预处理

##去除data中Group.name.3中的空值
data = data[-which(is.na(data$Group.name.3)),]

##将行名中的"-"," "变成”_“如果是"-"," "会影响后面做图 
colnames(data) <- gsub("-","_",colnames(data))
colnames(data) <- gsub(" ","_",colnames(data)

3. 统计检验+pvalue校正

##定义p_value校正的function 选择方法 默认为wilcox.test统计检验 与 FDR校正
get_adj_p <- function(data, .col, .grp = "Sample", comparisons = NULL,
                      method = "wilcox.test", p.adjust.method = "fdr", p.digits = 3L, ...)
 {
  # Compute p-values
  comparison.formula <- paste0(.col, "~", .grp) %>%
    as.form