Laboratorio1
Segunda_Parte
7/9/2021
UNIVERSIDAD DE EL SALVADOR
FACULTAD DE CIENCIAS ECONÓMICAS
ESCUELA DE ECONOMÍA
CICLO II-2021
“Laboratorio 1 Segunda parte”
Métodos para el Análisis Económico.
Docente:
Carlos Ademir Pérez Alas.
Integrantes:
Argueta Velasco Rocío Azucena. (AV14026)
Henríquez Urias Gerson Stanley. (HU17001)
Sanchez Reyes Paola Rosmery. (SR16050)
Velásquez Gonzáles Jóse Roberto. (VG16027)
Grupo Teorico:
03
Ciudad Universitaria, Miércoles 8 de septiembre de 2021.
Cargar las funciones
load("C:/Users/Admin/Desktop/MAE/Laboratorio 1/Parte 1/ base_comercio_std.RData")
nombre_archivo <- "C:/Users/Admin/Desktop/MAE/Laboratorio 1/nombres_iso_paises.xlsx"
source(file = "C:/Users/Admin/Desktop/MAE/Laboratorio 1/Parte 1/segunda_parte/funciones_comercio_exterior.R",echo=TRUE,max.deparse.length = 100,encoding = "utf-8")##
## > options(scipen = 999)
##
## > library(dplyr)
##
## > library(readxl)
##
## > library(stringr)
##
## > load("C:/Users/Admin/Desktop/MAE/Laboratorio 1/Parte 1/ tabla_comercio_exterior_update.RData")
##
## > nombre_archivo <- "C:/Users/Admin/Desktop/MAE/Laboratorio 1/nombres_iso_paises.xlsx"
##
## > nombres_iso_paises <- read_excel(nombre_archivo)
##
## > data_comercio_exterior_update <- data_comercio_exterior_update %>%
## + left_join(nombres_iso_pa .... [TRUNCATED]
##
## > herramientas_mostrar_paises <- function() {
## + paises <- nombres_iso_paises %>% select(nom_pais .... [TRUNCATED]
##
## > herramientas_mostrar_regiones <- function() {
## + nombres_iso_paises %>% group_by(region, cod_re .... [TRUNCATED]
##
## > herramientas_mostrar_paises_en_tabla <- function() {
## + data_comercio_exterior_update %>% group .... [TRUNCATED]
##
## > indicadores_IVCR_capitulo <- function(codigo_pais,
## + capitulo, anio) {
## + codigo_pais <- e .... [TRUNCATED]
##
## > indicadores_IVCR_partida <- function(codigo_pais,
## + partida, anio) {
## + codigo_pais <- enq .... [TRUNCATED]
##
## > indicadores_Balassa_capitulo <- function(codigo_pais,
## + capitulo, anio, normalizado = TRUE) { .... [TRUNCATED]
##
## > indicadores_Balassa_partida <- function(codigo_pais,
## + partida, anio, normalizado = TRUE) {
## + .... [TRUNCATED]
##
## > indicadores_IHH_Herfindahl_Hirschmann_anual_texto <- function(codigo_pais,
## + anio, normalizad .... [TRUNCATED]
##
## > indicadores_IHH_Herfindahl_Hirschmann_anual <- function(codigo_pais,
## + anio, normalizado = TR .... [TRUNCATED]Numeral 1
Herramienta de consulta de paises
# Cargamos libreria
library(kableExtra)
# Creamos herramienta
herramientas_mostrar_paises() %>% head() %>% kable() %>% kable_material(html_font = "sans-serif")| nom_pais_esp | codigo_pais | region | cod_region | region_intermedia | cod_region_intermedia |
|---|---|---|---|---|---|
| Afganistan | 4 | Asia | 142 | NA | NA |
| Albania | 8 | Europa | 150 | NA | NA |
| Alemania | 276 | Europa | 150 | NA | NA |
| Andorra | 20 | Europa | 150 | NA | NA |
| Angola | 24 | África | 2 | África Central | 17 |
| Anguila (R.U.) | 660 | Américas | 19 | Caribe | 29 |
Herramienta de consulta de regiones
# Cargamos libreria
library(kableExtra)
# Creamos herramienta
herramientas_mostrar_regiones() %>% head() %>% kable() %>% kable_material(html_font = "sans-serif")| region | cod_region | region_intermedia | cod_region_intermedia | # de Paises incluidos |
|---|---|---|---|---|
| África | 2 | África Central | 17 | 8 |
| África | 2 | Africa Del Sur | 18 | 5 |
| África | 2 | África Occidental | 11 | 16 |
| África | 2 | África Oriental | 14 | 17 |
| África | 2 | NA | NA | 6 |
| Américas | 19 | Caribe | 29 | 26 |
Herramientas de Consulta de Paises con transacciones en la tabla
# Cargamos libreria
library(kableExtra)
# Creamos herramienta
herramientas_mostrar_paises_en_tabla() %>% head() %>% kable() %>% kable_material(html_font = "sans-serif")| pais | anio | # Transacciones |
|---|---|---|
| Afganistan | 2018 | 7 |
| Afganistan | 2019 | 10 |
| Afganistan | 2020 | 29 |
| Afganistan | 2021 | 6 |
| Albania | 2018 | 27 |
| Albania | 2019 | 46 |
Reproduzca los ejemplos visto en clases sobre el uso de la librería dplyr.
Uso de Indicador de Ventaja Comparativa Revelada
Ejemplo 1.1
En este ejemplo se analizará el IVCR con Costa Rica, para el periodo 2018-2021, para al capitulo “02”
Resultados_IVCR<-data.frame("años"=2018:2021,
"IVCR"=sapply(X=2018:2021,FUN = indicadores_IVCR_capitulo,codigo_pais=188,capitulo="02"))
Resultados_IVCR %>% kable(align = "l") %>%
kable_material(html_font = "sans-serif")| años | IVCR |
|---|---|
| 2018 | -0.0041732 |
| 2019 | -0.0009633 |
| 2020 | -0.0036519 |
| 2021 | -0.0096628 |
Ejemplo 1.2
En este ejemplo se analizará el IVCR con Costa Rica, para el periodo 2018-2021, para el producto “0105110000”
Resultados_IVCR_sac<-data.frame("años"=2018:2021,
"IVCR"=sapply(X=2018:2021,FUN = indicadores_IVCR_partida,codigo_pais =188,partida ="0105110000"))
Resultados_IVCR_sac %>% kable(align = "l") %>%
kable_material(html_font = "sans-serif")| años | IVCR |
|---|---|
| 2018 | -0.0093796 |
| 2019 | -0.0755637 |
| 2020 | -0.1306645 |
| 2021 | -0.1699538 |
Uso de Indicador de Balassa, se usará a Costa Rica como pais
Ejemplo 1.3
En este ejemplo se analizará el indice de Balassa con Costa Rica, para el periodo 2018-2021, para al capitulo “02”
Resultados_Balassa<-data.frame("años"=2018:2021,
"IB"=sapply(X=2018:2021,FUN = indicadores_Balassa_capitulo,codigo_pais=188,capitulo="02"))
Resultados_Balassa %>% kable(align = "l") %>%
kable_material(html_font = "sans-serif")| años | IB |
|---|---|
| 2018 | -0.6496485 |
| 2019 | -1.0000000 |
| 2020 | -0.3696722 |
| 2021 | -1.0000000 |
Ejemplo 1.4
En este ejemplo se analizará el IB con Costa Rica, para el periodo 2017-2020, para el producto “0105110000”
Resultados_IB_sac<-data.frame("años"=2018:2021,
"IB"=sapply(X=2018:2021,
FUN = indicadores_Balassa_partida,
codigo_pais=188,
partida="0105110000",
normalizado=TRUE))
Resultados_IB_sac %>% kable(align = "l") %>%
kable_material(html_font = "sans-serif")| años | IB |
|---|---|
| 2018 | 0.5560822 |
| 2019 | 0.5561221 |
| 2020 | 0.6081481 |
| 2021 | 0.6866135 |
Uso de Indicador de Herfindahl Hirschmann anual
Ambos ejemplos para Costa Rica para el período 2018-2021
Ejemplo 1.5 Uso de Indicador de Herfindahl Hirschmann anual
Resultados_IHH<-data.frame("años"=2018:2021,"IHH"=sapply(X = 2018:2021,
FUN = indicadores_IHH_Herfindahl_Hirschmann_anual,codigo_pais=188,normalizado=TRUE))
Resultados_IHH %>% kable(align = "l",digits = 6) %>%
kable_material(html_font = "sans-serif")| años | IHH |
|---|---|
| 2018 | 0.002505 |
| 2019 | 0.002879 |
| 2020 | 0.004631 |
| 2021 | 0.007052 |
Ejemplo 1.6 Uso de Indicador de Herfindahl Hirschmann anual propuesto en 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=188,normalizado=TRUE))
Resultados_IHH_texto %>% kable(align = "l",digits = 6) %>%
kable_material(html_font = "sans-serif")| años | IHH |
|---|---|
| 2018 | -0.000171 |
| 2019 | -0.000168 |
| 2020 | -0.000178 |
| 2021 | -0.000241 |
Numeral 2
Generé una función personalizada que obtenga un ranking anual, de los “k” principales socios comerciales de El Salvador, para el período año inicio - año final
Función personalizada
library(tidyr)
library(kableExtra)
Ranking<-function(anio,k){
anios_ranking<-2018:2021
data_comercio_exterior_update %>%
filter(anio %in% anios_ranking) ->data_ranking
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() -> insumo_reporte
insumo_reporte %>%
pivot_wider(names_from = rank,values_from = data)%>%
filter(anio==!!anio) %>%
kable(caption = "Top k socios comerciales de El Salvador para el periodo año inicio -año final, datos en porcentaje de las exportaciones totales",align = "c")%>%
add_footnote(label = "Elaboración propia con base en datos del BCR",
notation = "symbol") %>% kable_styling()}Acontinuación se presentán dos ejemplos para verificar nuestra función personalizada.
Ejemplo 2.1
Se presentará el Top 10 socios comerciales de El Salvador, para el periodo 2018 - 2020
library(tidyr)
library(kableExtra)
Ranking(anio = 2018:2020,k=10)| anio | 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 |
| * Elaboración propia con base en datos del BCR |
Ejemplo 2.2
Se presentará el Top 7 socios comerciales de El Salvador, para el periodo 2018 - 2021
library(tidyr)
library(kableExtra)
Ranking(anio = 2018:2021,k=7)| anio | 1 | 2 | 3 | 4 | 5 | 6 | 7 |
|---|---|---|---|---|---|---|---|
| 2018 | USA|44.07 | HND|15.34 | GTM|14.36 | NIC|6.87 | CRI|4.39 | MEX|2.44 | PAN|2.09 |
| 2019 | USA|42.03 | GTM|15.92 | HND|15.9 | NIC|6.66 | CRI|4.46 | MEX|2.54 | PAN|2.03 |
| 2020 | USA|39.18 | GTM|16.89 | HND|15.44 | NIC|7.26 | CRI|4.65 | MEX|2.46 | PAN|1.89 |
| 2021 | USA|40.59 | GTM|16.98 | HND|15.82 | NIC|6.99 | CRI|4.04 | MEX|2.15 | KOR|1.98 |
| * Elaboración propia con base en datos del BCR |