生物信息学与算法科研信息学生物信息学从零开始学

R语言获取四分位分组\高低表达分组

2019-04-24  本文已影响3人  落寞的橙子

如题

data<-read.csv()#读入矩阵表达数据,row为gene names ,col为表达值
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(data, 2, Quantile)
new_median<-apply(data, 2, Median)

自由发挥自由扩散

上一篇 下一篇

猜你喜欢

热点阅读