Working With Eurostat Data

# Not run
not_included_directory <- file.path("..", "not_included")
if (!dir.exists(not_included_directory)) dir.create(not_included_directory)
# The contents of the 'not_included' directory can be found on GitHub,
# but they are not released and distributed with the package.

naio_10_cp1700 <- iotables_download(
  "naio_10_cp1700", # SIOT
  data_directory = not_included_directory
)

# For inclusion in the package, the files must be smaller.
# Reducing the size of the bulk files will not affect
# the demonstration.

naio_10_cp1700 <- naio_10_cp1700 %>%
  dplyr::filter(geo %in% c("CZ", "SK")) %>%
  dplyr::filter(year %in% c(2010, 2015))

# Conforming employment data both sexes from 15 years old, year 2015.
# prod_na vocabulary for product x product conformity
emp_cz <- employment_get(
  geo = "CZ", year = "2015", sex = "Total",
  age = "Y_GE15", labelling = "prod_na",
  data_directory = not_included_directory,
  force_download = TRUE
)

# Conforming employment data #both sexes from 15 years old, year 2017.
emp_sk <- employment_get(
  geo = "SK",
  year = "2017", sex = "Total",
  age = "Y_GE15", labelling = "prod_na",
  data_directory = not_included_directory,
  force_download = TRUE
)

save(naio_10_cp1700, emp_sk, emp_cz,
  file = file.path(
    "..", "inst", "extdata",
    "naio_10_product_x_product.rda"
  )
)

# load from pre-saved file to increase running speed
load(system.file("extdata",
  "naio_10_product_x_product.rda",
  package = "iotables"
))

cz_io <- iotable_get(
  source = "naio_10_cp1700",
  geo = "CZ",
  year = 2015,
  unit = "MIO_NAC",
  stk_flow = "TOTAL",
  labelling = "short"
)

sk_io <- iotable_get(
  source = "naio_10_cp1700",
  geo = "SK",
  year = 2015,
  unit = "MIO_EUR",
  stk_flow = "TOTAL",
  labelling = "short"
)

cz_input_flow <- input_flow_get(data_table = cz_io)

sk_input_flow <- input_flow_get(data_table = sk_io)

cz_output <- output_get(data_table = cz_io)
sk_output <- output_get(data_table = sk_io)

primary_inputs_cz <- coefficient_matrix_create(
  data_table = cz_io,
  total = "output",
  return = "primary_inputs"
)

primary_inputs_sk <- coefficient_matrix_create(
  data_table = sk_io,
  total = "output",
  return = "primary_inputs"
)
#> Columns and rows of CPA_L68A, CPA_T98, CPA_U are all zeros and will be removed.

direct_cz <- direct_effects_create(
  input_requirements = primary_inputs_cz,
  inverse = I_cz,
  digits = 4
)
direct_sk <- direct_effects_create(primary_inputs_sk, I_sk)

knitr::kable(head(direct_cz[, 1:8]), digits = 4)

cz_bw <- backward_linkages(I_cz)
sk_bw <- backward_linkages(I_sk)

knitr::kable(head(cz_bw[, 1:8]), digits = 4)

output_coeff_cz <- output_coefficient_matrix_create(
  data_table = cz_io, total = "TFU", digits = 4
)

output_coeff_sk <- output_coefficient_matrix_create(
  data_table = sk_io, total = "TFU", digits = 4
)
#> Columns and rows of CPA_L68A, CPA_T98, CPA_U are all zeros and will be removed.

knitr::kable(head(output_coeff_cz[, 1:8]))

cz_fw <- forward_linkages(output_coeff_cz)
sk_fw <- forward_linkages(output_coeff_sk)

knitr::kable(head(cz_fw), digits = 4)

require(xlsx)
cz_file_name <- file.path("..", "not_included", "CzechRep_test.xlsx")
# Czech Republic data
xlsx::write.xlsx(cz_io,
  file = cz_file_name, sheetName = "io_table",
  col.names = TRUE, row.names = TRUE, append = FALSE
)
xlsx::write.xlsx(cz_output,
  file = cz_file_name, sheetName = "cz_output",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(input_coeff_matrix_cz,
  file = cz_file_name,
  sheetName = "input_coeff_matrix_cz",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(L_cz,
  file = cz_file_name, sheetName = "L_cz",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(I_cz,
  file = cz_file_name, sheetName = "I_cz",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(direct_cz,
  file = cz_file_name,
  sheetName = "direct_cz",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(multipliers_cz,
  file = cz_file_name,
  sheetName = "multipliers_cz",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(emp_effect_cz,
  file = cz_file_name,
  sheetName = "emp_effect_cz_2015",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(emp_indicator_cz,
  file = cz_file_name,
  sheetName = "emp_indicator_cz",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(emp_multiplier_cz,
  file = cz_file_name,
  sheetName = "emp_multiplier_cz",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(cz_bw,
  file = cz_file_name,
  sheetName = "cz_backward_linkages",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(output_coeff_cz,
  file = cz_file_name,
  sheetName = "output_coeff_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(cz_fw,
  file = cz_file_name,
  sheetName = "cz_forward_linkages",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(output_multipliers_cz,
  file = cz_file_name,
  sheetName = "output_multipliers_cz",
  col.names = TRUE, row.names = TRUE, append = TRUE
)


sk_file_name <- file.path("..", "not_included", "SlovakRep_test.xlsx")
# Czech Republic data
xlsx::write.xlsx(sk_io,
  file = sk_file_name, sheetName = "io_table",
  col.names = TRUE, row.names = TRUE, append = FALSE
)
xlsx::write.xlsx(sk_output,
  file = sk_file_name, sheetName = "sk_output",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(input_coeff_matrix_sk,
  file = sk_file_name,
  sheetName = "input_coeff_matrix_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(L_sk,
  file = sk_file_name, sheetName = "L_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(I_sk,
  file = sk_file_name, sheetName = "I_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(direct_sk,
  file = sk_file_name,
  sheetName = "direct_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(multipliers_sk,
  file = sk_file_name,
  sheetName = "multipliers_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(emp_effect_sk,
  file = sk_file_name,
  sheetName = "emp_effect_sk_2015",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(emp_indicator_sk,
  file = sk_file_name,
  sheetName = "emp_indicator_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(emp_multiplier_sk,
  file = sk_file_name,
  sheetName = "employment_multipliers_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(sk_bw,
  file = sk_file_name,
  sheetName = "sk_backward_linkages",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(output_coeff_sk,
  file = sk_file_name,
  sheetName = "output_coeff_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(sk_fw,
  file = sk_file_name,
  sheetName = "sk_forward_linkages",
  col.names = TRUE, row.names = TRUE, append = TRUE
)
xlsx::write.xlsx(output_multipliers_sk,
  file = sk_file_name,
  sheetName = "output_multipliers_sk",
  col.names = TRUE, row.names = TRUE, append = TRUE
)