Δήλωση ομάδας εργαστηρίου EYO

Ώρα & Μέρα προσέλευσης EYO 25-29/11/2024 (προετοιμασία 4η προόδου)

ΑΕΠ και ανεργία από την ΕΛΣΤΑΤ 


library(tidyverse)
library(readxl)
library(lubridate)
library(scales)
#library(modelsummary)
#library(gridExtra)
#library(ggpubr)

u1 <- "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(u1, destfile = "../data/GDP.xls")

u2 <- "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=116021&_documents_WAR_publicationsportlet_INSTANCE_Mr0GiQJSgPHd_locale=el"
download.file(u2, destfile = "../data/une.xls")


DF_g <- read_excel("../data/GDP.xls", skip = 4) %>% 
    select(gdp_adj = 3) %>% 
    filter(is.na(gdp_adj) == FALSE) %>% 
    mutate(Time = seq.Date(from = as.Date('1995-01-01'), by = "3 months", length.out = nrow(.))) %>% 
    select(Time, gdp_adj) %>% 
    mutate(gdp = 100 * (gdp_adj - lag(gdp_adj, 1)) / lag(gdp_adj, 1)) %>% 
    select(Time, gdp)


DF_u <- read_excel("../data/une.xls", skip = 2) %>% 
    select(une = 9) %>% 
    mutate(une = as.numeric(une)) %>% 
    filter(is.na(une) == FALSE) %>% 
    mutate(Date = seq.Date(from = as.Date('2004-01-01'), by = "1 months", length.out = nrow(.))) %>% 
    mutate(Time = floor_date(Date, unit = "quarter")) %>% 
    group_by(Time) %>% 
    summarise(une = mean(une), .groups = "drop") %>% 
    mutate(une = c(NA, diff(une))) %>% 
    filter(is.na(une) == FALSE)

DF <- inner_join(DF_g, DF_u, by = "Time") %>% 
    filter(is.na(gdp) == FALSE) %>% 
    filter(is.na(une) == FALSE)

lm.1 <- lm(gdp ~ une, data = DF)
lm.r <- broom::tidy(lm.1) %>% 
    mutate_if(is.numeric,  list(~round(., 3)))

DF %>% 
    ggplot(aes(x = une, y = gdp, colour = Time)) +
    geom_point(size = 3) +
    geom_smooth(method = "lm") +
    scale_x_continuous(breaks = pretty_breaks(10)) +
    scale_y_continuous(breaks = pretty_breaks(10)) +
    labs(x = "Μεταβολή στο ποσοστό ανεργίας",
         y = "% Μεταβολή μεταβολή του ΑΕΠ ", 
         title = "Νόμος Okun",
         subtitle = "Τριμήνια στοιχεία με εποχική και ημερολογιακή διόρθωση",
         caption = "Πηγή: ΕΛΣΤΑΤ\nΕπεξεργασία: Α. Σταυρακούδης",
         colour = "Χρόνος") +
    theme_minimal(base_size = 16) +
    #theme(legend.position = "none") +
    theme(legend.key.height = unit(1, 'cm')) +
    #annotation_custom(tableGrob(lm.r), xmin = -1, ymin = -10) +
    stat_regline_equation(label.x = 1, label.y = 4, size = 6) +
    NULL

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

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


Athanassios Stavrakoudis, Sat 23 Nov 2024