I'm having trouble plotting a geom_point() layer on top of a geom_boxplot() layer in ggplot2, and have done some research and there don't seem to be any reported issues precisely of this nature. There are 3 factors in my data set: name, genotype, and region, and my response variable is volume. I have working code to produce a plot with both layers. The problem is that the points ignore the fill factor for geom_point(), but not for geom_boxplot(). The result is that the points are all plotted in the middle of a set of boxplots for each value of name. Here is my code for constructing the plot.
meansPlot = ggplot(data=meansData,aes(x=factor(name), y=volume, fill=factor(genotype)))
meansPlot = meansPlot +
geom_boxplot() +
geom_point() +
facet_wrap( ~ region, scales='free')
My apologies for not creating a reproducible data set -- I am not really well versed in simulating data quite yet. If there isn't an easy answer (which, I expect there is, and I'm probably just missing something), I will add simulated data to help answer the question.
Thanks!