apikey <- "AIzaSyBPxhDW9_bUJ1juxXEGvY_QcNIZNJPUy7A"

library(ggplot2)
library(ggmap)

register_google(apikey)

qmap(c(-76.5,3.43), zoom=12, maptype="hybrid")

paradasmio <- read.delim("C:/Users/Carolina/OneDrive/Carolina/Maestria/Visualizacion de datos/Mapa/MIO_STOPS.txt")

mapa <- ggplot(data=paradasmio, aes(x=DECIMALLONGITUDE, y=DECIMALLATITUDE))
mapa <- mapa + geom_point()
mapa

city<-get_map(c(-76.5,3.43), zoom=12)
mapa <- ggmap(city)
mapa <- mapa + geom_point(data=paradasmio,aes(x=DECIMALLONGITUDE,y=DECIMALLATITUDE),color=paradasmio$TYPE_BUS,size=paradasmio$PASSANGERS_NUM*0.02)
mapa

Adicionar columna de agrupamiento

k<-7
nIter<-100

grupos <- kmeans(paradasmio[,c(7,8)], k, iter.max=nIter, nstart=1)

paradasmio$grupos<-grupos$cluster

city<-get_map(c(-76.5,3.43), zoom=12)
mapa <- ggmap(city)
mapa <- mapa + geom_point(data=paradasmio,aes(x=DECIMALLONGITUDE,y=DECIMALLATITUDE),color=paradasmio$grupos,size=paradasmio$PASSANGERS_NUM*0.02)
mapa