Lise Vaudor
21 mars 2017
On peut tracer un nuage de points à l’aide de la fonction plot().
Pour tracer les valeurs d’une variable x (Temp):
plot(Temp)
Pour tracer une variable y (Ozone) en fonction d’une variable x (Temp)
plot(Temp,Ozone)
La fonction plot() a de nombreux arguments, permettant entre autres de modifier
plot(Temp,
Ozone,
col="deepskyblue",
pch=20,
main="Ozone en fonction de Température",
xlab="Température(°C)",
ylab="Ozone(ppm)")
La fonction plot() a de nombreux arguments, permettant entre autres de modifier
plot(Temp,
type="h",
col="hotpink",
lwd=2,
main="Evolution de la température",
xlab="Mesures",
ylab="Température")
La fonction plot() a de nombreux arguments, permettant entre autres de modifier
plot(Temp,
type="l",
col="purple",
lwd=3,
lty=3,
xlab="Mesures",
ylab="Conc. en ozone (ppm)")
On peut superposer plusieurs graphiques, par exemple en appelant d’abord une fonction graphique “de base” puis en appelant une des fonctions suivantes:
points()
pour superposer un nuage de pointsabline()
pour rajouter une ligne, de pente “a” et d’ordonnée à l’origine “b” ou une ligne verticale d’abscisse “v” ou une ligne horizontale d’ordonnée “h”text()
pour rajouter du texteplot(Ozone,Temp)
abline(h=mean(Temp),
col="hotpink")
text(150,
80,
labels="y=temp. moyenne",
col="hotpink")
Les boîtes à moustaches (ou boxplots, en anglais) sont un moyen extrêmement classique de représenter les distributions.
La “boîte” est délimitée par les 1er et 3ème quartiles. Une barre horizontale marque la médiane.
Soit \(H\) la hauteur de la boîte. Les moustaches s’étendent:
On peut modifier la valeur de \(\tau\) (argument range) de la fonction boxplot():
layout(matrix(1:3,nrow=1))
boxplot(Ozone, col="pink")
# pour que les moustaches soient plus courtes:
boxplot(Ozone, col="pink", range=0.5)
# pour que les moustaches s'étendent jusqu'aux extrêmes:
boxplot(Ozone, col="pink", range=0)
Le plus souvent, les boîtes à moustaches sont utilisées pour représenter simultanément les distributions des observations pour plusieurs groupes, par exemple:
boxplot(Temp~gWind, col="pink")
resume <- summarise(group_by(air,gWind),mTemp=mean(Temp))
resume
## # A tibble: 5 × 2
## gWind mTemp
## <fctr> <dbl>
## 1 (1.7,7.4] 83.85714
## 2 (7.4,9.7] 77.55263
## 3 (9.7,11.5] 77.41176
## 4 (11.5,20.7] 72.07895
## 5 NA 76.00000
boxplot(Temp~gWind, col="pink")
points(resume$gWind,resume$mTemp,
cex=2,pch=20,col="purple")
Les histogrammes offrent un mode de représentation des distributions plus précis que les boîtes à moustaches.
hist(Ozone, col="light blue")
L’argument breaks permet de spécifier comment former les classes pour le tracé de l’histogramme:
layout(matrix(1:2,nrow=2))
hist(Ozone, col="khaki1", breaks=seq(0,200,25))
hist(Ozone, col="pink", breaks=20)