ddply dans une fonction

[Bigot Anthony]

Bonjour,

j'essaie de faire appel à ddply dans une fonction.

La fonction utilisée dans ddply est là pour calculer la valeur du 3è quartile (ou 75è centile pour le coup).

Code : Tout sélectionner

f<-function(df, colGrp="group", colVal="values", s=0.75){
require(plyr)

print(colVal)

res<-ddply(df, .(colGrp), transform, seuil=quantile(colVal, eval(s)))

return(res)
}
f(dat)


Pour colGrp, pas de soucis l'évaluation se fait bien.
Pour s, obligé de passer par eval(s).
Pour colVal par contre je ne comprend pas, j'ai essayé avec eval et get mais j'ai toujours la même erreur ; pourtant le print se fait bien.

"colVal" seul :

Code : Tout sélectionner

Erreur dans eval(expr, envir, enclos) : objet 'colVal' introuvable

eval(colVal) ou get(colVal)

Code : Tout sélectionner

Erreur dans eval(colVal) : objet 'colVal' introuvable
Erreur dans get(colVal) : objet 'colVal' introuvable

Si je comprend bien, ddply ne trouve pas colVal dans l'environnement (d'ailleurs si je crée un variable colVal qui prend pour valeur le nom de la colonne "values", tout fonctionne).

Je ne vois vraiment pas comment faire.

[Bigot Anthony]

Un jeu pour tester (avec mtcars):

Code : Tout sélectionner

# selection de la cylindrée et du poids
dat1=dt=mtcars[, c("cyl","wt")]

# fonction pour calculé le centile voulu du poids en fonction de la cylindée
f<-function(df, colGrp="group", colVal="values", s=0.75){
require(plyr)
print(colVal)
res<-ddply(df, .(colGrp), transform, seuil=quantile(colVal, eval(s)))
return(res)
}

# appel de la fonction
f(dat1, colGrp="cyl", colVal="wt")

# resultats
# [1] "wt"
# Erreur dans quantile(colVal, eval(s)) : objet 'colVal' introuvable
 


si on le fait directement avec ddply ; sans les guillemets ça fonctionne ; avec les guillemets, il faut

Code : Tout sélectionner

ddply(df, c("cyl")...)

au lieu de

Code : Tout sélectionner

ddply(df, .("cyl")...)

et pour "wt" pas moyen de faire fonctionner...

[Didier Plat]

Bonjour,
Une suggestion :

Code : Tout sélectionner

dat1=dt=mtcars[, c("cyl","wt")]
require(plyr)

# fonction pour calculer le centile voulu du poids en fonction de la cylindée
f<-function(df, colGrp="group", colVal="values", s=0.75){
res<-ddply(df, colGrp, .fun=function(x) quantile(x[,colVal], probs=s))
return(res)
}

# appel de la fonction
f(dat1, colGrp="cyl", colVal="wt")




[dicko ahmadou]

La fonction here du package plyr a été écrite dans ce but précis (capturer l'environnement), faudra juste rajouter un appel à get (ou eval(parse(text = ...))

Code : Tout sélectionner

f <- function(df, colGrp = "group", colVal = "values", s = 0.75) {
require(plyr)
ddply(df, colGrp, here(transform), seuil = quantile(get(colVal), s))

}