Estou trabalhando com o ggmap e gostaria de plotar a taxa de uma doença usando o cols[findInterval...]. Isso é possível ou só posso usar o geom_point?

Meu script é:

stores <- data.frame(name=c("Jaboatão dos Guararapes", "Olinda","Paulista", "Recife", "Maceió"),

longitude=c(-35.00329,-34.86707,-34.88437,-34.93294, -35.71009),

latitude=c(-8.152749,-7.993280,-7.927770,-8.039684,-9.523036),

disease=c(143,226,24,142,14))

location = c(-45.6,-10)

nordeste = get_map(location = location, source = "osm", zoom=5)

nordesteMap = ggmap(nordeste)

Obrigado