创建具有下三角形的二次矩阵作为热图和上三角字符数据

时间:2015-07-25 10:09:07

标签: r panel heatmap

我有以下(示例)数据,我想绘制: matrix(c("Variable1",0.4,0.5,"string1","Variable2",0.3, "string2","string3","Variable3"),ncol=3)

我想用矩阵的数字0.4,0.5,0.3创建热图,并使用下三角形:

  numeric<- matrix(c(1,0.4,0.5,0.4,1,0.3,0.5,0.3,1),ncol=3)
  require(graphics)
  require(grDevices)
  library(colorRamps)

  library(lattice)
  x <- seq(pi/4, 5*pi, length.out=10)
  y <- seq(pi/4, 5*pi, length.out=10)
  grid <- expand.grid(X=x, Y=y)

   myPanel <- function(x, y, z, ...) { 
   panel.levelplot(x,y,z,...)
   panel.text(x, y, round(z,2),cex=0.578,...)}
  levelplot(numeric,panel=myPanel)

在包含热图的对角线的上三角形上,我希望在给定矩阵的同一网格位置上有值"string1", "string2", "string3", "Variable1", "Variable2", "Variable3"的字符!

我想到了以某种方式创建一个网格并使用图形包的image()函数为数字着色,但这对我不起作用。

你有一个简单而优雅的功能,你只需要给定的矩阵作为输入吗?

1 个答案:

答案 0 :(得分:1)

您可以尝试使用ComplexHeatmap包来自定义热图体:

首先创建一个对称矩阵:

set.seed(123)
mat = matrix(rnorm(100), 10)
colnames(mat) = letters[1:10]

mat2 = cor(mat)

进行群集并调整行和列的顺序:

od = hclust(dist(mat2))$order
mat2 = mat2[od, od]

现在我们可以自定义热图。在以下代码中,col控制如何将值映射到颜色。 由于已应用群集,因此我们将cluster_rowscluster_columns设置为FALSE。 在type中将none设置为gpar()时,会初始化热图,但不会添加任何内容, 然后我们可以使用自定义函数cell_fun将图形添加到其中。

cell_fun将应用于热图中的每个小网格。 在cell_fun中,有七个参数:

  • jmat2
  • 中的列索引
  • imat2
  • 中的行索引
  • x:情节上的位置
  • y:情节上的位置
  • w:当前网格的宽度
  • h:当前网格的高度
  • col:当前网格的填充颜色

使用Heatmap()生成最终热图:

Heatmap(mat2, col = colorRamp2(c(-1, 0, 1), c("green", "white", "red")),
    cluster_rows = FALSE, cluster_columns = FALSE,
    heatmap_legend_title = "cor", rect_gp = gpar(type = "none"), 
    cell_fun = function(j, i, x, y, w, h, col) {
        if(i > j) {
            grid.rect(x, y, w, h, gp = gpar(fill = col))
        } else if(j == i) {
            grid.text(labels[i], x, y)
        } else {
            grid.text(sprintf("%.2f", mat2[i, j]), x, y)
        }
        grid.rect(x, y, w, h, gp = gpar(fill = NA, col = "black"))
    })

enter image description here