Soy un R-Principiante, pero quiero usarlo para algunos de estadística espacial. Quiero comprobar la relación entre la ocurrencia de las personas enfermas y de la existencia de agua en un área de interés. Para dar una idea más específica de casos de uso, he creado algunos datos de ejemplo:
En primer lugar, quiero trazar en un mapa, incluyendo un mapa base. Esto ya funciona!
persons_dataframe = read.csv("persons.csv", header = TRUE)
city_dataframe = read.csv("city.csv", header = TRUE)
water_dataframe = read.csv("water.csv", header = TRUE)
# plot them on a map
# load used packages
library(RgoogleMaps)
library(ggplot2)
library(ggmap)
library(sp)
persons_ggplot2 <- persons_dataframe
city_ggplot2 <- city_dataframe
water_ggplot2 <- water_dataframe
gc <- geocode('new york, usa')
center <- as.numeric(gc)
G <- ggmap(get_googlemap(center = center, color = 'bw', scale = 1, zoom = 10, maptype = "terrain", frame=T), extent="device")
G1 <- G + geom_point(aes(x=POINT_X, y=POINT_Y ),data=city_dataframe, shape = 22, color="black", fill = "yellow", size = 4) + geom_point(aes(x=POINT_X, y=POINT_Y ),data=persons_dataframe, shape = 8, color="red", size=2.5) + geom_point(aes(x=POINT_X, y=POINT_Y ),data=water_ggplot2, color="blue", size=1)
plot(G1)
Esto tiene el siguiente resultado:
No quiero crear algo así como un clúster de mapa. Ya he buscado algunos paquetes útiles y he terminado con los paquetes SpatialEpi y DCluser. Ambos proporcionan algunas funciones, que pueden ser útiles.
Pero en este punto no sé, cómo proceder. Algunas preguntas abiertas:
- Es lo justo que necesito polígonos para agrupar el conjunto de datos? Me lo he leído en el SpatialEpi pdf pero solo tengo estos de punto conjuntos de datos.
- Además no sé si el hecho de que todos los puntos se encuentran en sólo tres lugares podrían ser un problema?!
¿Tiene sugerencias de cómo continuar con mi trabajo?
