2

我想使用包中的heatmap.2函数创建一个热图gplots。这是一个最小的例子。

require(gplots)

# create symmetric matrix
x = matrix(rnorm(100), nrow=10)
diag(x)         <- 1
x[upper.tri(x)] <-  t(x)[upper.tri(x)]
colnames(x)     <- rownames(x) <- letters[1:nrow(x)]

# create side colours
varcols = setNames(rainbow(nrow(x)), rownames(x))

# create heatmap
heatmap.2(x,
          symm = TRUE,             
          trace = "none",          
          revC=TRUE, # <-- THIS IS THE PROBLEM
          ColSideColors = varcols,
          RowSideColors = varcols
          )

问题是sidecolors。x是一个对称矩阵，因此列和行应该具有相同的侧色。只要revC = FALSE. 但是，当我使用revC = TRUE颜色的顺序时就搞砸了。有时 - 在小例子中 - 它有助于扭转ColSideColors，但这并不总是有效。

我做错了什么还是这是一个gplots错误？

4

1 回答 1

4

对于遇到此问题的其他任何人，这就是我解决它的方法：

thing = heatmap.2(my_matrix,...RowSideColors=row_cols, revC=F)
ordinary_order = thing$rowInd
reversal = cbind(ordinary_order, rev(ordinary_order))
rev_col = row_cols[reversal[,2]]; rev_col = rev_col[order(reversal[,1])];

heatmap.2(my_matrix, RowSideColors=rev_col, revC=T)

于 2016-07-11T15:19:23.413 回答