Discussion :

[Wallentyne]

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 :

Sélectionner tout - Visualiser dans une fenêtre à part

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

[Wallentyne]

J'ai trouvé c'était tout simplement à cause des NA. J'ai donc rajouté la commande en gras et ça fonctionne.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
       a <- tapply(data$painmedspk1, data$group,mean, na.rm=TRUE)
       points(a,pch=3)