Colors for heatmap in R

Question

I'm making a heatmap in R, but being very new to R, I have some questions:

My data is a big matrix 21 columns and 89 rows, containing numbers from 0 to 16. I would like to get the heatmap colored in a heatmappy way from 0 (white) to 16 (dark red - or any color for that sake). Or maybe even fancier, have a color palette going from 0 to <10 (so that the points having above 10 "hits" get the same color).

Can anyone help me with this ? Thanks alot!

My code:

library(ggplot2)
library("RColorBrewer")

AS <- read.csv("L:/Pseudoalteromonas/Heatmap antismash/HM_phyl.csv", sep=";")

row.names(AS) <- AS$Strain

AS <- AS[,2:21]

## The colors you specify.
my_palette <- colorRampPalette(c("white", "yellow","orange", "red"))(n = 299)

AS_matrix <- data.matrix(AS)

AS_heatmap <- heatmap(AS_matrix, Rowv=NA, Colv=NA, col = my_palette, scale="row", margins=c(5,10))

My data looks like this:

tail(HM)
   Sideophore Bacteriocin Aryl.polyene Nrps T1pks T2pks T3pks T1pks.Nrps Lantipeptide Terpene Hserlactone Transatpks

S4048          0           2            0    2     0     0     0          0            1       0           0          1
S3655          1           2            2    0     0     0     0          0            0       0           0          0
S4060          0           2            0    7     0     1     1          2            1       0           0          1
S2607          0           2            0   10     1     1     1          4            1       0           0          1
S4054          0           2            1    3     0     0     0          4            1       0           1          1
S4047          0           2            1    7     0     0     0          4            1       0           1          1
  Butyrolactone Indole Thiopeptide Ladderane Pufa Resorcinol Otherks Other
S4048             0      0           0         0    0          0       0     0
S3655             0      0           0         0    0          1       0     0
S4060             0      1           0         0    0          0       0     2
S2607             0      1           0         0    0          0       0     2
S4054             0      1           0         1    0          0       0     0
S4047             0      1           0         1    0          0       0     2

Answer 1

You could keep the data as a data.frame and use ggplot2 (it looks like that is what you are intending since you called ggplot2?)

library(ggplot2)
library(RColorBrewer)
library(tidyverse)

set.seed(12343)
# create matrix with 21 columns and 89 rows
# with numbers between 0 - 16

AS <- runif(n= 1869, min = 0, max = 16) %>%
  matrix(., nrow = 89)

colnames(AS) <- LETTERS[1:21]

AS <- as.data.frame(AS)
AS$train <- 1:89

AS <- gather(AS, A:U, key = "colname", value="value")

ggplot(AS, aes(x = colname, y=train)) + geom_tile(aes(fill = value), 
                                                  colour = "white") +
  scale_fill_distiller(palette = "Reds", limits = c(0,10), na.value = "#de2d26",
                       direction = 1, labels = c(0.0, 2.5, 5.0, 7.5, "> 10.0"))

So using your code, maybe something like this:

library(ggplot2)
library(RColorBrewer)
library(tidyverse)


AS <- read.csv("L:/Pseudoalteromonas/Heatmap antismash/HM_phyl.csv", sep=";")

# assume "train" is the row indicator, so we will use
# gather with -train argument to gather all columns but "train"

AS <- gather(AS, -train, key = "colname", value="value")


ggplot(AS, aes(x = colname, y=train)) + geom_tile(aes(fill = value), 
                                                  colour = "white") +
  scale_fill_distiller(palette = "Reds", limits = c(0,10), na.value = "#de2d26",
                       direction = 1, labels = c(0.0, 2.5, 5.0, 7.5, "> 10.0"))

# na.values (values >10) take maximum red color in Reds colorbrewer