Correction des exercices sur les courbes (logiciel R)

tablo = read.table("exercice5.txt",sep="\t",header=T)

colnames(tablo)

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")

# Calculer les valeurs moyennes pour chaque jour et chaque condition grâce à means

result_O = by(tablo$Fatigue[tablo$Traitement=="O"],tablo$Temps[tablo$Traitement=="O"], mean)

result_O

result_N = by(tablo$Fatigue[tablo$Traitement=="N"],tablo$Temps[tablo$Traitement=="N"],mean)

result_N

# Tracer le graphique

plot(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")

# Charger les librairies qui permettront de tracer les intervalles de confiance

install.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))