Legenda com gradientes-valores em superfície resposta!

Boa tarde senhores! Gostaria de saber se é possível criar uma legenda com gradiente-valores em um gráfico de superfície resposta? Vejamos o seguinte CMR. par(bg = "white") x <- seq(-1.95, 1.95, length = 30) y <- seq(-1.95, 1.95, length = 35) z <- outer(x, y, function(a,b) a*b^2) nrz <- nrow(z) ncz <- ncol(z) # Create a function interpolating colors in the range of specified colors jet.colors <- colorRampPalette( c("blue", "green") ) # Generate the desired number of colors from this palette nbcol <- 100 color <- jet.colors(nbcol) # Compute the z-value at the facet centres zfacet <- z[-1, -1] + z[-1, -ncz] + z[-nrz, -1] + z[-nrz, -ncz] # Recode facet z-values into color indices facetcol <- cut(zfacet, nbcol) persp(x, y, z, col=color[facetcol], phi=30, theta=-30) Aqui é criado um gráfico de superfície resposta, com cores que vão do azul ao verde. Vejamos outro CMR, para chegarmos onde eu quero. library(fields) loc <- make.surface.grid(list( seq(1,dim(z)[1],length=500), seq(1,dim(z)[2],length=500) )) r <- interp.surface( list(x=1:dim(z)[1], y=1:dim(z)[2], z=z), loc ) op <- par(mfrow=c(2,1)) image.plot(z, main="Raw data") image.plot(as.surface(loc,r), main="Linear interpolation") par(op) Percebam que aqui, é criado um gráfico com uma legenda indicando cores que correspondem a determinados valores. O que eu quero é justamente esta legenda no gráfico do primeiro CMR. É possível? Desde já agradeço 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, Usa o lattice::wireframe() ao invés da persp() que a legenda em escala de cores sai de graça. À 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