Ακαθάριστο Εθνικό Προϊόν από την ΕΛΣΤΑΤ

## Ακαθάριστο Εθνικό Προϊόν
## https://www.statistics.gr/el/statistics/-/publication/SEL84/-

library(tidyverse)
library(lubridate)
library(readxl)

u <- "https://www.statistics.gr/el/statistics?p_p_id=documents_WAR_publicationsportlet_INSTANCE_Mr0GiQJSgPHd&p_p_lifecycle=2&p_p_state=normal&p_p_mode=view&p_p_cacheability=cacheLevelPage&p_p_col_id=column-2&p_p_col_count=4&p_p_col_pos=2&_documents_WAR_publicationsportlet_INSTANCE_Mr0GiQJSgPHd_javax.faces.resource=document&_documents_WAR_publicationsportlet_INSTANCE_Mr0GiQJSgPHd_ln=downloadResources&_documents_WAR_publicationsportlet_INSTANCE_Mr0GiQJSgPHd_documentID=115388&_documents_WAR_publicationsportlet_INSTANCE_Mr0GiQJSgPHd_locale=el"
download.file(u, destfile = "../data/GDP.xls")

df <- read_excel("../data/GDP.xls", skip = 4) %>% 
    select(gdp_adj = 3, gdp_cp = 6)

dates <- seq.Date(from = as.Date('1995-01-01'), by = "3 months", length.out = nrow(df))

df %>% 
    mutate(Time = dates) %>% 
    select(Time, gdp_adj, gdp_cp)

DF <- read_excel("../data/GDP.xls", skip = 4) %>% 
    select(gdp_adj = 3, gdp_cp = 6) %>% 
    filter(is.na(gdp_adj) == FALSE) %>% 
    mutate(Time = seq.Date(from = as.Date('1995-01-01'), by = "3 months", length.out = nrow(.))) %>% 
    mutate(gdp_adj = 1e-3*gdp_adj, gdp_cp = 1e-3*gdp_cp) %>% 
    select(Time, gdp_adj, gdp_cp) %>% 
    mutate(qoq_adj = 100 * (gdp_adj - lag(gdp_adj, 1)) / lag(gdp_adj, 1)) %>% 
    mutate(yoy_adj = 100 * (gdp_adj - lag(gdp_adj, 4)) / lag(gdp_adj, 4)) %>% 
    mutate(qoq_cp = 100 * (gdp_cp - lag(gdp_cp, 1)) / lag(gdp_cp, 1)) %>% 
    mutate(yoy_cp = 100 * (gdp_cp - lag(gdp_cp, 4)) / lag(gdp_cp, 4)) %>% 
    pivot_longer(-Time, names_to = "key", values_to = "value") %>% 
    identity()


DF %>% 
    filter(key == 'gdp_adj') %>% 
    ggplot(aes(x = Time, y = value)) +
    geom_line()


DF %>% 
    filter(key == 'yoy_adj') %>% 
    ggplot(aes(x = Time, y = value)) +
    geom_line()


DF %>% 
    filter(key == 'yoy_adj') %>% 
    slice_tail(n = 8) %>% 
    ggplot(aes(x = Time, y = value)) +
    geom_col()

DF %>% 
    filter(key == 'yoy_adj') %>%
    mutate(type = ifelse(value < 0, "t1", "t2")) %>% 
    slice_tail(n = 12) %>% 
    ggplot(aes(x = Time, y = value, fill = type)) +
    geom_col() +
    scale_y_continuous(breaks = scales::pretty_breaks(10)) +
    labs(x = "", y = "% Μεταβολή", 
         title = "Μεταβολή του ΑΕΠ (y-o-y)",
         subtitle = "Τριμήνια στοιχεία με εποχική και ημερολογιακή διόρθωση",
         caption = "Πηγή: ΕΛΣΤΑΤ  https://www.statistics.gr/el/statistics/-/publication/SEL84/-\nΕπεξεργασία: Α. Σταυρακούδης") +
    theme_minimal(base_size = 16) +
    theme(legend.position = "none") 




Συνδεθείτε για περισσότερες δυνατότητες αλληλεπίδρασης,
σχολιασμοί, εξωτερικοί σύνδεσμοι, βοήθεια, ψηφοφορίες, αρχεία, κτλ.

Creative Commons License
Εκπαιδευτικό υλικό από τον Αθανάσιο Σταυρακούδη σας παρέχετε κάτω από την άδεια Creative Commons Attribution-NonCommercial-ShareAlike 4.0 License.
Σας παρακαλώ να ενημερωθείτε για κάποιους επιπλέον περιορισμούς
http://stavrakoudis.econ.uoi.gr/stavrakoudis/?iid=401.