Correction des exercices sur les courbes (logiciel R)
- Correction - Exercice 2 - Courbe avec intervalles de confiance
Consignes
1- Ouvrir le tableau "exercice5.txt".
tablo = read.table("exercice5.txt",sep="\t",header=T)colnames(tablo)2- Tracer un graphique exprimant l'Acuité en fonction du Temps pour tous les individus testés.
plot(tablo$Temps[tablo$Traitement=="O"],tablo$Fatigue[tablo$Traitement=="O"], col="blue", xlab="Temps",ylab="Acuité")points(tablo$Temps[tablo$Traitement=="N"],tablo$Fatigue[tablo$Traitement=="N"], col="red")3- Tracer cette fois-ci l'évolution de l'acuité moyenne en fonction du temps selon l'exemple donné (cf. exemple ci-dessous)
# Calculer les valeurs moyennes pour chaque jour et chaque condition grâce à meansresult_O = by(tablo$Fatigue[tablo$Traitement=="O"],tablo$Temps[tablo$Traitement=="O"], mean)result_Oresult_N = by(tablo$Fatigue[tablo$Traitement=="N"],tablo$Temps[tablo$Traitement=="N"],mean)result_N# Tracer le graphiqueplot(names(result_O),result_O,col="blue",gap=0,pch=16,ylim=c(0,10),type="o")points(names(result_N),result_N,col="red",gap=0,pch=16,type="o")4- Refaire le même graphique en superposant cette fois-ci les intervalles de confiance...
# Charger les librairies qui permettront de tracer les intervalles de confianceinstall.packages("gplots")library("gplots")# Tracer les graphiques en remplaçant plot par plotCI (qui ajouter les intervalles de confiance)plotCI(result_O$parametres,result_O$moyennes, uiw=result_O$ic,col="blue",gap=0,pch=16,ylim=c(0,10),type="o")plotCI(result_N$parametres,result_N$moyennes, uiw=result_N$ic,col="red",gap=0,pch=16,type="o",add=T)# le paramètre add permet de superposer les courbes# Tracer une légende (cf. aide légende)legend(x="topright", legend=c("Avec traitement","Sans traitement"), col=c("blue","red"), pch=c(16,16))