I managed to solve the problem myself using: q=sort(unique(foram$stage[which(foram$length>0)]))
boxplot(foram$logl~foram$stage, data=foram, outline=F, at=q, axes=TRUE, add=TRUE, col="gray82", medlwd=1) points(foram$stage,foram$logl, cex=.1) However, for future reference since this was my first post and I was trying to give code like you requested, what did I do wrong? -- View this message in context: http://r.789695.n4.nabble.com/Boxplots-over-a-Scatterplot-tp2281526p2281690.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.