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@ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================