De la couleur sur un graph 3D pour indiquer valeur de Z

Bonjour,

Voici un script 'hello world' pour afficher un graphique en 3D :

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

x <- seq(1, 200, by=2) y <- x z<-array(NA,c(length(x),length(y)))   for (i in 1:length(x)) for (j in 1:length(y)) z[i,j] <- i*j - sin(i+j)*(i+j)   contour(x, y, z) x11() persp(x,y,z,theta = -30, phi = 40)

---

Savez-vous s'il est possible - avec les fonctions contour et persp - de mettre des couleurs pour indiquer la valeur de "z"? (Un truc comme "plus c'est rouge, plus z est grand")

Merci d'avance :)

Ok, j'ai trouvé la solution ici .

Pour information, voici l'exemple donné sur le site:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17

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)

---