R help - Dec 2004 - Arrow and text on a 3d plot

I am currently working on surface plots (with R 2.0.1) and I want to
improve the aspect of them. Having a typical 'cow-boy hat' surface
[source: Spector (1994), p. 206]:

pts <- seq(from = -8, to = 8, length = 50)
cow <- function(x, y)
{
z <- sqrt(x^2+y^2)
sin(z)/z
}
out <- outer(pts, pts, cow)
sur <- persp(x=pts, y=pts, z=out, theta=45, phi=30)

What I exactly wanted is:

a) how to add an arrow pointing the maximum point in the surface, for
example: a short arrow starting on the maximum point, with slope 45??
and finishing near the top-edge of the plot.

b) how to add a typical text at the end of the arrow, as "this is the
maximum point", as described in the above reference.

N.B.: the author used the S function "perspp", but R does not have it.

Could anyone help me? 

Thank you very much in advance.

Arturo De Zan
PhD candidate
UPC-Universitat Politecnica de Catalunya
Barcelona, Spain.

Hi

the best way to do it is to make a function like the trans3d() function 
in help(persp)
and use it in your own function.  trans3d()
translates from (x,y,z) coords to plottable positions.  The following  
works for me:



  xyz <- matrix(rnorm(60),20,3)
xyz1 <- xyz+rnorm(60)/10

p3dpairs(xyz,xyz1,col="red",theta=40,r=1e9)

[ignore the warnings]

with

"p3dpairs" <-
function(x,x1, xlim=NULL,ylim=NULL,zlim=NULL,col=par("col"), 
pch=par("pch"), cex=par("cex"), ...){
   if(is.matrix(x)){
     z <- x[,3]
     y <- x[,2]
     x <- x[,1]
   }

     if(is.matrix(x1)){
     z1 <- x1[,3]
     y1 <- x1[,2]
     x1 <- x1[,1]
   }

   if(missing(zlim)) {
     z.grid <- matrix(range(z),2,2)
   } else {
     z.grid <- matrix(zlim,2,2)
   }

   if(missing(xlim)){ xlim <- range(x) }
   if(missing(ylim)){ ylim <- range(y) }

   persp(xlim, ylim, z.grid, col = NA, border=NA, ...) -> res

   trans3d <- function(x,y,z, pmat) {
     tr <- cbind(x,y,z,1) %*% pmat
     list(x = tr[,1]/tr[,4], y= tr[,2]/tr[,4])
   }

   out <- trans3d(x,y,z,pm=res)
   out1 <- trans3d(x1,y1,z1,pm=res)
   points(out, col=col, pch=pch, cex=cex, ...)

   for(i in 1:length(out$x)){
     lines(c(out$x[i],out1$x[i]),c(out$y[i],out1$y[i]), col="gray",
...)
   }
   return(invisible(out))
}


On Dec 2, 2004, at 11:14 am, Arturo T. De Zan wrote:
> I am currently working on surface plots (with R 2.0.1) and I want to
> improve the aspect of them. Having a typical 'cow-boy hat' surface
> [source: Spector (1994), p. 206]:
>
> pts <- seq(from = -8, to = 8, length = 50)
> cow <- function(x, y)
> {
> z <- sqrt(x^2+y^2)
> sin(z)/z
> }
> out <- outer(pts, pts, cow)
> sur <- persp(x=pts, y=pts, z=out, theta=45, phi=30)
>
> What I exactly wanted is:
>
> a) how to add an arrow pointing the maximum point in the surface, for
> example: a short arrow starting on the maximum point, with slope 45??
> and finishing near the top-edge of the plot.
>
> b) how to add a typical text at the end of the arrow, as "this is the
> maximum point", as described in the above reference.
>
> N.B.: the author used the S function "perspp", but R does not
have it.
>
> Could anyone help me?
>
> Thank you very much in advance.
>
> Arturo De Zan
> PhD candidate
> UPC-Universitat Politecnica de Catalunya
> Barcelona, Spain.
>
> ______________________________________________
> R-help at 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.html
>