Spatial Data Analysis Using R - Collaborative Project: Nearshoring
Team 5: Gabriela Monforte, David Saucedo, Samuel Montenegro, and Raúl Cantú
2023-01-11
Professor: Dohyeong Kim, Ph.D
Project Overview
According to an analysis by Morgan Stanley, the trend of relocating companies near their target markets, known as nearshoring, will prop up Mexican exports with an estimated 94.2 billion additional dollars in five years
The objective of this project is to perform an Exploratory Spatial Data Analysis (ESDA) to
Descriptive statistics of variables of interest
Importing data base (db)
# file.choose()
db <- read.csv("F:\\TEC\\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\\Proyecto\\database_nearshoring\\nearshoring.csv")Understanding db
summary (db)## id state region land_area
## Min. : 888 Length:32 Length:32 Min. : 1499
## 1st Qu.:1047 Class :character Class :character 1st Qu.: 23861
## Median :1081 Mode :character Mode :character Median : 58978
## Mean :1219 Mean : 60038
## 3rd Qu.:1118 3rd Qu.: 73730
## Max. :2357 Max. :247087
## border_distance fdi_2020 fdi_2021 fdi_2022
## Min. : 8.83 Min. : 85.0 Min. : 69.00 Min. : 61.00
## 1st Qu.: 613.26 1st Qu.: 112.8 1st Qu.: 97.75 1st Qu.: 78.25
## Median : 751.64 Median : 182.5 Median : 159.50 Median : 134.00
## Mean : 704.92 Mean : 264.7 Mean : 252.47 Mean : 197.50
## 3rd Qu.: 875.76 3rd Qu.: 324.5 3rd Qu.: 300.50 3rd Qu.: 253.00
## Max. :1252.66 Max. :1429.0 Max. :1542.00 Max. :1020.00
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## Min. : 944 Min. : 10.36 Min. : 10.69 Min. : 4.994
## 1st Qu.: 1295 1st Qu.: 37.49 1st Qu.: 37.92 1st Qu.: 18.063
## Median : 2142 Median : 64.91 Median : 66.36 Median : 26.793
## Mean : 3544 Mean : 306.43 Mean : 308.31 Mean : 143.559
## 3rd Qu.: 3454 3rd Qu.: 173.08 3rd Qu.: 174.65 3rd Qu.: 68.429
## Max. :30639 Max. :6135.57 Max. :6145.37 Max. :2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## Min. : 5.478 Min. :0.00100 Min. :0.00100 Min. :20.52
## 1st Qu.: 18.936 1st Qu.:0.00875 1st Qu.:0.00775 1st Qu.:20.52
## Median : 28.520 Median :0.01700 Median :0.01600 Median :20.52
## Mean : 145.644 Mean :0.02497 Mean :0.02284 Mean :20.52
## 3rd Qu.: 69.803 3rd Qu.:0.03050 3rd Qu.:0.02650 3rd Qu.:20.52
## Max. :3020.262 Max. :0.09100 Max. :0.08800 Max. :20.52
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## Min. :19.18 Min. :253.7 Min. :250.0 Min. : 0.00
## 1st Qu.:19.18 1st Qu.:291.7 1st Qu.:286.3 1st Qu.: 0.00
## Median :19.18 Median :312.3 Median :310.3 Median : 0.00
## Mean :19.18 Mean :318.0 Mean :312.4 Mean : 24.53
## 3rd Qu.:19.18 3rd Qu.:336.1 3rd Qu.:333.5 3rd Qu.: 16.75
## Max. :19.18 Max. :441.9 Max. :428.4 Max. :204.00
## industrial_parks crime_rate_2018 crime_rate_2019
## Min. : 0.000 Min. : 1.99 Min. : 3.08
## 1st Qu.: 0.000 1st Qu.:13.46 1st Qu.:10.90
## Median : 0.000 Median :18.98 Median :19.25
## Mean : 8.781 Mean :29.96 Mean :29.95
## 3rd Qu.: 9.500 3rd Qu.:42.97 3rd Qu.:43.02
## Max. :51.000 Max. :97.19 Max. :91.69library(dplyr)##
## Attaching package: 'dplyr'## The following objects are masked from 'package:stats':
##
## filter, lag## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, unionhead(db)## id state region land_area border_distance fdi_2020
## 1 1057 Aguascalientes Centro-Norte 5589 625.59 204
## 2 2304 Baja California Norte 70113 8.83 564
## 3 2327 Baja California Sur Norte-Occidente 73677 800.32 154
## 4 1086 Campeche Sur 51833 978.33 85
## 5 1182 Chiapas Sur 73887 1111.82 107
## 6 888 Chihuahua Norte 247087 350.02 375
## fdi_2021 fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019
## 1 166 151 1820 246.541 253.853
## 2 553 494 5752 53.469 53.716
## 3 153 103 4577 10.363 10.685
## 4 69 67 1142 17.042 17.709
## 5 93 74 1213 72.414 74.052
## 6 384 311 3390 14.913 15.059
## jobdensity_2018 jobdensity_2019 patentsrate_2016 patentsrate_2017
## 1 103.343 104.754 0.073 0.036
## 2 24.594 25.108 0.017 0.020
## 3 4.994 5.478 0.005 0.010
## 4 7.467 7.954 0.017 0.016
## 5 26.168 28.252 0.007 0.002
## 6 6.705 6.733 0.019 0.016
## exchangerate_2018 exchangerate_2019 realwages_2016 realwages_2017
## 1 20.521 19.181 313.673 314.075
## 2 20.521 19.181 333.210 333.984
## 3 20.521 19.181 317.766 311.398
## 4 20.521 19.181 423.567 396.244
## 5 20.521 19.181 307.359 296.384
## 6 20.521 19.181 318.540 320.650
## automotive_industry industrial_parks crime_rate_2018 crime_rate_2019
## 1 14 0 6.77 5.85
## 2 7 51 73.34 77.30
## 3 0 0 12.07 8.61
## 4 0 0 8.61 7.06
## 5 0 0 11.33 9.60
## 6 19 32 71.65 91.69count(db, region, sort = TRUE)## region n
## 1 Centro 8
## 2 Sur 8
## 3 Norte 6
## 4 Centro-Norte 5
## 5 Norte-Occidente 5hist(db$land_area)
boxplot(db$land_area, horizontal = TRUE)
boxplot(db$land_area,plot=FALSE)$out## [1] 247087 151571 184934db %>% filter_all(any_vars(. %in% c(247087)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 888 Chihuahua Norte 247087 350.02 375 384 311
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 3390 14.913 15.059 6.705 6.733
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.019 0.016 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 318.54 320.65 19 32
## crime_rate_2018 crime_rate_2019
## 1 71.65 91.69db %>% filter_all(any_vars(. %in% c(151571)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 933 Coahuila Centro-Norte 151571 367.26 308 273
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 245 2547 20.166 20.612 8.603
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 8.929 0.032 0.037 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 335.838 333.337 55
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 31 8.93 7.27db %>% filter_all(any_vars(. %in% c(184934)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 2303 Sonora Norte 184934 251.32 246 220 187
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 2758 15.719 15.906 7.141 7.351
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.015 0.021 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 294.761 293.158 1 0
## crime_rate_2018 crime_rate_2019
## 1 46.47 52.8hist(db$border_distance)
boxplot(db$border_distance, horizontal = TRUE)
boxplot(db$border_distance,plot=FALSE)$out## [1] 8.83db %>% filter_all(any_vars(. %in% c(8.83)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 2304 Baja California Norte 70113 8.83 564 553
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 494 5752 53.469 53.716 24.594
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 25.108 0.017 0.02 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 333.21 333.984 7
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 51 73.34 77.3hist(db$fdi_2020)
boxplot(db$fdi_2020, horizontal = TRUE)
boxplot(db$fdi_2020,plot=FALSE)$out## [1] 1429db %>% filter_all(any_vars(. %in% c(1429)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44hist(db$fdi_2021)
boxplot(db$fdi_2021, horizontal = TRUE)
boxplot(db$fdi_2021,plot=FALSE)$out## [1] 1542db %>% filter_all(any_vars(. %in% c(1542)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44hist(db$fdi_1999_2022)
boxplot(db$fdi_1999_2022, horizontal = TRUE)
boxplot(db$fdi_1999_2022,plot=FALSE)$out## [1] 30639db %>% filter_all(any_vars(. %in% c(30639)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44hist(db$popdensity_2018)
boxplot(db$popdensity_2018, horizontal = TRUE)
boxplot(db$popdensity_2018,plot=FALSE)$out## [1] 6135.566 551.082 397.481db %>% filter_all(any_vars(. %in% c(6135.566)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44db %>% filter_all(any_vars(. %in% c(551.082)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1067 Estado de Mexico Centro 30589 647.48 428 399
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 323 6091 551.082 554.086 85.763
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 87.566 0.062 0.069 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 285.788 286.906 77
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 40 19.02 17.92db %>% filter_all(any_vars(. %in% c(397.481)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 1117 Morelos Centro 4941 795.6 126 99 95
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 1403 397.481 399.635 166.558 169.288
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.03 0.026 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 342.864 334.042 0 0
## crime_rate_2018 crime_rate_2019
## 1 53.18 49.87hist(db$popdensity_2019)
boxplot(db$popdensity_2019, horizontal = TRUE)
boxplot(db$popdensity_2019,plot=FALSE)$out## [1] 6145.374 554.086 399.635db %>% filter_all(any_vars(. %in% c(6145.374)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44db %>% filter_all(any_vars(. %in% c(554.086)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1067 Estado de Mexico Centro 30589 647.48 428 399
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 323 6091 551.082 554.086 85.763
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 87.566 0.062 0.069 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 285.788 286.906 77
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 40 19.02 17.92db %>% filter_all(any_vars(. %in% c(399.635)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 1117 Morelos Centro 4941 795.6 126 99 95
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 1403 397.481 399.635 166.558 169.288
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.03 0.026 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 342.864 334.042 0 0
## crime_rate_2018 crime_rate_2019
## 1 53.18 49.87hist(db$jobdensity_2018)
boxplot(db$jobdensity_2018, horizontal = TRUE)
boxplot(db$jobdensity_2018,plot=FALSE)$out## [1] 2990.353 340.456 166.558db %>% filter_all(any_vars(. %in% c(2990.353)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44db %>% filter_all(any_vars(. %in% c(340.456)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 1098 Jalisco Centro 21461 767.1 484 470 352
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 5919 383.722 389.537 340.456 340.224
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.018 0.016 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 337.018 328.755 10 19
## crime_rate_2018 crime_rate_2019
## 1 29.68 25.76db %>% filter_all(any_vars(. %in% c(166.558)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 1117 Morelos Centro 4941 795.6 126 99 95
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 1403 397.481 399.635 166.558 169.288
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.03 0.026 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 342.864 334.042 0 0
## crime_rate_2018 crime_rate_2019
## 1 53.18 49.87hist(db$jobdensity_2019)
boxplot(db$jobdensity_2019, horizontal = TRUE)
boxplot(db$jobdensity_2019,plot=FALSE)$out## [1] 3020.262 340.224 169.288db %>% filter_all(any_vars(. %in% c(3020.262)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44db %>% filter_all(any_vars(. %in% c(340.224)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 1098 Jalisco Centro 21461 767.1 484 470 352
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 5919 383.722 389.537 340.456 340.224
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.018 0.016 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 337.018 328.755 10 19
## crime_rate_2018 crime_rate_2019
## 1 29.68 25.76db %>% filter_all(any_vars(. %in% c(169.288)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 1117 Morelos Centro 4941 795.6 126 99 95
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 1403 397.481 399.635 166.558 169.288
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.03 0.026 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 342.864 334.042 0 0
## crime_rate_2018 crime_rate_2019
## 1 53.18 49.87hist(db$patentsrate_2016)
boxplot(db$patentsrate_2016, horizontal = TRUE)
boxplot(db$patentsrate_2016,plot=FALSE)$out## [1] 0.073 0.091 0.065db %>% filter_all(any_vars(. %in% c(0.073)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1057 Aguascalientes Centro-Norte 5589 625.59 204 166
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 151 1820 246.541 253.853 103.343
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 104.754 0.073 0.036 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 313.673 314.075 14
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 6.77 5.85db %>% filter_all(any_vars(. %in% c(0.091)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44db %>% filter_all(any_vars(. %in% c(0.065)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 976 Nuevo Leon Norte 64555 216.33 579 591 472
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 6511 85.521 88.557 38.326 40.043
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.065 0.064 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 381.792 375.977 165 49
## crime_rate_2018 crime_rate_2019
## 1 17.16 15.62hist(db$patentsrate_2017)
boxplot(db$patentsrate_2017, horizontal = TRUE)
boxplot(db$patentsrate_2017,plot=FALSE)$out## [1] 0.088 0.069 0.069 0.064 0.061db %>% filter_all(any_vars(. %in% c(0.088)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44db %>% filter_all(any_vars(. %in% c(0.069)))## id state region land_area border_distance fdi_2020
## 1 1067 Estado de Mexico Centro 30589 647.48 428
## 2 1051 Hidalgo Centro-Norte 80137 798.67 106
## fdi_2021 fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019
## 1 399 323 6091 551.082 554.086
## 2 97 79 1124 37.991 38.296
## jobdensity_2018 jobdensity_2019 patentsrate_2016 patentsrate_2017
## 1 85.763 87.566 0.062 0.069
## 2 45.976 47.676 0.060 0.069
## exchangerate_2018 exchangerate_2019 realwages_2016 realwages_2017
## 1 20.521 19.181 285.788 286.906
## 2 20.521 19.181 340.683 334.774
## automotive_industry industrial_parks crime_rate_2018 crime_rate_2019
## 1 77 40 19.02 17.92
## 2 0 0 15.58 11.76db %>% filter_all(any_vars(. %in% c(0.064)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 976 Nuevo Leon Norte 64555 216.33 579 591 472
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 6511 85.521 88.557 38.326 40.043
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.065 0.064 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 381.792 375.977 165 49
## crime_rate_2018 crime_rate_2019
## 1 17.16 15.62db %>% filter_all(any_vars(. %in% c(0.061)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 1069 Queretaro Centro 11769 652.31 391 333 273
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 3644 195.397 201.383 77.325 78.519
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.061 0.061 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 381.937 377.635 20 0
## crime_rate_2018 crime_rate_2019
## 1 9.32 9.42hist(db$realwages_2016)
boxplot(db$realwages_2016, horizontal = TRUE)
boxplot(db$realwages_2016,plot=FALSE)$out## [1] 423.567 441.882db %>% filter_all(any_vars(. %in% c(423.567)))## id state region land_area border_distance fdi_2020 fdi_2021 fdi_2022
## 1 1086 Campeche Sur 51833 978.33 85 69 67
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018 jobdensity_2019
## 1 1142 17.042 17.709 7.467 7.954
## patentsrate_2016 patentsrate_2017 exchangerate_2018 exchangerate_2019
## 1 0.017 0.016 20.521 19.181
## realwages_2016 realwages_2017 automotive_industry industrial_parks
## 1 423.567 396.244 0 0
## crime_rate_2018 crime_rate_2019
## 1 8.61 7.06db %>% filter_all(any_vars(. %in% c(441.882)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44hist(db$realwages_2017)
boxplot(db$realwages_2017, horizontal = TRUE)
boxplot(db$realwages_2017,plot=FALSE)$out## [1] 428.362db %>% filter_all(any_vars(. %in% c(428.362)))## id state region land_area border_distance fdi_2020 fdi_2021
## 1 1114 Ciudad de Mexico Norte 1499 738.2 1429 1542
## fdi_2022 fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## 1 1020 30639 6135.566 6145.374 2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## 1 3020.262 0.091 0.088 20.521
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## 1 19.181 441.882 428.362 95
## industrial_parks crime_rate_2018 crime_rate_2019
## 1 0 13.93 14.44Findings
Mexico City deviate markedly from other observations.
Choropleth Maps
library("mxmaps")
summary(df_mxstate_2020)## region state_name state_name_official state_abbr
## Length:32 Length:32 Length:32 Length:32
## Class :character Class :character Class :character Class :character
## Mode :character Mode :character Mode :character Mode :character
##
##
##
## state_abbr_official year pop pop_male
## Length:32 Min. :2020 Min. : 731391 Min. : 360622
## Class :character 1st Qu.:2020 1st Qu.: 1851651 1st Qu.: 928801
## Mode :character Median :2020 Median : 3054892 Median :1488097
## Mean :2020 Mean : 3937938 Mean :1921043
## 3rd Qu.:2020 3rd Qu.: 4947592 3rd Qu.:2406242
## Max. :2020 Max. :16992418 Max. :8251295
## pop_female afromexican indigenous_language
## Min. : 370769 Min. : 10416 Min. : 2539
## 1st Qu.: 926140 1st Qu.: 34747 1st Qu.: 30331
## Median :1557615 Median : 50479 Median : 73507
## Mean :2016895 Mean : 80507 Mean : 230145
## 3rd Qu.:2541349 3rd Qu.:100404 3rd Qu.: 264067
## Max. :8741123 Max. :303923 Max. :1459648df_mxstate_2020$value <- db$automotive_industry
mxstate_choropleth(df_mxstate_2020, title = "Automotive Supply Chain Companies Distribution in Mexico") 
Spatial Data Construction and Mapping
Install packages for spatial analysis
require(proj4) || install.packages("proj4", dependencies = T)## Loading required package: proj4## [1] TRUErequire(rgdal) || install.packages("rgdal", dependencies = T)## Loading required package: rgdal## Loading required package: sp## Please note that rgdal will be retired by the end of 2023,
## plan transition to sf/stars/terra functions using GDAL and PROJ
## at your earliest convenience.
##
## rgdal: version: 1.5-32, (SVN revision 1176)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 3.4.3, released 2022/04/22
## Path to GDAL shared files: C:/Users/raulc/AppData/Local/R/win-library/4.2/rgdal/gdal
## GDAL binary built with GEOS: TRUE
## Loaded PROJ runtime: Rel. 7.2.1, January 1st, 2021, [PJ_VERSION: 721]
## Path to PROJ shared files: C:/Users/raulc/AppData/Local/R/win-library/4.2/rgdal/proj
## PROJ CDN enabled: FALSE
## Linking to sp version:1.5-0
## To mute warnings of possible GDAL/OSR exportToProj4() degradation,
## use options("rgdal_show_exportToProj4_warnings"="none") before loading sp or rgdal.##
## Attaching package: 'rgdal'## The following object is masked from 'package:proj4':
##
## project## [1] TRUErequire(sp) || install.packages("sp", dependencies = T)## [1] TRUErequire(spdep) || install.packages("spdep", dependencies = T)## Loading required package: spdep## Loading required package: spData## To access larger datasets in this package, install the spDataLarge
## package with: `install.packages('spDataLarge',
## repos='https://nowosad.github.io/drat/', type='source')`## Loading required package: sf## Linking to GEOS 3.9.3, GDAL 3.5.2, PROJ 8.2.1; sf_use_s2() is TRUE## [1] TRUErequire(ape) || install.packages("ape", dependencies = T)## Loading required package: ape## Registered S3 method overwritten by 'ape':
## method from
## plot.mst spdep## [1] TRUErequire(gstat) || install.packages("gstat", dependencies = T)## Loading required package: gstat## [1] TRUErequire(MASS) || install.packages("MASS", dependencies = T)## Loading required package: MASS##
## Attaching package: 'MASS'## The following object is masked from 'package:dplyr':
##
## select## [1] TRUErequire(devtools) || install.packages("devtools")## Loading required package: devtools## Loading required package: usethis## [1] TRUErequire(spgwr) || install.packages("spgwr", dependencies = T)## Loading required package: spgwr## NOTE: This package does not constitute approval of GWR
## as a method of spatial analysis; see example(gwr)## [1] TRUErequire(spData) || install.packages("spData", dependencies = T)## [1] TRUErequire(Matrix) || install.packages("Matrix", dependencies = T)## Loading required package: Matrix## [1] TRUErequire(psych) || install.packages("psych", dependencies = T)## Loading required package: psych## [1] TRUErequire(tmap) || install.packages("tmap")## Loading required package: tmap## [1] TRUErequire(readr) || install.packages("readr")## Loading required package: readr## [1] TRUErequire(stargazer) || install.packages("stargazer")## Loading required package: stargazer##
## Please cite as:## Hlavac, Marek (2022). stargazer: Well-Formatted Regression and Summary Statistics Tables.## R package version 5.2.3. https://CRAN.R-project.org/package=stargazer## [1] TRUErequire(automap) || install.packages("automap")## Loading required package: automap## [1] TRUErequire(spatialEco) || install.packages("spatialEco")## Loading required package: spatialEco##
## Attaching package: 'spatialEco'## The following object is masked from 'package:dplyr':
##
## combine## [1] TRUErequire(SpatialKDE) || install.packages("SpatialKDE")## Loading required package: SpatialKDE## [1] TRUErequire(smacpod) || install.packages("smacpod")## Loading required package: smacpod## [1] TRUErequire(lme4) || install.packages("lme4")## Loading required package: lme4## [1] TRUEPrepare the data
setwd("F:\\TEC\\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\\Proyecto\\database_nearshoring")
getwd()## [1] "F:/TEC/2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R/Proyecto/database_nearshoring"library(maptools)## Checking rgeos availability: TRUE
## Please note that 'maptools' will be retired by the end of 2023,
## plan transition at your earliest convenience;
## some functionality will be moved to 'sp'.library(sp)
library(readr)
nearshoring_map2 <- readShapePoly("nearshoring_map.shp") ## Warning: readShapePoly is deprecated; use rgdal::readOGR or sf::st_readnearshoring.csv <- read_csv("nearshoring.csv") ## Rows: 32 Columns: 23## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (2): state, region
## dbl (21): id, land_area, border_distance, fdi_2020, fdi_2021, fdi_2022, fdi_...
##
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.nearshoring_map.sp <- merge(nearshoring_map2,nearshoring.csv)
save(nearshoring_map.sp, file="nearshoring_map.RData")
library(rgdal)
st <- readOGR("nearshoring_map.shp") ## OGR data source with driver: ESRI Shapefile
## Source: "F:\TEC\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\Proyecto\database_nearshoring\nearshoring_map.shp", layer: "nearshoring_map"
## with 32 features
## It has 30 fields
## Integer64 fields read as strings: OBJECTID id land_area fdi_2020 fdi_2021 fdi_2022 fdi_1999_2st1 <- data.frame(st)
library(sf)
cc = st_read("nearshoring_map.shp") ## Reading layer `nearshoring_map' from data source
## `F:\TEC\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\Proyecto\database_nearshoring\nearshoring_map.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 32 features and 30 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: -118.4042 ymin: 14.55055 xmax: -86.73862 ymax: 32.71846
## Geodetic CRS: WGS 84summary(db) ## id state region land_area
## Min. : 888 Length:32 Length:32 Min. : 1499
## 1st Qu.:1047 Class :character Class :character 1st Qu.: 23861
## Median :1081 Mode :character Mode :character Median : 58978
## Mean :1219 Mean : 60038
## 3rd Qu.:1118 3rd Qu.: 73730
## Max. :2357 Max. :247087
## border_distance fdi_2020 fdi_2021 fdi_2022
## Min. : 8.83 Min. : 85.0 Min. : 69.00 Min. : 61.00
## 1st Qu.: 613.26 1st Qu.: 112.8 1st Qu.: 97.75 1st Qu.: 78.25
## Median : 751.64 Median : 182.5 Median : 159.50 Median : 134.00
## Mean : 704.92 Mean : 264.7 Mean : 252.47 Mean : 197.50
## 3rd Qu.: 875.76 3rd Qu.: 324.5 3rd Qu.: 300.50 3rd Qu.: 253.00
## Max. :1252.66 Max. :1429.0 Max. :1542.00 Max. :1020.00
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## Min. : 944 Min. : 10.36 Min. : 10.69 Min. : 4.994
## 1st Qu.: 1295 1st Qu.: 37.49 1st Qu.: 37.92 1st Qu.: 18.063
## Median : 2142 Median : 64.91 Median : 66.36 Median : 26.793
## Mean : 3544 Mean : 306.43 Mean : 308.31 Mean : 143.559
## 3rd Qu.: 3454 3rd Qu.: 173.08 3rd Qu.: 174.65 3rd Qu.: 68.429
## Max. :30639 Max. :6135.57 Max. :6145.37 Max. :2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## Min. : 5.478 Min. :0.00100 Min. :0.00100 Min. :20.52
## 1st Qu.: 18.936 1st Qu.:0.00875 1st Qu.:0.00775 1st Qu.:20.52
## Median : 28.520 Median :0.01700 Median :0.01600 Median :20.52
## Mean : 145.644 Mean :0.02497 Mean :0.02284 Mean :20.52
## 3rd Qu.: 69.803 3rd Qu.:0.03050 3rd Qu.:0.02650 3rd Qu.:20.52
## Max. :3020.262 Max. :0.09100 Max. :0.08800 Max. :20.52
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## Min. :19.18 Min. :253.7 Min. :250.0 Min. : 0.00
## 1st Qu.:19.18 1st Qu.:291.7 1st Qu.:286.3 1st Qu.: 0.00
## Median :19.18 Median :312.3 Median :310.3 Median : 0.00
## Mean :19.18 Mean :318.0 Mean :312.4 Mean : 24.53
## 3rd Qu.:19.18 3rd Qu.:336.1 3rd Qu.:333.5 3rd Qu.: 16.75
## Max. :19.18 Max. :441.9 Max. :428.4 Max. :204.00
## industrial_parks crime_rate_2018 crime_rate_2019
## Min. : 0.000 Min. : 1.99 Min. : 3.08
## 1st Qu.: 0.000 1st Qu.:13.46 1st Qu.:10.90
## Median : 0.000 Median :18.98 Median :19.25
## Mean : 8.781 Mean :29.96 Mean :29.95
## 3rd Qu.: 9.500 3rd Qu.:42.97 3rd Qu.:43.02
## Max. :51.000 Max. :97.19 Max. :91.69summary(nearshoring_map.sp) ## Object of class SpatialPolygonsDataFrame
## Coordinates:
## min max
## x -118.40417 -86.73862
## y 14.55055 32.71846
## Is projected: NA
## proj4string : [NA]
## Data attributes:
## id state region land_area
## Min. : 888 Aguascalientes : 1 Centro :8 Min. : 1499
## 1st Qu.:1047 Baja California : 1 Centro-Norte :5 1st Qu.: 23861
## Median :1081 Baja California Sur: 1 Norte :6 Median : 58978
## Mean :1219 Campeche : 1 Norte-Occidente:5 Mean : 60038
## 3rd Qu.:1118 Chiapas : 1 Sur :8 3rd Qu.: 73730
## Max. :2357 Chihuahua : 1 Max. :247087
## (Other) :26
## fdi_2020 fdi_2021 fdi_2022 OBJECTID
## Min. : 85.0 Min. : 69.00 Min. : 61.00 Min. : 888
## 1st Qu.: 112.8 1st Qu.: 97.75 1st Qu.: 78.25 1st Qu.:1047
## Median : 182.5 Median : 159.50 Median : 134.00 Median :1081
## Mean : 264.7 Mean : 252.47 Mean : 197.50 Mean :1219
## 3rd Qu.: 324.5 3rd Qu.: 300.50 3rd Qu.: 253.00 3rd Qu.:1118
## Max. :1429.0 Max. :1542.00 Max. :1020.00 Max. :2357
##
## GMI_ADMIN SQKM SQMI COLOR_MAP Shape_Leng
## MEX-AGS: 1 Min. : 1343 Min. : 518.4 3 : 4 Min. : 1.313
## MEX-BCN: 1 1st Qu.: 23484 1st Qu.: 9067.4 7 : 4 1st Qu.: 8.544
## MEX-BCS: 1 Median : 58628 Median :22636.4 8 : 4 Median :13.341
## MEX-CDZ: 1 Mean : 61289 Mean :23663.8 2 : 3 Mean :13.385
## MEX-CHH: 1 3rd Qu.: 74078 3rd Qu.:28601.3 6 : 3 3rd Qu.:18.266
## MEX-CHP: 1 Max. :247935 Max. :95727.7 1 : 2 Max. :26.273
## (Other):26 (Other):12
## Shape_Area longitude latitude border_dis
## Min. : 0.115 Min. :-115.10 Min. :16.50 Min. : 8.83
## 1st Qu.: 2.000 1st Qu.:-103.70 1st Qu.:19.07 1st Qu.: 613.26
## Median : 5.126 Median : -99.87 Median :20.69 Median : 751.64
## Mean : 5.441 Mean :-100.50 Mean :21.76 Mean : 704.92
## 3rd Qu.: 6.662 3rd Qu.: -98.08 3rd Qu.:24.46 3rd Qu.: 875.76
## Max. :22.891 Max. : -88.27 Max. :30.58 Max. :1252.66
##
## fdi_1999_2 popdensity popdensi_1 jobdensity
## Min. : 944 Min. : 10.36 Min. : 10.69 Min. : 4.99
## 1st Qu.: 1295 1st Qu.: 37.49 1st Qu.: 37.92 1st Qu.: 18.07
## Median : 2142 Median : 64.91 Median : 66.36 Median : 26.80
## Mean : 3544 Mean : 306.43 Mean : 308.31 Mean : 143.56
## 3rd Qu.: 3454 3rd Qu.: 173.08 3rd Qu.: 174.65 3rd Qu.: 68.43
## Max. :30639 Max. :6135.57 Max. :6145.37 Max. :2990.35
##
## jobdensi_1 patentsrat patentsr_1 exchangera
## Min. : 5.48 Min. :0.00000 Min. :0.00000 Min. :20.52
## 1st Qu.: 18.93 1st Qu.:0.01000 1st Qu.:0.01000 1st Qu.:20.52
## Median : 28.52 Median :0.02000 Median :0.02000 Median :20.52
## Mean : 145.64 Mean :0.02625 Mean :0.02406 Mean :20.52
## 3rd Qu.: 69.81 3rd Qu.:0.03000 3rd Qu.:0.03000 3rd Qu.:20.52
## Max. :3020.26 Max. :0.09000 Max. :0.09000 Max. :20.52
##
## exchange_1 realwages_ realwages1 crime_rate
## Min. :19.18 Min. :253.7 Min. :250.0 Min. : 1.99
## 1st Qu.:19.18 1st Qu.:291.7 1st Qu.:286.3 1st Qu.:13.46
## Median :19.18 Median :312.2 Median :310.3 Median :18.98
## Mean :19.18 Mean :318.0 Mean :312.4 Mean :29.96
## 3rd Qu.:19.18 3rd Qu.:336.1 3rd Qu.:333.5 3rd Qu.:42.97
## Max. :19.18 Max. :441.9 Max. :428.4 Max. :97.19
##
## crime_ra_1 border_distance fdi_1999_2022 popdensity_2018
## Min. : 3.08 Min. : 8.83 Min. : 944 Min. : 10.36
## 1st Qu.:10.90 1st Qu.: 613.26 1st Qu.: 1295 1st Qu.: 37.49
## Median :19.25 Median : 751.64 Median : 2142 Median : 64.91
## Mean :29.95 Mean : 704.92 Mean : 3544 Mean : 306.43
## 3rd Qu.:43.02 3rd Qu.: 875.76 3rd Qu.: 3454 3rd Qu.: 173.08
## Max. :91.69 Max. :1252.66 Max. :30639 Max. :6135.57
##
## popdensity_2019 jobdensity_2018 jobdensity_2019 patentsrate_2016
## Min. : 10.69 Min. : 4.994 Min. : 5.478 Min. :0.00100
## 1st Qu.: 37.92 1st Qu.: 18.063 1st Qu.: 18.936 1st Qu.:0.00875
## Median : 66.36 Median : 26.793 Median : 28.520 Median :0.01700
## Mean : 308.31 Mean : 143.559 Mean : 145.644 Mean :0.02497
## 3rd Qu.: 174.65 3rd Qu.: 68.429 3rd Qu.: 69.803 3rd Qu.:0.03050
## Max. :6145.37 Max. :2990.353 Max. :3020.262 Max. :0.09100
##
## patentsrate_2017 exchangerate_2018 exchangerate_2019 realwages_2016
## Min. :0.00100 Min. :20.52 Min. :19.18 Min. :253.7
## 1st Qu.:0.00775 1st Qu.:20.52 1st Qu.:19.18 1st Qu.:291.7
## Median :0.01600 Median :20.52 Median :19.18 Median :312.3
## Mean :0.02284 Mean :20.52 Mean :19.18 Mean :318.0
## 3rd Qu.:0.02650 3rd Qu.:20.52 3rd Qu.:19.18 3rd Qu.:336.1
## Max. :0.08800 Max. :20.52 Max. :19.18 Max. :441.9
##
## realwages_2017 automotive_industry industrial_parks crime_rate_2018
## Min. :250.0 Min. : 0.00 Min. : 0.000 Min. : 1.99
## 1st Qu.:286.3 1st Qu.: 0.00 1st Qu.: 0.000 1st Qu.:13.46
## Median :310.3 Median : 0.00 Median : 0.000 Median :18.98
## Mean :312.4 Mean : 24.53 Mean : 8.781 Mean :29.96
## 3rd Qu.:333.5 3rd Qu.: 16.75 3rd Qu.: 9.500 3rd Qu.:42.97
## Max. :428.4 Max. :204.00 Max. :51.000 Max. :97.19
##
## crime_rate_2019
## Min. : 3.08
## 1st Qu.:10.90
## Median :19.25
## Mean :29.95
## 3rd Qu.:43.02
## Max. :91.69
## summary(cc)## OBJECTID GMI_ADMIN SQKM SQMI
## Min. : 888 Length:32 Min. : 1343 Min. : 518.4
## 1st Qu.:1047 Class :character 1st Qu.: 23484 1st Qu.: 9067.4
## Median :1081 Mode :character Median : 58628 Median :22636.4
## Mean :1219 Mean : 61289 Mean :23663.8
## 3rd Qu.:1118 3rd Qu.: 74078 3rd Qu.:28601.3
## Max. :2357 Max. :247935 Max. :95727.7
## COLOR_MAP Shape_Leng Shape_Area longitude
## Length:32 Min. : 1.313 Min. : 0.115 Min. :-115.10
## Class :character 1st Qu.: 8.544 1st Qu.: 2.000 1st Qu.:-103.70
## Mode :character Median :13.341 Median : 5.126 Median : -99.87
## Mean :13.385 Mean : 5.441 Mean :-100.50
## 3rd Qu.:18.266 3rd Qu.: 6.662 3rd Qu.: -98.08
## Max. :26.273 Max. :22.891 Max. : -88.27
## latitude id state region
## Min. :16.50 Min. : 888 Length:32 Length:32
## 1st Qu.:19.07 1st Qu.:1047 Class :character Class :character
## Median :20.69 Median :1081 Mode :character Mode :character
## Mean :21.76 Mean :1219
## 3rd Qu.:24.46 3rd Qu.:1118
## Max. :30.58 Max. :2357
## land_area border_dis fdi_2020 fdi_2021
## Min. : 1499 Min. : 8.83 Min. : 85.0 Min. : 69.00
## 1st Qu.: 23861 1st Qu.: 613.26 1st Qu.: 112.8 1st Qu.: 97.75
## Median : 58978 Median : 751.64 Median : 182.5 Median : 159.50
## Mean : 60038 Mean : 704.92 Mean : 264.7 Mean : 252.47
## 3rd Qu.: 73730 3rd Qu.: 875.76 3rd Qu.: 324.5 3rd Qu.: 300.50
## Max. :247087 Max. :1252.66 Max. :1429.0 Max. :1542.00
## fdi_2022 fdi_1999_2 popdensity popdensi_1
## Min. : 61.00 Min. : 944 Min. : 10.36 Min. : 10.69
## 1st Qu.: 78.25 1st Qu.: 1295 1st Qu.: 37.49 1st Qu.: 37.92
## Median : 134.00 Median : 2142 Median : 64.91 Median : 66.36
## Mean : 197.50 Mean : 3544 Mean : 306.43 Mean : 308.31
## 3rd Qu.: 253.00 3rd Qu.: 3454 3rd Qu.: 173.08 3rd Qu.: 174.65
## Max. :1020.00 Max. :30639 Max. :6135.57 Max. :6145.37
## jobdensity jobdensi_1 patentsrat patentsr_1
## Min. : 4.99 Min. : 5.48 Min. :0.00000 Min. :0.00000
## 1st Qu.: 18.07 1st Qu.: 18.93 1st Qu.:0.01000 1st Qu.:0.01000
## Median : 26.80 Median : 28.52 Median :0.02000 Median :0.02000
## Mean : 143.56 Mean : 145.64 Mean :0.02625 Mean :0.02406
## 3rd Qu.: 68.43 3rd Qu.: 69.81 3rd Qu.:0.03000 3rd Qu.:0.03000
## Max. :2990.35 Max. :3020.26 Max. :0.09000 Max. :0.09000
## exchangera exchange_1 realwages_ realwages1
## Min. :20.52 Min. :19.18 Min. :253.7 Min. :250.0
## 1st Qu.:20.52 1st Qu.:19.18 1st Qu.:291.7 1st Qu.:286.3
## Median :20.52 Median :19.18 Median :312.2 Median :310.3
## Mean :20.52 Mean :19.18 Mean :318.0 Mean :312.4
## 3rd Qu.:20.52 3rd Qu.:19.18 3rd Qu.:336.1 3rd Qu.:333.5
## Max. :20.52 Max. :19.18 Max. :441.9 Max. :428.4
## crime_rate crime_ra_1 geometry
## Min. : 1.99 Min. : 3.08 MULTIPOLYGON :32
## 1st Qu.:13.46 1st Qu.:10.90 epsg:4326 : 0
## Median :18.98 Median :19.25 +proj=long...: 0
## Mean :29.96 Mean :29.95
## 3rd Qu.:42.97 3rd Qu.:43.02
## Max. :97.19 Max. :91.69db_data <- data.frame(db)
nearshoring_map_data <- data.frame(nearshoring_map.sp)Create map
plot(nearshoring_map2, pbg="blue", axes = T)
title("Spatial polygon plot") 
Mapping for points for the same area
file.choose()## [1] "F:\\TEC\\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\\Proyecto\\database_nearshoring\\assembly_plants.csv"#assembly_plants <- read.csv(F:\\TEC\\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\\Proyecto\\assembly_plants.csv)
#plot(assembly_plants, axes=T, lwd = 1, pch=16, cex=0.1)
#title("Spatial Points Plot") Geocoding using Open API
Install package to generate geocoordinates using Google Map API
library(dplyr) # for data manipulation
library(ggplot2) # for visualiziation##
## Attaching package: 'ggplot2'## The following objects are masked from 'package:psych':
##
## %+%, alphalibrary(maptools) # for mapping in R
library(ggmap) # for using Google Map with API key## Google's Terms of Service: https://cloud.google.com/maps-platform/terms/.## Please cite ggmap if you use it! See citation("ggmap") for details.Google’s Geocoding API
register_google(key = "AIzaSyDJcqoz9ff23VlhazQyZthFpruiOpIVTHE")Convert address of Audi Assembly Plant in Mexico into coordinates for mapping
audi_puebla.location <- data.frame(addr = c("75012 San José Chiapa, Pue."), stringsAsFactors = FALSE)
audi_puebla.latlon <- audi_puebla.location %>% mutate_geocode(addr)## Source : https://maps.googleapis.com/maps/api/geocode/json?address=75012+San+Jos%C3%A9+Chiapa,+Pue.&key=xxxaudi_puebla.latlon## addr lon lat
## 1 75012 San José Chiapa, Pue. -97.76884 19.21945bmw_slp.location <- data.frame(addr = c("Boulevard BMW No. 655 Parque Industrial Desarrollo Logistik II, 79526 S.L.P."), stringsAsFactors = FALSE)
bmw_slp.latlon <- bmw_slp.location %>% mutate_geocode(addr)## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Boulevard+BMW+No.+655+Parque+Industrial+Desarrollo+Logistik+II,+79526+S.L.P.&key=xxx## "Boulevard BMW No...." not uniquely geocoded, using "av. america 102, 79525 laguna de san vicente, s.l.p., mexico"bmw_slp.latlon## addr
## 1 Boulevard BMW No. 655 Parque Industrial Desarrollo Logistik II, 79526 S.L.P.
## lon lat
## 1 -100.8647 21.97944Mark the location on the map using coordinates
audi_puebla.cen <- c(mean(audi_puebla.latlon$lon), mean(audi_puebla.latlon$lat))
audi_puebla <- audi_puebla.latlon[,2:3]
bmw_slp.cen <- c(mean(bmw_slp.latlon$lon), mean(bmw_slp.latlon$lat))
bmw_slp <- bmw_slp.latlon[,2:3]Create a road map of Audi Assembly Plant in Mexico
audi_puebla.road_map <- get_googlemap(center = audi_puebla.cen, maptype = "roadmap", zoom=12, markers = audi_puebla) ## Source : https://maps.googleapis.com/maps/api/staticmap?center=19.219454,-97.768836&zoom=12&size=640x640&scale=2&maptype=roadmap&markers=19.219454,-97.768836&key=xxxggmap(audi_puebla.road_map)
Convert a number of addresses of Assembly Plants in Mexico into coordinates using csv file
setwd("F:\\TEC\\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\\Proyecto\\database_nearshoring")
getwd## function ()
## .Internal(getwd())
## <bytecode: 0x0000026db10b8ec8>
## <environment: namespace:base>Load Assembly Plants in Mexico Address Data
assembly_plants.data <- read.csv('assembly_plants.csv', header=T, encoding = 'utf-8' )
names(assembly_plants.data) ## [1] "id" "addr" "name"Create data frame for geocoding
assembly_plants.data$addr <- enc2utf8(assembly_plants.data$addr)Geocode the Assembly Plants addresses in Mexico
assembly_plants.latlon <- assembly_plants.data %>% mutate_geocode(addr) ## Source : https://maps.googleapis.com/maps/api/geocode/json?address=75012+San+Jos%C3%A9+Chiapa,+Pue.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Boulevard+BMW+No.+655+Parque+Industrial+Desarrollo+Logistik+II,+79526+S.L.P.&key=xxx## "Boulevard BMW No...." not uniquely geocoded, using "av. america 102, 79525 laguna de san vicente, s.l.p., mexico"## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Boulevard+Miguel+Ramos+Arizpe+No+124+Sur,+Saltillo,+Zona+Industrial,+25900+Ramos+Arizpe,+Coah.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Km+60.5,+Carr+Toluca+-+M%C3%A9xico,+Delegaci%C3%B3n+Sta+Ana+Tlapaltitl%C3%A1n,+50160+Toluca+de+Lerdo,+M%C3%A9x.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=KM+4.5,+Carr.+a+la+Colorada,+Parque+Industrial,+83299+Hermosillo,+Son.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Quer%C3%A9taro+%E2%80%93+Mexico+Highway+Km+36.5,+Lomas+del+Salitre,+54730+Cuautitl%C3%A1n+Izcalli&key=xxx## "Querétaro – Mexic..." not uniquely geocoded, using "54730 cuautitlán izcalli, state of mexico, mexico"## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Carretera+a+Ciudad+Juarez+Km.+11.5,+Complejo+Industrial+Chihuahua,+31000+Chihuahua&key=xxx## "Carretera a Ciuda..." not uniquely geocoded, using "complejo industrial chihuahua, chihuahua, chih., mexico"## Source : https://maps.googleapis.com/maps/api/geocode/json?address=R%C3%ADo+San+Lorenzo+1475,+Parque+Industrial,+36814+Castro+del+R%C3%ADo,+Gto.&key=xxx## "Río San Lorenzo 1..." not uniquely geocoded, using "guanajuato, mexico"## Source : https://maps.googleapis.com/maps/api/geocode/json?address=25934+Ramos+Arizpe,+Coah.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Carr.+Guanajuato+-+Silao+Kil%C3%B3metro+3.8,+36100+Silao,+Gto.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=M%C3%A9xico+120,+Central,+78390+San+Luis,+S.L.P.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=C.+Industria+Minera+700,+Delegaci%C3%B3n+Sta+Ana+Tlapaltitl%C3%A1n,+50160+Toluca+de+Lerdo,+M%C3%A9x.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Parque+Industrial,+Carretera+a+El+Castillo+7250,+45680+El+Salto,+Jal.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=La+Luz+Sur,+38140+Celaya,+Gto.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Kia+Motors+Avenue,+66664+N.L.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=36875+Salamanca,+Gto.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Carretera+Federal+Lagos+de+Moreno+Km+75,+20290+Aguascalientes,+Ags.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Carretera+Federal+Cuernavaca+-+Cuatla+Km.+4.5,+Civac,+62578+Jiutepec,+Mor.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Carretera+Libre,+Tijuana+-+Tecate+33143,+El+Realito,+22550+Tecate,+B.C.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Quer%C3%A9taro+-+Celaya+KM+45,+38160+Apaseo+el+Grande,+Gto.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Av.+San+Lorenzo+Almecatla+16,+Sanctorum,+72730+Sanctorum,+Pue.&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=Ave.+Mineral+de+Valenciana+Puerto,+Interior+611,+36275+Silao,+Gto.&key=xxxhead(assembly_plants.latlon)## id
## 1 1
## 2 2
## 3 3
## 4 4
## 5 5
## 6 6
## addr
## 1 75012 San José Chiapa, Pue.
## 2 Boulevard BMW No. 655 Parque Industrial Desarrollo Logistik II, 79526 S.L.P.
## 3 Boulevard Miguel Ramos Arizpe No 124 Sur, Saltillo, Zona Industrial, 25900 Ramos Arizpe, Coah.
## 4 Km 60.5, Carr Toluca - México, Delegación Sta Ana Tlapaltitlán, 50160 Toluca de Lerdo, Méx.
## 5 KM 4.5, Carr. a la Colorada, Parque Industrial, 83299 Hermosillo, Son.
## 6 Querétaro – Mexico Highway Km 36.5, Lomas del Salitre, 54730 Cuautitlán Izcalli
## name lon lat
## 1 audi_pue -97.76884 19.21945
## 2 bmw_slp -100.86473 21.97944
## 3 fca_coa -100.94397 25.53833
## 4 fca_edm -99.62117 19.28812
## 5 ford_son -110.91624 29.00784
## 6 ford_edm -99.19544 19.66013Create a base map
basemap <- qmap('mexico', zoom = 4, maptype = "roadmap")## Source : https://maps.googleapis.com/maps/api/staticmap?center=mexico&zoom=4&size=640x640&scale=2&maptype=roadmap&language=en-EN&key=xxx## Source : https://maps.googleapis.com/maps/api/geocode/json?address=mexico&key=xxxMark the locations of Assembly Plants on the base map
basemap + geom_point(data = assembly_plants.latlon, aes(lon, lat), size = 1.5, colour= "red")
Identify spatial distribution of variables of interest
Loading Libraries
library(ggplot2)
library(spdep)
library(MASS)
library(spmoran)
library(spatialreg)
library(sphet)
library(maptools)
library(sp)
library(sf)
library(maps)
library(rgeos)
library(ggmap)
library(mapproj)
library(RColorBrewer)
library(ggsn)
### importing spatial data
map<-readShapePoly("F:\\TEC\\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\\Proyecto\\database_nearshoring\\nearshoring_map.shp",IDvar="id",proj4string=CRS("+proj=longlat"))
lmat<-coordinates(map)
map.centroid<-coordinates(map)
summary(map)## Object of class SpatialPolygonsDataFrame
## Coordinates:
## min max
## x -118.40417 -86.73862
## y 14.55055 32.71846
## Is projected: FALSE
## proj4string : [+proj=longlat +datum=WGS84 +no_defs]
## Data attributes:
## OBJECTID GMI_ADMIN SQKM SQMI COLOR_MAP
## Min. : 888 MEX-AGS: 1 Min. : 1343 Min. : 518.4 3 : 4
## 1st Qu.:1047 MEX-BCN: 1 1st Qu.: 23484 1st Qu.: 9067.4 7 : 4
## Median :1081 MEX-BCS: 1 Median : 58628 Median :22636.4 8 : 4
## Mean :1219 MEX-CDZ: 1 Mean : 61289 Mean :23663.8 2 : 3
## 3rd Qu.:1118 MEX-CHH: 1 3rd Qu.: 74078 3rd Qu.:28601.3 6 : 3
## Max. :2357 MEX-CHP: 1 Max. :247935 Max. :95727.7 1 : 2
## (Other):26 (Other):12
## Shape_Leng Shape_Area longitude latitude
## Min. : 1.313 Min. : 0.115 Min. :-115.10 Min. :16.50
## 1st Qu.: 8.544 1st Qu.: 2.000 1st Qu.:-103.70 1st Qu.:19.07
## Median :13.341 Median : 5.126 Median : -99.87 Median :20.69
## Mean :13.385 Mean : 5.441 Mean :-100.50 Mean :21.76
## 3rd Qu.:18.266 3rd Qu.: 6.662 3rd Qu.: -98.08 3rd Qu.:24.46
## Max. :26.273 Max. :22.891 Max. : -88.27 Max. :30.58
##
## id state region land_area
## Min. : 888 Aguascalientes : 1 Centro :8 Min. : 1499
## 1st Qu.:1047 Baja California : 1 Centro-Norte :5 1st Qu.: 23861
## Median :1081 Baja California Sur: 1 Norte :6 Median : 58978
## Mean :1219 Campeche : 1 Norte-Occidente:5 Mean : 60038
## 3rd Qu.:1118 Chiapas : 1 Sur :8 3rd Qu.: 73730
## Max. :2357 Chihuahua : 1 Max. :247087
## (Other) :26
## border_dis fdi_2020 fdi_2021 fdi_2022
## Min. : 8.83 Min. : 85.0 Min. : 69.00 Min. : 61.00
## 1st Qu.: 613.26 1st Qu.: 112.8 1st Qu.: 97.75 1st Qu.: 78.25
## Median : 751.64 Median : 182.5 Median : 159.50 Median : 134.00
## Mean : 704.92 Mean : 264.7 Mean : 252.47 Mean : 197.50
## 3rd Qu.: 875.76 3rd Qu.: 324.5 3rd Qu.: 300.50 3rd Qu.: 253.00
## Max. :1252.66 Max. :1429.0 Max. :1542.00 Max. :1020.00
##
## fdi_1999_2 popdensity popdensi_1 jobdensity
## Min. : 944 Min. : 10.36 Min. : 10.69 Min. : 4.99
## 1st Qu.: 1295 1st Qu.: 37.49 1st Qu.: 37.92 1st Qu.: 18.07
## Median : 2142 Median : 64.91 Median : 66.36 Median : 26.80
## Mean : 3544 Mean : 306.43 Mean : 308.31 Mean : 143.56
## 3rd Qu.: 3454 3rd Qu.: 173.08 3rd Qu.: 174.65 3rd Qu.: 68.43
## Max. :30639 Max. :6135.57 Max. :6145.37 Max. :2990.35
##
## jobdensi_1 patentsrat patentsr_1 exchangera
## Min. : 5.48 Min. :0.00000 Min. :0.00000 Min. :20.52
## 1st Qu.: 18.93 1st Qu.:0.01000 1st Qu.:0.01000 1st Qu.:20.52
## Median : 28.52 Median :0.02000 Median :0.02000 Median :20.52
## Mean : 145.64 Mean :0.02625 Mean :0.02406 Mean :20.52
## 3rd Qu.: 69.81 3rd Qu.:0.03000 3rd Qu.:0.03000 3rd Qu.:20.52
## Max. :3020.26 Max. :0.09000 Max. :0.09000 Max. :20.52
##
## exchange_1 realwages_ realwages1 crime_rate
## Min. :19.18 Min. :253.7 Min. :250.0 Min. : 1.99
## 1st Qu.:19.18 1st Qu.:291.7 1st Qu.:286.3 1st Qu.:13.46
## Median :19.18 Median :312.2 Median :310.3 Median :18.98
## Mean :19.18 Mean :318.0 Mean :312.4 Mean :29.96
## 3rd Qu.:19.18 3rd Qu.:336.1 3rd Qu.:333.5 3rd Qu.:42.97
## Max. :19.18 Max. :441.9 Max. :428.4 Max. :97.19
##
## crime_ra_1
## Min. : 3.08
## 1st Qu.:10.90
## Median :19.25
## Mean :29.95
## 3rd Qu.:43.02
## Max. :91.69
## ### create spatial connectivity matrix (swm)
map.link<-poly2nb(map,queen=T) ### set SWM to queen
map.linkW<-nb2listw(map.link, style="W") ### SWM queen standardized
plot(map,border="blue",axes=FALSE,las=1)
plot(map,col="grey",border=grey(0.9),axes=T,add=T)
plot(map.linkW,coords=map.centroid,pch=19,cex=0.1,col="red",add=T) ### plot SWM queen standarized
title("Spatial Connectivity Matrix - Contiguity Case")
### Testing for Spatial Autocorrelation using Moran's I
# detecting spatial autocorrelation in main variable of interest
moran.mc(map$fdi_2020,map.linkW,nsim=9999) ### no statistically significant ##
## Monte-Carlo simulation of Moran I
##
## data: map$fdi_2020
## weights: map.linkW
## number of simulations + 1: 10000
##
## statistic = 0.047224, observed rank = 8177, p-value = 0.1823
## alternative hypothesis: greatermoran.mc(map$fdi_2021,map.linkW,nsim=9999) ### no statistically significant ##
## Monte-Carlo simulation of Moran I
##
## data: map$fdi_2021
## weights: map.linkW
## number of simulations + 1: 10000
##
## statistic = 0.032588, observed rank = 7836, p-value = 0.2164
## alternative hypothesis: greatermoran.mc(map$fdi_2022,map.linkW,nsim=9999) ### no statistically significant ##
## Monte-Carlo simulation of Moran I
##
## data: map$fdi_2022
## weights: map.linkW
## number of simulations + 1: 10000
##
## statistic = 0.037805, observed rank = 7787, p-value = 0.2213
## alternative hypothesis: greatermoran.mc(map$fdi_1999_2,map.linkW,nsim=9999) ### no statistically significant ##
## Monte-Carlo simulation of Moran I
##
## data: map$fdi_1999_2
## weights: map.linkW
## number of simulations + 1: 10000
##
## statistic = 0.032409, observed rank = 8648, p-value = 0.1352
## alternative hypothesis: greater# detecting spatial autocorrelation in control variable - job density
moran.mc(map$jobdensity,map.linkW,nsim=9999) ### statistically significant ##
## Monte-Carlo simulation of Moran I
##
## data: map$jobdensity
## weights: map.linkW
## number of simulations + 1: 10000
##
## statistic = 0.08376, observed rank = 9984, p-value = 0.0016
## alternative hypothesis: greatermoran.mc(map$jobdensi_1,map.linkW,nsim=9999) ### statistically significant ##
## Monte-Carlo simulation of Moran I
##
## data: map$jobdensi_1
## weights: map.linkW
## number of simulations + 1: 10000
##
## statistic = 0.083154, observed rank = 9980, p-value = 0.002
## alternative hypothesis: greaterdetecting spatial autocorrelation in control variable - crime rate
moran.mc(map$crime_rate,map.linkW,nsim=9999) ### statistically significant ##
## Monte-Carlo simulation of Moran I
##
## data: map$crime_rate
## weights: map.linkW
## number of simulations + 1: 10000
##
## statistic = -0.018213, observed rank = 5792, p-value = 0.4208
## alternative hypothesis: greatermoran.mc(map$crime_ra_1,map.linkW,nsim=9999) ### statistically significant ##
## Monte-Carlo simulation of Moran I
##
## data: map$crime_ra_1
## weights: map.linkW
## number of simulations + 1: 10000
##
## statistic = -0.070396, observed rank = 4016, p-value = 0.5984
## alternative hypothesis: greaterTesting for Spatial Autocorrelation using Moran’s I
# Lets use GeoDa to conduct cluster analysis
# install.packages("rgeoda")
library(rgeoda)## Loading required package: digest##
## Attaching package: 'rgeoda'## The following object is masked from 'package:spdep':
##
## skatermap_shapefile<-st_read("F:\\TEC\\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\\Proyecto\\database_nearshoring\\nearshoring_map.shp") ## Reading layer `nearshoring_map' from data source
## `F:\TEC\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\Proyecto\database_nearshoring\nearshoring_map.shp'
## using driver `ESRI Shapefile'
## Simple feature collection with 32 features and 30 fields
## Geometry type: MULTIPOLYGON
## Dimension: XY
## Bounding box: xmin: -118.4042 ymin: 14.55055 xmax: -86.73862 ymax: 32.71846
## Geodetic CRS: WGS 84# create spatial weights matrix object (contiguity case)
queen_w<-queen_weights(map_shapefile)
# calculate LISA as per GEODA
lisa_dv<-local_moran(queen_w, map_shapefile["fdi_2022"])
lisa_cv1<-local_moran(queen_w, map_shapefile["jobdensity"])
lisa_cv2<-local_moran(queen_w, map_shapefile["patentsrat"])
lisa_cv3<-local_moran(queen_w, map_shapefile["crime_rate"])
# add cluster calculations to dataset
map_shapefile$cluster_dv<-as.factor(lisa_dv$GetClusterIndicators())
levels(map_shapefile$cluster_dv)<-lisa_dv$GetLabels() ### dependent variable
map_shapefile$cluster_cv1<-as.factor(lisa_cv1$GetClusterIndicators())
levels(map_shapefile$cluster_cv1)<-lisa_cv1$GetLabels() ### job density
map_shapefile$cluster_cv2<-as.factor(lisa_cv2$GetClusterIndicators())
levels(map_shapefile$cluster_cv2)<-lisa_cv2$GetLabels() ### patents rate
map_shapefile$cluster_cv3<-as.factor(lisa_cv3$GetClusterIndicators())
levels(map_shapefile$cluster_cv3)<-lisa_cv3$GetLabels() ### crime rate
# visualize clusters using GeoDa view
ggplot(data=map_shapefile) +
geom_sf(aes(fill=cluster_dv)) +
ggtitle(label = "Foreign Direct Investment 2022", subtitle = "Number of Companies - Recipients of FDI's Inflows")
ggplot(data=map_shapefile) +
geom_sf(aes(fill=cluster_cv1)) +
ggtitle(label = "Job Density 2018", subtitle = "Population per km2")
ggplot(data=map_shapefile) +
geom_sf(aes(fill=cluster_cv2)) +
ggtitle(label = "R&D Activity 2016", subtitle = "Patents per 100,000 Population")
ggplot(data=map_shapefile) +
geom_sf(aes(fill=cluster_cv3)) +
ggtitle(label = "Crime 2018", subtitle = "Homicide Rate per 100,000 Population")
### Bivariate Moran's I
lisa_bi <- local_bimoran(queen_w,map_shapefile[c('fdi_2022','jobdensi_1')])
bivariate_m <- lisa_values(lisa_bi)
bivariate_m## [1] -0.149741647 -0.060446836 -0.357083500 -0.061212625 0.120402093
## [6] 0.158797531 0.129265178 0.027608417 0.113414498 0.054051392
## [11] 0.053969583 -0.133997214 -0.027204696 0.099445131 0.056132536
## [16] 0.032506321 0.151446874 0.002192283 0.040700508 0.472159801
## [21] -0.024360918 0.135208601 0.889435209 -0.733999145 -0.008395307
## [26] 0.114141307 0.115150229 0.136668836 0.013246490 -0.409137378
## [31] 0.112890258 0.166408594map_shapefile$bivariate<-as.factor(lisa_bi$GetClusterIndicators())
levels(map_shapefile$bivariate)<-lisa_bi$GetLabels() ### bivariate moran
# visualize Bivariate Moran's i using GeoDa view
ggplot(data=map_shapefile) +
geom_sf(aes(fill=bivariate)) +
ggtitle(label = "Foreign Direct Investment & Job Density", subtitle = "Bivariate Local's Moran's I")
Suggest hypotheses about potential relationships among the variables of interest
Hypothesis 1: There is a positive relationship between fdi inflows and job density
Hypothesis 2: There is a positive relationship between fdi inflows and R&D activity
Hypothesis 3: There is a negative relationship between fdi inflows and crime rate
Spatial regression analysis and interpretation of estimated results
# lets import dataset
data=read.csv("F:\\TEC\\2023-01-09 Lunes a 13 Viernes Curso Análisis Espacial de Datos usando R\\Proyecto\\database_nearshoring\\nearshoring.csv")
summary(data)## id state region land_area
## Min. : 888 Length:32 Length:32 Min. : 1499
## 1st Qu.:1047 Class :character Class :character 1st Qu.: 23861
## Median :1081 Mode :character Mode :character Median : 58978
## Mean :1219 Mean : 60038
## 3rd Qu.:1118 3rd Qu.: 73730
## Max. :2357 Max. :247087
## border_distance fdi_2020 fdi_2021 fdi_2022
## Min. : 8.83 Min. : 85.0 Min. : 69.00 Min. : 61.00
## 1st Qu.: 613.26 1st Qu.: 112.8 1st Qu.: 97.75 1st Qu.: 78.25
## Median : 751.64 Median : 182.5 Median : 159.50 Median : 134.00
## Mean : 704.92 Mean : 264.7 Mean : 252.47 Mean : 197.50
## 3rd Qu.: 875.76 3rd Qu.: 324.5 3rd Qu.: 300.50 3rd Qu.: 253.00
## Max. :1252.66 Max. :1429.0 Max. :1542.00 Max. :1020.00
## fdi_1999_2022 popdensity_2018 popdensity_2019 jobdensity_2018
## Min. : 944 Min. : 10.36 Min. : 10.69 Min. : 4.994
## 1st Qu.: 1295 1st Qu.: 37.49 1st Qu.: 37.92 1st Qu.: 18.063
## Median : 2142 Median : 64.91 Median : 66.36 Median : 26.793
## Mean : 3544 Mean : 306.43 Mean : 308.31 Mean : 143.559
## 3rd Qu.: 3454 3rd Qu.: 173.08 3rd Qu.: 174.65 3rd Qu.: 68.429
## Max. :30639 Max. :6135.57 Max. :6145.37 Max. :2990.353
## jobdensity_2019 patentsrate_2016 patentsrate_2017 exchangerate_2018
## Min. : 5.478 Min. :0.00100 Min. :0.00100 Min. :20.52
## 1st Qu.: 18.936 1st Qu.:0.00875 1st Qu.:0.00775 1st Qu.:20.52
## Median : 28.520 Median :0.01700 Median :0.01600 Median :20.52
## Mean : 145.644 Mean :0.02497 Mean :0.02284 Mean :20.52
## 3rd Qu.: 69.803 3rd Qu.:0.03050 3rd Qu.:0.02650 3rd Qu.:20.52
## Max. :3020.262 Max. :0.09100 Max. :0.08800 Max. :20.52
## exchangerate_2019 realwages_2016 realwages_2017 automotive_industry
## Min. :19.18 Min. :253.7 Min. :250.0 Min. : 0.00
## 1st Qu.:19.18 1st Qu.:291.7 1st Qu.:286.3 1st Qu.: 0.00
## Median :19.18 Median :312.3 Median :310.3 Median : 0.00
## Mean :19.18 Mean :318.0 Mean :312.4 Mean : 24.53
## 3rd Qu.:19.18 3rd Qu.:336.1 3rd Qu.:333.5 3rd Qu.: 16.75
## Max. :19.18 Max. :441.9 Max. :428.4 Max. :204.00
## industrial_parks crime_rate_2018 crime_rate_2019
## Min. : 0.000 Min. : 1.99 Min. : 3.08
## 1st Qu.: 0.000 1st Qu.:13.46 1st Qu.:10.90
## Median : 0.000 Median :18.98 Median :19.25
## Mean : 8.781 Mean :29.96 Mean :29.95
## 3rd Qu.: 9.500 3rd Qu.:42.97 3rd Qu.:43.02
## Max. :51.000 Max. :97.19 Max. :91.69# lag_model
model1<-lagsarlm(log(fdi_2022)~log(jobdensity_2019)+patentsrate_2017+log(crime_rate_2019)+log(realwages_2017)+log(border_distance),
data=data,listw=map.linkW, zero.policy=TRUE, na.action=na.omit)
# summary(model1)
# error_model
model2<-errorsarlm(log(fdi_2022)~log(jobdensity_2019)+patentsrate_2017+log(crime_rate_2019)+log(realwages_2017)+log(border_distance),
data=data,listw=map.linkW, zero.policy=TRUE, na.action=na.omit)
# summary(model2)
# durbin_model
model3<-lagsarlm(log(fdi_2022)~log(jobdensity_2019)+patentsrate_2017+log(crime_rate_2019)+log(realwages_2017)+log(border_distance), data=data, listw=map.linkW, Durbin=TRUE)
# summary(model3)
# compare spatial regression results
# install.packages("huxtable")
library(jtools)
jtools::export_summs(model1, model2, model3)## Registered S3 methods overwritten by 'broom':
## method from
## tidy.glht jtools
## tidy.summary.glht jtools| Model 1 | Model 2 | Model 3 | |
|---|---|---|---|
| rho | -0.27 | -0.60 ** | |
| (0.19) | (0.22) | ||
| (Intercept) | 2.32 | 3.27 | 6.37 |
| (5.14) | (4.94) | (13.09) | |
| log(jobdensity_2019) | 0.07 | 0.03 | 0.07 |
| (0.08) | (0.07) | (0.07) | |
| patentsrate_2017 | 12.99 * | 14.03 ** | 11.38 * |
| (5.05) | (4.97) | (5.79) | |
| log(crime_rate_2019) | 0.13 | 0.11 | 0.11 |
| (0.11) | (0.09) | (0.10) | |
| log(realwages_2017) | 0.90 | 0.58 | 0.77 |
| (0.85) | (0.82) | (0.88) | |
| log(border_distance) | -0.33 ** | -0.37 *** | -0.35 *** |
| (0.11) | (0.10) | (0.10) | |
| lambda | -0.61 ** | ||
| (0.22) | |||
| lag.log(jobdensity_2019) | -0.08 | ||
| (0.16) | |||
| lag.patentsrate_2017 | 7.28 | ||
| (10.10) | |||
| lag.log(crime_rate_2019) | 0.16 | ||
| (0.23) | |||
| lag.log(realwages_2017) | 0.19 | ||
| (1.69) | |||
| lag.log(border_distance) | -0.45 | ||
| (0.29) | |||
| nobs | 32 | 32 | 32 |
| r.squared | 0.59 | 0.65 | 0.68 |
| AIC | 56.94 | 54.81 | 61.49 |
| BIC | 68.66 | 66.54 | 80.55 |
| deviance | 6.62 | 5.74 | 5.20 |
| logLik | -20.47 | -19.41 | -17.75 |
| nobs.1 | 32.00 | 32.00 | 32.00 |
| *** p < 0.001; ** p < 0.01; * p < 0.05. | |||
Create interactive map
Conclusions
LS0tDQp0aXRsZTogIlNwYXRpYWwgRGF0YSBBbmFseXNpcyBVc2luZyBSIC0gQ29sbGFib3JhdGl2ZSBQcm9qZWN0OiBOZWFyc2hvcmluZyINCmF1dGhvcjogIlRlYW0gNTogR2FicmllbGEgTW9uZm9ydGUsIERhdmlkIFNhdWNlZG8sIFNhbXVlbCBNb250ZW5lZ3JvLCBhbmQgUmHDumwgQ2FudMO6Ig0KZGF0ZTogIjIwMjMtMDEtMTEiDQpvdXRwdXQ6IA0KIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiB0cnVlDQogICAgdG9jX2Zsb2F0OiB0cnVlDQogICAgY29kZV9kb3dubG9hZDogdHJ1ZQ0KLS0tDQoNClByb2Zlc3NvcjogRG9oeWVvbmcgS2ltLCBQaC5EDQoNCjxpbWcgc3JjPSJGOlxcVEVDXFwyMDIzLTAxLTA5IEx1bmVzIGEgMTMgVmllcm5lcyBDdXJzbyBBbsOhbGlzaXMgRXNwYWNpYWwgZGUgRGF0b3MgdXNhbmRvIFJcXFByb3llY3RvXFxkYXRhYmFzZV9uZWFyc2hvcmluZ1xcbmVhcnNob3JpbmcucG5nIj4NCg0KIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBibHVlOyI+UHJvamVjdCBPdmVydmlldzwvc3Bhbj4gDQoNCkFjY29yZGluZyB0byBhbiBhbmFseXNpcyBieSBNb3JnYW4gU3RhbmxleSwgdGhlIHRyZW5kIG9mIHJlbG9jYXRpbmcgY29tcGFuaWVzIG5lYXIgdGhlaXIgdGFyZ2V0IG1hcmtldHMsIGtub3duIGFzIG5lYXJzaG9yaW5nLCB3aWxsIHByb3AgdXAgTWV4aWNhbiBleHBvcnRzIHdpdGggYW4gZXN0aW1hdGVkICoqOTQuMiBiaWxsaW9uIGFkZGl0aW9uYWwgZG9sbGFycyBpbiBmaXZlIHllYXJzKioNCg0KVGhlIG9iamVjdGl2ZSBvZiB0aGlzIHByb2plY3QgaXMgdG8gcGVyZm9ybSBhbiBFeHBsb3JhdG9yeSBTcGF0aWFsIERhdGEgQW5hbHlzaXMgKEVTREEpIHRvDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmx1ZTsiPkRlc2NyaXB0aXZlIHN0YXRpc3RpY3Mgb2YgdmFyaWFibGVzIG9mIGludGVyZXN0PC9zcGFuPiANCg0KIyMjIEltcG9ydGluZyBkYXRhIGJhc2UgKGRiKQ0KYGBge3J9DQojIGZpbGUuY2hvb3NlKCkNCmRiIDwtIHJlYWQuY3N2KCJGOlxcVEVDXFwyMDIzLTAxLTA5IEx1bmVzIGEgMTMgVmllcm5lcyBDdXJzbyBBbsOhbGlzaXMgRXNwYWNpYWwgZGUgRGF0b3MgdXNhbmRvIFJcXFByb3llY3RvXFxkYXRhYmFzZV9uZWFyc2hvcmluZ1xcbmVhcnNob3JpbmcuY3N2IikNCmBgYA0KIyMjIFVuZGVyc3RhbmRpbmcgZGINCmBgYHtyfQ0Kc3VtbWFyeSAoZGIpDQpsaWJyYXJ5KGRwbHlyKQ0KaGVhZChkYikNCmNvdW50KGRiLCByZWdpb24sIHNvcnQgPSBUUlVFKQ0KDQpoaXN0KGRiJGxhbmRfYXJlYSkNCmJveHBsb3QoZGIkbGFuZF9hcmVhLCBob3Jpem9udGFsID0gVFJVRSkNCmJveHBsb3QoZGIkbGFuZF9hcmVhLHBsb3Q9RkFMU0UpJG91dA0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoMjQ3MDg3KSkpDQpkYiAlPiUgZmlsdGVyX2FsbChhbnlfdmFycyguICVpbiUgYygxNTE1NzEpKSkNCmRiICU+JSBmaWx0ZXJfYWxsKGFueV92YXJzKC4gJWluJSBjKDE4NDkzNCkpKQ0KDQpoaXN0KGRiJGJvcmRlcl9kaXN0YW5jZSkNCmJveHBsb3QoZGIkYm9yZGVyX2Rpc3RhbmNlLCBob3Jpem9udGFsID0gVFJVRSkNCmJveHBsb3QoZGIkYm9yZGVyX2Rpc3RhbmNlLHBsb3Q9RkFMU0UpJG91dA0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoOC44MykpKQ0KDQpoaXN0KGRiJGZkaV8yMDIwKQ0KYm94cGxvdChkYiRmZGlfMjAyMCwgaG9yaXpvbnRhbCA9IFRSVUUpDQpib3hwbG90KGRiJGZkaV8yMDIwLHBsb3Q9RkFMU0UpJG91dA0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoMTQyOSkpKQ0KDQpoaXN0KGRiJGZkaV8yMDIxKQ0KYm94cGxvdChkYiRmZGlfMjAyMSwgaG9yaXpvbnRhbCA9IFRSVUUpDQpib3hwbG90KGRiJGZkaV8yMDIxLHBsb3Q9RkFMU0UpJG91dA0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoMTU0MikpKQ0KDQpoaXN0KGRiJGZkaV8xOTk5XzIwMjIpDQpib3hwbG90KGRiJGZkaV8xOTk5XzIwMjIsIGhvcml6b250YWwgPSBUUlVFKQ0KYm94cGxvdChkYiRmZGlfMTk5OV8yMDIyLHBsb3Q9RkFMU0UpJG91dA0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoMzA2MzkpKSkNCg0KaGlzdChkYiRwb3BkZW5zaXR5XzIwMTgpDQpib3hwbG90KGRiJHBvcGRlbnNpdHlfMjAxOCwgaG9yaXpvbnRhbCA9IFRSVUUpDQpib3hwbG90KGRiJHBvcGRlbnNpdHlfMjAxOCxwbG90PUZBTFNFKSRvdXQNCmRiICU+JSBmaWx0ZXJfYWxsKGFueV92YXJzKC4gJWluJSBjKDYxMzUuNTY2KSkpDQpkYiAlPiUgZmlsdGVyX2FsbChhbnlfdmFycyguICVpbiUgYyg1NTEuMDgyKSkpDQpkYiAlPiUgZmlsdGVyX2FsbChhbnlfdmFycyguICVpbiUgYygzOTcuNDgxKSkpDQoNCmhpc3QoZGIkcG9wZGVuc2l0eV8yMDE5KQ0KYm94cGxvdChkYiRwb3BkZW5zaXR5XzIwMTksIGhvcml6b250YWwgPSBUUlVFKQ0KYm94cGxvdChkYiRwb3BkZW5zaXR5XzIwMTkscGxvdD1GQUxTRSkkb3V0DQpkYiAlPiUgZmlsdGVyX2FsbChhbnlfdmFycyguICVpbiUgYyg2MTQ1LjM3NCkpKQ0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoNTU0LjA4NikpKQ0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoMzk5LjYzNSkpKQ0KDQpoaXN0KGRiJGpvYmRlbnNpdHlfMjAxOCkNCmJveHBsb3QoZGIkam9iZGVuc2l0eV8yMDE4LCBob3Jpem9udGFsID0gVFJVRSkNCmJveHBsb3QoZGIkam9iZGVuc2l0eV8yMDE4LHBsb3Q9RkFMU0UpJG91dA0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoMjk5MC4zNTMpKSkNCmRiICU+JSBmaWx0ZXJfYWxsKGFueV92YXJzKC4gJWluJSBjKDM0MC40NTYpKSkNCmRiICU+JSBmaWx0ZXJfYWxsKGFueV92YXJzKC4gJWluJSBjKDE2Ni41NTgpKSkNCg0KaGlzdChkYiRqb2JkZW5zaXR5XzIwMTkpDQpib3hwbG90KGRiJGpvYmRlbnNpdHlfMjAxOSwgaG9yaXpvbnRhbCA9IFRSVUUpDQpib3hwbG90KGRiJGpvYmRlbnNpdHlfMjAxOSxwbG90PUZBTFNFKSRvdXQNCmRiICU+JSBmaWx0ZXJfYWxsKGFueV92YXJzKC4gJWluJSBjKDMwMjAuMjYyKSkpDQpkYiAlPiUgZmlsdGVyX2FsbChhbnlfdmFycyguICVpbiUgYygzNDAuMjI0KSkpDQpkYiAlPiUgZmlsdGVyX2FsbChhbnlfdmFycyguICVpbiUgYygxNjkuMjg4KSkpDQoNCmhpc3QoZGIkcGF0ZW50c3JhdGVfMjAxNikNCmJveHBsb3QoZGIkcGF0ZW50c3JhdGVfMjAxNiwgaG9yaXpvbnRhbCA9IFRSVUUpDQpib3hwbG90KGRiJHBhdGVudHNyYXRlXzIwMTYscGxvdD1GQUxTRSkkb3V0DQpkYiAlPiUgZmlsdGVyX2FsbChhbnlfdmFycyguICVpbiUgYygwLjA3MykpKQ0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoMC4wOTEpKSkNCmRiICU+JSBmaWx0ZXJfYWxsKGFueV92YXJzKC4gJWluJSBjKDAuMDY1KSkpDQoNCmhpc3QoZGIkcGF0ZW50c3JhdGVfMjAxNykNCmJveHBsb3QoZGIkcGF0ZW50c3JhdGVfMjAxNywgaG9yaXpvbnRhbCA9IFRSVUUpDQpib3hwbG90KGRiJHBhdGVudHNyYXRlXzIwMTcscGxvdD1GQUxTRSkkb3V0DQpkYiAlPiUgZmlsdGVyX2FsbChhbnlfdmFycyguICVpbiUgYygwLjA4OCkpKQ0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoMC4wNjkpKSkNCmRiICU+JSBmaWx0ZXJfYWxsKGFueV92YXJzKC4gJWluJSBjKDAuMDY0KSkpDQpkYiAlPiUgZmlsdGVyX2FsbChhbnlfdmFycyguICVpbiUgYygwLjA2MSkpKQ0KDQpoaXN0KGRiJHJlYWx3YWdlc18yMDE2KQ0KYm94cGxvdChkYiRyZWFsd2FnZXNfMjAxNiwgaG9yaXpvbnRhbCA9IFRSVUUpDQpib3hwbG90KGRiJHJlYWx3YWdlc18yMDE2LHBsb3Q9RkFMU0UpJG91dA0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoNDIzLjU2NykpKQ0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoNDQxLjg4MikpKQ0KDQpoaXN0KGRiJHJlYWx3YWdlc18yMDE3KQ0KYm94cGxvdChkYiRyZWFsd2FnZXNfMjAxNywgaG9yaXpvbnRhbCA9IFRSVUUpDQpib3hwbG90KGRiJHJlYWx3YWdlc18yMDE3LHBsb3Q9RkFMU0UpJG91dA0KZGIgJT4lIGZpbHRlcl9hbGwoYW55X3ZhcnMoLiAlaW4lIGMoNDI4LjM2MikpKQ0KDQpgYGANCg0KIyMjIEZpbmRpbmdzDQpNZXhpY28gQ2l0eSBkZXZpYXRlIG1hcmtlZGx5IGZyb20gb3RoZXIgb2JzZXJ2YXRpb25zLg0KDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmx1ZTsiPkNob3JvcGxldGggTWFwczwvc3Bhbj4gDQpgYGB7cn0NCmxpYnJhcnkoIm14bWFwcyIpDQpzdW1tYXJ5KGRmX214c3RhdGVfMjAyMCkNCmRmX214c3RhdGVfMjAyMCR2YWx1ZSA8LSBkYiRhdXRvbW90aXZlX2luZHVzdHJ5DQpteHN0YXRlX2Nob3JvcGxldGgoZGZfbXhzdGF0ZV8yMDIwLCB0aXRsZSA9ICJBdXRvbW90aXZlIFN1cHBseSBDaGFpbiBDb21wYW5pZXMgRGlzdHJpYnV0aW9uIGluIE1leGljbyIpIA0KYGBgDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmx1ZTsiPlNwYXRpYWwgRGF0YSBDb25zdHJ1Y3Rpb24gYW5kIE1hcHBpbmc8L3NwYW4+IA0KDQojIyMgSW5zdGFsbCBwYWNrYWdlcyBmb3Igc3BhdGlhbCBhbmFseXNpcw0KYGBge3J9DQpyZXF1aXJlKHByb2o0KSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJwcm9qNCIsIGRlcGVuZGVuY2llcyA9IFQpDQpyZXF1aXJlKHJnZGFsKSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJyZ2RhbCIsIGRlcGVuZGVuY2llcyA9IFQpDQpyZXF1aXJlKHNwKSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJzcCIsIGRlcGVuZGVuY2llcyA9IFQpDQpyZXF1aXJlKHNwZGVwKSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJzcGRlcCIsIGRlcGVuZGVuY2llcyA9IFQpDQpyZXF1aXJlKGFwZSkgfHwgaW5zdGFsbC5wYWNrYWdlcygiYXBlIiwgZGVwZW5kZW5jaWVzID0gVCkNCnJlcXVpcmUoZ3N0YXQpIHx8IGluc3RhbGwucGFja2FnZXMoImdzdGF0IiwgZGVwZW5kZW5jaWVzID0gVCkNCnJlcXVpcmUoTUFTUykgfHwgaW5zdGFsbC5wYWNrYWdlcygiTUFTUyIsIGRlcGVuZGVuY2llcyA9IFQpDQpyZXF1aXJlKGRldnRvb2xzKSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJkZXZ0b29scyIpDQpyZXF1aXJlKHNwZ3dyKSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJzcGd3ciIsIGRlcGVuZGVuY2llcyA9IFQpDQpyZXF1aXJlKHNwRGF0YSkgfHwgaW5zdGFsbC5wYWNrYWdlcygic3BEYXRhIiwgZGVwZW5kZW5jaWVzID0gVCkNCnJlcXVpcmUoTWF0cml4KSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJNYXRyaXgiLCBkZXBlbmRlbmNpZXMgPSBUKQ0KcmVxdWlyZShwc3ljaCkgfHwgaW5zdGFsbC5wYWNrYWdlcygicHN5Y2giLCBkZXBlbmRlbmNpZXMgPSBUKQ0KcmVxdWlyZSh0bWFwKSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJ0bWFwIikNCnJlcXVpcmUocmVhZHIpIHx8IGluc3RhbGwucGFja2FnZXMoInJlYWRyIikNCnJlcXVpcmUoc3RhcmdhemVyKSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJzdGFyZ2F6ZXIiKQ0KcmVxdWlyZShhdXRvbWFwKSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJhdXRvbWFwIikNCnJlcXVpcmUoc3BhdGlhbEVjbykgfHwgaW5zdGFsbC5wYWNrYWdlcygic3BhdGlhbEVjbyIpDQpyZXF1aXJlKFNwYXRpYWxLREUpIHx8IGluc3RhbGwucGFja2FnZXMoIlNwYXRpYWxLREUiKQ0KcmVxdWlyZShzbWFjcG9kKSB8fCBpbnN0YWxsLnBhY2thZ2VzKCJzbWFjcG9kIikNCnJlcXVpcmUobG1lNCkgfHwgaW5zdGFsbC5wYWNrYWdlcygibG1lNCIpDQoNCmBgYA0KDQojIyMgUHJlcGFyZSB0aGUgZGF0YQ0KYGBge3J9DQpzZXR3ZCgiRjpcXFRFQ1xcMjAyMy0wMS0wOSBMdW5lcyBhIDEzIFZpZXJuZXMgQ3Vyc28gQW7DoWxpc2lzIEVzcGFjaWFsIGRlIERhdG9zIHVzYW5kbyBSXFxQcm95ZWN0b1xcZGF0YWJhc2VfbmVhcnNob3JpbmciKQ0KZ2V0d2QoKQ0KbGlicmFyeShtYXB0b29scykNCmxpYnJhcnkoc3ApDQpsaWJyYXJ5KHJlYWRyKQ0KbmVhcnNob3JpbmdfbWFwMiA8LSByZWFkU2hhcGVQb2x5KCJuZWFyc2hvcmluZ19tYXAuc2hwIikgIA0KbmVhcnNob3JpbmcuY3N2IDwtIHJlYWRfY3N2KCJuZWFyc2hvcmluZy5jc3YiKSAgDQpuZWFyc2hvcmluZ19tYXAuc3AgPC0gbWVyZ2UobmVhcnNob3JpbmdfbWFwMixuZWFyc2hvcmluZy5jc3YpICAgIA0Kc2F2ZShuZWFyc2hvcmluZ19tYXAuc3AsIGZpbGU9Im5lYXJzaG9yaW5nX21hcC5SRGF0YSIpICANCmxpYnJhcnkocmdkYWwpDQpzdCA8LSByZWFkT0dSKCJuZWFyc2hvcmluZ19tYXAuc2hwIikgICANCnN0MSA8LSBkYXRhLmZyYW1lKHN0KSAgICAgDQpsaWJyYXJ5KHNmKQ0KY2MgPSBzdF9yZWFkKCJuZWFyc2hvcmluZ19tYXAuc2hwIikgDQpzdW1tYXJ5KGRiKSAgIA0Kc3VtbWFyeShuZWFyc2hvcmluZ19tYXAuc3ApIA0Kc3VtbWFyeShjYykNCmRiX2RhdGEgPC0gZGF0YS5mcmFtZShkYikNCm5lYXJzaG9yaW5nX21hcF9kYXRhIDwtIGRhdGEuZnJhbWUobmVhcnNob3JpbmdfbWFwLnNwKQ0KYGBgDQoNCiMjIyBDcmVhdGUgbWFwDQpgYGB7cn0NCnBsb3QobmVhcnNob3JpbmdfbWFwMiwgcGJnPSJibHVlIiwgYXhlcyA9IFQpICANCnRpdGxlKCJTcGF0aWFsIHBvbHlnb24gcGxvdCIpICAgDQpgYGANCg0KIyMjIE1hcHBpbmcgZm9yIHBvaW50cyBmb3IgdGhlIHNhbWUgYXJlYQ0KYGBge3J9DQpmaWxlLmNob29zZSgpDQojYXNzZW1ibHlfcGxhbnRzIDwtIHJlYWQuY3N2KEY6XFxURUNcXDIwMjMtMDEtMDkgTHVuZXMgYSAxMyBWaWVybmVzIEN1cnNvIEFuw6FsaXNpcyBFc3BhY2lhbCBkZSBEYXRvcyB1c2FuZG8gUlxcUHJveWVjdG9cXGFzc2VtYmx5X3BsYW50cy5jc3YpDQojcGxvdChhc3NlbWJseV9wbGFudHMsIGF4ZXM9VCwgbHdkID0gMSwgcGNoPTE2LCBjZXg9MC4xKSAgDQojdGl0bGUoIlNwYXRpYWwgUG9pbnRzIFBsb3QiKSANCmBgYA0KDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmx1ZTsiPkdlb2NvZGluZyB1c2luZyBPcGVuIEFQSTwvc3Bhbj4gDQoNCiMjIyBJbnN0YWxsIHBhY2thZ2UgdG8gZ2VuZXJhdGUgZ2VvY29vcmRpbmF0ZXMgdXNpbmcgR29vZ2xlIE1hcCBBUEkNCmBgYHtyfQ0KbGlicmFyeShkcGx5cikgIyBmb3IgZGF0YSBtYW5pcHVsYXRpb24NCmxpYnJhcnkoZ2dwbG90MikgIyBmb3IgdmlzdWFsaXppYXRpb24NCmxpYnJhcnkobWFwdG9vbHMpICMgZm9yIG1hcHBpbmcgaW4gUg0KbGlicmFyeShnZ21hcCkgIyBmb3IgdXNpbmcgR29vZ2xlIE1hcCB3aXRoIEFQSSBrZXkNCmBgYA0KDQojIyMgR29vZ2xlJ3MgR2VvY29kaW5nIEFQSQ0KYGBge3J9DQpyZWdpc3Rlcl9nb29nbGUoa2V5ID0gIkFJemFTeURKY3FvejlmZjIzVmxoYXpReVp0aEZwcnVpT3BJVlRIRSIpDQpgYGANCg0KIyMjIENvbnZlcnQgYWRkcmVzcyBvZiBBdWRpIEFzc2VtYmx5IFBsYW50IGluIE1leGljbyBpbnRvIGNvb3JkaW5hdGVzIGZvciBtYXBwaW5nDQpgYGB7cn0NCmF1ZGlfcHVlYmxhLmxvY2F0aW9uIDwtIGRhdGEuZnJhbWUoYWRkciA9IGMoIjc1MDEyIFNhbiBKb3PDqSBDaGlhcGEsIFB1ZS4iKSwgc3RyaW5nc0FzRmFjdG9ycyA9IEZBTFNFKQ0KYXVkaV9wdWVibGEubGF0bG9uIDwtIGF1ZGlfcHVlYmxhLmxvY2F0aW9uICU+JSBtdXRhdGVfZ2VvY29kZShhZGRyKQ0KYXVkaV9wdWVibGEubGF0bG9uDQoNCmJtd19zbHAubG9jYXRpb24gPC0gZGF0YS5mcmFtZShhZGRyID0gYygiQm91bGV2YXJkIEJNVyBOby4gNjU1IFBhcnF1ZSBJbmR1c3RyaWFsIERlc2Fycm9sbG8gTG9naXN0aWsgSUksIDc5NTI2IFMuTC5QLiIpLCBzdHJpbmdzQXNGYWN0b3JzID0gRkFMU0UpDQpibXdfc2xwLmxhdGxvbiA8LSBibXdfc2xwLmxvY2F0aW9uICU+JSBtdXRhdGVfZ2VvY29kZShhZGRyKQ0KYm13X3NscC5sYXRsb24NCg0KYGBgDQoNCiMjIyBNYXJrIHRoZSBsb2NhdGlvbiBvbiB0aGUgbWFwIHVzaW5nIGNvb3JkaW5hdGVzDQpgYGB7cn0NCmF1ZGlfcHVlYmxhLmNlbiA8LSBjKG1lYW4oYXVkaV9wdWVibGEubGF0bG9uJGxvbiksIG1lYW4oYXVkaV9wdWVibGEubGF0bG9uJGxhdCkpDQphdWRpX3B1ZWJsYSA8LSBhdWRpX3B1ZWJsYS5sYXRsb25bLDI6M10NCg0KYm13X3NscC5jZW4gPC0gYyhtZWFuKGJtd19zbHAubGF0bG9uJGxvbiksIG1lYW4oYm13X3NscC5sYXRsb24kbGF0KSkNCmJtd19zbHAgPC0gYm13X3NscC5sYXRsb25bLDI6M10NCmBgYA0KDQojIyMgQ3JlYXRlIGEgcm9hZCBtYXAgb2YgQXVkaSBBc3NlbWJseSBQbGFudCBpbiBNZXhpY28NCmBgYHtyfQ0KYXVkaV9wdWVibGEucm9hZF9tYXAgPC0gZ2V0X2dvb2dsZW1hcChjZW50ZXIgPSBhdWRpX3B1ZWJsYS5jZW4sIG1hcHR5cGUgPSAicm9hZG1hcCIsIHpvb209MTIsIG1hcmtlcnMgPSBhdWRpX3B1ZWJsYSkgDQpnZ21hcChhdWRpX3B1ZWJsYS5yb2FkX21hcCkNCmBgYA0KDQoNCiMjIyBDb252ZXJ0IGEgbnVtYmVyIG9mIGFkZHJlc3NlcyBvZiBBc3NlbWJseSBQbGFudHMgaW4gTWV4aWNvIGludG8gY29vcmRpbmF0ZXMgdXNpbmcgY3N2IGZpbGUNCmBgYHtyfQ0Kc2V0d2QoIkY6XFxURUNcXDIwMjMtMDEtMDkgTHVuZXMgYSAxMyBWaWVybmVzIEN1cnNvIEFuw6FsaXNpcyBFc3BhY2lhbCBkZSBEYXRvcyB1c2FuZG8gUlxcUHJveWVjdG9cXGRhdGFiYXNlX25lYXJzaG9yaW5nIikNCmdldHdkDQpgYGANCg0KIyMjIExvYWQgQXNzZW1ibHkgUGxhbnRzIGluIE1leGljbyBBZGRyZXNzIERhdGENCmBgYHtyfQ0KYXNzZW1ibHlfcGxhbnRzLmRhdGEgPC0gcmVhZC5jc3YoJ2Fzc2VtYmx5X3BsYW50cy5jc3YnLCAgaGVhZGVyPVQsIGVuY29kaW5nID0gJ3V0Zi04JyApICAgIA0KbmFtZXMoYXNzZW1ibHlfcGxhbnRzLmRhdGEpIA0KYGBgDQoNCiMjIyBDcmVhdGUgZGF0YSBmcmFtZSBmb3IgZ2VvY29kaW5nDQpgYGB7cn0NCmFzc2VtYmx5X3BsYW50cy5kYXRhJGFkZHIgPC0gZW5jMnV0ZjgoYXNzZW1ibHlfcGxhbnRzLmRhdGEkYWRkcikNCmBgYA0KDQojIyMgR2VvY29kZSB0aGUgQXNzZW1ibHkgUGxhbnRzIGFkZHJlc3NlcyBpbiBNZXhpY28gDQpgYGB7cn0NCmFzc2VtYmx5X3BsYW50cy5sYXRsb24gPC0gYXNzZW1ibHlfcGxhbnRzLmRhdGEgJT4lIG11dGF0ZV9nZW9jb2RlKGFkZHIpIA0KaGVhZChhc3NlbWJseV9wbGFudHMubGF0bG9uKQ0KYGBgDQoNCiMjIyBDcmVhdGUgYSBiYXNlIG1hcA0KYGBge3J9DQpiYXNlbWFwIDwtIHFtYXAoJ21leGljbycsIHpvb20gPSA0LCBtYXB0eXBlID0gInJvYWRtYXAiKQ0KYGBgDQoNCiMjIyBNYXJrIHRoZSBsb2NhdGlvbnMgb2YgQXNzZW1ibHkgUGxhbnRzIG9uIHRoZSBiYXNlIG1hcA0KYGBge3J9DQpiYXNlbWFwICsgZ2VvbV9wb2ludChkYXRhID0gYXNzZW1ibHlfcGxhbnRzLmxhdGxvbiwgYWVzKGxvbiwgbGF0KSwgc2l6ZSA9IDEuNSwgY29sb3VyPSAicmVkIikNCmBgYA0KDQojIyA8c3BhbiBzdHlsZT0iY29sb3I6IGJsdWU7Ij5JZGVudGlmeSBzcGF0aWFsIGRpc3RyaWJ1dGlvbiBvZiB2YXJpYWJsZXMgb2YgaW50ZXJlc3Q8L3NwYW4+IA0KDQojIyMgTG9hZGluZyBMaWJyYXJpZXMgDQpgYGB7ciBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQ0KbGlicmFyeShnZ3Bsb3QyKQ0KbGlicmFyeShzcGRlcCkNCmxpYnJhcnkoTUFTUykNCmxpYnJhcnkoc3Btb3JhbikNCmxpYnJhcnkoc3BhdGlhbHJlZykNCmxpYnJhcnkoc3BoZXQpDQpsaWJyYXJ5KG1hcHRvb2xzKQ0KbGlicmFyeShzcCkNCmxpYnJhcnkoc2YpDQpsaWJyYXJ5KG1hcHMpDQpsaWJyYXJ5KHJnZW9zKQ0KbGlicmFyeShnZ21hcCkNCmxpYnJhcnkobWFwcHJvaikNCmxpYnJhcnkoUkNvbG9yQnJld2VyKQ0KbGlicmFyeShnZ3NuKQ0KDQojIyMgaW1wb3J0aW5nIHNwYXRpYWwgZGF0YQ0KbWFwPC1yZWFkU2hhcGVQb2x5KCJGOlxcVEVDXFwyMDIzLTAxLTA5IEx1bmVzIGEgMTMgVmllcm5lcyBDdXJzbyBBbsOhbGlzaXMgRXNwYWNpYWwgZGUgRGF0b3MgdXNhbmRvIFJcXFByb3llY3RvXFxkYXRhYmFzZV9uZWFyc2hvcmluZ1xcbmVhcnNob3JpbmdfbWFwLnNocCIsSUR2YXI9ImlkIixwcm9qNHN0cmluZz1DUlMoIitwcm9qPWxvbmdsYXQiKSkNCmxtYXQ8LWNvb3JkaW5hdGVzKG1hcCkNCm1hcC5jZW50cm9pZDwtY29vcmRpbmF0ZXMobWFwKSAgDQpzdW1tYXJ5KG1hcCkNCg0KIyMjIGNyZWF0ZSBzcGF0aWFsIGNvbm5lY3Rpdml0eSBtYXRyaXggKHN3bSkNCm1hcC5saW5rPC1wb2x5Mm5iKG1hcCxxdWVlbj1UKSAgICAgICAgICAgICAgIyMjIHNldCBTV00gdG8gcXVlZW4gDQptYXAubGlua1c8LW5iMmxpc3R3KG1hcC5saW5rLCBzdHlsZT0iVyIpICAgICAjIyMgU1dNIHF1ZWVuIHN0YW5kYXJkaXplZCANCnBsb3QobWFwLGJvcmRlcj0iYmx1ZSIsYXhlcz1GQUxTRSxsYXM9MSkNCnBsb3QobWFwLGNvbD0iZ3JleSIsYm9yZGVyPWdyZXkoMC45KSxheGVzPVQsYWRkPVQpIA0KcGxvdChtYXAubGlua1csY29vcmRzPW1hcC5jZW50cm9pZCxwY2g9MTksY2V4PTAuMSxjb2w9InJlZCIsYWRkPVQpICAjIyMgcGxvdCBTV00gcXVlZW4gc3RhbmRhcml6ZWQgDQp0aXRsZSgiU3BhdGlhbCBDb25uZWN0aXZpdHkgTWF0cml4IC0gQ29udGlndWl0eSBDYXNlIikNCg0KIyMjIFRlc3RpbmcgZm9yIFNwYXRpYWwgQXV0b2NvcnJlbGF0aW9uIHVzaW5nIE1vcmFuJ3MgSSANCiMgZGV0ZWN0aW5nIHNwYXRpYWwgYXV0b2NvcnJlbGF0aW9uIGluIG1haW4gdmFyaWFibGUgb2YgaW50ZXJlc3QgDQptb3Jhbi5tYyhtYXAkZmRpXzIwMjAsbWFwLmxpbmtXLG5zaW09OTk5OSkgIyMjIG5vIHN0YXRpc3RpY2FsbHkgc2lnbmlmaWNhbnQgDQptb3Jhbi5tYyhtYXAkZmRpXzIwMjEsbWFwLmxpbmtXLG5zaW09OTk5OSkgIyMjIG5vIHN0YXRpc3RpY2FsbHkgc2lnbmlmaWNhbnQgDQptb3Jhbi5tYyhtYXAkZmRpXzIwMjIsbWFwLmxpbmtXLG5zaW09OTk5OSkgIyMjIG5vIHN0YXRpc3RpY2FsbHkgc2lnbmlmaWNhbnQgDQptb3Jhbi5tYyhtYXAkZmRpXzE5OTlfMixtYXAubGlua1csbnNpbT05OTk5KSAjIyMgbm8gc3RhdGlzdGljYWxseSBzaWduaWZpY2FudCANCg0KIyBkZXRlY3Rpbmcgc3BhdGlhbCBhdXRvY29ycmVsYXRpb24gaW4gY29udHJvbCB2YXJpYWJsZSAtIGpvYiBkZW5zaXR5IA0KbW9yYW4ubWMobWFwJGpvYmRlbnNpdHksbWFwLmxpbmtXLG5zaW09OTk5OSkgIyMjIHN0YXRpc3RpY2FsbHkgc2lnbmlmaWNhbnQgDQptb3Jhbi5tYyhtYXAkam9iZGVuc2lfMSxtYXAubGlua1csbnNpbT05OTk5KSAjIyMgc3RhdGlzdGljYWxseSBzaWduaWZpY2FudCANCmBgYA0KDQojIyMgZGV0ZWN0aW5nIHNwYXRpYWwgYXV0b2NvcnJlbGF0aW9uIGluIGNvbnRyb2wgdmFyaWFibGUgLSBjcmltZSByYXRlDQpgYGB7cn0NCm1vcmFuLm1jKG1hcCRjcmltZV9yYXRlLG1hcC5saW5rVyxuc2ltPTk5OTkpICMjIyBzdGF0aXN0aWNhbGx5IHNpZ25pZmljYW50IA0KbW9yYW4ubWMobWFwJGNyaW1lX3JhXzEsbWFwLmxpbmtXLG5zaW09OTk5OSkgIyMjIHN0YXRpc3RpY2FsbHkgc2lnbmlmaWNhbnQgDQpgYGANCg0KDQojIyA8c3BhbiBzdHlsZT0iY29sb3I6IGJsdWU7Ij5UZXN0aW5nIGZvciBTcGF0aWFsIEF1dG9jb3JyZWxhdGlvbiB1c2luZyBNb3JhbuKAmXMgSTwvc3Bhbj4gDQpgYGB7cn0NCiMgTGV0cyB1c2UgR2VvRGEgdG8gY29uZHVjdCBjbHVzdGVyIGFuYWx5c2lzIA0KIyBpbnN0YWxsLnBhY2thZ2VzKCJyZ2VvZGEiKQ0KbGlicmFyeShyZ2VvZGEpDQptYXBfc2hhcGVmaWxlPC1zdF9yZWFkKCJGOlxcVEVDXFwyMDIzLTAxLTA5IEx1bmVzIGEgMTMgVmllcm5lcyBDdXJzbyBBbsOhbGlzaXMgRXNwYWNpYWwgZGUgRGF0b3MgdXNhbmRvIFJcXFByb3llY3RvXFxkYXRhYmFzZV9uZWFyc2hvcmluZ1xcbmVhcnNob3JpbmdfbWFwLnNocCIpIA0KDQojIGNyZWF0ZSBzcGF0aWFsIHdlaWdodHMgbWF0cml4IG9iamVjdCAoY29udGlndWl0eSBjYXNlKQ0KcXVlZW5fdzwtcXVlZW5fd2VpZ2h0cyhtYXBfc2hhcGVmaWxlKQ0KDQojIGNhbGN1bGF0ZSBMSVNBIGFzIHBlciBHRU9EQSANCmxpc2FfZHY8LWxvY2FsX21vcmFuKHF1ZWVuX3csIG1hcF9zaGFwZWZpbGVbImZkaV8yMDIyIl0pIA0KbGlzYV9jdjE8LWxvY2FsX21vcmFuKHF1ZWVuX3csIG1hcF9zaGFwZWZpbGVbImpvYmRlbnNpdHkiXSkgDQpsaXNhX2N2MjwtbG9jYWxfbW9yYW4ocXVlZW5fdywgbWFwX3NoYXBlZmlsZVsicGF0ZW50c3JhdCJdKSANCmxpc2FfY3YzPC1sb2NhbF9tb3JhbihxdWVlbl93LCBtYXBfc2hhcGVmaWxlWyJjcmltZV9yYXRlIl0pIA0KDQojIGFkZCBjbHVzdGVyIGNhbGN1bGF0aW9ucyB0byBkYXRhc2V0IA0KbWFwX3NoYXBlZmlsZSRjbHVzdGVyX2R2PC1hcy5mYWN0b3IobGlzYV9kdiRHZXRDbHVzdGVySW5kaWNhdG9ycygpKQ0KbGV2ZWxzKG1hcF9zaGFwZWZpbGUkY2x1c3Rlcl9kdik8LWxpc2FfZHYkR2V0TGFiZWxzKCkgIyMjIGRlcGVuZGVudCB2YXJpYWJsZQ0KDQptYXBfc2hhcGVmaWxlJGNsdXN0ZXJfY3YxPC1hcy5mYWN0b3IobGlzYV9jdjEkR2V0Q2x1c3RlckluZGljYXRvcnMoKSkNCmxldmVscyhtYXBfc2hhcGVmaWxlJGNsdXN0ZXJfY3YxKTwtbGlzYV9jdjEkR2V0TGFiZWxzKCkgIyMjIGpvYiBkZW5zaXR5DQoNCm1hcF9zaGFwZWZpbGUkY2x1c3Rlcl9jdjI8LWFzLmZhY3RvcihsaXNhX2N2MiRHZXRDbHVzdGVySW5kaWNhdG9ycygpKQ0KbGV2ZWxzKG1hcF9zaGFwZWZpbGUkY2x1c3Rlcl9jdjIpPC1saXNhX2N2MiRHZXRMYWJlbHMoKSAjIyMgcGF0ZW50cyByYXRlDQoNCm1hcF9zaGFwZWZpbGUkY2x1c3Rlcl9jdjM8LWFzLmZhY3RvcihsaXNhX2N2MyRHZXRDbHVzdGVySW5kaWNhdG9ycygpKQ0KbGV2ZWxzKG1hcF9zaGFwZWZpbGUkY2x1c3Rlcl9jdjMpPC1saXNhX2N2MyRHZXRMYWJlbHMoKSAjIyMgY3JpbWUgcmF0ZSANCg0KIyB2aXN1YWxpemUgY2x1c3RlcnMgdXNpbmcgR2VvRGEgdmlldyANCmdncGxvdChkYXRhPW1hcF9zaGFwZWZpbGUpICsNCiAgZ2VvbV9zZihhZXMoZmlsbD1jbHVzdGVyX2R2KSkgKyANCiAgZ2d0aXRsZShsYWJlbCA9ICJGb3JlaWduIERpcmVjdCBJbnZlc3RtZW50IDIwMjIiLCBzdWJ0aXRsZSA9ICJOdW1iZXIgb2YgQ29tcGFuaWVzIC0gUmVjaXBpZW50cyBvZiBGREkncyBJbmZsb3dzIikNCg0KZ2dwbG90KGRhdGE9bWFwX3NoYXBlZmlsZSkgKw0KICBnZW9tX3NmKGFlcyhmaWxsPWNsdXN0ZXJfY3YxKSkgKyANCiAgZ2d0aXRsZShsYWJlbCA9ICJKb2IgRGVuc2l0eSAyMDE4Iiwgc3VidGl0bGUgPSAiUG9wdWxhdGlvbiBwZXIga20yIikNCg0KZ2dwbG90KGRhdGE9bWFwX3NoYXBlZmlsZSkgKw0KICBnZW9tX3NmKGFlcyhmaWxsPWNsdXN0ZXJfY3YyKSkgKyANCiAgZ2d0aXRsZShsYWJlbCA9ICJSJkQgQWN0aXZpdHkgMjAxNiIsIHN1YnRpdGxlID0gIlBhdGVudHMgcGVyIDEwMCwwMDAgUG9wdWxhdGlvbiIpDQoNCmdncGxvdChkYXRhPW1hcF9zaGFwZWZpbGUpICsNCiAgZ2VvbV9zZihhZXMoZmlsbD1jbHVzdGVyX2N2MykpICsgDQogIGdndGl0bGUobGFiZWwgPSAiQ3JpbWUgMjAxOCIsIHN1YnRpdGxlID0gIkhvbWljaWRlIFJhdGUgcGVyIDEwMCwwMDAgUG9wdWxhdGlvbiIpDQpgYGANCg0KYGBge3J9DQojIyMgQml2YXJpYXRlIE1vcmFuJ3MgSSANCmxpc2FfYmkgPC0gbG9jYWxfYmltb3JhbihxdWVlbl93LG1hcF9zaGFwZWZpbGVbYygnZmRpXzIwMjInLCdqb2JkZW5zaV8xJyldKQ0KYml2YXJpYXRlX20gPC0gbGlzYV92YWx1ZXMobGlzYV9iaSkNCmJpdmFyaWF0ZV9tDQoNCm1hcF9zaGFwZWZpbGUkYml2YXJpYXRlPC1hcy5mYWN0b3IobGlzYV9iaSRHZXRDbHVzdGVySW5kaWNhdG9ycygpKQ0KbGV2ZWxzKG1hcF9zaGFwZWZpbGUkYml2YXJpYXRlKTwtbGlzYV9iaSRHZXRMYWJlbHMoKSAjIyMgYml2YXJpYXRlIG1vcmFuDQoNCiMgdmlzdWFsaXplIEJpdmFyaWF0ZSBNb3JhbidzIGkgdXNpbmcgR2VvRGEgdmlldyANCmdncGxvdChkYXRhPW1hcF9zaGFwZWZpbGUpICsNCiAgZ2VvbV9zZihhZXMoZmlsbD1iaXZhcmlhdGUpKSArIA0KICBnZ3RpdGxlKGxhYmVsID0gIkZvcmVpZ24gRGlyZWN0IEludmVzdG1lbnQgJiBKb2IgRGVuc2l0eSIsIHN1YnRpdGxlID0gIkJpdmFyaWF0ZSBMb2NhbCdzIE1vcmFuJ3MgSSIpDQpgYGANCg0KDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmx1ZTsiPlN1Z2dlc3QgaHlwb3RoZXNlcyBhYm91dCBwb3RlbnRpYWwgcmVsYXRpb25zaGlwcyBhbW9uZyB0aGUgdmFyaWFibGVzIG9mIGludGVyZXN0PC9zcGFuPiANCg0KIyMjIEh5cG90aGVzaXMgMTogVGhlcmUgaXMgYSBwb3NpdGl2ZSByZWxhdGlvbnNoaXAgYmV0d2VlbiBmZGkgaW5mbG93cyBhbmQgam9iIGRlbnNpdHkgDQojIyMgSHlwb3RoZXNpcyAyOiBUaGVyZSBpcyBhIHBvc2l0aXZlIHJlbGF0aW9uc2hpcCBiZXR3ZWVuIGZkaSBpbmZsb3dzIGFuZCBSJkQgYWN0aXZpdHkgDQojIyMgSHlwb3RoZXNpcyAzOiBUaGVyZSBpcyBhIG5lZ2F0aXZlIHJlbGF0aW9uc2hpcCBiZXR3ZWVuIGZkaSBpbmZsb3dzIGFuZCBjcmltZSByYXRlIA0KDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmx1ZTsiPlNwYXRpYWwgcmVncmVzc2lvbiBhbmFseXNpcyBhbmQgaW50ZXJwcmV0YXRpb24gb2YgZXN0aW1hdGVkIHJlc3VsdHM8L3NwYW4+DQpgYGB7cn0NCiMgbGV0cyBpbXBvcnQgZGF0YXNldCANCmRhdGE9cmVhZC5jc3YoIkY6XFxURUNcXDIwMjMtMDEtMDkgTHVuZXMgYSAxMyBWaWVybmVzIEN1cnNvIEFuw6FsaXNpcyBFc3BhY2lhbCBkZSBEYXRvcyB1c2FuZG8gUlxcUHJveWVjdG9cXGRhdGFiYXNlX25lYXJzaG9yaW5nXFxuZWFyc2hvcmluZy5jc3YiKQ0Kc3VtbWFyeShkYXRhKQ0KDQojIGxhZ19tb2RlbA0KbW9kZWwxPC1sYWdzYXJsbShsb2coZmRpXzIwMjIpfmxvZyhqb2JkZW5zaXR5XzIwMTkpK3BhdGVudHNyYXRlXzIwMTcrbG9nKGNyaW1lX3JhdGVfMjAxOSkrbG9nKHJlYWx3YWdlc18yMDE3KStsb2coYm9yZGVyX2Rpc3RhbmNlKSwNCiAgICAgICAgICAgICAgICAgICBkYXRhPWRhdGEsbGlzdHc9bWFwLmxpbmtXLCB6ZXJvLnBvbGljeT1UUlVFLCBuYS5hY3Rpb249bmEub21pdCkNCiMgc3VtbWFyeShtb2RlbDEpDQoNCiMgZXJyb3JfbW9kZWwNCm1vZGVsMjwtZXJyb3JzYXJsbShsb2coZmRpXzIwMjIpfmxvZyhqb2JkZW5zaXR5XzIwMTkpK3BhdGVudHNyYXRlXzIwMTcrbG9nKGNyaW1lX3JhdGVfMjAxOSkrbG9nKHJlYWx3YWdlc18yMDE3KStsb2coYm9yZGVyX2Rpc3RhbmNlKSwNCiAgICAgICAgICAgZGF0YT1kYXRhLGxpc3R3PW1hcC5saW5rVywgemVyby5wb2xpY3k9VFJVRSwgbmEuYWN0aW9uPW5hLm9taXQpDQojIHN1bW1hcnkobW9kZWwyKQ0KDQojIGR1cmJpbl9tb2RlbA0KbW9kZWwzPC1sYWdzYXJsbShsb2coZmRpXzIwMjIpfmxvZyhqb2JkZW5zaXR5XzIwMTkpK3BhdGVudHNyYXRlXzIwMTcrbG9nKGNyaW1lX3JhdGVfMjAxOSkrbG9nKHJlYWx3YWdlc18yMDE3KStsb2coYm9yZGVyX2Rpc3RhbmNlKSwgZGF0YT1kYXRhLCBsaXN0dz1tYXAubGlua1csIER1cmJpbj1UUlVFKQ0KIyBzdW1tYXJ5KG1vZGVsMykNCg0KIyBjb21wYXJlIHNwYXRpYWwgcmVncmVzc2lvbiByZXN1bHRzIA0KIyAgaW5zdGFsbC5wYWNrYWdlcygiaHV4dGFibGUiKQ0KbGlicmFyeShqdG9vbHMpDQpqdG9vbHM6OmV4cG9ydF9zdW1tcyhtb2RlbDEsIG1vZGVsMiwgbW9kZWwzKQ0KYGBgDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmx1ZTsiPkNyZWF0ZSBpbnRlcmFjdGl2ZSBtYXA8L3NwYW4+IA0KDQoNCg0KIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBibHVlOyI+Q29uY2x1c2lvbnM8L3NwYW4+IA0KDQo=