Considérons un exemple
la masse des larves de papillons varie avec la température et de la nourriture (plante)
Pour tracer un barplot avec les barres d’erreur de ces deux facteurs :
attach(data)
mean.mass<-tapply(mass,list(T,H),mean) # calcul de la moyenne pour chaque groupe
sd<-tapply(mass,list(T,H),sd) # ou T est le facteur Température et H le facteur plante
length<-tapply(DTp,list(T,H),length)
sqrt<-sqrt(length)
length<-as.numeric(length)
se<-sd/sqrt #calcul de l'erreur standard pour chaque groupe
graphemass<-barplot(mean.mass,beside=TRUE,density=c(0,30),col=1,ylim=c(0,15.2),font=3,ylab="masse (mg)",)
library(nmle) #nécessaire pour appeler la fonction "segments"
segments(graphemass, mean.mass -se,graphemass, mean.mass +se)
segments(graphemass-0.1, mean.mass -se,graphemass+0.1, mean.mass -se)
segments(graphemass-0.1, mean.mass +se,graphemass+0.1, mean.mass +se)]