Discussion :

[audrey13930]

Bonjour,
je réalise un projet d'études météorologiques; j'ai obtenu un certain nombre de données
selon différents postes météorologiques; voici un head afin que vous visualisiez ce jeu.

Code :

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

1
2
3
4
5
6
7
8
head(h2016) 
     POSTE       DATE   TX FXI UX 
1 13001009 2016061508 20.6 3.8 57 
2 13001009 2016061509 21.1 3.2 51 
3 13001009 2016061510 22.2 3.4 47 
4 13001009 2016061511 23.4 5.3 47 
5 13001009 2016061512 24.5 5.2 48 
6 13001009 2016061513 25.0 6.9 43

J'ai en réalité 4 postes différents, comment extraires les données relatives à un poste seulement?
Voici ce que j'ai essayé, cependant cela n'a pas fonctionné: j'ai toutes les données pour tous les postes...

Code :

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

1
2
3
4
subset(h2016, POSTE==13001009)->aixh2016 # Je sélectionne uniquement les donnees pour Aix 
subset(h2016, POSTE==13004003)->arlesh2016 
subset(h2016, POSTE==13103004)->salonh2016 
subset(h2016, POSTE==13054001)->marih2016 # Pour Marignagne

[mgdondon]

Bonjour,

Votre code fonctionne correctement (j'ai modifié le poste sur 2 lignes pour tester) :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
> h2016 <- data.frame(POSTE=c(13001009,13001009,13001009,13004003,13001009,13004003),
+               DATE=c(2016061508,2016061509,2016061510,2016061511,2016061512,2016061513),
+               TX=c(20.6,21.1,22.2,23.4,24.5,25.0),
+               FXI=c(3.8,3.2,3.4,5.3,5.2,6.9),
+               UX=c(57,51,47,47,48,43))
> h2016
     POSTE       DATE   TX FXI UX
1 13001009 2016061508 20.6 3.8 57
2 13001009 2016061509 21.1 3.2 51
3 13001009 2016061510 22.2 3.4 47
4 13004003 2016061511 23.4 5.3 47
5 13001009 2016061512 24.5 5.2 48
6 13004003 2016061513 25.0 6.9 43

Code :

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

1
2
3
4
5
6
7
8
9
10
> subset(h2016, POSTE==13001009)
     POSTE       DATE   TX FXI UX
1 13001009 2016061508 20.6 3.8 57
2 13001009 2016061509 21.1 3.2 51
3 13001009 2016061510 22.2 3.4 47
5 13001009 2016061512 24.5 5.2 48
> subset(h2016, POSTE==13004003)
     POSTE       DATE   TX FXI UX
4 13004003 2016061511 23.4 5.3 47
6 13004003 2016061513 25.0 6.9 43

Cordialement,