[R-br] Plotar dois tipos de gráficos em uma mesma janela gráfica com a lattice.

Ivan Bezerra Allaman ivanalaman em yahoo.com.br
Quarta Junho 15 09:55:28 BRT 2011


Bom dia senhores!


Estou tentando plotar dois tipos de gráficos em uma mesma janela gráfica, mais tá difícil. Não sei até onde vale a pena o 'tempo/benefício' em se usar a lattice. Segue abaixo o CMR.

require(plotrix)

jet.colors <- colorRampPalette( c("blue", "green") )  

x <- seq(-1.95, 1.95, length=30)
y <- seq(-1.95, 1.95, length=35)
da <- expand.grid(x=x, y=y)
da$z <- with(da, x*y^2)

require(lattice)
panel.3d.contour <-
  function(x, y, z, rot.mat, distance,
           nlevels = 20, zlim.scaled, ...)
  {
    add.line <- trellis.par.get("add.line")
    panel.3dwire(x, y, z, rot.mat, distance,
                 zlim.scaled = zlim.scaled, ...)
    clines <-
      contourLines(x, y, matrix(z, nrow = length(x), byrow = TRUE),
                   nlevels = nlevels)
    for (ll in clines) {
      m <- ltransform3dto3d(rbind(ll$x, ll$y, zlim.scaled[1]),
                            rot.mat, distance)
      panel.lines(m[1,], m[2,], col = add.line$col,
             
     lty = add.line$lty, lwd = add.line$lwd)
    }
  }

g1 <- wireframe(z~x+y, da, drape=TRUE, scales=list(arrows=FALSE),col.regions=jet.colors(100),panel.3d.wireframe="panel.3d.contour")

g2 <- contourplot(z ~ x+y,da) 

library(grid)
pushViewport(viewport(layout=grid.layout(2,1)))
print(g1,vp=viewport(layout.pos.row=1,layout.pos.col=1))
print(g2,vp=viewport(layout.pos.row=2,layout.pos.col=1))

Este último comando não plota os dois gráficos em uma mesma janela gráfica. Eu gostaria de um comando equivalente ao par(mfrow). 

Desde já agradeço novamente a atenção de todos.

Allaman
(S,f,P)
 
M.Sc Ivan Bezerra Allaman 
Zootecnista
Doutorando em Produção Animal/Aquicultura - UFLA 
email e msn - ivanalaman em yahoo.com.br 
Tel: (35)3826-6608/9900-2924
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110615/3f7e9eef/attachment.html>


Mais detalhes sobre a lista de discussão R-br