R help - Dec 2011 - How to color a region in a contour plot with the contour being the boundary?

Dear expeRts,

I would like to color a certain region in a levelplot. The region for z <=
0.02 should have a dark gray color. Below is a minimal example. It almost does
what I want, but The region between z=0.02 and z=1 is also colored in dark gray
(instead of just the region for z <= 0.02).
How can I solve this?

Cheers,

Marius


## z values for given x and y
f <- function(x) 4*((1-x)^(-1/2)-1)
eps <- 1e-12
x <- y <- seq(eps, 1-eps, length.out=500)
z <- outer(x, y, FUN=function(x, y) pmax(f(x) + f(y) - 10, 0))

## determine colors
zcols <- c(gray(0.2), gray(seq(0.5, 1, length.out=50)),
rep("#FFFFFF", max(z)-51)) # colors with dark gray in the beginning
for z <= 0.02 and many whites for z > 50

## trial 1
pdf(file="levelplot1.pdf", width=6, height=6)
image(x, y, z, xaxs="r", yaxs="r", col=zcols)
contour(x, y, z, levels=c(0.02, 1, 5, 10, 20, 50, 500), add=TRUE)
dev.off()

## trial 2
pdf(file="levelplot2.pdf", width=6, height=6)
image(x, y, z, xaxs="r", yaxs="r", col=zcols, oldstyle=TRUE)
contour(x, y, z, levels=c(0.02, 1, 5, 10, 20, 50, 500), add=TRUE)
dev.off()

Try the levelplot() function in package lattice.

library(lattice)

dat <- expand.grid(x=x, y=y)
dat$z <- pmax(f(dat$x) + f(dat$y) - 10, 0)

levelplot(z ~ x * y, dat, 
        at=c(-1, 0.02, 1, 5, 10, 20, 50, 500, 9000000), 
        labels=TRUE, contour=TRUE, colorkey=FALSE, 
        col.regions=gray(c(0.2, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1)))

Jean


Marius Hofert wrote on 12/31/2011 05:02:22 AM:
> Dear expeRts,
> 
> I would like to color a certain region in a levelplot. The region 
> for z <= 0.02 should have a dark gray color. Below is a minimal 
> example. It almost does what I want, but The region between z=0.02 
> and z=1 is also colored in dark gray (instead of just the region for
> z <= 0.02).
> How can I solve this?
> 
> Cheers,
> 
> Marius
> 
> 
> ## z values for given x and y
> f <- function(x) 4*((1-x)^(-1/2)-1)
> eps <- 1e-12
> x <- y <- seq(eps, 1-eps, length.out=500)
> z <- outer(x, y, FUN=function(x, y) pmax(f(x) + f(y) - 10, 0))
> 
> ## determine colors
> zcols <- c(gray(0.2), gray(seq(0.5, 1, length.out=50)), rep
> ("#FFFFFF", max(z)-51)) # colors with dark gray in the beginning
for
> z <= 0.02 and many whites for z > 50
> 
> ## trial 1
> pdf(file="levelplot1.pdf", width=6, height=6)
> image(x, y, z, xaxs="r", yaxs="r", col=zcols)
> contour(x, y, z, levels=c(0.02, 1, 5, 10, 20, 50, 500), add=TRUE)
> dev.off()
> 
> ## trial 2
> pdf(file="levelplot2.pdf", width=6, height=6)
> image(x, y, z, xaxs="r", yaxs="r", col=zcols,
oldstyle=TRUE)
> contour(x, y, z, levels=c(0.02, 1, 5, 10, 20, 50, 500), add=TRUE)
> dev.off()
	[[alternative HTML version deleted]]

I was hoping for a solution without lattice. 
But it solves the problem, thanks Jean.

Cheers,

Marius


On 2012-01-03, at 19:06 , Jean V Adams wrote:
> 
> Try the levelplot() function in package lattice. 
> 
> library(lattice) 
> 
> dat <- expand.grid(x=x, y=y) 
> dat$z <- pmax(f(dat$x) + f(dat$y) - 10, 0) 
> 
> levelplot(z ~ x * y, dat, 
>         at=c(-1, 0.02, 1, 5, 10, 20, 50, 500, 9000000), 
>         labels=TRUE, contour=TRUE, colorkey=FALSE, 
>         col.regions=gray(c(0.2, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1))) 
> 
> Jean 
> 
> 
> Marius Hofert wrote on 12/31/2011 05:02:22 AM:
> 
> > Dear expeRts,
> > 
> > I would like to color a certain region in a levelplot. The region 
> > for z <= 0.02 should have a dark gray color. Below is a minimal 
> > example. It almost does what I want, but The region between z=0.02 
> > and z=1 is also colored in dark gray (instead of just the region for
> > z <= 0.02).
> > How can I solve this?
> > 
> > Cheers,
> > 
> > Marius
> > 
> > 
> > ## z values for given x and y
> > f <- function(x) 4*((1-x)^(-1/2)-1)
> > eps <- 1e-12
> > x <- y <- seq(eps, 1-eps, length.out=500)
> > z <- outer(x, y, FUN=function(x, y) pmax(f(x) + f(y) - 10, 0))
> > 
> > ## determine colors
> > zcols <- c(gray(0.2), gray(seq(0.5, 1, length.out=50)), rep
> > ("#FFFFFF", max(z)-51)) # colors with dark gray in the
beginning for
> > z <= 0.02 and many whites for z > 50
> > 
> > ## trial 1
> > pdf(file="levelplot1.pdf", width=6, height=6)
> > image(x, y, z, xaxs="r", yaxs="r", col=zcols)
> > contour(x, y, z, levels=c(0.02, 1, 5, 10, 20, 50, 500), add=TRUE)
> > dev.off()
> > 
> > ## trial 2
> > pdf(file="levelplot2.pdf", width=6, height=6)
> > image(x, y, z, xaxs="r", yaxs="r", col=zcols,
oldstyle=TRUE)
> > contour(x, y, z, levels=c(0.02, 1, 5, 10, 20, 50, 500), add=TRUE)
> > dev.off()