Discussion :

[louisonl]

Bonjour,

J'ai un problème avec un jeu de données.

J'ai voulu faire des médianes sur les différents groupes de mon jeu de données. Le résultat est donc logiquement un jeu de données avec moins de lignes, puisqu'elles ont été regroupées selon leur groupe.
J'ai donc écrit la ligne de commande suivante:

Code :

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

racine2<-aggregate(racine[ ,"valeurSPAD",drop=FALSE],by=racine[,"groupe",drop=FALSE],median)

Ceci me donne bien ce que je veux mais enlève les autres colonnes de mon jeu de données, or j'aurais besoin de garder les informations données dans les autres colonnes pour la suite de mes calculs avec conditions.
Voici un extrait de ma base de données initiale :

Code :

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

1
2
3
4
5
6
7
8
9
10
   facteurherbicide           type facteurespece facteurbloc groupe valeurSPAD
31            traite    gesse/colza         colza           1    121       31.6
32            traite    gesse/colza         colza           1    121       32.7
33            traite    gesse/colza         colza           1    121       30.3
34         nontraite         temoin         gesse           1    122       40.6
35         nontraite         temoin         gesse           1    122       45.0
36         nontraite         temoin         gesse           1    122       50.4
37            traite feverole/colza      feverole           1    123       46.9
38            traite feverole/colza      feverole           1    123       49.4
39            traite feverole/colza      feverole           1    123       42.8

et voici ce que j'obtiens suite à ma ligne de commande:

Code :

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

1
2
3
4
 groupe valeurSPAD
    121       31.6
    122       45.0
    123       46.9

Merci d'avance pour votre aide,

Louison

[mamounMob]

D'abord il faut récupérer le nombre de mesures par groupe :

Code :

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

1
2
 
nb_mes_groupe <- tapply(racine[, 'groupe'], racine[, 'groupe'], length)

Ensuite il faut récupérer les indices dans ta base

Code :

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

1
2
 
ind_groupe <- cumsum(nb_mes_groupe)

Enfin il te suffit juste concaténer 'ind_groupe' avec 'racine2', le type et le facteur ne changent pas en fonction du groupe donc ça marche.
Rq : Tu peux décaler 'ind_groupe' par '-2' car il y a 3 observations du groupe 121.

[louisonl]

Bonjour,

Merci beaucoup pour la rapidité de votre réponse. Je suis débutante sur R et c'est peut être pour cela que j'ai mal compris votre message.

Cependant quand je reprend vos lignes de commandes, mon tableau final ne contient toujours pas les colonnes "facteur herbicide", "type","facteur espece" et "facteur bloc" :/

Merci d'avance,

Louison

[mamounMob]

Code :

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

1
2
 
cbind(racine[ind_groupe, ], racine2)

C'est la bonne commande pour retrouver le tableau que tu veux.
Rq : 'ind_groupe' contient les indices du tableau racine.