Linkedin: Mauricio Márquez Goa
Venezuela, un país ubicado en América del Sur es actualmente un territorio donde los índices delictivos van en aumento producto de las precarizadas condiciones económicas, sociales y políticas predominantes.
Para el análisis delictivo se obtuvieron los delitos ocurridos de forma extraoficial, los cuales quedan tipificados de la siguiente forma:
Homicidio. Violación. Secuestro. Robo Vehículo. Enfrentamiento. Robo Genérico. Hurto Vehículo. Lesiones Personales. Hurto Genérico.
# Procedemos a llamar las librerías para conectarnos con las base de datos Postgres + PostGIS y acceder los datos
library(tidyverse)
library(sf)
library(ggplot2)
library(lubridate)
library(osmdata)
library(ggmap)
library(leaflet)
library(ggspatial)
#Capa de Delitos
delitos <- st_read("C:/Users/20958699825/Desktop/QGIS Blog/criminalidad/R Analisis/delitos.geojson")
## Reading layer `delitos' from data source `C:\Users\20958699825\Desktop\QGIS Blog\criminalidad\R Analisis\delitos.geojson' using driver `GeoJSON'
## Simple feature collection with 58956 features and 69 fields
## geometry type: MULTIPOINT
## dimension: XY
## bbox: xmin: -72.66719 ymin: 4.577148 xmax: -60.90665 ymax: 12.05013
## epsg (SRID): 4326
## proj4string: +proj=longlat +datum=WGS84 +no_defs
#Capa de Estados
estado <- st_read("C:/Users/20958699825/Desktop/QGIS Blog/criminalidad/R Analisis/estados.geojson")
## Reading layer `estados' from data source `C:\Users\20958699825\Desktop\QGIS Blog\criminalidad\R Analisis\estados.geojson' using driver `GeoJSON'
## Simple feature collection with 24 features and 4 fields
## geometry type: MULTIPOLYGON
## dimension: XY
## bbox: xmin: -73.37767 ymin: 0.647146 xmax: -59.80029 ymax: 12.20101
## epsg (SRID): 4326
## proj4string: +proj=longlat +datum=WGS84 +no_defs
#Capa de Municipios
municipios <- st_read("C:/Users/20958699825/Desktop/QGIS Blog/criminalidad/R Analisis/municipios.geojson")
## Reading layer `municipios' from data source `C:\Users\20958699825\Desktop\QGIS Blog\criminalidad\R Analisis\municipios.geojson' using driver `GeoJSON'
## Simple feature collection with 336 features and 9 fields
## geometry type: MULTIPOLYGON
## dimension: XY
## bbox: xmin: -73.37941 ymin: 0.634551 xmax: -59.74104 ymax: 12.19504
## epsg (SRID): 4326
## proj4string: +proj=longlat +datum=WGS84 +no_defs
#Capa de Parroquias
parroquias <- st_read("C:/Users/20958699825/Desktop/QGIS Blog/criminalidad/R Analisis/parroquias.geojson")
## Reading layer `parroquias' from data source `C:\Users\20958699825\Desktop\QGIS Blog\criminalidad\R Analisis\parroquias.geojson' using driver `GeoJSON'
## Simple feature collection with 1125 features and 8 fields
## geometry type: MULTIPOLYGON
## dimension: XY
## bbox: xmin: -73.37967 ymin: 0.6439212 xmax: -59.79583 ymax: 12.19785
## epsg (SRID): 4326
## proj4string: +proj=longlat +datum=WGS84 +no_defs
dim(delitos)
## [1] 58956 70
La cantidad de delitos denunciados es de 58956 con 70 columnas
summary(delitos)
## id estado municipio
## Min. : 156 MIRANDA : 8740 LIBERTADOR: 8493
## 1st Qu.:16558 DISTRITO CAPITAL: 7449 SUCRE : 3976
## Median :32972 CARABOBO : 5416 VALENCIA : 2282
## Mean :33030 ZULIA : 4335 IRIBARREN : 2222
## 3rd Qu.:49506 BOLÍVAR : 3915 CARONÍ : 2216
## Max. :65950 ARAGUA : 3604 MARACAIBO : 1508
## (Other) :25497 (Other) :38259
## parroquia capital ciudad
## SUCRE : 1583 CARACAS : 4426 CARACAS : 1964
## PETARE : 1348 VALENCIA : 2282 GUARENAS : 234
## CHACAO : 1178 CIUDAD GUAYANA: 2216 HIGUEROTE : 194
## URBANA MIGUEL PEÑA: 874 BARQUISIMETO : 2203 LOS TEQUES : 186
## CATEDRAL : 850 MARACAIBO : 1470 CARACAS-PETARE: 133
## EL CARMEN : 833 PETARE : 1348 (Other) : 5150
## (Other) :52290 (Other) :45011 NA's :51095
## sector cuadrante_delito organo.actuante
## CENTRO : 687 P-01 : 7672 CICPC :58954
## CHACAO : 504 ZONA PERIURBANA: 6578 GAES 41 CARABOBO: 1
## CATIA : 365 P-02 : 6152 GAES SUCRE : 1
## MATURÍN : 344 P-03 : 4763
## SECTOR SAN JUAN: 328 P-04 : 4192
## (Other) :56661 P-06 : 2900
## NA's : 67 (Other) :26699
## número delito
## Min. : 1.00 HURTO :23247
## 1st Qu.: 17.39 ROBO : 9733
## Median : 33.81 LESIONES PERSONALES : 7863
## Mean : 40.42 ROBO DE VEHÍCULO : 5796
## 3rd Qu.: 50.35 HOMICIDIO : 4714
## Max. :998.00 RESISTENCIA A LA AUTORIDAD: 3605
## (Other) : 3998
## móvil día mes año
## AJUSTE DE CUENTA: 2846 Min. : 1.00 Min. :1.000 Min. :2019
## ROBO : 632 1st Qu.: 8.00 1st Qu.:2.000 1st Qu.:2019
## POR DETERMINAR : 578 Median :15.00 Median :5.000 Median :2019
## RIÑA : 200 Mean :15.49 Mean :4.752 Mean :2019
## PASIONAL : 116 3rd Qu.:23.00 3rd Qu.:7.000 3rd Qu.:2019
## (Other) : 351 Max. :31.00 Max. :9.000 Max. :2019
## NA's :54233
## hora
## 5:00:00 p. m. : 2909
## 6:00:00 p. m. : 2547
## 11:00:00 a. m.: 2291
## 12:00:00 p. m.: 2253
## 1:00:00 p. m. : 2123
## (Other) :46825
## NA's : 8
## dirección.exacta
## SECTOR ALTAMIRA, DISTRIBUIDOR DE ALTAMIRA, VIA PUBLICA, PARROQUIA CHACAO, MUNICIPIO CHACAO, ESTADO MIRANDA : 51
## SECTOR AEROPUERTO, LAS INSTALACIONES DEL CUERPO DE ATENCIÓN AL CIUDADANO 623, DE LA GUARDIA NACIONAL BOLIVARIANA, UBICADO EN LA POBLACIÓN DE SANTA ELENA DE UAIREN, PARROQUIA SANTA ELENA DE UAIREN, MUNICIPIO GRAN SABANA, ESTADO BOLÍVAR.: 43
## SECTOR CASCO CENTRAL , CALLE PRINCIPAL VIA PÚBLICA, POBLACIÓN DE SANTA ELENA DE UAIREN, MUNICIPIO GRAN SABANA, ESTADO BOLÍVAR : 34
## SECTOR CENTRO, ENTRE CALLE 15 Y 16 CARRERA 5TA, GUANARE, MUNICIPIO PAEZ, PARROQUIA GUANARE, ESTADO PORTUGUESA : 23
## SECTOR KUMARAKAPAY, TRONCAL 10, COMUNIDAD INDÍGENA SAN FRANCISCO DE YURUARIN DE KUMARAKAPAY, PARROQUIA GRAN SABANA MUNICIPIO GRAN SABANA, ESTADO BOLÍVAR, : 15
## KILOMETRO 8 DE LA CARRETERA PANAMERICANA, ENTRADA DE LA URBANIZACIÓN FELIPE ACOSTA, VÍA PÚBLICA, PARROQUIA COCHE, MUNICIPIO BOLIVARIANO LIBERTADOR : 14
## (Other) :58776
## categorización.de.referencias
## VÍA PÚBLICA :22071
## BARRIO O CASERÍO :15173
## URBANIZACIÓN O CONJUNTO RESIDENCIAL: 4043
## INSTALACIONES DEL ESTADO VENEZOLANO: 3238
## COMERCIO : 3187
## (Other) :11234
## NA's : 10
## número.de.cadáveres
## : 29
## : 1
## 0 : 735
## 1 : 7663
## NA's :50528
##
##
## arma.utilizada heridas.producidas
## ARMA DE FUEGO :22544 MÚLTIPLES : 3927
## FUERZA FÍSICA : 6182 UNA : 752
## ARMA BLANCA : 1986 MULTIPLES : 13
## OBJETO CONTUNDENTE: 766 GRAVE : 3
## OTRA : 497 NO APLICA : 2
## (Other) : 223 (Other) : 5
## NA's :26758 NA's :54254
## despacho tipo.de.lesión
## SUB DELEGACIÓN VALENCIA : 2130 LEVE : 6784
## SUB DELEGACIÓN CIUDAD GUAYANA: 1443 GRAVE : 1060
## SUB DELEGACIÓN CHACAO : 1418 LEVE : 8
## SUB DELEGACIÓN El LLANITO : 1400 N/A : 7
## SUB DELEGACIÓN MATURÍN : 1318 ESTABLE: 1
## (Other) :51246 (Other): 3
## NA's : 1 NA's :51093
## nombre.y.apellido sexo.victima
## POR DETERMINAR : 2201 MASCULINO :38813
## EL ESTADO VENEZOLANO: 322 FEMENINO :18334
## POR IDENTIFICAR : 167 POR DETERMINAR: 1064
## POR DETERMINAR : 144 NO APLICA : 633
## ESTADO VENEZOLANO : 46 MASCULINO : 45
## DE OFICIO : 43 (Other) : 65
## (Other) :56033 NA's : 2
## edad nacionalidad prontuario.policial_1
## POR DETERMINAR: 3509 VENEZOLANA :54038 POR DETERMINAR :52116
## 38 : 1596 POR DETERMINAR: 2671 NO POSEE : 2675
## 40 : 1570 NO APLICA : 723 POSEE : 463
## 35 : 1563 VENEZOLANA : 629 POR DETERMINAR : 33
## 39 : 1516 COLOMBIA : 271 NINGUNO : 17
## (Other) :49199 (Other) : 621 (Other) : 44
## NA's : 3 NA's : 3 NA's : 3608
## número.de.victimarios comisión.policial hubo.víctima tipo.de.víctima
## 1 : 3609 PNB / FAES: 788 NO : 243 CIVIL : 3423
## 2 : 143 PNB : 114 SI : 3438 FUNCIONARIO: 19
## 0 : 118 PNB/FAES : 111 SI : 3 NA's :55514
## 4 : 74 GNB : 97 NA's:55272
## 3 : 58 CICPC : 77
## (Other): 134 (Other) : 2501
## NA's :54820 NA's :55268
## condición.de.la.víctima número.de.victimas.civiles número.de.detenidos
## OCCISO : 7636 : 2 0 : 3587
## ESTABLE : 6552 0 : 9 1 : 61
## DETENIDO: 780 1 : 3413 2 : 14
## ILESO : 582 NA's:55532 9 : 5
## GRAVE : 486 3 : 3
## (Other) : 216 (Other): 7
## NA's :42704 NA's :55279
## patrimonio.afectado objeto.hurtado.robado
## PERSONA :12298 UN TELEFONO CELULAR : 559
## VIVIENDA : 8493 ELECTRODOMESTICOS : 283
## ESTADO VENEZOLANO: 5767 UN (01) TELEFONO CELULAR: 223
## ENTIDAD PRIVADA : 3191 TELEFONO CELULAR : 213
## COMERCIO : 3063 UN TELEFONO CELULAR : 202
## (Other) : 159 (Other) :31532
## NA's :25985 NA's :25944
## monto.aprox..de.lo.hurtado.robado cant..veh clasedevehículo
## 1000000: 1975 1 : 9207 MOTOCICLETA: 3389
## 2000000: 1825 2 : 35 AUTOMÓVIL : 2971
## 500000 : 1586 3 : 11 CAMIONETA : 1704
## 3000000: 1439 : 3 MOTO : 560
## 5000000: 1308 | : 2 CAMIÓN : 536
## (Other):34104 (Other): 7 (Other) : 111
## NA's :16719 NA's :49691 NA's :49685
## marca modelo año.del.vehículo
## CHEVROLET : 1551 HORSE KW-150: 229 2013 : 1163
## EMPIRE KEEWAY : 1173 BR-150 : 222 2012 : 1009
## TOYOTA : 842 HORSE : 221 2011 : 805
## FORD : 817 COROLLA : 216 2008 : 676
## BERA : 659 AVEO : 209 2007 : 674
## (Other) : 4229 (Other) : 8174 (Other): 4940
## NA's :49685 NA's :49685 NA's :49689
## color.del.veh. tipo.de.veh. uso.del.vehículo
## AZUL : 1664 PASEO : 3651 PARTICULAR : 8536
## NEGRO : 1545 SEDAN : 2600 CARGA : 535
## BLANCO : 1416 SPORT WAGON: 986 OFICIAL : 94
## ROJO : 1333 PICK UP : 374 TRANSPORTE PÚBLICO : 58
## GRIS : 836 PASEO : 208 PARTICULAR / PARTICULAR: 9
## (Other): 2477 (Other) : 1449 (Other) : 39
## NA's :49685 NA's :49688 NA's :49685
## hubo.pago vehículo.incriminado características.del.vehículo
## NO : 43 POR DETERMINAR:29304 POR DETERMINAR: 316
## SI : 36 NO :26345 COLOR NEGRO : 43
## NA's:58877 MOTO : 1437 COLOR ROJO : 32
## AUTOMÓVIL : 1207 COLOR AZUL : 21
## CAMIONETA : 487 COLOR BLANCO : 17
## (Other) : 137 (Other) : 2822
## NA's : 39 NA's :55705
## observaciones
## NINGUNA : 57
## SE INCLUYE POSTERIOR AL MOMENTO DEL REPORTE DEBIDO A QUE EL REDIP NO LO REPORTO EN MISMO DÍA. / / VARIOS SUJETOS POR IDENTIFICAR Y APREHENDER, TOMARON POR ASALTO LAS INSTALACIONES DEL CUERPO DE ATENCIÓN AL CIUDADANO 623, DE LA GUARDIA NACIONAL BOLIVARIANA: 43
## CASOS DE MANIFESTACIONES Y GUARIMBAS / LOS OCCISOS Y LOS LESIONADOS PARA EL MOMENTO EN QUE SE ENCONTRABAN EN UNA MANIFESTACIÓN, LE FUERON CAUSADAS LESIONES POR ARMA DE FUEGO Y OBJETOS CONTUNDENTES. LOGRANDO EN HORAS DE LA TARDE DE HOY MEDIANTE INVESTIGACI: 23
## SE COLOCO EL SECTOR MAS CERCANO A LA DIRECCION EXACTA ANTES MENCIONADA : 19
## SE RECIBE LLAMADA TELEFÓNICA DE PARTE DE LA INSPECTORA AGREGADO DIOSELINA MÁRQUEZ, ADSCRITA A LA RED DE EMERGENCIAS TÁCHIRA VEN 911, INFORMANDO QUE UN GRAN GRUPO DE PERSONAS SE ENCONTRABAN CONCENTRADOS EN LA DIRECCIÓN ANTES MENCIONADA, CON EL FIN DE REALI: 16
## (Other) : 5124
## NA's :53674
## cant. ejes fecha
## Min. :0.0000 - :48825 18/2/2019: 394
## 1st Qu.:1.0000 ALTOS MIRANDINOS: 976 18/3/2019: 382
## Median :1.0000 BARLOVENTO : 828 28/3/2019: 350
## Mean :0.9533 GUARENAS-GUATIRE: 808 11/2/2019: 348
## 3rd Qu.:1.0000 METROPOLITANO : 4180 6/3/2019 : 343
## Max. :1.0000 VALLES DEL TUY : 1586 8/4/2019 : 341
## NA's : 1753 (Other) :56798
## sem nombre_mes nombre_dia clase.dias
## Min. : 1.00 ENERO : 7598 DOMINGO : 5403 ENTRE SEMANA :37969
## 1st Qu.: 9.00 FEBRERO: 7240 JUEVES : 9227 FIN DE SEMANA:20987
## Median :18.00 MAYO : 7063 LUNES :10176
## Mean :19.01 ABRIL : 6907 MARTES : 9392
## 3rd Qu.:28.00 MARZO : 6730 MIÉRCOLES: 9174
## Max. :40.00 JULIO : 6127 SÁBADO : 6589
## (Other):17291 VIERNES : 8995
## semana.calendario rango.de.horas rango.de.edades.de.las.víctimas
## 12 : 1981 #N/D : 10 35 a 39: 7441
## 8 : 1901 00:00 a 06:00: 7695 30 a 34: 6958
## 6 : 1861 06:01 a 12:00:15760 40 a 44: 6779
## 7 : 1790 12:01 a 18:00:26340 25 a 29: 6550
## 3 : 1768 18:01 a 24:00: 8364 45 a 49: 5719
## 11 : 1761 IMPRECISA : 787 20 a 24: 5102
## (Other):47894 (Other):20407
## categoría
## CONTRA LA COSA PÚBLICA : 3085
## CONTRA LA PROPIEDAD :36161
## CONTRA LAS BUENAS COSTUMBRES Y BUEN ORDEN DE LAS FAMILIAS: 370
## CONTRA LAS PERSONAS :10785
## NA's : 8555
##
##
## municipios.priorizados comision.policial.estandarizado
## - : 8638 _ : 8603
## PRIORIZADOS:27199 PNB : 1313
## NA's :23119 CICPC : 1057
## GNB : 522
## POLICIA ESTADAL : 465
## (Other) : 1176
## NA's :45820
## casos.de.resistencias.sin.victimas fecha_latina hora_anglo
## 0,25 : 156 18-2-2019: 394 17:00:00: 2909
## 0,5 : 129 18-3-2019: 382 18:00:00: 2547
## 0,75 : 115 28-3-2019: 350 11:00:00: 2291
## 0,708333333: 114 11-2-2019: 348 00:00:00: 2253
## 0,458333333: 108 6-3-2019 : 343 13:00:00: 2123
## (Other) : 3245 8-4-2019 : 341 (Other) :46024
## NA's :55089 (Other) :56798 NA's : 809
## fecha_hora_anglo long_geom lat_geom
## Min. :2019-01-01 00:00:00 Min. :-72.67 Min. : 4.577
## 1st Qu.:2019-02-28 13:00:00 1st Qu.:-69.28 1st Qu.: 9.636
## Median :2019-05-04 19:57:30 Median :-67.01 Median :10.213
## Mean :2019-05-09 05:41:32 Mean :-67.41 Mean : 9.882
## 3rd Qu.:2019-07-14 09:00:00 3rd Qu.:-66.59 3rd Qu.:10.481
## Max. :2019-09-30 23:30:00 Max. :-60.91 Max. :12.050
##
## geometry
## MULTIPOINT :58956
## epsg:4326 : 0
## +proj=long...: 0
##
##
##
##
Con la consulta a la capa de delitos podemos saber que:
Para hacer el análisis haremos especial énfasis en la columna “delito” la cual define el tipo de delito denunciado en 9 categorias.
Verificamos si posee registros nulos:
sum(is.na(delitos$delito))
## [1] 0
Vamos muy bien!, no hay datos nulos.
Verificamos ahora las fechas y horas
sum(is.na(delitos$hora_anglo))
## [1] 809
Existen 809 casos donde la hora que se comete el delito no es conocida (imprecisa)
sum(is.na(delitos$fecha_hora_anglo))
## [1] 0
Perfecto!, todos los delitos tienen fecha.
Veamos ahora como es la categorización y la cantidad de delitos.
#Creamos una tabla que cuente los delitos por tipos
contador_delitos<- delitos %>%
count(delito) %>%
arrange(desc(n))
Grafiquemos esos datos por el contador
ggplot(contador_delitos) +
geom_col(aes(x=reorder(delito,n), y=n, fill=delito)) +
theme(legend.position="none") +
geom_text(size=6,
aes(
x = delito,
y = n,
label = n))+
coord_flip() +
labs(title = "Total Delitos Denunciados",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
x = "Delitos",
y = "Cantidad")
En la previa ya sabíamos que el delito más denunciado es el hurto con 23247 casos en el período de enero a septiembre del 2019.
Veamos como se comportan por cada entidad en Venezuela los delitos
#Creo un vector con los tipos de delitos
delitos_tipo <- delitos$delito
#Creo una tabla tipo de delitos por Estados
delitos_estados <- delitos %>%
filter(delito %in% delitos_tipo) %>%
count(delito,estado) %>%
arrange(desc(n))
ggplot() +
geom_bar(data=delitos_estados,
aes(x=reorder(delito,n),
weight=n,
fill=factor(delito))) +
facet_wrap(ncol=6,~estado) +
labs(title = "Histograma de Delitos Denunciados",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
fill = "Tipo de Delito",
x = "Delitos",
y = "Cantidad")
El gráfico muestra claramente que los estados centrales del país Miranda, Aragua, Carabobo y Distrito Capital así como para el occidente Lara y Zulia y el Oriente-Sur Anzoátegui y Bolívar son los que poseen mayor cantidad de denuncias.
Veamos ahora como se comportan en el tiempo
delitos <- delitos %>% mutate(fecha = ymd_hms(fecha_hora_anglo))
delitos %>%
ggplot() +
geom_bar(aes(x = month(fecha_hora_anglo, label = TRUE))) +
labs(title = "Total Delitos Denunciados",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
x = "Delitos",
y = "Cantidad")
Existe una tendencia decreciente a la denuncia de delitos, se pronuncia desde el mes de mayo. Examinando los delitos por tipología se verificará los altibajos para cada uno.
delitos %>%
filter(delito %in% delitos_tipo) %>%
ggplot() +
geom_bar(aes(x = month(fecha_hora_anglo, label = TRUE), fill = delito)) +
labs(title = "Total Delitos Denunciados",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
fill = "Tipo de Delito",
x = "Meses",
y = "Cantidad")
A partir del mes de junio del 2019 se ve un decrecimiento en los delitos robo y robo de vehículos (delitos violentos), quizás esté asociado a la implementación de alguna política de seguridad la cual sea efectiva en la persuación de los delicuentes violentos o sea producto de la crisis económica en la obtención de municiones que afecta de manera directa a la delincuencia.
Denuncias Realizadas por mes mediante Poligono de Frecuencia
# Se realiza primeramente un conteo de denuncias por mes
conteo_mes <- delitos %>%
count(delito, mes = month(fecha_hora_anglo, label = TRUE))
# Luego se muestran las cantidades mensuales como polígono de frecuencia
ggplot(conteo_mes) +
geom_line(aes(x = mes, y = n, group = delito, color = delito), size = 2) +
facet_wrap(~delito)+
guides(color = FALSE) +
geom_text(size=3,
aes(
x = mes,
y = n,
label = n))+
labs(title = "Total Delitos Denunciados",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
fill = "Tipo de Delito",
x = "Meses",
y = "Cantidad") +
expand_limits(y = 0)
Ahora se detalla con mayor claridad que el robo de vehículos ha tenido un decrecimiento sostenido así como el hurto de vehículo en menor medida. De igual forma el delito del hurto genérico ha disminuído los últimos 4 meses a comparación de los primeros 5 meses del año. El resto de los delitos han mantenido un valor relativamente constante en el tiempo. Quizás la crisis haya generado la desincorporación de vehículos por mal estado.
Denuncias Realizadas por día de semana mediante grafico de radar
conteo_dia <- delitos %>%
count(delito, diasemana = wday(fecha_hora_anglo, label = TRUE))
ggplot(conteo_dia, aes(x = diasemana, y = n, col = n, group = delito)) +
guides(color = "NA") +
geom_polygon(fill = "NA", size = 2, color = "red") +
coord_polar() +
facet_wrap(~delito)+
labs(title = "Total Delitos Denunciados",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
fill = "Tipo de Delito") +
expand_limits(y = -2000) +
xlab("")+
ylab("")+
geom_text(size=2,
aes(label = n),
color="black")
Mientras que los delitos hurto, robo de vehículos y robo genérico descansan los fines de semana con respeto a los días lunes a viernes, el resto de los delitos tienen un comportamiento similar todos los días de la semana.
Denuncias Realizadas por hora mediante Poligono de Frecuencia
# Se realiza primeramente un conteo de denuncias por mes
conteo_dia <- delitos %>%
count(delito, hora = hour(fecha_hora_anglo))
# Luego se muestran las cantidades mensuales como polígono de frecuencia
ggplot(conteo_dia) +
geom_line(aes(x = hora, y = n, group = delito, color = delito), size = 2) +
facet_wrap(~delito)+
guides(color = FALSE) +
geom_text(size=2,
aes(
x = hora,
y = n,
label = n))+
labs(title = "Total Delitos Denunciados",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
fill = "Tipo de Delito",
x = "Horas",
y = "Cantidad") +
expand_limits(y = 0)
A pesar que los missing quedaron registrados a las 00:00:00 y obviando ese horario se puede ver claramente que los delitos de hurto, hurto de vehículo, robo y lesiones personales tienden a aumentar en el horario comprendido de las 13:00:00 hasta las 20:00:00 mientras que el resto de los delitos tienen un comportamiento bastante similar durante todas las horas del día.
Mapa Base
bbox <- c(min(delitos$long_geom, na.rm = TRUE),
min(delitos$lat_geom, na.rm = TRUE),
max(delitos$long_geom, na.rm = TRUE),
max(delitos$lat_geom, na.rm = TRUE))
venezuela <- get_stamenmap(bbox, maptype = "toner-lite")
Mapa de Venezuela con Georeferenciación Delictiva por Tipo
ggplot() +
geom_sf(data = estado, color = "gray39", fill = "white") +
geom_point(data = filter(delitos, delito %in% delitos_tipo), aes(x = long_geom, y = lat_geom, color = delito), size = 0.1, alpha = 0.1) +
guides(color = guide_legend(override.aes = list(size=2, alpha = 1))) +
scale_fill_distiller(palette = "Spectral") +
labs(title = "Ubicacion de los 58956 Delitos Denunciados por Tipo",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
x = "Longitud",
y = "Latitud") +
annotation_scale() +
annotation_north_arrow(location='tr')
Mapa de Venezuela con Georeferenciación Delictiva por Categoria
ggplot() +
geom_sf(data = estado, color = "gray39", fill = "white") +
geom_point(data = filter(delitos, delito %in% delitos_tipo), aes(x = long_geom, y = lat_geom, color = delito), size = 0.1, alpha = 0.1) +
guides(color = guide_legend(override.aes = list(size=2, alpha = 1))) +
scale_fill_distiller(palette = "Spectral") +
labs(title = "Ubicacion de los 58956 Delitos Denunciados por Categoria",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
x = "Longitud",
y = "Latitud") +
facet_wrap(~delito, nrow = 3)
Combinando Horario y Tipo de Delito (Homicidios por Horas del Día)
delitos <- delitos %>%
mutate(hora_base = hour(ymd_hms(fecha_hora_anglo)))
ggplot() +
geom_sf(data = estado, color = "gray39", fill = "white") +
geom_point(data = filter(delitos, delito == "HOMICIDIO"), aes(x = long_geom, y = lat_geom, color = delito), size = 0.1, alpha = 0.1) +
guides(color = guide_legend(override.aes = list(size=2, alpha = 1))) +
scale_fill_distiller(palette = "Spectral") +
labs(title = "Ubicacion de los Homicidios por Horas del Día",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
x = "Longitud",
y = "Latitud") +
facet_wrap(~hora_base, nrow = 4)
Mapa de Densidad de Homicidios por Horas del Día
ggplot() +
geom_sf(data = estado, color = "gray39", fill = "white") +
stat_density_2d(data = filter(delitos, delito == "HOMICIDIO"), aes(x = long_geom, y = lat_geom, fill = stat(level)), alpha = 0.65, geom = "polygon") +
scale_fill_distiller(palette = "Spectral")+
facet_wrap(~hora_base, nrow = 4) +
labs(title = "Densidad de los Homicidios por Horas del Día",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
x = "Longitud",
y = "Latitud",
fill = "Densidad Delictiva")
Mapa de Calor de Delitos
ggmap(venezuela) +
geom_bin2d(data = delitos, aes(x = long_geom, y = lat_geom), bins = 90, alpha = 0.8) +
scale_fill_distiller(palette = "Spectral") +
labs(title = "Mapa de Calor Delitos Denunciados",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
x = "Longitud",
y = "Latitud",
fill = "Concentración Delictiva")
## Warning: Removed 2 rows containing missing values (geom_tile).
Mapa de Densidad de Delitos
ggplot() +
geom_sf(data = estado, color = "gray39", fill = "white") +
stat_density_2d(data = delitos, aes(x = long_geom, y = lat_geom, fill = stat(level)), alpha = 0.65, geom = "polygon") +
scale_fill_distiller(palette = "Spectral")+
labs(title = "Densidad de los Delitos Denunciados",
subtitle = "Enero - Septiembre 2019",
caption = "Fuente: Extraoficial - Venezuela",
x = "Longitud",
y = "Latitud",
fill = "Densidad Delictiva") +
annotation_scale() +
annotation_north_arrow(location='tr')
El mapa muestra lo que previamente en el histograma se había analizado que los lugares (las capitales) en los estados centrales del país Miranda, Aragua, Carabobo y Distrito Capital así como para el occidente Lara y Zulia y el Oriente-Sur Anzoátegui y Bolívar son los que poseen mayor cantidad de denuncias.