Discussion :

[sam381]

Bonjour,
je commence tout juste à utiliser R, et je ne suis pas très à l'aise.

J'ai un tableau avec des noms de joueurs et leur temps et j'aimerais réaliser quelques opérations sur chacun des joueurs. Par exemple, j'aimerais calculer la médiane des temps pour chacun des joueurs et réduire l'intervalle des données que je prends (que j'ai nommé x1 et x2).
J'ai fait le code suivant mais il réalise tous les calculs pour chaque ligne et non pas pour le bloc de lignes appartenant au même joueur :

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
for(v in unique(data$joueur)){
    sel<-data$joueur==v
    extraction<-data[sel, ]
    M<-median(extraction$temps)
    MAD<-absolue(extraction$temps,M)
    x1<-M-3*MAD
    x2<-M+3*MAD
    print(sprintf("intervalle  pour procede %s = %f et %f",v,x1,x2))}
 
[1] "intervalle  pour procede sophie = 2.000000 et 2.000000" 
[2] "intervalle  pour procede sophie = -2.440000 et 6.440000"
[3] "intervalle  pour procede sophie = 2.000000 et 2.000000" 
[1] "intervalle  pour procede m‚m‚ = 3.000000 et 3.000000" 
[2] "intervalle  pour procede m‚m‚ = -1.440000 et 7.440000"
[3] "intervalle  pour procede m‚m‚ = -1.440000 et 7.440000"
[1] "intervalle  pour procede Ben = 2.000000 et 2.000000"
[1] "intervalle  pour procede Benjamin = 0.280000 et 4.720000"
[2] "intervalle  pour procede Benjamin = 0.280000 et 4.720000"
[1] "intervalle  pour procede h‚loise = 1.000000 et 1.000000"
[2] "intervalle  pour procede h‚loise = 1.000000 et 1.000000"

J'aimerais qu'il me fasse les calculs que je lui dit pour chaque "bloc" de joueur... une idée ? De plus, une fois que l'intervalle pour chaque joueur sera bien déterminé, comment sélectionner que les lignes répondant à cette intervalle ?

Merci pour vos réponses.

[Theta]

Ta boucle for fait bien les calculs "par bloc", c'est ton sprintf qui affiche plusieurs lignes car une ou plusieurs variables que tu lui donne n'est pas de longueur 1.

[sam381]

Pourtant quand je ne fais pas de sprintf et que je lui demande juste de me retourner le xa par exemple, il m'en affiche autant qu'avec le sprintf..

[Theta]

C'est surement justement cette variable qui n'est pas de longueur 1.

Enlève ta boucle for et initialise la variable v, puis execute le contenu de la boucle ligne par ligne en regardant le contenu de tes variables, le problème devrais te sauter aux yeux.

Je pense que c'est la fonction absolue (?) qui renvoie pas ce que tu pense.

[juliatheric]

Bonjour,
Comme vous l'a dit Theta, le souci se trouve certainement dans la fonction absolue.
Par ailleurs, ceci est un exemple typique où les itérations explicites devraient être évitées: la famille des fonctions *apply et apparentées (split, by, etc.) est taillée sur mesure pour pareils cas (vectorisation).

Code R :

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
eg <- list(joueur = rep(letters[seq_len(5)], 2), temps = sample.int(10))
 
données <- as.data.frame(eg, stringsAsFactors = FALSE)
 
faireTraitements <- function(x) {
    y <- median(x$temps)
    #Mettre ici les instructions de la fonction "absolue"
    x1 <- y - 3 * mad(x$temps)
    x2 <- y + 3 * mad(x$temps)
    print(sprintf("(%f, %f)", x1, x2))
    return(list(x1, x2))
}
 
résultats <- by(données, données$joueur, faireTraitements)
 
print(résultats)

Pourriez-vous poster le code de cette fonction absolue ?
Récemment, il y a eu une discussion sur comment "simplifier" les itérations.

[sam381]

Bonjour,

D'accord ! merci beaucoup !

La fonction absolue était la suivant

Code :

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

1
2
3
absolue<- function(x,y)
MAD=abs(x-y)*1.48
return(MAD)

[juliatheric]

Bonjour,
L'usage auquel vous réservez cette fonction "absolue" peut être remplacé par une méthode native en R, voir :

Code R :

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

help("mad")

En somme, il semble que vous essayez de faire quelque chose qui ressemble à ceci :

Code R :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
faireTraitements <- function(x) {
    y <- median(x$temps)
    tmp <- 3 * mad(x$temps)
    x1 <- y - tmp
    x2 <- y + tmp
    cat(x$joueur, " = (", x1, ", ", x2, ")\n", sep = "")
    return(list(x1 = x1, x2 = x2))
}
 
eg <- list(joueur = rep(letters[seq_len(5)], 2), temps = sample.int(10))
données <- as.data.frame(eg, stringsAsFactors = FALSE)
résultats <- by(données, données$joueur, faireTraitements)
print(résultats)

[sam381]

C'est parfait ça marche super bien ! Merci beaucoup !