Mettre les écart types sur un graph suite à une ANOVA

Bonsoir,

Je viens d'effectuer une ANOVA et une SNK sur une plage de donnée contenant 15 modalités * 4 répétitions. Maintenant j'aimerai synthétiser le tout sur un graphique, j'ai fait un premier jet dans lequel j'indique les groupes statistiques, mais je n'arrive absolument pas à mettre les écart types, quelqu'un aurait une idée?

voici mon code:

Code:

Y <- ble$Leaf.1
Produits <- ble$Produits
par(mfrow=c(1,2))
stripchart(Y~Produits, xlab="Produits", ylab="Leaf.1", vertical=T)
boxplot(Y~Produits, xlab="Produits", ylab="Leaf.1")
Res <- lm(Y~Produits)
summary(Res)
stripchart(Res$residuals~Produits, xlab="Produits", ylab="Résidus", vertical=T)
plot(Res$residuals~Res$fitted,xlab="Rendements attendus",ylab="Résidus")
hist(Res$residuals,freq=F)
curve(dnorm(x,mean(Res$res),sd(Res$res)),add=T)
ks.test(Res$residuals,"pnorm",mean(Res$residuals),sd(Res$residuals))
shapiro.test(Res$residuals)
bartlett.test(Res$residuals,Produits)
library(car)
leveneTest(Res)
anova(Res)
model <- aov(Leaf.1 ~ Produits, data=ble)
library(agricolae)
SNK.test(model,"Produits")
comp <- SNK.test(model,"Produits")
comp
par(mfrow=c(1,1))
comp <- SNK.test(model,"Produits")$groups
g <- barplot(comp$means,names.arg=comp$trt,ylim=c(0,100),col = "lightblue",border = "blue", main = "% Efficacité", ylab = "Efficacité des produits par rapport au témoin (%)", xlab = "Identifiant de la modalité")
text(g,comp$means+3,comp$M)

Et voici mes données:

Code:

Produits Leaf 1
B -5,38
B 15,07
B 41,75
B 17,45
C 56,77
C 83,14
C 73,81
C 81,40
D 58,82
D 85,02
D 83,20
D 81,12
E 68,17
E 81,70
E 77,22
E 77,25
F 75,38
F 78,57
F 79,18
F 76,11
G 80,43
G 82,69
G 74,02
G 85,55
H 75,27
H 81,52
H 75,36
H 71,53
I 73,55
I 80,09
I 71,86
I 62,23
J 80,75
J 74,08
J 80,31
J 74,54
K 65,70
K 72,02
K 57,84
K 77,40
L 80,43
L 86,10
L 82,58
L 82,12
M 2,69
M 52,91
M 44,02
M 35,05
N 29,57
N 45,74
N 35,05
N 27,04
O 52,69
O 71,84
O 77,63
O 75,97

Merci

Bonjour,

Cette fonction existe dans différents packages. Au hasard, en faisant une recherche google :
http://www.inside-r.org/packages/cran/plotrix/docs/plotCI

Merci pour votre réponse, elle m'a bien aidé.