Traitement et Analyse Informatiques de Données Biologiques (STAT0077)
TP 3: Graphique et data.frame
Prof. Patrick E. Meyer
Version 2.0
A l'aide des fonctions data(),pairs(), plot(), lines(),
title(), pdf() et dev.off()
Charger le dataset CO2 et
- Visualisez graphiquement les relations entre toutes les paires de variables parmi
$Type, $Treatment,$conc et $uptake.
- Dans le graphique obtenu, choisissez la paire de variables qui semble la plus pertinente et re-dessinez ce graphique avec la fonction
plot.
- Ajoutez, en pointillé vert, la diagonale du graphique, ainsi qu'un titre (en noir) ayant du sens.
- Sauvez le graphique dans un fichier .pdf et vérifiez le résultat obtenu dans le fichier spécifié.
A l'aide des fonctions write.table(), read.table(), intersect(), which() et which.max()
Charger le dataset CO2 et
- Écrire un code R qui écrit le contenu du dataset
CO2, avec des points virgules comme séparateur et des virgules pour les décimales, dans un fichier nommé mondataset.txt.
- Charger le dataset
mondataset.txt dans une variable mydata avec les noms de lignes et colonnes adéquatement formatés (identique à CO2).
- Écrire un code R qui compte le nombre de plantes provenant du Quebec qui ont, une absorption en
supérieur à 35 (
), et qui n'ont pas été refroidie dans la nuit qui précède la mesure.
- Écrire un code R qui retourne la plante (l'indice de la ligne) du dataset
CO2 qui absorbe le plus de
lorsque la concentration en
dans l'atmosphère est la plus basse.
- Écrire un code R qui crée un vecteur contenant le maximum d'absorption atteint par chaque plante du dataset
CO2 (utilisez une boucle implicite).
- Faites un
barplot avec les minima et maxima d'absorption pour chaque type de plantes, mettez y des couleurs et titres appropriés.
- Créez une
heatmap de la matrice construite avec les valeurs de $uptake, avec les lignes correspondantes aux différentes concentrations ($conc) et les colonnes correspondantes aux différentes plantes ($Plant).
- Combinez les graphiques en mettant côte à côte le nuage de points réalisé en début de séance (cf EX1.2) et le barplot (utilisez
par(mfrow=)).
<--- BACK