UNIVERSIDAD DE EL SALVADOR
FACULTAD DE CIENCIAS ECONOMICAS
ESCUELA DE ECONOMIA
Materia: Métodos para el análisis económico
Docente: Carlos Ademir Pérez Alas
Grupo teórico: 3
Laboratorio 1: Segunda parte
Integrantes:
Apellidos
Nombres
DUE
Escolán Chávez
Diana Alejandra
EC19003
García Cabrera
Bryan Alexander
GC17011
Ciclo: II-2021
Fecha: Vienes 10 de septiembre, 2021
Ciudad Universitaria, San Salvador, El Salvador
Carga de la Base de Comercio Exterior
library(dplyr)
library(kableExtra)
load("C:/Archivos R/MAE118/base_comercio.RData")
data_comercio_exterior <- base_comercio
1) Reproduciendo Funciones de Comercio Exterior
Indicación: Reproduzca los ejemplos vistos en clase sobre el uso de la librería dplyr, usando la información actualizada recopilada en la tarea 1, y usando los nombres estandarizados que se han agregado para los países, utilice países diferentes a los que se emplearon en los ejemplos (queda a elección del grupo los paises).
# Carga de funciones
source(file = "C:/Archivos R/MAE118/funciones_comercio_exterior.R",echo = FALSE,max.deparse.length = 100,encoding = "utf-8")
Consulta de países
library(kableExtra)
herramientas_mostrar_paises() %>% tail() %>%
kable(caption = "Países en la Base de Comercio Exterior 2018-2021 (Julio)",
align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
Países en la Base de Comercio Exterior 2018-2021 (Julio)
pais
codigo_pais
region
cod_region
region_intermedia
cod_region_intermedia
Zimbabwe
716
África
2
África Oriental
14
Zimbabwe
716
África
2
África Oriental
14
Zimbabwe
716
África
2
África Oriental
14
Zimbabwe
716
África
2
África Oriental
14
Zimbabwe
716
África
2
África Oriental
14
Zimbabwe
716
África
2
África Oriental
14
a Fuente: Elaboración propia, con base en información del BCR.
Consulta del número de transacciones por país
herramientas_mostrar_paises_en_tabla() %>%
tail() %>% kable(caption = "Países con transacciones 2018-2021 (Julio)",
align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
Países con transacciones 2018-2021 (Julio)
pais
anio
# Transacciones
Yugoeslavia -Nueva Organizacio
2020
1
Zambia
2021
4
Zimbabwe
2018
1
Zimbabwe
2019
8
Zimbabwe
2020
3
Zimbabwe
2021
1
a Fuente: Elaboración propia, con base en información del BCR.
Consulta de regiones en tabla
herramientas_mostrar_regiones() %>%
tail() %>% kable(caption = "Regiones en la Base de Comercio Exterior 2018-2021 (Julio)",
align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
Regiones en la Base de Comercio Exterior 2018-2021 (Julio)
region
cod_region
region_intermedia
cod_region_intermedia
# de Paises incluidos
Américas
19
NA
NA
132225
Asia
142
NA
NA
250262
Europa
150
Islas Del Canal
830
7
Europa
150
NA
NA
170958
Oceanía
9
NA
NA
4139
NA
NA
NA
NA
229
a Fuente: Elaboración propia, con base en información del BCR.
Indicador de Balassa
Por capítulo
Se analiza el IB con Alemania DEU, para el periodo 2018-2021, para el producto “1905200000”
Resultados_IB_sac <- data.frame("Años"=2018:2021,
"IB"=sapply(X=2018:2021,
FUN = indicadores_Balassa_partida,
codigo_pais=276,
partida="1905200000",
normalizado=TRUE))
Resultados_IB_sac %>% kable(caption = "IB capítulo: Alemania 2018-2021 (Julio)", align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
IB capítulo: Alemania 2018-2021 (Julio)
Años
IB
2018
0.5560822
2019
0.5561221
2020
0.6081481
2021
0.6866135
a Fuente: Elaboración propia, con base en información del BCR.
Por partida
Se analizará el IB con Alemania DEU, para el periodo 2018-2021, para el capítulo “02”
Resultados_Balassa <- data.frame("Años"=2018:2021,
"IB"=sapply(X=2018:2021,
FUN = indicadores_Balassa_capitulo,
codigo_pais=276,
capitulo="02"))
Resultados_IB_sac %>% kable(caption = "IB capítulo: Alemania 2018-2021 (Julio)",
align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia,con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
IB capítulo: Alemania 2018-2021 (Julio)
Años
IB
2018
0.5560822
2019
0.5561221
2020
0.6081481
2021
0.6866135
a Fuente: Elaboración propia,con base en información del BCR.
Uso de Indicador de Ventaja Comparativa Revelada
Por capítulo
Se analiza el IVCR con Argentina ARG, para el periodo 2018-2021, para al capitulo “01”
options(scipen = 99999)
Resultados_IVCR<-data.frame("Años"=2018:2021,
"IVCR"=sapply(X=2018:2021,
FUN = indicadores_IVCR_capitulo,
codigo_pais=32,
capitulo="01"))
Resultados_IVCR %>% kable(caption = "Indicador de Ventaja Comparativa Revelada (capítulo): Argentina 2018-2021 (Julio)",
align = "c") %>% add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>%
kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
Indicador de Ventaja Comparativa Revelada (capítulo): Argentina 2018-2021 (Julio)
Años
IVCR
2018
-0.0008839
2019
-0.0008268
2020
-0.0027043
2021
0.0000000
a Fuente: Elaboración propia, con base en información del BCR.
Por partida
Se analiza el IVCR con Argentina ARG, para el periodo 2018-2021, para el producto “2204100000”
options(scipen = 99999)
Resultados_IVCR<-data.frame("Años"=2018:2021,
"IVCR"=sapply(X=2018:2021,
FUN = indicadores_IVCR_partida,
codigo_pais=32,
partida="2204100000"))
Resultados_IVCR %>% kable(caption = "Indicador de Ventaja Comparativa Revelada (partida): Argentina 2018-2021 (Julio)", align = "c") %>% add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>%
kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
Indicador de Ventaja Comparativa Revelada (partida): Argentina 2018-2021 (Julio)
Años
IVCR
2018
-0.9467259
2019
-0.9966814
2020
-0.9360780
2021
-0.9944464
a Fuente: Elaboración propia, con base en información del BCR.
Indicador de Herfindahl Hirschmann anual
Indicador por default
Ambos casos para Canadá CAN 2018 - 2021.
Resultados_IHH <- data.frame("Años"=2018:2021,
"IHH"=sapply(X = 2018:2021,
FUN = indicadores_IHH_Herfindahl_Hirschmann_anual,
codigo_pais=124,
normalizado=TRUE))
Resultados_IHH %>%
kable(caption = "Indicador de Herfindahl Hirschmann anual: Canadá 2018-2021 (Julio)", align = "c", digits = 6) %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
Indicador de Herfindahl Hirschmann anual: Canadá 2018-2021 (Julio)
Años
IHH
2018
0.008444
2019
0.119374
2020
0.028965
2021
0.016459
a Fuente: Elaboración propia, con base en información del BCR.
Indicador propuesto por el texto
Resultados_IHH_texto<-data.frame("Años"=2018:2021,
"IHH"=sapply(X = 2018:2021,
FUN = indicadores_IHH_Herfindahl_Hirschmann_anual_texto,
codigo_pais=124,
normalizado=TRUE))
Resultados_IHH_texto %>%
kable(caption = "Indicador de Herfindahl Hirschmann (texto): Canadá 2018-2021 (Julio)", align = "c", digits = 6) %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
Indicador de Herfindahl Hirschmann (texto): Canadá 2018-2021 (Julio)
Años
IHH
2018
-0.001764
2019
-0.001827
2020
-0.002074
2021
-0.003663
a Fuente: Elaboración propia, con base en información del BCR.
2) Ranking de socios comerciales
Indicación: Genere una función personalizada que Obtenga un ranking, anual, de los “K” principales socios comerciales de El Salvador, para el periodo año inicio-año final. Presente sus resultados en el siguiente formato: Top K Socios comerciales, periodo año inicio – año final, datos en porcentaje de las exportaciones totales.
Nota: ISO país | %Expor. Significa que debe aparecer el código iso_3 del país y el porcentaje de exportaciones totales de el Salvador para ese país, en el año en cuestión, las columnas representan el lugar en el ranking, por ejemplo: GTM | 25.7
Creando la función
library(tidyr)
library(kableExtra)
# Seleccionar Años
anios_ranking <- 2018:2021
# Filtro
data_comercio_exterior %>%
filter(anio %in% anios_ranking) -> data_ranking
# Determinar el valor de k
# Función
insumo_reporte <- function(k){
data_ranking %>%
group_by(anio, iso_3) %>%
summarise(total = sum(valor_fob)) %>%
mutate(percent = round(prop.table(total) * 100, 2)) %>%
slice_max(n=k, order_by = total) %>%
as.data.frame() %>%
group_by(anio) %>%
mutate(rank = row_number(), data = paste(iso_3, " | ", percent, sep = "")) %>%
select(anio, data, rank) %>%
as.data.frame()
}
Mostrando los resultados
insumo_reporte(10) %>%
pivot_wider(names_from = rank,values_from = data)->tabla
tabla %>%
rename(Año = `anio`) %>%
kable(caption = "Top k Socios Comerciales 2018 - 2021", align = "c") %>%
footnote(general = "Fuente: Elaboración propia con base en datos del BCR",
number = c("K selecciona los primeros países del ranking")) %>%
kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")
Top k Socios Comerciales 2018 - 2021
Año
1
2
3
4
5
6
7
8
9
10
2018
USA | 44.07
HND | 15.34
GTM | 14.36
NIC | 6.87
CRI | 4.39
MEX | 2.44
PAN | 2.09
DOM | 1.65
CHN | 1.45
ESP | 0.82
2019
USA | 42.03
GTM | 15.92
HND | 15.9
NIC | 6.66
CRI | 4.46
MEX | 2.54
PAN | 2.03
DOM | 1.75
CAN | 1.09
CHN | 0.88
2020
USA | 39.18
GTM | 16.89
HND | 15.44
NIC | 7.26
CRI | 4.65
MEX | 2.46
PAN | 1.89
DOM | 1.78
CHN | 1.54
KOR | 1.21
2021
USA | 40.59
GTM | 16.98
HND | 15.82
NIC | 6.99
CRI | 4.04
MEX | 2.15
KOR | 1.98
PAN | 1.9
DOM | 1.6
DEU | 0.9
Note:
Fuente: Elaboración propia con base en datos del BCR
1 K selecciona los primeros países del ranking
UNIVERSIDAD DE EL SALVADOR
FACULTAD DE CIENCIAS ECONOMICAS
ESCUELA DE ECONOMIA
Materia: Métodos para el análisis económico
Docente: Carlos Ademir Pérez Alas
Grupo teórico: 3
Laboratorio 1: Segunda parte
Integrantes:
| Apellidos | Nombres | DUE |
|---|---|---|
| Escolán Chávez | Diana Alejandra | EC19003 |
| García Cabrera | Bryan Alexander | GC17011 |
Ciclo: II-2021
Fecha: Vienes 10 de septiembre, 2021
Ciudad Universitaria, San Salvador, El Salvador
Carga de la Base de Comercio Exterior
library(dplyr)
library(kableExtra)
load("C:/Archivos R/MAE118/base_comercio.RData")
data_comercio_exterior <- base_comercio1) Reproduciendo Funciones de Comercio Exterior
Indicación: Reproduzca los ejemplos vistos en clase sobre el uso de la librería dplyr, usando la información actualizada recopilada en la tarea 1, y usando los nombres estandarizados que se han agregado para los países, utilice países diferentes a los que se emplearon en los ejemplos (queda a elección del grupo los paises).
# Carga de funciones
source(file = "C:/Archivos R/MAE118/funciones_comercio_exterior.R",echo = FALSE,max.deparse.length = 100,encoding = "utf-8")Consulta de países
library(kableExtra)
herramientas_mostrar_paises() %>% tail() %>%
kable(caption = "Países en la Base de Comercio Exterior 2018-2021 (Julio)",
align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| pais | codigo_pais | region | cod_region | region_intermedia | cod_region_intermedia |
|---|---|---|---|---|---|
| Zimbabwe | 716 | África | 2 | África Oriental | 14 |
| Zimbabwe | 716 | África | 2 | África Oriental | 14 |
| Zimbabwe | 716 | África | 2 | África Oriental | 14 |
| Zimbabwe | 716 | África | 2 | África Oriental | 14 |
| Zimbabwe | 716 | África | 2 | África Oriental | 14 |
| Zimbabwe | 716 | África | 2 | África Oriental | 14 |
| a Fuente: Elaboración propia, con base en información del BCR. |
Consulta del número de transacciones por país
herramientas_mostrar_paises_en_tabla() %>%
tail() %>% kable(caption = "Países con transacciones 2018-2021 (Julio)",
align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| pais | anio | # Transacciones |
|---|---|---|
| Yugoeslavia -Nueva Organizacio | 2020 | 1 |
| Zambia | 2021 | 4 |
| Zimbabwe | 2018 | 1 |
| Zimbabwe | 2019 | 8 |
| Zimbabwe | 2020 | 3 |
| Zimbabwe | 2021 | 1 |
| a Fuente: Elaboración propia, con base en información del BCR. |
Consulta de regiones en tabla
herramientas_mostrar_regiones() %>%
tail() %>% kable(caption = "Regiones en la Base de Comercio Exterior 2018-2021 (Julio)",
align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| region | cod_region | region_intermedia | cod_region_intermedia | # de Paises incluidos |
|---|---|---|---|---|
| Américas | 19 | NA | NA | 132225 |
| Asia | 142 | NA | NA | 250262 |
| Europa | 150 | Islas Del Canal | 830 | 7 |
| Europa | 150 | NA | NA | 170958 |
| Oceanía | 9 | NA | NA | 4139 |
| NA | NA | NA | NA | 229 |
| a Fuente: Elaboración propia, con base en información del BCR. |
Indicador de Balassa
Por capítulo
Se analiza el IB con Alemania DEU, para el periodo 2018-2021, para el producto “1905200000”
Resultados_IB_sac <- data.frame("Años"=2018:2021,
"IB"=sapply(X=2018:2021,
FUN = indicadores_Balassa_partida,
codigo_pais=276,
partida="1905200000",
normalizado=TRUE))
Resultados_IB_sac %>% kable(caption = "IB capítulo: Alemania 2018-2021 (Julio)", align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| Años | IB |
|---|---|
| 2018 | 0.5560822 |
| 2019 | 0.5561221 |
| 2020 | 0.6081481 |
| 2021 | 0.6866135 |
| a Fuente: Elaboración propia, con base en información del BCR. |
Por partida
Se analizará el IB con Alemania DEU, para el periodo 2018-2021, para el capítulo “02”
Resultados_Balassa <- data.frame("Años"=2018:2021,
"IB"=sapply(X=2018:2021,
FUN = indicadores_Balassa_capitulo,
codigo_pais=276,
capitulo="02"))
Resultados_IB_sac %>% kable(caption = "IB capítulo: Alemania 2018-2021 (Julio)",
align = "c") %>%
add_footnote(label = "Fuente: Elaboración propia,con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| Años | IB |
|---|---|
| 2018 | 0.5560822 |
| 2019 | 0.5561221 |
| 2020 | 0.6081481 |
| 2021 | 0.6866135 |
| a Fuente: Elaboración propia,con base en información del BCR. |
Uso de Indicador de Ventaja Comparativa Revelada
Por capítulo
Se analiza el IVCR con Argentina ARG, para el periodo 2018-2021, para al capitulo “01”
options(scipen = 99999)
Resultados_IVCR<-data.frame("Años"=2018:2021,
"IVCR"=sapply(X=2018:2021,
FUN = indicadores_IVCR_capitulo,
codigo_pais=32,
capitulo="01"))
Resultados_IVCR %>% kable(caption = "Indicador de Ventaja Comparativa Revelada (capítulo): Argentina 2018-2021 (Julio)",
align = "c") %>% add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>%
kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| Años | IVCR |
|---|---|
| 2018 | -0.0008839 |
| 2019 | -0.0008268 |
| 2020 | -0.0027043 |
| 2021 | 0.0000000 |
| a Fuente: Elaboración propia, con base en información del BCR. |
Por partida
Se analiza el IVCR con Argentina ARG, para el periodo 2018-2021, para el producto “2204100000”
options(scipen = 99999)
Resultados_IVCR<-data.frame("Años"=2018:2021,
"IVCR"=sapply(X=2018:2021,
FUN = indicadores_IVCR_partida,
codigo_pais=32,
partida="2204100000"))
Resultados_IVCR %>% kable(caption = "Indicador de Ventaja Comparativa Revelada (partida): Argentina 2018-2021 (Julio)", align = "c") %>% add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>%
kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| Años | IVCR |
|---|---|
| 2018 | -0.9467259 |
| 2019 | -0.9966814 |
| 2020 | -0.9360780 |
| 2021 | -0.9944464 |
| a Fuente: Elaboración propia, con base en información del BCR. |
Indicador de Herfindahl Hirschmann anual
Indicador por default
Ambos casos para Canadá CAN 2018 - 2021.
Resultados_IHH <- data.frame("Años"=2018:2021,
"IHH"=sapply(X = 2018:2021,
FUN = indicadores_IHH_Herfindahl_Hirschmann_anual,
codigo_pais=124,
normalizado=TRUE))
Resultados_IHH %>%
kable(caption = "Indicador de Herfindahl Hirschmann anual: Canadá 2018-2021 (Julio)", align = "c", digits = 6) %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| Años | IHH |
|---|---|
| 2018 | 0.008444 |
| 2019 | 0.119374 |
| 2020 | 0.028965 |
| 2021 | 0.016459 |
| a Fuente: Elaboración propia, con base en información del BCR. |
Indicador propuesto por el texto
Resultados_IHH_texto<-data.frame("Años"=2018:2021,
"IHH"=sapply(X = 2018:2021,
FUN = indicadores_IHH_Herfindahl_Hirschmann_anual_texto,
codigo_pais=124,
normalizado=TRUE))
Resultados_IHH_texto %>%
kable(caption = "Indicador de Herfindahl Hirschmann (texto): Canadá 2018-2021 (Julio)", align = "c", digits = 6) %>%
add_footnote(label = "Fuente: Elaboración propia, con base en información del BCR.") %>% kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| Años | IHH |
|---|---|
| 2018 | -0.001764 |
| 2019 | -0.001827 |
| 2020 | -0.002074 |
| 2021 | -0.003663 |
| a Fuente: Elaboración propia, con base en información del BCR. |
2) Ranking de socios comerciales
Indicación: Genere una función personalizada que Obtenga un ranking, anual, de los “K” principales socios comerciales de El Salvador, para el periodo año inicio-año final. Presente sus resultados en el siguiente formato: Top K Socios comerciales, periodo año inicio – año final, datos en porcentaje de las exportaciones totales.
Nota: ISO país | %Expor. Significa que debe aparecer el código iso_3 del país y el porcentaje de exportaciones totales de el Salvador para ese país, en el año en cuestión, las columnas representan el lugar en el ranking, por ejemplo: GTM | 25.7
Creando la función
library(tidyr)
library(kableExtra)
# Seleccionar Años
anios_ranking <- 2018:2021
# Filtro
data_comercio_exterior %>%
filter(anio %in% anios_ranking) -> data_ranking
# Determinar el valor de k
# Función
insumo_reporte <- function(k){
data_ranking %>%
group_by(anio, iso_3) %>%
summarise(total = sum(valor_fob)) %>%
mutate(percent = round(prop.table(total) * 100, 2)) %>%
slice_max(n=k, order_by = total) %>%
as.data.frame() %>%
group_by(anio) %>%
mutate(rank = row_number(), data = paste(iso_3, " | ", percent, sep = "")) %>%
select(anio, data, rank) %>%
as.data.frame()
}Mostrando los resultados
insumo_reporte(10) %>%
pivot_wider(names_from = rank,values_from = data)->tabla
tabla %>%
rename(Año = `anio`) %>%
kable(caption = "Top k Socios Comerciales 2018 - 2021", align = "c") %>%
footnote(general = "Fuente: Elaboración propia con base en datos del BCR",
number = c("K selecciona los primeros países del ranking")) %>%
kable_styling() %>%
row_spec(0, bold = T, color = "white", background = "#D7261E")| Año | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 |
|---|---|---|---|---|---|---|---|---|---|---|
| 2018 | USA | 44.07 | HND | 15.34 | GTM | 14.36 | NIC | 6.87 | CRI | 4.39 | MEX | 2.44 | PAN | 2.09 | DOM | 1.65 | CHN | 1.45 | ESP | 0.82 |
| 2019 | USA | 42.03 | GTM | 15.92 | HND | 15.9 | NIC | 6.66 | CRI | 4.46 | MEX | 2.54 | PAN | 2.03 | DOM | 1.75 | CAN | 1.09 | CHN | 0.88 |
| 2020 | USA | 39.18 | GTM | 16.89 | HND | 15.44 | NIC | 7.26 | CRI | 4.65 | MEX | 2.46 | PAN | 1.89 | DOM | 1.78 | CHN | 1.54 | KOR | 1.21 |
| 2021 | USA | 40.59 | GTM | 16.98 | HND | 15.82 | NIC | 6.99 | CRI | 4.04 | MEX | 2.15 | KOR | 1.98 | PAN | 1.9 | DOM | 1.6 | DEU | 0.9 |
| Note: | ||||||||||
| Fuente: Elaboración propia con base en datos del BCR | ||||||||||
| 1 K selecciona los primeros países del ranking |