Problème de représentation graphique

Bonjour

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

Code:

---

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:

---

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:

---

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:

---

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

---

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

Cdt

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:

---

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

Bonjour,

Citation:

---

Envoyé par Hydro999

Code:

---

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 !

Ok merci pour vos réponses.

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

Code:

---

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:

---

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

---

en pensant que par défaut

Code:

---

ozone[selection,]

---

R allait utiliser la colonne d'index.

cdt

Citation:

---

Envoyé par Hydro999

...en pensant que par défaut

Code:

---

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:

---

row.names(ozone)

---

ou bien

Code:

---

as.integer(row.names(ozone))

---

Ok merci pour les indications

Cdt