Moyenne sur les boxplot

Bonjour je souhaite afficher la moyenne sur les boxplot. Pour cela j'ai mis le code qui est en gras mais cela ne fonctionne que pour le premier. Pour les 2 autres rien ne s'affiche et je ne sais pas pourquoi.
Savez-vous pourquoi ou y a t-il une autre solution ?

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
        
        
        par(mfrow=c(1,3))
        boxplot(data$painmedspk1~data$group, xlab = "Groupes à l'inclusion", 
                ylab = "Nombre d'antalgiques consommés",
                ylim = c(0,560),
                col = c("lightblue", "plum"),
                names = c("Groupe soin usuel", "Groupe acupuncture"))
       a <- tapply(data$painmedspk1, data$group,mean)
       points(a,pch=3)
        
        boxplot(data$painmedspk2~data$group, xlab = "Groupes à 3 mois", 
                ylab = "Nombre d'antalgiques consommés",
                ylim = c(0,560),
                col = c("lightblue", "plum"),
                names = c("Groupe soin usuel", "Groupe acupuncture"))
       b <- tapply(data$painmedspk2, data$group,mean)
       points(b,pch=3)
        
        boxplot(data$painmedspk5~data$group, xlab = "Groupes à 1 an", 
                ylab = "Nombre d'antalgiques consommés",
                ylim = c(0,560),
                col = c("lightblue", "plum"),
                names = c("Groupe soin usuel", "Groupe acupuncture"))
       c <- tapply(data$painmedspk1, data$group,mean)
       points(c,pch=3)

Merci pour votre aide