Extraire des valeurs d'un fichier, help me !

Bonjour tout le monde. :)

Je viens de m'inscrire sur ce forum ! Je suis thésard en hydrologie à Toulouse, et j'ai décidé de me mettre à R pour traiter des données. Je suis encore bien novice, et je risque sans doute de poser des questions bêtes parfois... J'espère dans tous les cas que certains d'entre vous pourront répondre à mes interrogations. :)

Ma question est donc la suivante :
J'aimerais extraire des valeurs d'un fichier X constitué de 365 valeurs. Il faudrait dans un premier temps que je divise ce fichier en paquets de 5 pour obtenir un fichier de 73 paquets de 5.

------
V1
V2
V3
V4
V5
--------
V6
V7
V8
V9
V10
--------
etc...

Ensuite j'aimerais sélectionner la plus grande valeur de chaque paquet. Je voudrais donc obtenir au final un nouveau fichier Y constituée de 73 valeurs (étant les valeurs les plus élévées de chaque paquet).

Comment formuler ça sous R... ? :roll:

Merci d'avance pour votre aide !!

Vincent

Bonjour et bienvenu sur le forum.

Quand vous dites que votre fichier X est constitué de 365 valeurs, ce sont bien 365 lignes et une seule colonne qui constituent ce fichier ?
Peut-on voir un exemple de fichier à traiter ?

Bonjour et merci d'avoir lu mon message !

Ce sont bien 365 lignes oui, et une colonne.
Sur R, le fichier se présente sous forme de deux colonnes : une qui représente le numéro de ligne (de 1 à 365), l'autre la valeur attribuée.

Code:

Code:

---

1 2 3

> X=read.table("/home/vincent/modcou/pluies_obs") > attach(X) > X

---

V1
1 51.27
2 45.34
3 40.43
4 36.49
5 33.43
6 31.13
7 29.41
etc...

Pour obtenir une valeur max d'un paquet de 5, je pensais utiliser la fonction "max", mais ce que j'aimerais pouvoir faire c'est une boucle sur les 73 paquets de 5 en fait.

Citation:

---

Envoyé par ced

Bonjour et bienvenu sur le forum.

Quand vous dites que votre fichier X est constitué de 365 valeurs, ce sont bien 365 lignes et une seule colonne qui constituent ce fichier ?
Peut-on voir un exemple de fichier à traiter ?

---

Pas besoin de boucles pour tout ça...
Personnellement, pour ce genre de choses, j'ai tendance à utiliser le package data.table.

Voici comment je ferais :

Code:

---

1 2 3 4 5 6 7 8 9 10

X <- data.frame(V1=rnorm(365)) X$gp <- as.integer((1:365-0.1)/5)   library(data.table) X2 <- data.table(X, key="V1")   X3 <- X2[, max(V1), by=gp] setkeyv(X3, "gp")   X3

---

Je crée un data frame de 365 valeurs aléatoires.
Pour obtenir des groupes de 5 lignes, je rajoute une colonne à mon data frame dans laquelle je divise le numéro de la ligne (le - 0.1 permet que le premier groupe ait bien 5 lignes) par 5, pour numéroter mes groupes de 0 à 72.

Ensuite, je transforme mon data frame en data.table (variable X2).
Puis je sélectionne, dans ce data.table, le max de chaque groupe de 5 lignes, et je mets le tout dans une variable X3.
En définissant ma colonne gp (groupe) comme clé de ce data.table, ça le trie automatiquement selon gp.
Ne reste plus qu'à afficher X3 pour vérifier qu'il est correct, et le tour est joué.

On peut faire aussi sans passer par data.table, mais ce package est super puissant pour ce genre de choses.

Au passage, si vraiment vous souhaitez récupérer un dataframe final avec une seule colonne, comme celui de départ, alors il n'y a qu'à faire :

Code:

---

as.data.frame(X3)[, "V1", drop = FALSE]

---

Je n'ai qu'une réponse à apporter ...

Balaise !!!

Votre technique marche très bien, je n'avais pas encore installé le package data.table, bonne chose de faite ! :)

ça fait franchement plaisir d'être accueilli de la sorte sur ce forum. Non seulement vous m'avez donné le code, mais en plus vous m'avez expliqué comment vous avez procédé, c'est vraiment sympa !! :)
R est un outil que j'utilise souvent dans la thèse que je débute, et vous ne pouvez pas vous imaginer comme votre aide à cette question m'a été précieuse ! Il est je trouve parfois difficile de trouver de la documentation pour des questions spécifiques de la sorte, sachant en plus que l'informatique n'est pas mon fort. Enfin, j'espère progresser grâce à ce forum et à des gens comme vous. :)
A très vite sur le forum donc !

Bien amicalement,
Vincent