je fais des boites à moustache par groupe et j'aimerais ajouter sur chaque panneau :
- la moyenne du groupe
- la moyenne général
Code : Tout sélectionner
# données test
test<-data.frame(classe=c(rep(1:3, 9)), ind=rep(1:9, 3), values=rnorm(n=27, 0, 1))
# aperçu
test[with(test, order(classe, ind)),]
# classe ind values
# 1 1 1 -0.982562538
# 10 1 1 0.671096635
# 19 1 1 -2.144761133
# 4 1 4 1.042936218
# 13 1 4 0.019593993
# 22 1 4 1.636128502
# 7 1 7 0.172440543
# 16 1 7 1.377762072
# 25 1 7 1.226224218
# 2 2 2 0.857910359
# 11 2 2 0.923269540
# 20 2 2 0.268502722
# 5 2 5 -0.372159562
# 14 2 5 0.367394526
# 23 2 5 1.801222973
# 8 2 8 -1.573680102
# 17 2 8 0.169825303
# 26 2 8 -0.008336949
# 3 3 3 0.591973841
# 12 3 3 0.569334749
# 21 3 3 0.272558998
# 6 3 6 0.233889287
# 15 3 6 1.390091984
# 24 3 6 1.466123905
# 9 3 9 1.054160179
# 18 3 9 0.525190927
# 27 3 9 -0.109793719
# boxplot par groupe
ggplot(test, aes(factor(ind), values, fill=factor(classe), group=ind)) + geom_boxplot() + facet_grid(. ~ classe, scales="free_x")
J'ai essayé avec
Code : Tout sélectionner
+ stat_summary(fun.y=mean ,geom="line",color="red", mapping=aes(group=1))
mais j'ai un warning comme quoi il n'y a qu'un élément par groupe ; en remplaçant "line" par "point" j'obtiens la moyenne de chaque "ind" représentée par un point...
j'ai également essayer stat_smooth sans succès.
est-ce qu'il faut passer par un nouveau layer? ou est-ce que c'est faisable comme ça?
Merci