Densité 2D

Si on a 2 variables avec un très grand nombre de points, difficile de voir quelque chose sur un plot 2D. Alternative est de représenter une densité :

x <- runif(500, 0, 2)
y <- rnorm(500, 2, 1)
image(kde2d(x, y, n = 100, lims = c(-1, 3, -1, 5)))

• n : nombre de mailles dans chaque direction (horizontale ou verticale), 25 par défaut.
• lims : xmin, xmax, ymin, ymax.

Traçage de la densité en 3D :

x <- runif(500, 0, 2)
y <- rnorm(500, 2, 1)
persp(kde2d(x, y, n = 20, lims = c(-1, 3, -1, 5)),
      phi = 30, theta = 20, d = 5, xlab = "x", ylab = "y", zlab = "z val",
      xlim = c(-1, 3), ylim = c(-1, 5), col = "skyblue",
      main = "in 3 d", ticktype = "detailed")

• theta : azimuth de visualisation.
• phi : colatitude de visualisation.
• d : valeur qui accentue plus ou moins l'effet de la perspective (plus c'est petit, plus l'effet de perspective est fort).
• ticktype : "simple" pour absence de graduation (défaut), "detailed" pour graduation comme en 2D.
• main : titre.
• xlab, ylab, zlab : les étiquettes des axes.
• xlim, ylim, zlim : les valeurs extrèmes selon les axes.