## Warning: package 'readxl' was built under R version 4.4.3
library(dplyr)
## 
## Adjuntando el paquete: 'dplyr'
## The following objects are masked from 'package:stats':
## 
##     filter, lag
## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union
fechas_adhesion <- data.frame(
  geo = c(
    "Austria", "Belgium", "Bulgaria", "Croatia", "Cyprus", "Czechia", "Denmark",
    "Estonia", "Finland", "France", "Germany", "Greece", "Hungary", "Ireland",
    "Italy", "Latvia", "Lithuania", "Luxembourg", "Malta", "Netherlands",
    "Poland", "Portugal", "Romania", "Slovakia", "Slovenia", "Spain", "Sweden"
  ),
  año_entrada_ue = c(
    1995, 1958, 2007, 2013, 2004, 2004, 1973,
    2004, 1995, 1958, 1958, 1981, 2004, 1973,
    1958, 2004, 2004, 1958, 2004, 1958,
    2004, 1986, 2007, 2004, 2004, 1986, 1995
  )
)


IE3_con_ue <- IE3 %>%
  left_join(fechas_adhesion, by = "geo") %>%
  mutate(
    es_miembro_ue = ifelse(!is.na(año_entrada_ue) & TIME_PERIOD >= año_entrada_ue, 1, 0),
    es_ue_desconocido = is.na(año_entrada_ue)
  )

IE_export <- IE3_con_ue %>%
  filter(na_item == "Exports of goods and services", !es_ue_desconocido)

library(ggplot2)

ggplot(IE_export, aes(x = TIME_PERIOD, y = OBS_VALUE, color = as.factor(es_miembro_ue), group = geo)) +
  geom_line() +
  facet_wrap(~geo, scales = "free_y") +
  labs(title = "Evolucion de exportaciones segun pertenencia a la UE",
       x = "Anyo", y = "Exportaciones",
       color = "Miembro de la UE?") +
  scale_x_continuous(breaks = seq(1980, 2020, by = 10)) +
  theme_minimal(base_size = 12) +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

comparacion_export <- IE_export %>%
  group_by(geo, es_miembro_ue) %>%
  summarise(media_export = mean(OBS_VALUE, na.rm = TRUE), .groups = "drop")

ggplot(comparacion_export, aes(x = as.factor(es_miembro_ue), y = media_export, fill = as.factor(es_miembro_ue))) +
  geom_bar(stat = "identity") +
  facet_wrap(~geo, scales = "free_y") +
  labs(title = "Exportaciones antes y despues de entrar en la UE",
       x = "Miembro de la UE?", y = "Media de exportaciones") +
  theme_minimal(base_size = 12) +
  theme(axis.text.x = element_text(angle = 0, hjust = 0.5))

En el primer grafico observamos que en cuanto a las exportaciones, la mayoria de los países estaban al alza antes de ser miembros, tambien observamos que después de unirse sus exportaciones fueron incrementando de forma general. Lo podemos observar de forma mas resumida en la segunda imagen donde se colorea en rojo el promedio de exportaciones antes de unirse y en azul el promedio despues de unirse, donde se observa que todos los países tuvieron un incremento de exportaciones. Cabe destacar que no todos los países tienen un crecimiento similar por lo que nos podríamos plantear si podemos agrupar países segun su comportamiento.

IE_import <- IE3_con_ue %>%
  filter(na_item == "Imports of goods and services", !es_ue_desconocido)

# Gráfico 1: Evolución temporal de importaciones
ggplot(IE_import, aes(x = TIME_PERIOD, y = OBS_VALUE, color = as.factor(es_miembro_ue), group = geo)) +
  geom_line() +
  facet_wrap(~geo, scales = "free_y") +
  labs(title = "Evolucion de importaciones segun pertenencia a la UE",
       x = "Anyo", y = "Importaciones",
       color = "Miembro de la UE?") +
  scale_x_continuous(breaks = seq(1980, 2020, by = 10)) +
  theme_minimal(base_size = 12) +
  theme(axis.text.x = element_text(angle = 45, hjust = 1))

# Cálculo de medias antes/después
comparacion_import <- IE_import %>%
  group_by(geo, es_miembro_ue) %>%
  summarise(media_import = mean(OBS_VALUE, na.rm = TRUE), .groups = "drop")

# Gráfico 2: Comparación antes vs después
ggplot(comparacion_import, aes(x = as.factor(es_miembro_ue), y = media_import, fill = as.factor(es_miembro_ue))) +
  geom_bar(stat = "identity") +
  facet_wrap(~geo, scales = "free_y") +
  labs(title = "Importaciones antes y despues de entrar en la UE",
       x = "Miembro de la UE?", y = "Media de importaciones") +
  scale_x_discrete(labels = c("0" = "Antes de entrar", "1" = "Despues de entrar")) +
  scale_fill_manual(values = c("0" = "#F8766D", "1" = "#00BFC4")) +
  theme_minimal(base_size = 12) +
  theme(axis.text.x = element_text(angle = 0, hjust = 0.5))

En cuanto a las importaciones observamos lo que ya habiamos comentado anteriormente, que todos estaban al alza al unirse y que una vez se unieron sus importaciones incrementaron. Pero aqui es donde notamos que no todos tienen el mismo comportamiento, algunos tienen comportamientos similares mientas que un comportamiento como el de Grecia no se asemeja a primera vista a los otros.