Пользовательская легенда для пузырьковой диаграммы с использованием ggplot

У меня есть фрейм данных, который содержит пространственную переменную с смоделированными и наблюдаемыми значениями.

df <- data.frame(sim = sample((20:30),10),
                 obs = sample(25:40,10),
                 long = rnorm(10,10,8),
                 lat = rnorm(10,30,15))

Я построил пространственный пузырьковый график, используя следующий код.

ggplot() +
geom_polygon(data = shp, aes(x = long, y = lat, group = group),
             col = "black", lwd = 0.8, fill = "slategray1",
             alpha = 0.5) +
coord_cartesian() + 
geom_point(data = df, aes(x = long, y = lat, group = sim, size = sim), col = "red", alpha = 0.5) +
geom_point(data = df, aes(x = long, y = lat, group = obs, size = obs), col = "blue", alpha = 0.5) +
scale_size_continuous(range = c(10,20)) +
coord_map(xlim = c(-1, 22), ylim = c(10, 45))

Результат выполнения приведенного выше кода показан на картинке. На графике я хочу показать разницу между смоделированными и наблюдаемыми данными, которая сделана идеально. Теперь я хочу показать две разные легенды для двух переменных sim и obs, которые будут иметь красный и синий цвета соответственно. Также я хочу показать величину переменной в легенде, но не как пузырьки переменного размера, а как-то иначе. Что-то вроде данного здесь < / а>. Может ли кто-нибудь помочь мне в этом?