Graphe a 3 variables sur R?

Version imprimable

01/07/2015, 21h41
Posthume

2 pièce(s) jointe(s)

Graphe a 3 variables sur R?

Bonjour,

J'ai un large fichier de données avec 3 variables par donnée:
x = position
y= longueur
z = nombre

Voila un exemple:
x y z
12 40 8
13 40 8
14 40 8
15 40 5
16 24 9

J'aimerais representer un genre de heatmap ou de surface 3D avec pour abscisse x, pour ordonnée y et chaque point(x, y) avec une couleur ou une hauteur proportionnelle a la valeur z correspondante.
J'ai essayé levelplot mais on dirait qu'il ne prend pas en considération la valeur z de mon fichier.
J'ai aussi essayé persp3D et filled.contour mais a chaque fois j'ai cette erreur :
Error in contour.default(x, y, z) : increasing 'x' and 'y' values expected

Est-ce qu'avec le genre de données que j'ai, il est possible de créer une surface 3D ou un heatmap comme ci-dessous?

Pièce jointe 182040

Pièce jointe 182041

Je vous remercie d'avance pour vos conseils.

1 pièce(s) jointe(s)

J'ai avancé et j'ai réussi a produire un wireframe. Cela se rapproche de ce que je veux.
Code:

1 2 3 4 5 wireframe(z~x*y, data=dfr, scales=list(arrows=FALSE), zlab=list(zlab, rot=90), xlab=list(xlab, rot=24), ylab=list(ylab, rot=-37), col.regions=colr(100), drape=TRUE, screen=list(z=40, x=-70))
Pièce jointe 182048

Est-ce que quelqu'un saurait comment lisser les courbes et se débarrasser des arretes du cube a part celles correspondant au x, y, z?

3 pièce(s) jointe(s)

Bonjour,
Oui c'est possible de supprimer les cadres : utiliser les paramètres de lattice. Quant au lissage, il est faisable directement dans les arguments de wireframe ; l'interpolation améliorera encore plus le résultat du lissage (par exemple, voir la bibliothèque akima).

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
invisible(sapply(
    c("akima", "fields", "lattice"),
    function(xxx) {
        library(xxx, character.only = TRUE)
    }
))
 
interpolate <- function(component) {
    return(seq(
        from = min(component),
        to = max(component),
        length.out = 2 * length(component)
    ))
}
 
#smp1 : jeu de données ; liste à 3 composantes (x, y, z)
 
smp1Smooth <- with(smp1, interp(x = x, y = y, z = z ,
                                xo = interpolate(x), yo = interpolate(y))
)
 
image(smp1Smooth)  #bibliothèque "base"
 
image.plot(smp1Smooth) #bibliothèque "fields"
 
wireframe(z ~ x * y, data = smp1, shade = TRUE, drape = TRUE, colorkey = TRUE,
          pretty = TRUE, par.settings = list(
              box.3d = list(col=NA),
              axis.line = list(col = "transparent")
          )
)

Pièce jointe 182603 Pièce jointe 182604 Pièce jointe 182605