1 Fonctions graphiques de R

Attention cependant aux excès d’esthétisme, un beau graphique n’est pas

Fonctions graphiques de R nécessairement un graphique pertinent ; éviter par exemple les abus d’effets de
perspective sans signification statistique.
Résumé
2 Commandes de base
Un rapide aperçu des très riches fonctions graphiques de R pour
gérer fenêtres et fichiers, tracer des camemberts, diagrammes en 2.1 Aide en ligne
colonnes, mosaïc plot, histogrammes, diagrammes boîtes, images,
contours, vues en 3D, ajouter des composantes graphiques. Comme pour tout utilisation avancée d’une fonction de R, le recours à l’aide
Organisation des tutoriels R. en ligne est incontournable. Apprendre à y naviguer et comprendre sa logique
• Démarrer rapidement avec R font partie de l’apprentissage.
• Initiation à R help(plot) offre un premier accès alors que
• Fonctions graphiques de R
help(par) embrouille vite les cartes.
• Programmation en R
• MapReduce pour le statisticien 2.2 Gestion de la fenêtre graphique
Les aspect statistiques sont développés dans les différents scénarios
de Wikistat. Toute commande graphique ouvre une fenêtre adaptée mais
x11() ouvre une autre fenêtre pour éviter d’écraser le graphe précédent,
1 Introduction dev.off() ferme correctement la dernière fenêtre,
Plus un problème et les données afférentes sont complexes, plus est impor- split.screen(c(1,2)) partage la fenêtre en 2 de même que l’option
tant la pertinence des graphes qui permettront d’en appréhender la structure mfrow de la commande par.
et visualiser l’information utile. Aussi, quelque soit le logiciel utilisé, com-
mercial comme SAS ou libre comme R, les versions basées sur une interface 2.3 Exportation de fichiers graphiques
amicale : SAS/Insight ou Rcmdr, sont évidemment séduisantes, car faciles à
Voici quelques possibilités pour intégrer un graphique dans un traitement de
utiliser par des menus, mais drastiquement limitées par les choix des options
texte ou générer un fichier contenant le graphique avant intégration ou archi-
nécessairement très contraints parmi des possibilités considérables.
vage.
demo(graphics) donnera un aperçu de la palette des possibilités tandis • Sous Windows, clic droit dans la fenêtre pour copier le graphique dans le
qu’une visite à l’un de ces sites : gallerie 1, gallerie 2, gallerie 3, ouvre de presse papier avant de le coller dans un texte,
vastes perspectives sur les possibilités offertes. • L’obtention du format standard jpeg nécessite d’enchaîner les com-
C’est la principale raison, qui motive de cheminer par l’apprentissage du mandes :
lange de commande R plutôt que par celui d’une interface graphique aussi jpeg("fichier.jpeg") # ou bmp(), png(), pdf()
amicale soit-elle. Bien entendu ce tutoriel ne présente que quelques facettes # le graphique n’apparaît pas à l’écran :
des graphiques les plus utilisés tandis que beaucoup de ressources et sites, plot(1:100); text(20,80,"abcdef")
comme ceux ci-dessus, rendent accessibles des exemples complexes et ciblés dev.off() # fermeture indispensable du fichier
sur des besoins particuliers.
 2 Fonctions graphiques de R

qui est alors créé dans le répertoire courant. Réponses

• Utiliser le menu Fichier, rubrique Sauver sous de la fenêtre gra- 1. Les 2 commandes découpent la fenêtre graphique en 4 cellules. La dif-
phique. férence entre les 2 commandes est illustrée par le tableau ci-dessous
(par(mfrow=c(2,2)) - par(mfcol=c(2,2))) :
3 Description élémentaire 1-1 2-3
3-2 4-4
3.1 Variables qualitatives Avec par(mfrow=c(2,2)), les 4 graphiques à venir sont intégrés
ligne par ligne ; avec par(mfcol=c(2,2)), ils le sont colonne par
Fonctions pie, barplot, mosaicplot. colonne.
vec=c(12,10,7,13,26,16,4,12) 2. Par exemple, pour modifier les couleurs du camembert :
pie(vec);pie(vec,clockwise=T) pie(vect,col=rgb(0,0,seq(0,1,l=8)))
names(vec)=LETTERS[1:8]
# de très nombreuses options sont disponibles pie(vect,col=rainbow(8))
pie(vec) ...
barplot(vec) 3. Commentaire de la ligne de commande :
# pour découper la fenêtre en 1 ligne 2 colonnes • n=200 # La valeur 200 est affectée à n.
par(mfrow=c(1,2)) • rep(1,n) # Crée le vecteur de taille n ne contenant que des 1.
pie(vec) • pie(rep(1,n) # Produit un camembert pour n variables d’effec-
barplot(vec) tifs tous égaux à 1. Les secteurs sont donc de même angle.
par(mfrow=c(1,1)) • rainbow(n) # Renvoie un vecteur de longueur n contenant un dé-
mosaicplot(Titanic, color = T) gradé des couleurs de l’arc-en-ciel.
colors() • border=NA # Évite de tracer les bords des secteurs.
Questions • labels="" # Aucune étiquette n’identifie les secteurs.
1. Quelle est la différence entre par(mfrow=c(2,2)) et 3.2 Variables quantitatives
par(mfcol=c(2,2)) ?
2. Consulter l’aide en ligne pour tester certains arguments optionnels des Commandes boxplot, hist, plot, pairs.
fonctions pie() et barplot(). x=rnorm(50)
3. Commenter les commandes suivantes (extraites de l’aide en ligne de boxplot(x)
pie()) : hist(x)
n=200 stripchart(x)
pie(rep(1,n), labels="", col=rainbow(n), border=NA) Question
1. Représenter dans la même fenêtre graphique, le "stripchart", le
diagramme-boîte horizontal et l’histogramme correspondant l’un sous
l’autre.
 3 Fonctions graphiques de R

Réponse # images , nappes et contours

image(M)
par(mfrow=c(3,1))
x = seq(-10, 10, length= 30);y=x
stripchart(x)
f = function(x,y){r=sqrt(x^2+y^2); 10 * sin(r)/r}
boxplot(x,horizontal=T)
z = outer(x, y, f)
hist(x)
z[is.na(z)] = 1
Croisement de variables. persp(x,y,z)
persp(x, y, z, theta=30, phi=30,expand = 0.5,
data(iris) col="lightblue")
pairs(iris[,1:4]) image(x,y,z)
par(mfrow=c(2,2)) contour(x,y,z)
plot(iris[,1],iris[,2],xlab="Sepal Length", filled.contour(x,y,z)
ylab="Sepal Width",main="Iris data", image(x,y,z)
col="red",type="l") contour(x,y,z,add=T)
points(iris[,1],iris[,2],col="green",pch=21)
boxplot(iris[,1:4]) Pour mieux comprendre la fonction outer(), tester :
hist(iris[,1],xlab="Sepal Length",
x=y=1:5
main="Histogramme")
z=outer(x,y,"+");z
Ajouter des éléments graphiques.
x=seq(-10,10,l=50)
plot(x,sin(x))
plot(x,sin(x),type="l")
abline(v=0,col="blue",lwd=5,lty=3)
abline(h=sin(0.7),col=3)
text(-5,-0.5,"texte",font=3)

Gestion des libellés.

par(mfrow=c(1,2))
plot(x,sin(x),type="l",col=1,main="sinus")
plot(x,cos(x),type="b",col=3,xlab="Abscisses")

3.3 Vers la 3D
Commandes image, persp, contour.
M=matrix(1:100,nc=10)