Courbe moyenne avec barres d'erreur et intervalle de confiance
Comment réaliser simplement un graphique avec intervalles de confiance en calculant une moyenne pour chaque paramètre.
Exemple : Graphique avec barres d'erreur
1- Charger les données - ici on va en simuler 2 variables : poids, jour
jour <- rep(1:10,each=6) ; jour # 10 jours, 6 mesures par jour
poids <- rnorm(60,20,4)+rnorm(60,5,2)*jour # Des valeurs au pif qui montent avec le jour
plot(jour,poids,pch=16,cex=1.5)
Bien, nous avons 6 mesures par jour : il reste à faire une courbe moyenne et ajouter les barres d'erreur (intervalles de confiance)
2- Calculer les moyennes avec by()
mes_moyennes <- by(poids,jour,mean)
3- Calculer les intervalles de confiance avec by() et int.ech de {KefiR}
library(KefiR) # Pour installer KefiR, cliquer ici.
mes_IC <- by(poids, jour,int.ech)
4- Tracer la courbe avec plotCI de {gplots}
#install.packages('gplots') # Cette commande installe le package si jamais installé.
library(gplots)
plotCI(as.numeric(names(mes_moyennes)),
mes_moyennes,
uiw=mes_IC,
col="red",gap=0,type="o")
Mise à jour de la page septembre 2023