作者:生信札记
来源:简书
准备数据,最简单粗暴的方式,从Excel中复制,
此时,这个基因表达矩阵已经存放在电脑的剪切板中,回到R语言运行环境,输入(如果是直接复制下面这行代码,那么注意,复制代码的同时已经替换了剪贴板内容,此时黏贴代码之后不要回车,先回去Excel重新复制表达开矩阵)Exp
注意,所有字符都是外文的,如果出错,注意是不是 横杆 在复制代码的时候被转码成 中文的,自己输入就没问题。
此时你只需要第三行代码,就可以输出一张可以看的热图pheatmap(Exp)
一般情况下,我们需要相对基因进行标准化pheatmap(Exp,scale="row")
列的顺序是时间序列,我们不应该对列进行聚类pheatmap(Exp,scale="row",cluster_cols=F)
配色似乎也不好看,pheatmap(Exp,scale="row",cluster_cols=F,color = colorRampPalette(c("navy", "white", "firebrick3