Hi there,
I'd like to create a heatmap from my matrix with
a) a defined color range (lets say from yellow to red)
b) using striking colors above and below a certain threshold (above =
"green",
below = "blue")
Example matrix (there should be a few outliers generated...) + simple levelplot
without outliers marked:
library(lattice)
my.mat <- matrix(rnorm(800), nrow = 40)
threshold <- c(-1,1) # should be used for the extreme colors
colorFun <- colorRampPalette(c("yellow","red"))
levelplot(my.mat, col.regions = colorFun(50))
I don't know how to handle the extrem values...
Can anybody help?
Antje
I played a little bit around and got the following solution which works for
now, though it seems to be too complicated to me.
If anybody else know another solution - please let me know!!!
library(lattice)
my.mat <- matrix(rnorm(800), nrow = 40)
colorFun <- colorRampPalette(c("yellow","red"))
b <- boxplot(my.mat, plot = FALSE)
thr <- c(b$stats[1],b$stats[5])
col.bins <- 100
step <- abs(thr[2] - thr[1])/50
limit <- ifelse(min(my.mat) > thr[1] - step, min(my.mat) - step,
min(my.mat))
lp <- rev(seq(thr[1] - step, limit - step, -step))
mp <- seq(thr[1], thr[2], step)
limit <- ifelse(max(my.mat) < thr[2] + step, max(my.mat) + step,
max(my.mat))
up <- seq(thr[2] + step, limit + step, step)
my.at <- c(lp,mp,up)
my.col.regions <- c(rep("green", length(lp)), colorFun(length(mp)),
rep("blue",
length(up)) )
levelplot(my.mat, at = my.at, col.regions = my.col.regions)
Antje schrieb:> Hi there,
>
> I'd like to create a heatmap from my matrix with
> a) a defined color range (lets say from yellow to red)
> b) using striking colors above and below a certain threshold (above =
> "green", below = "blue")
>
> Example matrix (there should be a few outliers generated...) + simple
> levelplot without outliers marked:
>
> library(lattice)
> my.mat <- matrix(rnorm(800), nrow = 40)
> threshold <- c(-1,1) # should be used for the extreme colors
> colorFun <- colorRampPalette(c("yellow","red"))
> levelplot(my.mat, col.regions = colorFun(50))
>
>
> I don't know how to handle the extrem values...
>
> Can anybody help?
>
> Antje
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
I played a little bit around and got the following solution which works for
now, though it seems to be too complicated to me.
If anybody else know another solution - please let me know!!!
library(lattice)
my.mat <- matrix(rnorm(800), nrow = 40)
colorFun <- colorRampPalette(c("yellow","red"))
b <- boxplot(my.mat, plot = FALSE)
thr <- c(b$stats[1],b$stats[5])
col.bins <- 100
step <- abs(thr[2] - thr[1])/50
limit <- ifelse(min(my.mat) > thr[1] - step, min(my.mat) - step,
min(my.mat))
lp <- rev(seq(thr[1] - step, limit - step, -step))
mp <- seq(thr[1], thr[2], step)
limit <- ifelse(max(my.mat) < thr[2] + step, max(my.mat) + step,
max(my.mat))
up <- seq(thr[2] + step, limit + step, step)
my.at <- c(lp,mp,up)
my.col.regions <- c(rep("green", length(lp)), colorFun(length(mp)),
rep("blue",
length(up)) )
levelplot(my.mat, at = my.at, col.regions = my.col.regions)
Antje schrieb:> Hi there,
>
> I'd like to create a heatmap from my matrix with
> a) a defined color range (lets say from yellow to red)
> b) using striking colors above and below a certain threshold (above =
> "green", below = "blue")
>
> Example matrix (there should be a few outliers generated...) + simple
> levelplot without outliers marked:
>
> library(lattice)
> my.mat <- matrix(rnorm(800), nrow = 40)
> threshold <- c(-1,1) # should be used for the extreme colors
> colorFun <- colorRampPalette(c("yellow","red"))
> levelplot(my.mat, col.regions = colorFun(50))
>
>
> I don't know how to handle the extrem values...
>
> Can anybody help?
>
> Antje
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>