我正在使用 r 中的 pheatmap 生成热图。我正在绘制基于倍数变化上调或下调的错误发现率校正 p 值。我希望色阶从 -0.1 分级到 0。然后在 0 处切换颜色并再次从 0 分级到 0.1。如何使色标仅在一侧分级?
data=read.csv("metabolism_heatmap.csv")
rownames(data)=data$Gene
data$Category=as.character(data$Category)
data=data[order(data$Category),]
#rows 1-29 = one category
#cols 11-12 = signed FDR-value for 2 groups of samples
colors2<- colorRampPalette(c("white", "blue", "red", "white"))(n = 100)
pheatmap(data[1:29,c(11,12)],
color= colors2,
cluster_rows = FALSE,
breaks = seq(-0.1,0.1,length.out = 101))
这是情节。我想摆脱紫色并在零处有一个清晰的边界。 热图