数据库代码公开:四分位值分组获得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,]