Year Area s mean sd se
1 2004 Gootebank 9 0.2158556 0.1188472 0.03961573
2 2004 Thornton 4 1.9564700 1.9369257 0.96846283
3 2017 Gootebank 13 1.0664641 1.7131108 0.47513144
4 2017 Thornton 10 1.9384720 2.3308575 0.73708186
5 2018 Gootebank 13 1.4408562 1.6793808 0.46577643
6 2018 Thornton 10 1.8738613 3.0972972 0.97945136

I am trying to make a grouped boxplot that presents the mean values shown, in the table for each sandbank, for each year (which is 2004, 2017 and 2018 here). So should be a simple grouped boxplot. Im new with ggplot2 and can't find a way to make R separate each year in the plot

I'm not super confident about what you hope to see but perhaps this gets closer ?
Also could you gather the actual min and max values for your samples ? I fake them to the boxplot by assuming they are at 3 standard deviations

Thank you very much for your reply, yes that definately gets closer!

I was hoping to obtain a boxplot similar to this:

Where x-axis is the years 2014, 2017 and 2018, y axis the mean and the subgroups are the sandbanks Goote & Thornton. I thought using "fill= factor(Area)" in aes() would be sufficient for R to recognize the subgroups for each year.