我有以下矩阵:
row genome1 genome2 genome3
genome1 100 99.99 98.88
genome2 99.99 100 NA
genome3 98.88 NA 100
我想将其绘制为热图。
使用下面发布的代码,我在读取文件时遇到问题 - 没有保留行名和标题。另外,我希望这些值%
在热图或NA
.
ani <-read.table(file="matrix", header=T, sep="\t", dec=".", row.names=1)
anim<-as.matrix(as.data.frame(lapply(ani, as.numeric)))
library(gplots)
my_palette <- colorRampPalette(c("red", "yellow", "green"))(n = 299)
heatmap.2(anim, notecol="black", density.info="none", trace="none", col=my_palette, cellnote=anim)