Calculating and plotting standard errors using a dataset in long format

Question

I am making a plot for likeability of different ice cream flavours by school. I already have the code for making the plot, but I\'m stuck on calculating and adding the error bar