tengo un shapefile del area de Camden de UK el cual se puede descargar aquí (junto con los datos): https://drive.google.com/file/d/1Py9BrlB8i0Nz9P_tq92s7Q124iIW-1JQ/view?usp=sharing
hasta el momento cargo las base de datos del Shapefile y luego la base de datos de los puntos y proyecto ambas capas sobre el mismo grafico. esto es lo que llevo hasta ahora.
Census.Data <- read.csv("practical_data.csv")
if (!require("rgdal")) install.packages("rgdal")
library("rgdal")
if (!require("rgeos")) install.packages("rgeos")
library("rgeos")
Output.Areas<- readOGR("D:/Camden/shapefiles", "Camden_oa11")
OA.Census <- merge(Output.Areas, Census.Data, by.x="OA11CD", by.y="OA")
houses <- read.csv("camdenhousesales15.csv")
#solo necesito el UID, el precio y las cordenadas.
houses <- houses[,c(1,2,8,9)]
#lo anterior no tiene atributos espaciales, para poder hacerlo es necesario utilizar el paquete sp.
if (!require("sp")) install.packages("sp")
library("sp")
#Podemos consultar las referencias en https://spatialreference.org/
houses.points <- SpatialPointsDataFrame(coords = houses[,3:4],data = houses,proj4string = CRS("+init=EPSG:27700"))
# Creando un mapa base con la libreria "tmap".
if (!require("tmap")) install.packages("tmap")
library("tmap")
tm_shape(OA.Census) + tm_borders(alpha=.4) +
tm_shape(houses.points) + tm_dots(col = "Price", scale = 2.5, palette = "Reds", style = "quantile", title = "Price Paid (£)")
Con lo anterior obtengo el siguiente gráfico.
Lo que realmente quiero es contar cuantos puntos caen dentro de cada polígono del mapa de Camden y guardar esto como una variable mas en los datos de OA.Census.
he intentado lo siguiente pero sin éxito alguno
cuenta <- over(houses.points, OA.Census)
¿¿ count <- as.integer(table(cuenta$id)) ???