Group by a dataframe and iterate over subplots for each group