col分类表头注释文件:head_group
表达量文件:tpm
基因组信息文件: genome
library(ComplexHeatmap)
iw=read.table('head',sep='\t',header = T) ###col_ha,group信息
cs=iw[,1]
head(cs)
df<-read.table('gene_tpm1',sep='\t',header=F)
df[1:5,1:5]
#添加rownames的基因名
#iw=read.table('ttest_CV_LD.sort_non_top300.gene',sep = '\t')
rownames(df)=df[,1]
#iw=as.matrix(iw[,-c(1:2)])
#pp=read.table(