UNIVERSIDAD DE EL SALVADOR

FACULTAD DE CIENCIAS ECONÓMICAS

ESCUELA DE ECONOMÍA

MÉTODOS PARA EL ANÁLISIS ECONÓMICO


LABORATORIO 1. PRIMERA PARTE


DOCENTE:

MSF. CARLOS ADEMIR PÉREZ ALAS.

INTEGRANTES:

ARITA POLANCO, BRIAN GUADALUPE.      AP18007.

MEJÍA HERRERA, NORMAN JAVIER.       MH18013.

RAMÍREZ URQUIZA, ALEJANDRA EUNICE.      RU18002.

GRUPO:

3

CIUDAD UNIVERSITARIA, DOMINGO 5 DE SEPTIEMBRE DE 2021.

1. Compile la información disponible en la “Base de Datos de Comercio Exterior” del BCR, para los años 2018-2020.

1.1 Tabla Importaciones

library(dplyr)
library(readr)
library(kableExtra)

setwd("C:/Users/EunicePC/Downloads/Importaciones/")

nombres_archivos_importar<-list.files()

lista_importacion<-lapply(nombres_archivos_importar,
                          FUN = read_delim,"|",
                          col_names=c("pais","sac","anio","mes","valor_cif","kilogramos"),
col_types=cols(pais="c",sac="c", anio= "d", mes="d", valor_cif="d", kilogramos="d"), skip=1)

#Ensamblar los archivos en un solo dataframe
importaciones_mensuales<-bind_rows(lista_importacion)

#Mostrar la tabla con solo 6 elementos
importaciones_mensuales %>% head()%>% kable(caption="El Salvador: Importaciones Mensuales desde 2018 a 2020", align= "c") %>% add_footnote(label= "Elaboración propia con base en datos del BCR de El Salvador", notation= "symbol") %>% kable_minimal()
El Salvador: Importaciones Mensuales desde 2018 a 2020
pais sac anio mes valor_cif kilogramos
Canada 0404900000 2018 2 68040.00 18810.00
Canada 0404900000 2018 3 92725.74 25528.60
Canada 0404900000 2018 4 92546.45 25525.00
Canada 0404900000 2018 6 91531.00 25525.00
Canada 0405100000 2018 6 26.56 20.00
Canada 0511100000 2018 1 16365.04 23.65
* Elaboración propia con base en datos del BCR de El Salvador

1.2 Tabla Exportaciones

library(dplyr)
library(readr)
library(kableExtra)

setwd("C:/Users/EunicePC/Downloads/Exportaciones/")

nombres_archivos_importar2 <-list.files()


lista_exportacion2<-lapply(nombres_archivos_importar2,
                          FUN = read_delim,"|",
                          col_names=c("pais","sac","anio","mes","valor_fob","kilogramos"),
col_types=cols(pais="c",sac="c", anio= "d", mes="d", valor_fob="d", kilogramos="d"), skip=1)

#Ensamblar los archivos en un solo dataframe
exportaciones_mensuales<-bind_rows(lista_exportacion2)

#Mostrar la tabla con solo 6 elementos
exportaciones_mensuales %>% head()%>% kable(caption="El Salvador: Exportaciones Mensuales desde 2018 a 2020", align= "c") %>% add_footnote(label= "Elaboración propia con base en datos del BCR de El Salvador", notation= "symbol") %>% kable_minimal()
El Salvador: Exportaciones Mensuales desde 2018 a 2020
pais sac anio mes valor_fob kilogramos
Canada 0106900000 2018 2 4667.2 15.0
Canada 0106900000 2018 3 4237.1 14.0
Canada 0106900000 2018 4 3918.0 13.5
Canada 0106900000 2018 5 3155.4 10.0
Canada 0106900000 2018 6 3097.4 9.0
Canada 0305720000 2018 2 2702.3 280.0
* Elaboración propia con base en datos del BCR de El Salvador

1.3 Tabla Comercio Exterior

library(dplyr)
library(tidyr)
library(kableExtra)

#Unión de las tablas de Importaciones y Exportaciones
importaciones_mensuales %>% full_join(
  exportaciones_mensuales,
  by= c("pais", "sac", "anio", "mes"),
  suffix= c("_importaciones", "_exportaciones")
) -> data_comercio_exterior

#Sustituir los valores NULLS (NA) por 0
data_comercio_exterior %>% replace_na(
  list(
    valor_cif= 0,
    valor_fob= 0,
    kilogramos_importaciones= 0,
    kilogramos_exportaciones= 0
  )
) %>%
  arrange(pais, anio, mes, sac) -> data_comercio_exterior

#Mostrar la tabla con solo 6 elementos
data_comercio_exterior %>% head()%>% kable(caption="El Salvador: Comercio Exterior desde 2018 a 2020", align= "c") %>% add_footnote(label= "Elaboración propia con base en datos del BCR de El Salvador", notation= "symbol") %>% kable_minimal()
El Salvador: Comercio Exterior desde 2018 a 2020
pais sac anio mes valor_cif kilogramos_importaciones valor_fob kilogramos_exportaciones
Afganistan 0806200000 2018 6 6448.43 1463.92 0 0
Afganistan 6104220000 2018 10 3153.37 2407.61 0 0
Afganistan 6104620000 2018 10 946.01 722.28 0 0
Afganistan 6105100000 2018 10 9405.39 7181.03 0 0
Afganistan 6106100000 2018 10 1353.32 1725.55 0 0
Afganistan 6405900000 2018 10 2260.03 1725.55 0 0
* Elaboración propia con base en datos del BCR de El Salvador

1.4 Guardar Tabla Comercio Exterior

#Para guardar la Tabla de Comercio Exterior
save(data_comercio_exterior, file = "C:/Users/EunicePC/Downloads/data_comercio_exterior.RData")

2. Actualización de la tabla anterior con información disponible para 2021.

library(dplyr)
library(tidyr)
library(kableExtra)

#Agregar Importaciones hasta Julio 2021
setwd("C:/Users/EunicePC/Desktop/Importaciones 1/")

archivos_importar_ac <- list.files()

lista_importacion_ac <- lapply(
  archivos_importar_ac,
  FUN = read_delim,
  delim = "|",
  col_names = c("pais",
                "sac",
                "anio",
                "mes",
                "valor_cif",
                "kilogramos"),
  col_types = cols(
    pais = "c",
    sac = "c",
    anio = "d",
    mes = "d",
    valor_cif = "d",
    kilogramos = "d"
  ),
  skip = 1
)

#Ensamblar los archivos en un sólo dataframe
data_importaciones_update <- bind_rows(lista_importacion_ac)

#Agregar Exportaciones hasta Julio 2021
setwd("C:/Users/EunicePC/Desktop/Exportaciones 1/")

archivos_importar_acex <- list.files()

lista_importacion_acex <- lapply(
  archivos_importar_acex,
  FUN = read_delim,
  delim = "|",
  col_names = c("pais",
                "sac",
                "anio",
                "mes",
                "valor_fob",
                "kilogramos"),
  col_types = cols(
    pais = "c",
    sac = "c",
    anio = "d",
    mes = "d",
    valor_fob = "d",
    kilogramos = "d"
  ),
  skip = 1
)

#Ensamblar los archivos en un sólo dataframe
data_exportaciones_update <- bind_rows(lista_importacion_acex)

#Concatenar las tablas de data_importaciones y data_exportaciones
data_importaciones_update %>% full_join(
  data_exportaciones_update,
  by = c("pais", "sac", "anio", "mes"),
  suffix = c("_importaciones", "_exportaciones")
) -> data_comercio_exterior_update

#Sustituir los Valores Nulos "NA" y guardar en el mismo objeto
data_comercio_exterior_update %>% replace_na(
  list(
    valor_cif = 0,
    valor_fob = 0,
    kilogramos_importaciones = 0,
    kilogramos_exportaciones = 0
  )
) %>%
  arrange(pais, anio, mes, sac) -> data_comercio_exterior_update

#Unir la tabla de comercio exterior y su actualización
data_comercio_exterior %>% bind_rows(data_comercio_exterior_update) %>% arrange(pais, anio, mes, sac) ->
  data_comercio_exterior_update_unida

#Volviendo a guardar
save(data_comercio_exterior_update_unida, file = "C:/Users/EunicePC/Downloads/data_comercio_exterior_actualizada_unida.RData")

#Mostrar la tabla con solo 6 elementos
data_comercio_exterior_update_unida %>% filter(anio == 2021, mes == 07) %>% head() %>%
  kable(caption = "El Salvador: Comercio Exterior desde enero 2018 hasta julio 2021",
        align = "c") %>%
  add_footnote(label = "Elaboración propia con base en datos del BCR de El Salvador",
               notation = "symbol") %>%  kable_minimal()
El Salvador: Comercio Exterior desde enero 2018 hasta julio 2021
pais sac anio mes valor_cif kilogramos_importaciones valor_fob kilogramos_exportaciones
Albania 9801002000 2021 7 5115.02 632.00 0 0.0
Alemania 0406400000 2021 7 196.06 7.81 0 0.0
Alemania 0406909000 2021 7 4720.73 680.42 0 0.0
Alemania 0409000000 2021 7 0.00 0.00 85050 21994.0
Alemania 0602100000 2021 7 1024.40 5.00 0 0.0
Alemania 0901113000 2021 7 0.00 0.00 707275 169606.7
* Elaboración propia con base en datos del BCR de El Salvador

3. Agregue la información estandarizada de los países.

library(dplyr)
library(readxl)
library(kableExtra)
options(scipen = 9999999)

load("C:/Users/EunicePC/Downloads/data_comercio_exterior_actualizada_unida.RData")

nombre_archivo <- "C:/Users/EunicePC/Downloads/nombres_iso_paises.xlsx"

nombres_iso_paises <- read_excel(nombre_archivo)

data_comercio_exterior_update_unida %>% 
  left_join(nombres_iso_paises, 
            by= c("pais"="nom_pais_esp")) %>% 
  group_by(pais, sub_region, anio, iso_3) %>% 
  summarise(`Exportaciones Totales MM US$` = sum(valor_fob),
            `Importaciones Totales MM US$` = sum(valor_cif)) %>% 
  head()%>% kable(digits = 2) %>% kable_minimal() 
pais sub_region anio iso_3 Exportaciones Totales MM US$ Importaciones Totales MM US$
Afganistan Asia Meridional 2018 AFG 0.00 47146.22
Afganistan Asia Meridional 2019 AFG 764.26 26227.54
Afganistan Asia Meridional 2020 AFG 0.00 10853.48
Afganistan Asia Meridional 2021 AFG 37873.33 0.00
Albania Europa Del Sur 2018 ALB 0.00 36614.72
Albania Europa Del Sur 2019 ALB 102785.00 51320.32

4. Obtenga un ranking, anual de los 5 principales socios comerciales de El Salvador, para el periodo 2018-2021.

4.1 Exportaciones Totales 2018

library(dplyr)
library(readxl)
library(kableExtra)
options(scipen = 9999999)

load("C:/Users/EunicePC/Downloads/data_comercio_exterior_actualizada_unida.RData")

nombre_archivo_ran <- "C:/Users/EunicePC/Downloads/nombres_iso_paises.xlsx"

nombres_iso_paises_ran <- read_excel(nombre_archivo)

#Tabla exportaciones Totales 2018
data_comercio_exterior_update_unida %>% 
  left_join(nombres_iso_paises_ran, 
            by= c("pais"="nom_pais_esp")) %>% 
  group_by(pais,iso_3, anio) %>% 
  filter(anio==2018)%>%
  summarise(`Exportaciones` = sum(valor_fob)/11810822422*100)%>%
  arrange(desc(`Exportaciones`))%>%
  select("iso_3","Exportaciones")%>%
  head(n=5)-> tabla_exp_totales18

  tabla_exp_totales18 %>% kable(digits = 2) %>% kable_minimal()
pais iso_3 Exportaciones
Estados Unidos (U.S.A.) USA 44.07
Honduras HND 15.34
Guatemala GTM 14.36
Nicaragua NIC 6.87
Costa Rica CRI 4.39

4.1.1 Exportaciones Totales 2018 en porcentaje

expo2018<-data.frame("USA|44.1",
                     "GTM|15.3",
                     "HND|14.4",
                     "NIC|6.9",
                     "CRI|4.4")
anio<-2018
colum<-c("1","2","3","4","5")
row.names(expo2018)<-anio
colnames(expo2018)<-colum
print(expo2018)
        1        2        3       4       5

2018 USA|44.1 GTM|15.3 HND|14.4 NIC|6.9 CRI|4.4

4.2 Exportaciones Totales 2019

library(dplyr)
library(readxl)
library(kableExtra)
options(scipen = 9999999)

load("C:/Users/EunicePC/Downloads/data_comercio_exterior_actualizada_unida.RData")

nombre_archivo_ran <- "C:/Users/EunicePC/Downloads/nombres_iso_paises.xlsx"

nombres_iso_paises_ran <- read_excel(nombre_archivo)

#Tabla exportaciones Totales 2019
data_comercio_exterior_update_unida %>% 
  left_join(nombres_iso_paises_ran, 
            by= c("pais"="nom_pais_esp")) %>% 
  group_by(pais,iso_3, anio) %>% 
  filter(anio==2019)%>%
  summarise(`Exportaciones19` = sum(valor_fob)/11809459629*100) %>%
  arrange(desc(`Exportaciones19`))%>%
  select("iso_3","Exportaciones19")%>%
  head(n=5)-> tabla_exp_totales19

  tabla_exp_totales19 %>% kable(digits = 2) %>% kable_minimal()
pais iso_3 Exportaciones19
Estados Unidos (U.S.A.) USA 42.03
Guatemala GTM 15.92
Honduras HND 15.90
Nicaragua NIC 6.66
Costa Rica CRI 4.46

4.2.1 Exportaciones Totales 2019 en porcentaje

expo2019<-data.frame("USA|42.0",
                     "GTM|15.9",
                     "HND|15.9",
                     "NIC|6.7",
                     "CRI|4.5")
anio<-2019
colum<-c("1","2","3","4","5")
row.names(expo2019)<-anio
colnames(expo2019)<-colum
print(expo2019)
        1        2        3       4       5

2019 USA|42.0 GTM|15.9 HND|15.9 NIC|6.7 CRI|4.5

4.3 Exportaciones Totales 2020

library(dplyr)
library(readxl)
library(kableExtra)
options(scipen = 9999999)

load("C:/Users/EunicePC/Downloads/data_comercio_exterior_actualizada_unida.RData")

nombre_archivo_ran <- "C:/Users/EunicePC/Downloads/nombres_iso_paises.xlsx"

nombres_iso_paises_ran <- read_excel(nombre_archivo)

#Tabla exportaciones Totales 2020
data_comercio_exterior_update_unida %>% 
  left_join(nombres_iso_paises_ran, 
            by= c("pais"="nom_pais_esp")) %>% 
  group_by(pais,iso_3, anio) %>% 
  filter(anio==2020)%>%
  summarise(`Exportaciones20` = sum(valor_fob)/10088021186  *100)%>%
  arrange(desc(`Exportaciones20`))%>%
  select("iso_3","Exportaciones20")%>%
  head(n=5)-> tabla_exp_totales20

  tabla_exp_totales20 %>% kable(digits = 2) %>% kable_minimal()
pais iso_3 Exportaciones20
Estados Unidos (U.S.A.) USA 39.18
Guatemala GTM 16.89
Honduras HND 15.44
Nicaragua NIC 7.26
Costa Rica CRI 4.65

4.3.1 Exportaciones Totales 2020 en porcentaje

expo2020<-data.frame("USA|39.2",
                     "GTM|16.9",
                     "HND|15.4",
                     "NIC|7.3",
                     "CRI|4.7")
anio<-2020
colum<-c("1","2","3","4","5")
row.names(expo2020)<-anio
colnames(expo2020)<-colum
print(expo2020)
        1        2        3       4       5

2020 USA|39.2 GTM|16.9 HND|15.4 NIC|7.3 CRI|4.7

4.4 Exportaciones Totales 2021

library(dplyr)
library(readxl)
library(kableExtra)
options(scipen = 9999999)

load("C:/Users/EunicePC/Downloads/data_comercio_exterior_actualizada_unida.RData")

nombre_archivo_ran <- "C:/Users/EunicePC/Downloads/nombres_iso_paises.xlsx"

nombres_iso_paises_ran <- read_excel(nombre_archivo)

#Tabla exportaciones Totales 2021
data_comercio_exterior_update_unida %>% 
  left_join(nombres_iso_paises_ran, 
            by= c("pais"="nom_pais_esp")) %>% 
  group_by(pais,iso_3, anio) %>% 
  filter(anio==2021)%>%
  summarise(`Exportaciones21` = sum(valor_fob)/3846314150   *100)%>%
  arrange(desc(`Exportaciones21`))%>%
  head(n=5)%>%
  select("iso_3","Exportaciones21")-> tabla_exp_totales21

  tabla_exp_totales21 %>% kable(digits = 2) %>% kable_minimal()
pais iso_3 Exportaciones21
Estados Unidos (U.S.A.) USA 40.59
Guatemala GTM 16.98
Honduras HND 15.82
Nicaragua NIC 6.99
Costa Rica CRI 4.04

4.4.1 Exportaciones Totales 2021 en porcentaje

expo2021<-data.frame("USA|40.6",
                     "GTM|16.9",
                     "HND|15.8",
                     "NIC|6.9",
                     "CRI|4.0")
anio<-2021
colum<-c("1","2","3","4","5")
row.names(expo2021)<-anio
colnames(expo2021)<-colum
print(expo2021)
        1        2        3       4       5

2021 USA|40.6 GTM|16.9 HND|15.8 NIC|6.9 CRI|4.0

5. Presente sus resultados

library(kableExtra)

#Ensamblar los archivos en un solo dataframe
Top5_Socios<-bind_rows(expo2018,expo2019,expo2020,expo2021)

Top5_Socios %>% head()%>% kable(caption= "Top 5 Socios comerciales, periodo 2018-2021, datos en porcentaje de las exportaciones totales.", align= "c") %>% add_footnote(label= "Fuente: Elaboración propia con base en información del BCR de El Salvador", notation= "symbol") %>% kable_minimal()
Top 5 Socios comerciales, periodo 2018-2021, datos en porcentaje de las exportaciones totales.
1 2 3 4 5
2018 USA|44.1 GTM|15.3 HND|14.4 NIC|6.9 CRI|4.4
2019 USA|42.0 GTM|15.9 HND|15.9 NIC|6.7 CRI|4.5
2020 USA|39.2 GTM|16.9 HND|15.4 NIC|7.3 CRI|4.7
2021 USA|40.6 GTM|16.9 HND|15.8 NIC|6.9 CRI|4.0
* Fuente: Elaboración propia con base en información del BCR de El Salvador