I have create a plot with two aesthetics color and shape, but only legend of color is displayed not the shape.
sav %>% 
  ggplot(aes(x = Subject, y = Percentage, color = Subject, shape = Location)) +
  geom_point(data = filter(sav, Location == "IRELAND"), shape = "square",  size = 3) + 
  geom_point(data = filter(sav, Location == "OECD"), shape = "diamond",  size = 3) +
  geom_point(data = filter(sav, !Location %in% c("OECD", "IRELAND")), shape = "circle", 
             size = 2, position = "jitter") +
  coord_flip() +
  theme(legend.position = "top")+
  theme_light()
dput(sav):
structure(list(Location = c("CANADA", "CANADA", "FRANCE", "FRANCE", 
"GERMANY", "GERMANY", "IRELAND", "IRELAND", "ITALY", "ITALY", 
"NETHERLANDS", "NETHERLANDS", "SPAIN", "SPAIN", "UNITED STATES", 
"UNITED STATES", "OECD", "OECD", "CANADA", "ITALY", "SPAIN", 
"FRANCE", "IRELAND", "UNITED STATES", "NETHERLANDS", "OECD", 
"GERMANY"), Subject = c("Food", "Total", "Food", "Total", "Food", 
"Total", "Food", "Total", "Food", "Total", "Food", "Total", "Food", 
"Total", "Food", "Total", "Food", "Total", "Total_Minus_Food_Energy", 
"Total_Minus_Food_Energy", "Total_Minus_Food_Energy", "Total_Minus_Food_Energy", 
"Total_Minus_Food_Energy", "Total_Minus_Food_Energy", "Total_Minus_Food_Energy", 
"Total_Minus_Food_Energy", "Total_Minus_Food_Energy"), Frequency = c("Monthly", 
"Monthly", "Monthly", "Monthly", "Monthly", "Monthly", "Monthly", 
"Monthly", "Monthly", "Monthly", "Monthly", "Monthly", "Monthly", 
"Monthly", "Monthly", "Monthly", "Monthly", "Monthly", "Monthly", 
"Monthly", "Monthly", "Monthly", "Monthly", "Monthly", "Monthly", 
"Monthly", "Monthly"), Time = c("2022-12", "2022-12", "2022-12", 
"2022-12", "2022-12", "2022-12", "2022-12", "2022-12", "2022-12", 
"2022-12", "2022-12", "2022-12", "2022-12", "2022-12", "2022-12", 
"2022-12", "2022-12", "2022-12", "2022-12", "2022-12", "2022-12", 
"2022-12", "2022-12", "2022-12", "2022-12", "2022-12", "2022-12"
), Percentage = c(11.02015, 6.319445, 12.86678, 5.850718, 19.75631, 
8.550855, 11.74636, 8.224299, 13.14815, 11.63227, 16.7983, 9.586879, 
15.68565, 5.70769, 11.88275, 6.454401, 15.60381, 9.438622, 5.58275, 
4.469475, 4.442303, 3.36004, 4.999758, 5.707835, 6.034873, 7.221961, 
5.05511)), class = "data.frame", row.names = c(NA, -27L))

 
     
    

