所用R包:pheatmap
#pheatmap画热图,TP53基因三个探针的表达量
#重点在于整理数据:
#1,分组信息对应的样本编号annotation-col
#2,所需基因表达量对应的样本编号qq
library(pheatmap)
tp53=as.data.frame(probe_tp53)
colnames(tp53) <- "rowname"
qq=merge(tp53,exp_tab,by="rowname")
qq$rowname
qq=qq[,-1]
rownames(qq) <- c("1939_at","1974_s_at","31618_at" )
library(stringr)
group_list=ifelse(str_detect(aa$Disease,"progres"),"progres","stable")
annotation_col=data.frame(Disease=group_list)#分组信息
rownames(annotation_col)=colnames(qq)#让annotation_col的行名等于qq的列名
pheatmap(qq,
show_colnames =T,#F不显示列名,T为显示
show_rownames = T,#F不显示行名
annotation_col=annotation_col,
scale = "row")
dev.off()
纪念第一次自己整理数据做出来的热图...