R：在水平图中将中心red_to_blue调色板设置为0

Question

我正在创建一个levelplot，其中我的数据框的一个变量用于为单元格着色（fold.change），而另一个变量（map.signif）用于顶部。在这种情况下，我为重要的单元格写*和**。

这是我的MWE：

set.seed(150)
pv.df <- data.frame(compound=rep(LETTERS[1:8], each=3), comparison=rep(c("a/b","b/c","a/c"), 8), p.value=runif(24, 0, 0.2), fold.change=runif(24, -0.3, 0.9))
pv.df$map.signif <- ifelse(pv.df$p.value > 0.05, "", ifelse(pv.df$p.value > 0.01,"*", "**"))
pv.df
myPanel <- function(x, y, z, ...) {
  panel.levelplot(x, y, z, ...)
  panel.text(x, y, pv.df$map.signif, cex=3)
}
#install.packages("latticeExtra")
library(latticeExtra)
library(RColorBrewer)
cols <- colorRampPalette(brewer.pal(11, "RdBu"))(11)
png(filename="test.png", height=800, width=400)
print(
    levelplot(fold.change ~ comparison*compound, #p.value instead of p.adjust depending on map.signif
      pv.df,
      panel = myPanel,
      col.regions = cols,
      at = do.breaks(range(pv.df$fold.change), 11),
      colorkey = list(col = cols, 
                      at = do.breaks(range(pv.df$fold.change), 11)),
      xlab = "", ylab = "",              # remove axis titles
      scales = list(x = list(rot = 45),  # change rotation for x-axis text
                    cex = 0.8),          # change font size for x- & y-axis text
      main = list(label = "Test\nfold change color\n*pv<0.05\t**pv<0.01",
                  cex = 1.5))
)
dev.off()

产生：

我的问题是：由于fold.change包含负值和正值，如何在调色板中使0与白色重合，因此负值为红色，正值为蓝色？

对于胜利，是否有可能在细胞背景清晰时用黑色表示*，在背景暗时用白色表示？非常感谢！

Answer 1

您的范围不是simetrical。另一种选择是：

max_abs <- max(abs(pv.df$fold.change))
brk <- do.breaks(c(-max_abs, max_abs), 11)
levelplot(fold.change ~ comparison*compound, #p.value instead of p.adjust depending on map.signif
          pv.df,
          panel = myPanel,
          col.regions = cols,
          at = brk,
          colorkey = list(col = cols, 
                          at = brk),
          xlab = "", ylab = "",              # remove axis titles
          scales = list(x = list(rot = 45),  # change rotation for x-axis text
                        cex = 0.8),          # change font size for x- & y-axis text
          main = list(label = "Test\nfold change color\n*pv<0.05\t**pv<0.01",
                      cex = 1.5))

修改

如果你不想要额外的休息时间：

max_abs <- max(abs(pv.df$fold.change))
brk <- do.breaks(c(-max_abs, max_abs), 11)
first_true <- which.max(brk > min(pv.df$fold.change))
brk <- brk[(first_true -1):length(brk)]
cols <- cols[(first_true -1):length(cols)]
levelplot(fold.change ~ comparison*compound, #p.value instead of p.adjust depending on map.signif
          pv.df,
          panel = myPanel,
          col.regions = cols,
          at = brk,
          colorkey = list(col = cols, 
                          at = brk),
          xlab = "", ylab = "",              # remove axis titles
          scales = list(x = list(rot = 45),  # change rotation for x-axis text
                        cex = 0.8),          # change font size for x- & y-axis text
          main = list(label = "Test\nfold change color\n*pv<0.05\t**pv<0.01",
                      cex = 1.5))