206
Le logiciel R - Maˆ ıtriser le langage
5.4- Recr´ eez la matrice de taille 256 × 256 contenant les observations pr´ esentes
dans valeurs.
5.5- Utilisez la fonction image() sur la matrice X ainsi cr´ e´ ee. Vous utiliserez
un d´ egrad´ e de couleurs en niveaux de gris avec une centaine de teintes,
obtenu grˆ ace ` a la fonction gray().
5.6- Notez que c’est exactement ce qui est fait dans le package R AnalyzeFMRI.
Apr` es avoir t´ el´ echarg´ e les deux fichiers http://www.biostatisticien.
eu/springeR/anat.img et http://www.biostatisticien.eu/springeR/
anat.hdr, vous auriez donc pu obtenir la mˆ eme chose que pr´ ec´ edemment
(apr` es installation du package AnalyzeFMRI) en tapant :
require("AnalyzeFMRI")
Y <- f.read.volume("/chemin/vers/anat.img") # Chemin à
# remplacer.
image(X,col=gray(0:1000 / 1000))
E- Dessiner la carte d’une r´ egion fran¸ caise
Le package maps contient de nombreuses cartes de divers pays. Nous allons
l’utiliser pour repr´ esenter les contours d’un d´ epartement fran¸ cais.
5.1- Installez, puis chargez les packages maps et mapdata.
5.2- Tracez la carte de la France : map("france").
5.3- R´ ecup´ erez les donn´ ees des contours des r´ egions fran¸ caises :
france <- map("france",plot=FALSE)
5.4- Affichez le contenu de l’objet france et assurez-vous de bien comprendre
son organisation. Vous noterez par exemple que les donn´ ees de
latitude/longitude sont organis´ ees dans france$x/france$y pour chaque
d´ epartement dans france$names (jusqu’au prochain NA).
5.5- Cr´ eez un vecteur indNA contenant le vecteur des indices des valeurs manquantes.
5.6- Cr´ eez un objet contenant le nom de votre d´ epartement (par exemple
nomdept <- "Gard").
5.7- Cr´ eez un objet nomm´ e inddept qui contient l’indice du d´ epartement
nomdept dans le vecteur france$names.
5.8- Tracez la carte de votre d´ epartement.
5.9- Rajoutez le point d’une localit´ e (ou d’une adresse) sur la carte. Vous
pouvez r´ ecup´ erer les coordonn´ ees (latitude/longitude) de cette localit´ e
en vous rendant par exemple sur le site http://www.gpsvisualizer.
com/geocode.
Le logiciel R - Maˆ ıtriser le langage
5.4- Recr´ eez la matrice de taille 256 × 256 contenant les observations pr´ esentes
dans valeurs.
5.5- Utilisez la fonction image() sur la matrice X ainsi cr´ e´ ee. Vous utiliserez
un d´ egrad´ e de couleurs en niveaux de gris avec une centaine de teintes,
obtenu grˆ ace ` a la fonction gray().
5.6- Notez que c’est exactement ce qui est fait dans le package R AnalyzeFMRI.
Apr` es avoir t´ el´ echarg´ e les deux fichiers http://www.biostatisticien.
eu/springeR/anat.img et http://www.biostatisticien.eu/springeR/
anat.hdr, vous auriez donc pu obtenir la mˆ eme chose que pr´ ec´ edemment
(apr` es installation du package AnalyzeFMRI) en tapant :
require("AnalyzeFMRI")
Y <- f.read.volume("/chemin/vers/anat.img") # Chemin à
# remplacer.
image(X,col=gray(0:1000 / 1000))
E- Dessiner la carte d’une r´ egion fran¸ caise
Le package maps contient de nombreuses cartes de divers pays. Nous allons
l’utiliser pour repr´ esenter les contours d’un d´ epartement fran¸ cais.
5.1- Installez, puis chargez les packages maps et mapdata.
5.2- Tracez la carte de la France : map("france").
5.3- R´ ecup´ erez les donn´ ees des contours des r´ egions fran¸ caises :
france <- map("france",plot=FALSE)
5.4- Affichez le contenu de l’objet france et assurez-vous de bien comprendre
son organisation. Vous noterez par exemple que les donn´ ees de
latitude/longitude sont organis´ ees dans france$x/france$y pour chaque
d´ epartement dans france$names (jusqu’au prochain NA).
5.5- Cr´ eez un vecteur indNA contenant le vecteur des indices des valeurs manquantes.
5.6- Cr´ eez un objet contenant le nom de votre d´ epartement (par exemple
nomdept <- "Gard").
5.7- Cr´ eez un objet nomm´ e inddept qui contient l’indice du d´ epartement
nomdept dans le vecteur france$names.
5.8- Tracez la carte de votre d´ epartement.
5.9- Rajoutez le point d’une localit´ e (ou d’une adresse) sur la carte. Vous
pouvez r´ ecup´ erer les coordonn´ ees (latitude/longitude) de cette localit´ e
en vous rendant par exemple sur le site http://www.gpsvisualizer.
com/geocode.