Discussion :

[Hydro999]

Bonjour

Je regarde les exemples du livre : Statistiques avec R
Voici les données que j'utilise:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
          T12 maxO3  vent pluie    Vx12
20010601 18.5    87  Nord   Sec -1.7101
20010602 18.4    82  Nord   Sec -4.0000
20010603 17.6    92   Est   Sec  1.8794
20010604 19.7   114  Nord   Sec  0.3473
20010605 20.5    94 Ouest   Sec -2.9544
20010606 19.8    80 Ouest Pluie -5.0000
20010607 15.6    79 Ouest   Sec -1.8794
20010610 17.5    79  Nord   Sec -1.0419
20010611 19.6   101  Nord   Sec -1.0261
20010612 21.9   106 Ouest   Sec -2.2981

pour utiliser un symbole spécial marquant les jours de pluie on utilise le code suivant:

Code :

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

1
2
3
4
plot(maxO3~T12,data=ozone,type="p",pch=3,cex=.75)
text(ozone[,"T12"],ozone[,"maxO3"],substr(rownames(ozone),5,8),cex=.75,pos=3,offset=.3)
selection <- ozone[,"pluie"]=="Pluie"
points(ozone[selection,"T12"],ozone[selection,"maxO3"],pch=21,bg="grey70",cex=.75)

maintenant je voudrai faire la même chose mais en abscisse au lieu d'utiliser "T12" je souhaite uniquement utiliser l'index c'est à dire les dates
voici ce que j'écris

Code :

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

1
2
> plot(ozone[,"maxO3"],type="l")
> points(ozone[selection,],ozone[selection,"maxO3"],pch=21,bg="grey70",cex=.75)

Et j'ai l'erreur suivante

Code :

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

Erreur dans xy.coords(x, y) : 'x' and 'y' lengths differ

Comment faire pour marquer sur le graph les jours de pluie.

Cdt

[manoir]

Bonjour Hydro999,

Si j'ai bien compris, tu veux afficher les dates sur la courbe et aussi distinguer le jour de pluie: je te suggère de tourner ce code

Code :

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

1
2
3
4
5
6
 
toto<-cbind(ozone,indic=c(1:10))  # ajouter une colonne indice
plot(toto[,"maxO3"],type="l")
text(toto[,"indic"],toto[,"maxO3"],substr(rownames(toto),5,8),cex=.75,pos=3,offset=.3)
selection <- ozone[,"pluie"]=="Pluie"
points(toto[selection,"indic"],toto[selection,"maxO3"],pch=21,bg="grey70",cex=.75)

J'espère que ce code répondrait à ta question

Manoir

[pitipoisson]

Bonjour,

Code :

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

1
2
> plot(ozone[,"maxO3"],type="l")
> points(ozone[selection,],ozone[selection,"maxO3"],pch=21,bg="grey70",cex=.75)

En tout cas, normal que tu ais un problème : pour les y, tu sélectionnes une colonne, mais pas pour les x !

[Hydro999]

Ok merci pour vos réponses.

Je voulais éviter de créer une autre colonne et utiliser directement la colonne d'indice

Code :

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

1
2
3
4
5
6
7
8
9
10
20010601 
20010602 
20010603 
20010604 
20010605 
20010606 
20010607 
20010610 
20010611 
20010612

Mais ça ne semble pas possible.
C'est pour cela que j'ai écrit cette ligne de code

Code :

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

 points(ozone[selection,],ozone[selection,"maxO3"],pch=21,bg="grey70",cex=.75)

en pensant que par défaut

Code :

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

ozone[selection,]

R allait utiliser la colonne d'index.

cdt

[pitipoisson]

...en pensant que par défaut

Code :

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

ozone[selection,]

R allait utiliser la colonne d'index.

Ah bin non !
Selon que tu veuilles l'utiliser comme des noms (dans l'ordre où ils apparaissent... faut que la méthode plot existe pour ce type de donnée) ou comme entier, ce sera respectivement

Code :

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

row.names(ozone)

ou bien

Code :

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

as.integer(row.names(ozone))

[Hydro999]

Ok merci pour les indications

Cdt