myprojects-hg-reborn

Fork

rm(list=ls())

library(tidyverse)
library(janitor)
library(openxlsx)


source("/home/lorenzo/myprojects-hg/R-codes/stat_lib.R")

uk_ini <- read_excel("Copy_of_UK_measures.xlsx")

uk_series <- uk_ini |>
    summarise(aid_element_mio_eur=sum(AID_ELEMENT_EUR, na.rm=T),
              aid_element_deflated_mio_eur=sum(AID_ELEMENT_EUR_ADJ, na.rm=T),
              .by=c(EXPENDITURE_YEAR)) |>
    arrange(desc(EXPENDITURE_YEAR))

save_excel(uk_series, "UK-time-series.xlsx")

uk_series1 <- uk_series |>
    select(-aid_element_deflated_mio_eur) |>
    pivot_wider(names_from=EXPENDITURE_YEAR,
                values_from=aid_element_mio_eur,
                names_sort=T)


uk_series2 <- uk_series |>
    select(-aid_element_mio_eur) |>
    pivot_wider(names_from=EXPENDITURE_YEAR,
                values_from=aid_element_deflated_mio_eur,
                names_sort=T)


uk_series3 <- bind_rows(uk_series1, uk_series2)

extra <- tibble(year=c("aid element (MIO EUR) current prices",
                       "aid element (MIO EUR) constant prices"))

uk_series4 <- bind_cols(extra, uk_series3)

save_excel(uk_series4, "UK-series-wide.xlsx")

pan_ini <- read_excel("Copy_of_pan_european_guarantee_fund_measures.xlsx") 


pan_agg <- pan_ini |>
        summarise(aid_element_mio_eur=sum(AID_ELEMENT_EUR, na.rm=T),
              aid_element_deflated_mio_eur=sum(AID_ELEMENT_EUR_ADJ, na.rm=T),
              .by=c(EXPENDITURE_YEAR, MEMBER_STATE)) |>
    arrange(MEMBER_STATE,desc(EXPENDITURE_YEAR)) |>
    filter(aid_element_mio_eur!=0 & aid_element_deflated_mio_eur!=0)

save_excel(pan_agg, "pan-european_fund.xlsx")

pan_agg_wide_current <- pan_agg |>
    clean_names() |>
    select(-aid_element_deflated_mio_eur) |> 
    pivot_wider(names_from=expenditure_year,
                values_from=aid_element_mio_eur, values_fill=0,
                names_sort=T)


save_excel(pan_agg_wide_current, "pan_european_wide.xlsx", "pan_European_fund_MIO_EUR")

pan_agg_wide_constant <- pan_agg |>
    clean_names() |>
    select(-aid_element_mio_eur) |> 
    pivot_wider(names_from=expenditure_year,
                values_from=aid_element_deflated_mio_eur, values_fill=0,
                names_sort=T)


print("So far so good")