我有一个基因组数据文件,我正在尝试创建数据热图以及染色体信息的侧栏。为了制作热图,我从我的文件转换了数值,以创建一个我能够绘制的数据矩阵:
Heatmap <- heatmap.2(DataMatrix, trace="none", col=greenred, key=T, keysize=1.5, labRow=NA, density.info="none", Colv=FALSE)
但是,我想添加一个额外的颜色侧栏来指示所述数据的染色体位置。特别是我想只有两种颜色,如果它来自任何其他染色体的“chrX”或“黄色”,则说“黑色”。包含此信息的列不在我的数据矩阵中,我不确定如何将此信息包含在我的热图中。任何帮助将不胜感激!
答案 0 :(得分:4)
heatmap.2
函数有一个RowSideColors
参数用于此目的,正如您可能从文档中推断的那样:
RowSideColors :(可选)长度为'nrow(x)'的字符向量 包含可能的垂直边栏的颜色名称 用于注释'x'的行。
我最近还在NMF包中发现了aheatmap
函数,该函数在热图上制作了具有任意数量的注释行/列的精美热图。在该主页上,您会找到several links到heatmap examples。