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.

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}

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