Chapitre 5. Techniques pour tracer des courbes et des graphiques
205
1 4
1
1
1
2
2
4
4
1 5
1
1
1
1
1
1
1
1 6
1
1
1
1
1
1
1
Nous allons proposer une repr´ esentation graphique int´ eressante pour ce type
de donn´ ees.
5.1- Pour la semaine W7, calculez le vecteur ( f 1 , 1− f 1 , f 2 , 1− f 2 , f 3 , 1− f 3 , f 4 , 1− f 4 )
o` u f i est la fr´ equence de la modalit´ e i (1 ≤ i ≤ 4) observ´ ee pendant la
semaine W7 sur les seize sujets. (indice : utilisez la fonction tabulate(),
cbind(), t() et as.vector()).
5.2- Maintenant, utilisez la fonction apply() pour faire le mˆ eme calcul sur
toutes les semaines. Le r´ esultat sera stock´ e dans une matrice.
5.3- Utilisez la fonction barplot() et le param` etre col=c("black","white")
sur cette matrice. Le graphique que vous obtenez permet de rapidement se faire une id´ ee de l’´ evolution de la distribution de la variable
Sensation de brulure au cours du temps.
5.4- Modifiez le graphique pr´ ec´ edent afin de changer la couleur des barres
repr´ esentant les fr´ equences en rouge. Les num´ eros des semaines devront
ˆ etre ´ ecrits en bleu, en haut du graphique et non plus en bas. Les num´ eros
des modalit´ es devront ˆ etre ´ ecrits `
a gauche, en bleu. Un titre sera ajout´ e
au graphique.
D- Affichage d’images anatomiques du cerveau
Les donn´ ees acquises au cours d’un examen d’imagerie par r´ esonance magn´ etique (IRM) du cerveau humain sont en g´ en´ eral stock´ ees dans un fichier
binaire d’extension * .img. Nous allons voir comment lire et afficher ce type de
donn´ ees.
5.1- Importez le fichier http://www.biostatisticien.eu/springeR/anat.
img, qui contient l’image d’une seule coupe c´ er´ ebrale de 256 × 256 pixels,
au moyen de la fonction readBin(). Ces donn´ ees peuvent ˆ etre consid´ er´ ees
comme une suite de 256 × 256 paires d’octets (raw). Vous stockerez ces
donn´ ees dans un objet nomm´ e octets.
5.2- Lors de l’´ ecriture de ces donn´ ees, chaque paire d’octets a en fait ´ et´ e ´ ecrite
en sens inverse (c’est-` a-dire par exemple la paire d’octets 02 56 a ´ et´ e
´ ecrite 56 02). Vous devez donc permuter toutes les paires d’octets deux
` a deux. Vous stockerez le r´ esultat de cette op´ eration dans x.
5.3- Il faut maintenant transformer cette suite de paires d’octets en valeurs
num´ eriques que l’on pourra repr´ esenter graphiquement. Lorsque l’on a
deux octets (par exemple 02 56), il faut utiliser l’instruction
as.numeric("0x0256") pour obtenir la valeur d´ ecimale correspondant ` a
cette paire d’octets (dans ce cas on obtient 598). Transformez x en valeurs
d´ ecimales et stockez le r´ esultat dans un objet nomm´ e valeurs (indice :
utilisez les fonctions matrix(), apply() et paste()).
205
1 4
1
1
1
2
2
4
4
1 5
1
1
1
1
1
1
1
1 6
1
1
1
1
1
1
1
Nous allons proposer une repr´ esentation graphique int´ eressante pour ce type
de donn´ ees.
5.1- Pour la semaine W7, calculez le vecteur ( f 1 , 1− f 1 , f 2 , 1− f 2 , f 3 , 1− f 3 , f 4 , 1− f 4 )
o` u f i est la fr´ equence de la modalit´ e i (1 ≤ i ≤ 4) observ´ ee pendant la
semaine W7 sur les seize sujets. (indice : utilisez la fonction tabulate(),
cbind(), t() et as.vector()).
5.2- Maintenant, utilisez la fonction apply() pour faire le mˆ eme calcul sur
toutes les semaines. Le r´ esultat sera stock´ e dans une matrice.
5.3- Utilisez la fonction barplot() et le param` etre col=c("black","white")
sur cette matrice. Le graphique que vous obtenez permet de rapidement se faire une id´ ee de l’´ evolution de la distribution de la variable
Sensation de brulure au cours du temps.
5.4- Modifiez le graphique pr´ ec´ edent afin de changer la couleur des barres
repr´ esentant les fr´ equences en rouge. Les num´ eros des semaines devront
ˆ etre ´ ecrits en bleu, en haut du graphique et non plus en bas. Les num´ eros
des modalit´ es devront ˆ etre ´ ecrits `
a gauche, en bleu. Un titre sera ajout´ e
au graphique.
D- Affichage d’images anatomiques du cerveau
Les donn´ ees acquises au cours d’un examen d’imagerie par r´ esonance magn´ etique (IRM) du cerveau humain sont en g´ en´ eral stock´ ees dans un fichier
binaire d’extension * .img. Nous allons voir comment lire et afficher ce type de
donn´ ees.
5.1- Importez le fichier http://www.biostatisticien.eu/springeR/anat.
img, qui contient l’image d’une seule coupe c´ er´ ebrale de 256 × 256 pixels,
au moyen de la fonction readBin(). Ces donn´ ees peuvent ˆ etre consid´ er´ ees
comme une suite de 256 × 256 paires d’octets (raw). Vous stockerez ces
donn´ ees dans un objet nomm´ e octets.
5.2- Lors de l’´ ecriture de ces donn´ ees, chaque paire d’octets a en fait ´ et´ e ´ ecrite
en sens inverse (c’est-` a-dire par exemple la paire d’octets 02 56 a ´ et´ e
´ ecrite 56 02). Vous devez donc permuter toutes les paires d’octets deux
` a deux. Vous stockerez le r´ esultat de cette op´ eration dans x.
5.3- Il faut maintenant transformer cette suite de paires d’octets en valeurs
num´ eriques que l’on pourra repr´ esenter graphiquement. Lorsque l’on a
deux octets (par exemple 02 56), il faut utiliser l’instruction
as.numeric("0x0256") pour obtenir la valeur d´ ecimale correspondant ` a
cette paire d’octets (dans ce cas on obtient 598). Transformez x en valeurs
d´ ecimales et stockez le r´ esultat dans un objet nomm´ e valeurs (indice :
utilisez les fonctions matrix(), apply() et paste()).