数据库代码公开:四分位值分组获得Q1 Q2 Q3 Q4

作者:半步博导 时间:2024年6月30日 04:37 阅读量:158


#以PDCD1为例
geneName = "PDCD1"
geneNamedata = as.data.frame(data[geneName,])
colnames(geneNamedata) = geneName
#设置函数,计算四分位分组
Quantile <- function(x){
  ifelse(x>quantile(x,.75),"Q1",
         ifelse(x>quantile(x,.5),"Q2",
                ifelse(x>quantile(x,.25),"Q3",
                       "Q4")
                )
         )
}
#如采用中位值分组
#Median <- function(x){
#  ifelse(x > median(x),"High","Low")
#}
new_quantile = apply(geneNamedata, 2, Quantile)
colnames(new_quantile) = "ExpressionType"
ExpressionType = cbind(geneNamedata,new_quantile)
#获得排序后索引
idx = order(ExpressionType$ExpressionType)
#按索引重新提取以达到排序效果
ExpressionType = ExpressionType[idx,]