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

[Walmes Zeviani]

Ivan,

Modifiquei a função para sinalizar a coordenada do ponto de sela, veja

panel.3d.contour <- function(x, y, z, zlim.scaled,
                             rot.mat, distance, nlevels=40,
                             X, Y, Z, pcol=2, pcex=2, ppch="+", ...){
  panel.3dwire(x=x, y=y, z=z, rot.mat=rot.mat,
               distance=distance, zlim.scaled=zlim.scaled, ...)
  panel.3dscatter(x=X, y=Y, z=Z, dist=distance,
                  zlim.scaled=zlim.scaled, type="p", col=pcol,
                  cex=pcex, pch=ppch, ...)
  add.line <- trellis.par.get("add.line")
  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), nlevels=10,
          X=0, Y=0, Z=0, ppch=19, pcex=1.2, pcol="white",
          col.regions=jet.colors(100), panel.3d.wireframe=panel.3d.contour)

À disposição.
Walmes.


==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================


2011/6/11 Walmes Zeviani <walmeszeviani em gmail.com>

> Ivan,
>
> A última mensagem dessa seqüência de mensagens faz a adição de um ponto a
> superfície. Tente partir daí. O segundo link também pode ser útil.
>
>
> http://r.789695.n4.nabble.com/lattice-wireframe-quot-eats-up-quot-points-how-to-make-points-on-wireframe-visible-td3334836.html
>
>
> http://stackoverflow.com/questions/1406202/plotting-a-wireframe-and-a-cloud-with-lattice-in-r
>
> À disposição.
> Walmes.
>
> ==========================================================================
> Walmes Marques Zeviani
> LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
> Departamento de Estatística - Universidade Federal do Paraná
> fone: (+55) 41 3361 3573
> VoIP: (3361 3600) 1053 1173
> e-mail: walmes em ufpr.br
> twitter: @walmeszeviani
> homepage: http://www.leg.ufpr.br/~walmes
> linux user number: 531218
> ==========================================================================
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110611/79dce6aa/attachment.html>