I got a dataframe that looks like this. I want to plot for EACH "localidad" the year-mean "rinde" value for the three different "agua_inicial"