[R-es] Grafico en 3D

[Marcuzzi, Javier Rubén]

Estimado Jorge

Miro este ejemplo, que funciona, pero me intriga su problema que no puede 
resolver:

wireframe(volcano, shade = TRUE,
          aspect = c(61/87, 0.4),
          light.source = c(10,0,10))

g <- expand.grid(x = 1:10, y = 5:15, gr = 1:2)
g$z <- log((g$x^g$g + g$y^2) * g$gr)
wireframe(z ~ x * y, data = g, groups = gr,
          scales = list(arrows = FALSE),
          drape = TRUE, colorkey = TRUE,
          screen = list(z = 30, x = -60))

Pero veo dos cosas diferentes que honestamente no las se, el expand.grid y 
el groups

Javier Marcuzzi

-----Mensaje original----- 
From: Jorge I Velez
Sent: Wednesday, March 28, 2012 3:59 PM
To: R-help-es
Subject: [R-es] Grafico en 3D

Buenas tardes,

Podria por favor alguien ayudarme a descifrar por que no es posible ver la
grafica 3D cuando se utiliza lattice?  Ya se que algunas veces los duendes
hacen de las suyas, pero no estoy seguro de que este sea el caso.  He
reiniciado R, cambiado a RStudio y utilizado XP, pero no pasa "nada".

Mi sessionInfo() es:

R version 2.14.0 Patched (2011-11-12 r57642)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] copula_0.9-9         pspline_1.0-14       mvtnorm_0.9-9992
[4] scatterplot3d_0.3-33 locfit_1.5-6         akima_0.5-7
[7] lattice_0.20-6

loaded via a namespace (and not attached):
[1] grid_2.14.0  tools_2.14.0


Y este es el codigo que estoy usando:

# mi funcion
mypdf <- function(x, y){
    rho <- cor(x, y)
    one <- 1/(sqrt(1-rho^2))
    two <- exp((-1/2*(1-rho^2))*(x^2+y^2-2*rho*x*y))
    three <- exp(.5*(x^2+y^2))
    one*two*three
    }

# ejemplo
set.seed(123)
x <- rnorm(100)
y <- rnorm(100)
z <- mypdf(x, y)
dd <- data.frame(x = x, y = y, z = z)
head(dd)

# ventana grafica
par(mfrow = c(1, 2), pty = 's')

# marginales -- OK!
plot(x, z, las = 1, ylab = 'f(x, y)')
plot(y, z, las = 1, ylab = 'f(x, y)')

# con lattice -- no veo nada!
require(lattice)  # no plot
wireframe(z ~ x*y, data = dd, scales = list(arrows = FALSE),
         drape = TRUE, colorkey = TRUE, col = 1, zlim = c(0.9, 1.2))

# con scatterplot3d -- si aparece!
require(scatterplot3d)  # plot
par(mfrow = c(1,1), pty = 's')
scatterplot3d(x, y, z, highlight.3d=TRUE, col.axis=1,
      col.grid=1,pch=20)

# otro ejemplo con lattice  -- si funciona!
wireframe(volcano, shade = TRUE,
          aspect = c(61/87, 0.4), light.source = c(10,0,10))


Muchas gracias,
Jorge.-

[[alternative HTML version deleted]]

_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es