fonction aggregate et drop

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:

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:

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:

1 2 3 4 groupe valeurSPAD 121 31.6 122 45.0 123 46.9
Merci d'avance pour votre aide,

Louison :)

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

1 2 nb_mes_groupe <- tapply(racine[, 'groupe'], racine[, 'groupe'], length)
Ensuite il faut récupérer les indices dans ta base
Code:

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.

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

Code:

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.