1

我有这个数据

line_tab<-                

    elements   scaffold  start    end Lengths        A           B        C        D                      E                     F                 G                  
  BovB_1297  KQ417326.1    763   2186    1423 0.1565491 0.02291220 1.9857781  1.4346243          2.81754891            1.23339700         1.0394839
  BovB_3317  KQ416600.1 660334 661626    1292 0.0000000 0.02523534 0.0000000  0.1934799          0.11710299            0.06318396         0.0000000
  BovB_5130  KQ418760.1 327212 328873    1661 0.0000000 0.00000000 0.0265819  0.0000000          0.02277198            0.03686048         0.0000000
  BovB_6215  KQ417629.1 339534 340586    1052 0.0000000 0.00000000 0.8394018  1.2277020          0.43145549            0.34919346         0.2685752
  BovB_9578  KQ420742.1 157232 158657    1425 0.0000000 0.00000000 0.0000000  0.0000000          0.02654334            1.00251877         0.0000000
  BovB_9976  KQ423592.1  30989  32425    1436 0.0000000 0.02270478 0.4612034  0.3191430          0.13170009            0.24160383         0.1388863

然后我pheatmap用这个代码创建了一个热图

  pheatmap(log10(line_tab+1), color = col.pal.red, 
           fontsize = 8, cluster_cols = FALSE, cluster_rows= FALSE, cellwidth = 15, cellheight = 5, 
           fontsize_row = 5, main = "LINEs length")

但我想在我的热图中添加更多信息。我正在尝试在每一行附近添加一个图例，该图例显示每个“元素”的“长度”列。你能建议我怎么做吗？

4

1 回答 1

0

如果我理解正确，您想将标签添加到热图的每一行。然后您可以使用以下参数labels_row：pheatmap()

library(pheatmap)
library(RColorBrewer)

pheatmap(log10(line_tab[, -(1:2)] + 1), 
    color = colorRampPalette(rev(brewer.pal(n = 7, name ="Oranges")))(100),
    fontsize = 8, cluster_cols = FALSE, cluster_rows= FALSE, 
    cellwidth = 15, cellheight = 5, 
    fontsize_row = 5, main = "LINEs length",
    labels_row = line_tab$Lengths, 
    filename = "heat_map.png", heigh = 1.75, width = 2.8)

在此处输入图像描述

请让我知道这是否是您需要获得的。

请注意，这pheatmap(log10(line_tab+1),...)不起作用，因为pheatmap()需要一个数字矩阵。这就是为什么我必须排除前两列。

于 2019-07-03T14:36:33.383 回答