Discussion :

[hoccha]

Bonjour,

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[hoccha]

Ok, j'ai trouvé la solution ici .

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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)