Discussion :

[dachboden]

Bonjour,

Je souhaiterais savoir s'il existe une panel.function qui calcule des moyennes pondérées pour des trellis graphs (package lattice)
Avec le code suivant je superpose une moyenne des y pour chaque valeur de x au scatterplot de (x,y) mais ce n'est pas une moyenne pondérée.

Code :

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

1
2
3
4
5
6
 
xyplot(y~x|cat,data=mytab,groups=mygroup,layout=c(1,nbr.panel), 
       panel = function(x, y) {
           panel.xyplot(x, y, pch="*")
           panel.average(x,y,fun=mean,horizontal=FALSE,col.line="red")
})

Je souhaiterais pouvoir faire quelque chose comme ça

Code :

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

1
2
3
4
5
6
 
xyplot(y~x|cat,data=mytab,groups=mygroup,layout=c(1,nbr.panel), 
       panel = function(x, y) {
           panel.xyplot(x, y, pch="*")
           panel.average(x,y,weight,fun=mean_ponderee,horizontal=FALSE,col.line="red")
})

Merci bcp !

[dev_ggy]

Si j'ai bien compris le fonctionnement de xyplot, tu crées la fonction

Code :

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

1
2
 
mean_pond<-function(x){mean(x*weight)/mean(weight)}

dans la fonction xyplot, tu prends les paramètres notés ci-dessous.

Code :

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

1
2
3
4
5
6
xyplot(y~x|cat,data=mytab,groups=mygroup,layout=c(1,nbr.panel), 
       panel = function(x, y) {
           panel.xyplot(x, y, pch="*")
           panel.average(x,y,weight=c,fun=mean_pond,horizontal=FALSE,col.line="red")
})

Avec c ton vecteur poids.

et cela devrait fonctionner.

[dachboden]

Merci pour ton aide.

C'est une bonne idée d'essayer de définir sa propre fonction le problème c'est dans le passage de l'argument weight
Que ce soit au niveau de la fonction panel.average elle-même comme tu le
suggères ou bien dans les arguments de la fonction xyplot, je n'arrive pas à faire reconnaitre le paramètre weight à l’exécution.

Note qu'une solution alternative consiste à ajouter dans mon data.frame les moyennes pondérées déjà calculées et à les flagger avec une nouvelle valeur pour la variable mygroup.
Ensuite il ne me reste plus qu'a faire un

Code :

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

1
2
 
xyplot(y~x|cat,data=mytab,groups=mygroup,layout=c(1,nbr.panel))

sans même mentionner de fonction car par default

Code :

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

1
2
 
panel = panel.xyplot(x, y)

C'est cette solution que je voulais éviter en posant ma question.

J'aime ta suggestion mais je bute sur le passage d'argument !

[dev_ggy]

Si cela ne marche pas, tu peux réécrire la fonction panel.average

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
18
19
20
21
22
23
24
25
26
27
28
panel.average_weight<-function (x, y, fun = mean, horizontal = TRUE, lwd = reference.line$lwd, 
    lty = reference.line$lty, col, col.line = reference.line$col, 
    type = "l", ...) 
{
    x <- as.numeric(x)
    y <- as.numeric(y)
    reference.line = trellis.par.get("reference.line")
    if (!missing(col)) {
        if (missing(col.line)) 
            col.line <- col
    }
    if (horizontal) {
        vals <- unique(sort(y))
        yy <- seq_along(vals)
        xx <- numeric(length(yy))
        for (i in yy) xx[i] <- fun(x[y == vals[i]], weight)
        panel.lines(xx, vals[yy], col = col.line, lty = lty, 
            lwd = lwd, ...)
    }
    else {
        vals <- unique(sort(x))
        xx <- seq_along(vals)
        yy <- numeric(length(xx))
        for (i in xx) yy[i] <- fun(y[x == vals[i]], weight)
        panel.lines(vals[xx], yy, col = col.line, lty = lty, 
            lwd = lwd, ...)
    }
}

et là ta fonction panel.average_weight devrait fonctionner, toujours avec la fonction mean_pond.

C'est l'idée alternative que j'avais.

[dachboden]

Merci,

Je ne suis pas encore suffisamment expert pour écrire des fonctions, même si pour commencer je peux essayer de modifier à la marge le code de fonctions existantes.

J'ai essayé ta solution qui devrait marcher modulo qq ajustements. Pour le moment ça bugue encore.
Je vais creuser. Je donnerai un update quand ça marchera.