Dear R-users,

	I would like to place two 3D plots side-by-side in a rgl-setting. It would nice to have something like "par(mfrow=c(1,2))" for basic plots, or an array framework for wireframe(lattice) (see example below). I only managed to overlap two persp3d plots.
My final idea would be to animate both surfaces using play3d(rgl).

Thanks in advance for any help.

Best,
Carlo Giovanni Camarda


# simple datasets
x <- seq(-10, 10, length= 30)
y <- x
f <- function(x,y) { r <- sqrt(x^2+y^2); 10 * sin(r)/r }
z1 <- outer(x, y, f)
z1[is.na(z1)] <- 1
z2 <- matrix(outer(1:30, 1:30)/200, 30, 30)

# simple persp
par(mfrow=c(1,2))
persp(x, y, z1)
persp(x, y, z2)

# wireframe
library(lattice)
grid. <- expand.grid(list(x=x,y=y,name=c("Z1", "Z2")))
grid.$Z <- c(c(z1), c(z2))
wireframe(Z ~ x * y | name, data = grid.)

# single plot
library(rgl)
persp3d(x, y, z1)
persp3d(x, y, z2)

# overlapping plots
persp3d(x,y,z1, color="#CCCCFF", front="lines") 
persp3d(x,y,z2, col="red", add=TRUE) 
play3d(spin3d(axis=c(0,0,1), rpm=10), duration=5)



-----------------------------------------------------
Dr. Carlo Giovanni Camarda
Research Scientist
Max Planck Institute for Demographic Research
Laboratory of Statistical Demography
Konrad-Zuse-Straße 1
18057 Rostock - Germany
Phone: +49 (0)381 2081 172
Fax: +49 (0)381 2081 472
camarda@demogr.mpg.de
http://www.demogr.mpg.de/en/staff/camarda/default.htm
-----------------------------------------------------





----------
This mail has been sent through the MPI for Demographic ...{{dropped:10}}


