[R-br] Como adicionar ponto(s) de sela(s) em superfície resposta?

Ivan Bezerra Allaman ivanalaman em yahoo.com.br
Sexta Junho 10 22:27:57 BRT 2011


Boa noite senhores!!

Após ler quase todos os tópicos do nabble internacional de como adicionar pontos com "wireframe" e manipular de tudo qualquer jeito as funções panel.3d.INFERNO e outras milhões de funções da lattice para tentar fazer uma coisa simples, venho recorrer novamente a ajuda dos senhores. Tenho o seguinte 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)
    }
  }
wireframe(z~x+y, da, drape=TRUE, scales=list(arrows=FALSE),col.regions=jet.colors(100),panel.3d.wireframe="panel.3d.contour")
 Gostaria de indicar no gráfico o(s) ponto(s) (coordenadas x,y e z) de sela, para tornar o gráfico autoexplicativo.

Desde já agradeço mais uma vez 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/20110610/56a6281a/attachment.html>


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