library(ggplot2) library(plyr) d<-read.csv("boxplot_gamma.csv", header = TRUE, sep = ",", dec = ".") attach(d) names(d) mydata <- within(data.frame(x=Groupe), { y<-Median.Gamma}) ## Dot plot over boxplot gp = ggplot(mydata, aes(x, y, col =x)) + geom_boxplot(outlier.size=0, outlier.stroke=0) + labs(title= "Distribution of individual median 60-110Hz activity \n over parieto-occipital electrodes \n during meditation in each group", x="Groups", y="Median gamma power") gp+stat_boxplot(geom ='errorbar') + geom_point(size = 3, position = position_jitter(width = .3, height=0))+ labs(colour="")