R:热图的read.table

Question

我有以下矩阵：

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)

Answer 1

试试这个：

ani <- read.table(text='
row genome1 genome2 genome3
genome1 100 99.99 98.88
genome2 99.99 100 NA
genome3 98.88 NA 100
', header=T, row.names=1, quote="")

anim <- as.matrix(ani)
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, margins=c(10,10))

​

​