Lectura del mapa de municipios de Antioquia.
require(spdep)
## Loading required package: spdep
## Loading required package: sp
## Loading required package: Matrix
require(maps)
## Loading required package: maps
##
## # maps v3.1: updated 'world': all lakes moved to separate new #
## # 'lakes' database. Type '?world' or 'news(package="maps")'. #
require(maptools)
## Loading required package: maptools
## Checking rgeos availability: FALSE
## Note: when rgeos is not available, polygon geometry computations in maptools depend on gpclib,
## which has a restricted licence. It is disabled by default;
## to enable gpclib, type gpclibPermit()
require(classInt)
## Loading required package: classInt
require(RColorBrewer)
## Loading required package: RColorBrewer
require(data.table)
## Loading required package: data.table
require(rgdal)
## Loading required package: rgdal
## rgdal: version: 1.1-8, (SVN revision 616)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 1.11.2, released 2015/02/10
## Path to GDAL shared files: /usr/share/gdal/1.11
## Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
## Path to PROJ.4 shared files: (autodetected)
## Linking to sp version: 1.2-3
# Con el paquete "maptools"
muniantioq1 <- readShapeSpatial("muniantioq.shp")
# Otra forma de leer el mapa de municipios de Antioquia con el paqueta "rgdal""
muniantioq2 <- readOGR(".","muniantioq",encoding="latin1")
## OGR data source with driver: ESRI Shapefile
## Source: ".", layer: "muniantioq"
## with 125 features
## It has 12 fields
# Cualquiera de los dos mapas puede ser utilizado, pero se recomienda
# utilizar el "readOGR()" porque lee directamente el sistema de coordenadas
# si el mapa original tiene sistema de coordenadas.
muniantioq <- muniantioq1
Despliegue del mapa de municipios de Antioquia
plot(muniantioq)
Manejo de la estructura SpatialPolygonsDataFrame
summary(muniantioq)
## Object of class SpatialPolygonsDataFrame
## Coordinates:
## min max
## x 663198 1023924
## y 1096917 1475523
## Is projected: NA
## proj4string : [NA]
## Data attributes:
## AREA PERIMETER MUNICIP_ MUNICIP_ID
## Min. :1.622e+07 Min. : 17151 Min. :165.0 Min. :406.0
## 1st Qu.:1.491e+08 1st Qu.: 57851 1st Qu.:297.0 1st Qu.:619.0
## Median :2.676e+08 Median : 75085 Median :373.0 Median :690.0
## Mean :5.110e+08 Mean :102330 Mean :369.2 Mean :685.4
## 3rd Qu.:5.963e+08 3rd Qu.:122722 3rd Qu.:443.0 3rd Qu.:760.0
## Max. :3.121e+09 Max. :524575 Max. :563.0 Max. :837.0
##
## MUNICIPIO TIPOCAB NOMB_CAB DEPTO
## ABEJORRAL : 1 Min. :2601 ABEJORRAL : 1 ANTIOQUIA:125
## ABRIAQUI : 1 1st Qu.:2601 ABRIAQUI : 1
## ALEJANDRIA: 1 Median :2601 ALEJANDRIA: 1
## AMAGA : 1 Mean :2601 AMAGA : 1
## AMALFI : 1 3rd Qu.:2601 AMALFI : 1
## ANDES : 1 Max. :2605 ANDES : 1
## (Other) :119 (Other) :119
## CODIGO COD_DANE CUENACTUAL JURISDICCI
## Min. :5001 5001 : 1 Min. :0.00 Min. :1
## 1st Qu.:5147 5002 : 1 1st Qu.:1.00 1st Qu.:1
## Median :5376 5004 : 1 Median :1.00 Median :1
## Mean :5416 5021 : 1 Mean :0.88 Mean :1
## 3rd Qu.:5659 5030 : 1 3rd Qu.:1.00 3rd Qu.:1
## Max. :5895 5031 : 1 Max. :1.00 Max. :1
## (Other):119 NA's :120
class(muniantioq)
## [1] "SpatialPolygonsDataFrame"
## attr(,"package")
## [1] "sp"
str(muniantioq, max.level=2)
## Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots
## ..@ data :'data.frame': 125 obs. of 12 variables:
## .. ..- attr(*, "data_types")= chr [1:12] "N" "N" "N" "N" ...
## ..@ polygons :List of 125
## .. .. [list output truncated]
## ..@ plotOrder : int [1:125] 5 56 19 9 23 13 20 25 41 10 ...
## ..@ bbox : num [1:2, 1:2] 663198 1096917 1023924 1475523
## .. ..- attr(*, "dimnames")=List of 2
## ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
bbox(muniantioq)
## min max
## x 663198 1023924
## y 1096917 1475523
Sistema de coordenadas para el mapa de antioquia
proj4string(muniantioq)
## [1] NA
# Sitio web para obtener parámetros de la proyección
# http://spatialreference.org/
proj4string(muniantioq) <- CRS("+proj=tmerc +lat_0=4.599047222222222 +lon_0=-74.08091666666667 +k=1 +x_0=1000000 +y_0=1000000 +ellps=intl +towgs84=307,304,-318,0,0,0,0 +units=m +no_defs")
proj4string(muniantioq)
## [1] "+proj=tmerc +lat_0=4.599047222222222 +lon_0=-74.08091666666667 +k=1 +x_0=1000000 +y_0=1000000 +ellps=intl +towgs84=307,304,-318,0,0,0,0 +units=m +no_defs"
Manejo de la estrucutra espacial SpatialPolygonsDataFrame
summary(muniantioq)
## Object of class SpatialPolygonsDataFrame
## Coordinates:
## min max
## x 663198 1023924
## y 1096917 1475523
## Is projected: TRUE
## proj4string :
## [+proj=tmerc +lat_0=4.599047222222222 +lon_0=-74.08091666666667
## +k=1 +x_0=1000000 +y_0=1000000 +ellps=intl
## +towgs84=307,304,-318,0,0,0,0 +units=m +no_defs]
## Data attributes:
## AREA PERIMETER MUNICIP_ MUNICIP_ID
## Min. :1.622e+07 Min. : 17151 Min. :165.0 Min. :406.0
## 1st Qu.:1.491e+08 1st Qu.: 57851 1st Qu.:297.0 1st Qu.:619.0
## Median :2.676e+08 Median : 75085 Median :373.0 Median :690.0
## Mean :5.110e+08 Mean :102330 Mean :369.2 Mean :685.4
## 3rd Qu.:5.963e+08 3rd Qu.:122722 3rd Qu.:443.0 3rd Qu.:760.0
## Max. :3.121e+09 Max. :524575 Max. :563.0 Max. :837.0
##
## MUNICIPIO TIPOCAB NOMB_CAB DEPTO
## ABEJORRAL : 1 Min. :2601 ABEJORRAL : 1 ANTIOQUIA:125
## ABRIAQUI : 1 1st Qu.:2601 ABRIAQUI : 1
## ALEJANDRIA: 1 Median :2601 ALEJANDRIA: 1
## AMAGA : 1 Mean :2601 AMAGA : 1
## AMALFI : 1 3rd Qu.:2601 AMALFI : 1
## ANDES : 1 Max. :2605 ANDES : 1
## (Other) :119 (Other) :119
## CODIGO COD_DANE CUENACTUAL JURISDICCI
## Min. :5001 5001 : 1 Min. :0.00 Min. :1
## 1st Qu.:5147 5002 : 1 1st Qu.:1.00 1st Qu.:1
## Median :5376 5004 : 1 Median :1.00 Median :1
## Mean :5416 5021 : 1 Mean :0.88 Mean :1
## 3rd Qu.:5659 5030 : 1 3rd Qu.:1.00 3rd Qu.:1
## Max. :5895 5031 : 1 Max. :1.00 Max. :1
## (Other):119 NA's :120
slotNames(muniantioq)
## [1] "data" "polygons" "plotOrder" "bbox" "proj4string"
class(muniantioq@data)
## [1] "data.frame"
class(muniantioq@polygons)
## [1] "list"
slotNames(muniantioq@polygons[[1]])
## [1] "Polygons" "plotOrder" "labpt" "ID" "area"
datos_muniantioq <- data.table(muniantioq@data)
# Nombres de municipios en formato de windows
muniantioq@data$MUNICIPIO
## [1] ARBOLETES SAN JUAN DE URABA NECOCLI
## [4] SAN PEDRO DE URABA TURBO NECHI
## [7] APARTADO CAUCASIA CACERES
## [10] EL BAGRE CAREPA CHIGORODO
## [13] TARAZA ZARAGOZA MUTATA
## [16] ANORI SEGOVIA VALDIVIA
## [19] ITUANGO YONDO AMALFI
## [22] YARUMAL REMEDIOS BRICE\xd10
## [25] DABEIBA PEQUE CAMPAMENTO
## [28] TOLEDO SABANALARGA MURINDO
## [31] SAN ANDRES URAMITA CA\xd1ASGORDAS
## [34] ANGOSTURA BURITICA GUADALUPE
## [37] VEGACHI SAN JOSE DE LA MONTANA FRONTINO
## [40] SANTA ROSA DE OSOS VIGIA DEL FUERTE GOMEZ PLATA
## [43] YOLOMBO CAROLINA LIBORINA
## [46] YALI BELMIRA ABRIAQUI
## [49] PUERTO BERRIO GIRALDO ENTRERRIOS
## [52] SANTA FE DE ANTIOQUIA OLAYA MACEO
## [55] SOPETRAN URRAO CISNEROS
## [58] DON MATIAS SANTO DOMINGO SAN ROQUE
## [61] BARBOSA SAN PEDRO SAN JERONIMO
## [64] CAICEDO GIRARDOTA CONCEPCION
## [67] CARACOLI BELLO ALEJANDRIA
## [70] COPACABANA EBEJICO ANZA
## [73] SAN RAFAEL SAN VICENTE MEDELLIN
## [76] EL PE\xd1OL SAN CARLOS GUARNE
## [79] GUATAPE PUERTO NARE BETULIA
## [82] HELICONIA MARINILLA ARMENIA
## [85] RIONEGRO GRANADA ITAGUI
## [88] EL RETIRO ANGELOPOLIS ENVIGADO
## [91] LA ESTRELLA CONCORDIA SANTUARIO
## [94] SABANETA CARMEN DE VIBORAL SAN LUIS
## [97] CALDAS TITIRIBI COCORNA
## [100] SALGAR AMAGA LA UNION
## [103] LA CEJA PUERTO TRIUNFO VENECIA
## [106] SANTA BARBARA SAN FRANCISCO FREDONIA
## [109] MONTEBELLO ABEJORRAL TARSO
## [112] SONSON CIUDAD BOLIVAR PUEBLO RICO
## [115] JERICO HISPANIA LA PINTADA
## [118] TAMESIS BETANIA ANDES
## [121] ARGELIA VALPARAISO NARI\xd1O
## [124] JARDIN CARAMANTA
## 125 Levels: ABEJORRAL ABRIAQUI ALEJANDRIA AMAGA AMALFI ... ZARAGOZA
# Convertir a formato de linux
muniantioq@data$MUNICIPIO <- factor(iconv(datos_muniantioq$MUNICIPIO,
"latin1", "UTF-8"))
muniantioq@data$MUNICIPIO
## [1] ARBOLETES SAN JUAN DE URABA NECOCLI
## [4] SAN PEDRO DE URABA TURBO NECHI
## [7] APARTADO CAUCASIA CACERES
## [10] EL BAGRE CAREPA CHIGORODO
## [13] TARAZA ZARAGOZA MUTATA
## [16] ANORI SEGOVIA VALDIVIA
## [19] ITUANGO YONDO AMALFI
## [22] YARUMAL REMEDIOS BRICEÑ0
## [25] DABEIBA PEQUE CAMPAMENTO
## [28] TOLEDO SABANALARGA MURINDO
## [31] SAN ANDRES URAMITA CAÑASGORDAS
## [34] ANGOSTURA BURITICA GUADALUPE
## [37] VEGACHI SAN JOSE DE LA MONTANA FRONTINO
## [40] SANTA ROSA DE OSOS VIGIA DEL FUERTE GOMEZ PLATA
## [43] YOLOMBO CAROLINA LIBORINA
## [46] YALI BELMIRA ABRIAQUI
## [49] PUERTO BERRIO GIRALDO ENTRERRIOS
## [52] SANTA FE DE ANTIOQUIA OLAYA MACEO
## [55] SOPETRAN URRAO CISNEROS
## [58] DON MATIAS SANTO DOMINGO SAN ROQUE
## [61] BARBOSA SAN PEDRO SAN JERONIMO
## [64] CAICEDO GIRARDOTA CONCEPCION
## [67] CARACOLI BELLO ALEJANDRIA
## [70] COPACABANA EBEJICO ANZA
## [73] SAN RAFAEL SAN VICENTE MEDELLIN
## [76] EL PEÑOL SAN CARLOS GUARNE
## [79] GUATAPE PUERTO NARE BETULIA
## [82] HELICONIA MARINILLA ARMENIA
## [85] RIONEGRO GRANADA ITAGUI
## [88] EL RETIRO ANGELOPOLIS ENVIGADO
## [91] LA ESTRELLA CONCORDIA SANTUARIO
## [94] SABANETA CARMEN DE VIBORAL SAN LUIS
## [97] CALDAS TITIRIBI COCORNA
## [100] SALGAR AMAGA LA UNION
## [103] LA CEJA PUERTO TRIUNFO VENECIA
## [106] SANTA BARBARA SAN FRANCISCO FREDONIA
## [109] MONTEBELLO ABEJORRAL TARSO
## [112] SONSON CIUDAD BOLIVAR PUEBLO RICO
## [115] JERICO HISPANIA LA PINTADA
## [118] TAMESIS BETANIA ANDES
## [121] ARGELIA VALPARAISO NARIÑO
## [124] JARDIN CARAMANTA
## 125 Levels: ABEJORRAL ABRIAQUI ALEJANDRIA AMAGA AMALFI ... ZARAGOZA
muniantioq$rn <- row.names(muniantioq)
NBI <- data.table(read.csv2("NBI.csv"))
head(NBI)
## Cod.Depto Nombre.Departamento Cod.Mpio Nombre.Municipio
## 1: 5 ANTIOQUIA 1 MEDELLIN
## 2: 5 ANTIOQUIA 2 ABEJORRAL
## 3: 5 ANTIOQUIA 4 ABRIAQUI
## 4: 5 ANTIOQUIA 21 ALEJANDRIA
## 5: 5 ANTIOQUIA 30 AMAGA
## 6: 5 ANTIOQUIA 31 AMALFI
## Cod.MpioCompleto NBI Nota.Aclaratoria
## 1: 5001 12.42
## 2: 5002 33.85
## 3: 5004 28.54
## 4: 5021 33.14
## 5: 5030 27.07
## 6: 5031 45.29
setkey(NBI,Cod.MpioCompleto)
muniantioqBD <- data.table(muniantioq@data)
setkey(muniantioqBD, CODIGO)
temp1 <- NBI[muniantioqBD]
setkey(temp1, rn)
muniantioq@data <- temp1[row.names(muniantioq)]
Despligue de colores de acuerdo a la variable NBI
escalagrises <- muniantioq@data$NBI/100
plot(muniantioq, col=gray(escalagrises))
escalagrises <- 1-muniantioq@data$NBI/100
plot(muniantioq, col=gray(escalagrises))
require(RColorBrewer)
rojo_azul <- brewer.pal(11,"RdBu")
categorias <- as.integer(cut(1-muniantioq@data$NBI/100, seq(0,1,l=11)))
plot(muniantioq, col=rojo_azul[categorias])
display.brewer.all()
am_nar_roj <- brewer.pal(9,"YlOrRd")
categorias <- as.integer(cut(muniantioq@data$NBI/100, seq(0,1,l=9)))
plot(muniantioq, col=am_nar_roj[categorias])
azul_rojo <- colorRampPalette(c("blue", "red"))(100)
categorias <- as.integer(cut(muniantioq@data$NBI/100, seq(0,1,l=101)))
plot(muniantioq, col=azul_rojo[categorias])
categorias <- as.integer(cut(muniantioq@data$NBI/100, seq(0,1,l=101)))
plot(muniantioq, col=terrain.colors(100)[categorias])
require(sp)
spplot(muniantioq, c("NBI"))
l2 <- list("SpatialPolygonsRescale",
layout.north.arrow(),
offset = c(1000000,1120000),
scale = 50000)
spplot(muniantioq, c("NBI"), sp.layout=list(l2))
Creación de la matriz de adjacencia
vecinosAnt <- poly2nb(muniantioq)
vecinosAnt
## Neighbour list object:
## Number of regions: 125
## Number of nonzero links: 678
## Percentage nonzero weights: 4.3392
## Average number of links: 5.424
summary(vecinosAnt)
## Neighbour list object:
## Number of regions: 125
## Number of nonzero links: 678
## Percentage nonzero weights: 4.3392
## Average number of links: 5.424
## Link number distribution:
##
## 2 3 4 5 6 7 8 9 10 12 13
## 9 9 26 24 22 20 7 4 2 1 1
## 9 least connected regions:
## 1 3 5 6 19 56 120 122 124 with 2 links
## 1 most connected region:
## 74 with 13 links
str(vecinosAnt)
## List of 125
## $ : int [1:4] 2 3 4 5
## $ : int [1:2] 1 3
## $ : int [1:3] 1 2 5
## $ : int [1:2] 1 5
## $ : int [1:7] 1 3 4 7 11 12 15
## $ : int [1:2] 8 10
## $ : int [1:2] 5 11
## $ : int [1:4] 6 9 10 14
## $ : int [1:4] 8 13 14 16
## $ : int [1:4] 6 8 14 17
## $ : int [1:3] 5 7 12
## $ : int [1:3] 5 11 15
## $ : int [1:4] 9 16 18 19
## $ : int [1:6] 8 9 10 16 17 21
## $ : int [1:4] 5 12 19 25
## $ : int [1:10] 9 13 14 17 18 21 22 27 34 36
## $ : int [1:5] 10 14 16 21 23
## $ : int [1:5] 13 16 19 22 24
## $ : int [1:8] 13 15 18 24 25 26 28 29
## $ : int [1:2] 23 49
## $ : int [1:9] 14 16 17 23 36 37 42 43 46
## $ : int [1:8] 16 18 24 27 28 31 34 40
## $ : int [1:7] 17 20 21 37 43 46 49
## $ : int [1:5] 18 19 22 28 31
## $ : int [1:7] 15 19 26 30 32 33 39
## $ : int [1:6] 19 25 29 32 33 35
## $ : int [1:4] 16 22 34 36
## $ : int [1:5] 19 22 24 29 31
## $ : int [1:6] 19 26 28 31 35 45
## $ : int [1:3] 25 39 41
## $ : int [1:7] 22 24 28 29 38 40 45
## $ : int [1:4] 25 26 33 39
## $ : int [1:7] 25 26 32 35 39 48 50
## $ : int [1:6] 16 22 27 36 40 44
## $ : int [1:7] 26 29 33 45 50 52 53
## $ : int [1:6] 16 21 27 34 42 44
## $ : int [1:3] 21 23 46
## $ : int [1:4] 31 40 45 47
## $ : int [1:7] 25 30 32 33 41 48 56
## $ : int [1:12] 22 31 34 38 42 43 44 47 51 58 ...
## $ : int [1:3] 30 39 56
## $ : int [1:5] 21 36 40 43 44
## $ : int [1:10] 21 23 40 42 46 49 54 57 59 60
## $ : int [1:4] 34 36 40 42
## $ : int [1:7] 29 31 35 38 47 52 53
## $ : int [1:4] 21 23 37 43
## $ : int [1:8] 38 40 45 51 53 55 62 63
## $ : int [1:6] 33 39 50 52 56 64
## $ : int [1:6] 20 23 43 54 67 80
## $ : int [1:4] 33 35 48 52
## $ : int [1:4] 40 47 58 62
## $ : int [1:9] 35 45 48 50 53 55 64 71 72
## $ : int [1:5] 35 45 47 52 55
## $ : int [1:4] 43 49 60 67
## $ : int [1:5] 47 52 53 63 71
## $ : int [1:7] 39 41 48 64 72 81 100
## $ : int [1:2] 43 59
## $ : int [1:6] 40 51 59 61 62 65
## $ : int [1:8] 40 43 57 58 60 61 66 69
## $ : int [1:7] 43 54 59 67 69 73 77
## $ : int [1:5] 58 59 65 66 74
## $ : int [1:9] 40 47 51 58 63 65 68 70 75
## $ : int [1:6] 47 55 62 68 71 75
## $ : int [1:4] 48 52 56 72
## $ : int [1:6] 58 61 62 70 74 78
## $ : int [1:5] 59 61 69 74 76
## $ : int [1:5] 49 54 60 77 80
## $ : int [1:4] 62 63 70 75
## $ : int [1:6] 59 60 66 73 76 79
## $ : int [1:5] 62 65 68 75 78
## $ : int [1:8] 52 55 63 72 75 81 82 84
## $ : int [1:5] 52 56 64 71 81
## $ : int [1:5] 60 69 77 79 86
## $ : int [1:7] 61 65 66 76 78 83 85
## $ : int [1:13] 62 63 68 70 71 78 82 85 87 88 ...
## $ : int [1:7] 66 69 74 79 83 86 93
## $ : int [1:7] 60 67 73 79 80 86 96
## $ : int [1:5] 65 70 74 75 85
## $ : int [1:5] 69 73 76 77 86
## $ : int [1:5] 49 67 77 96 104
## $ : int [1:6] 56 71 72 84 92 100
## $ : int [1:4] 71 75 84 89
## $ : int [1:5] 74 76 85 93 95
## $ : int [1:6] 71 81 82 89 92 98
## $ : int [1:7] 74 75 78 83 88 95 103
## $ : int [1:7] 73 76 77 79 93 96 99
## $ : int [1:4] 75 90 91 94
## $ : int [1:7] 75 85 90 97 103 106 109
## $ : int [1:7] 75 82 84 91 97 98 101
## $ : int [1:6] 75 87 88 91 94 97
## $ : int [1:6] 75 87 89 90 94 97
## $ : int [1:5] 81 84 98 100 105
## $ : int [1:5] 76 83 86 95 99
## $ : int [1:4] 87 90 91 97
## $ : int [1:7] 83 85 93 99 102 103 112
## $ : int [1:7] 77 80 86 99 104 107 112
## $ : int [1:8] 88 89 90 91 94 101 106 108
## $ : int [1:5] 84 89 92 101 105
## $ : int [1:6] 86 93 95 96 107 112
## [list output truncated]
## - attr(*, "class")= chr "nb"
## - attr(*, "region.id")= chr [1:125] "0" "1" "2" "3" ...
## - attr(*, "call")= language poly2nb(pl = muniantioq)
## - attr(*, "type")= chr "queen"
## - attr(*, "sym")= logi TRUE
subset(muniantioq@data, muniantioq@data$rn == "1")
## Cod.Depto Nombre.Departamento Cod.Mpio Nombre.Municipio
## 1: 5 ANTIOQUIA 659 SAN JUAN DE URABA
## Cod.MpioCompleto NBI Nota.Aclaratoria AREA PERIMETER MUNICIP_
## 1: 5659 75.86 253573453 71762.33 169
## MUNICIP_ID MUNICIPIO TIPOCAB NOMB_CAB DEPTO
## 1: 815 SAN JUAN DE URABA 2601 SAN JUAN DE URABA ANTIOQUIA
## COD_DANE CUENACTUAL JURISDICCI rn
## 1: 5659 0 NA 1
subset(muniantioq@data, muniantioq@data$rn == "124")
## Cod.Depto Nombre.Departamento Cod.Mpio Nombre.Municipio
## 1: 5 ANTIOQUIA 145 CARAMANTA
## Cod.MpioCompleto NBI Nota.Aclaratoria AREA PERIMETER MUNICIP_
## 1: 5145 25 71974707 38707.89 563
## MUNICIP_ID MUNICIPIO TIPOCAB NOMB_CAB DEPTO COD_DANE CUENACTUAL
## 1: 698 CARAMANTA 2601 CARAMANTA ANTIOQUIA 5145 1
## JURISDICCI rn
## 1: NA 124
subset(muniantioq@data, muniantioq@data$rn == "74")
## Cod.Depto Nombre.Departamento Cod.Mpio Nombre.Municipio
## 1: 5 ANTIOQUIA 1 MEDELLIN
## Cod.MpioCompleto NBI Nota.Aclaratoria AREA PERIMETER MUNICIP_
## 1: 5001 12.42 382964071 92409.79 398
## MUNICIP_ID MUNICIPIO TIPOCAB NOMB_CAB DEPTO COD_DANE CUENACTUAL
## 1: 676 MEDELLIN 2605 MEDELLIN ANTIOQUIA 5001 1
## JURISDICCI rn
## 1: NA 74
matrizAdjAnt <- nb2mat(vecinosAnt, style="B")