From: Paul Murrell <p.murrell_at_auckland.ac.nz>

Date: Tue 31 May 2005 - 10:14:34 EST

z[is.na(z)] <- 1

op <- par(bg = "white")

persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue",

z2[xi, yi] <- z[xi, yi]

persp(x, y, z2, theta = 30, phi = 30, expand = 0.5,

Date: Tue 31 May 2005 - 10:14:34 EST

Hi

jjorgensen@fastmail.fm wrote:

> Hello R-sters,

*>
**> I'm trying to add several lines to a response surface that I've plotted
**> using persp(). I've tried lines() using the "trans3d" function but I've
**> been unsuccessful in getting it to work (R v2.0.1). Essentially, I'm
**> trying to highlight one or more of the surface wireframe lines in a
**> bolder (or different) color. Any tips from those of you who have some
**> experience with this would be greatly appreciated. [Would it be easier
**> using wireframe() in library(lattice) instead?]
*

Here's an example that just overlays two persp() plots ...

x <- seq(-10, 10, length= 30) y <- x f <- function(x,y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r } z <- outer(x, y, f)

z[is.na(z)] <- 1

op <- par(bg = "white")

persp(x, y, z, theta = 30, phi = 30, expand = 0.5, col = "lightblue",

zlim=range(z))

# overlay plot with just "highlighted" surface area par(new=TRUE)

z2 <- matrix(NA, ncol=30, nrow=30) xi <- 15:18 yi <- 13:14

z2[xi, yi] <- z[xi, yi]

persp(x, y, z2, theta = 30, phi = 30, expand = 0.5,

zlim=range(z), border="red", col="pink", box=FALSE, axes=FALSE)

> And, any suggestions on how to add text outside of the persp() plot next

*> to the highlighted line would be much appreciated.
*

Do you mean like a legend?

Paul

-- Dr Paul Murrell Department of Statistics The University of Auckland Private Bag 92019 Auckland New Zealand 64 9 3737599 x85392 paul@stat.auckland.ac.nz http://www.stat.auckland.ac.nz/~paul/ ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.htmlReceived on Tue May 31 10:21:04 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:32:16 EST
*