[R-br] Legenda com gradientes-valores em superfície resposta!

[Ivan Bezerra Allaman]

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 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/20110608/bf6eb086/attachment.html>