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]]
______________________________________________
[email protected] 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.
