Como adicionar ponto(s) de sela(s) em superfície resposta?
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@yahoo.com.br Tel: (35)3826-6608/9900-2924
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... http://stackoverflow.com/questions/1406202/plotting-a-wireframe-and-a-cloud-... À 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 ==========================================================================
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 ========================================================================== 2011/6/11 Walmes Zeviani <walmeszeviani@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...
http://stackoverflow.com/questions/1406202/plotting-a-wireframe-and-a-cloud-...
À 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 ==========================================================================
participantes (2)
-
Ivan Bezerra Allaman -
Walmes Zeviani