## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
library(iotables)
require(dplyr)

## ----download, eval=FALSE-----------------------------------------------------
# # 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---------------------------------------------------------------------
# load from pre-saved file to increase running speed
load(system.file("extdata",
  "naio_10_product_x_product.rda",
  package = "iotables"
))

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

sk_io <- iotable_get(
  labelled_io_data = naio_10_cp1700,
  source = "naio_10_cp1700", geo = "SK",
  year = 2010, 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)

## ----inputcoeff, results='asis'-----------------------------------------------
input_coeff_matrix_cz <- input_coefficient_matrix_create(
  data_table = cz_io
)

input_coeff_matrix_sk <- input_coefficient_matrix_create(
  data_table = sk_io
)

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

## ----leontieff, results='asis'------------------------------------------------
L_cz <- leontieff_matrix_create(input_coeff_matrix_cz)
I_cz <- leontieff_inverse_create(input_coeff_matrix_cz)

L_sk <- leontieff_matrix_create(input_coeff_matrix_sk)
I_sk <- leontieff_inverse_create(input_coeff_matrix_sk)

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

## ----direct, results='asis'---------------------------------------------------
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"
)

direct_cz <- direct_effects_create(primary_inputs_cz, I_cz)
direct_sk <- direct_effects_create(primary_inputs_sk, I_sk)

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

## ----total, results='asis'----------------------------------------------------
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"
)

multipliers_cz <- input_multipliers_create(primary_inputs_cz, I_cz)
multipliers_sk <- input_multipliers_create(primary_inputs_sk, I_sk)

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

## ----employmenteffect, results='asis', message=FALSE--------------------------
# New function is needed to add employment vector to SIOT
names(emp_sk)[1] <- "prod_na"
names(emp_cz)[1] <- "prod_na"

emp_indicator_sk <- rbind(
  sk_io[, 1:66],
  emp_sk
) %>%
  coefficient_matrix_create(.,
    return_part = "primary_inputs"
  ) %>%
  filter(prod_na == "employment_total")

emp_indicator_cz <- full_join(
  cz_io,
  emp_cz
) %>%
  coefficient_matrix_create(.,
    return_part = "primary_inputs"
  ) %>%
  filter(prod_na == "employment_total")


emp_effect_sk <- direct_effects_create(emp_indicator_sk, I_sk)
emp_effect_cz <- direct_effects_create(emp_indicator_cz, I_cz)

knitr::kable(emp_effect_cz[1:8], digits = 5)

## ----employmentindicator, results='asis'--------------------------------------
# New function is needed to add employment vector to SIOT

emp_multiplier_sk <- input_multipliers_create(emp_indicator_sk, I_sk)
emp_multiplier_cz <- input_multipliers_create(emp_indicator_cz, I_cz)

knitr::kable(emp_multiplier_cz[1:8], digits = 5)

## ----output_multipliers, results='asis'---------------------------------------
output_multipliers_cz <- output_multiplier_create(input_coeff_matrix_cz)
output_multipliers_sk <- output_multiplier_create(input_coeff_matrix_sk)

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

## ----backward, results='asis'-------------------------------------------------
cz_bw <- backward_linkages(I_cz)
sk_bw <- backward_linkages(I_sk)

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

## ----output_coeff, results='asis', eval=FALSE---------------------------------
# output_coeff_cz <- output_coefficient_matrix_create(
#   io_table = cz_io, total = "tfu", digits = 4
# )
# 
# output_coeff_sk <- output_coefficient_matrix_create(
#   io_table = sk_io, total = "tfu"
# )
# 
# knitr::kable(head(output_coeff_cz[, 1:8]))

## ----forward, results='asis', eval=FALSE--------------------------------------
# cz_fw <- forward_linkages(output_coeff_cz)
# sk_fw <- forward_linkages(output_coeff_sk)
# 
# knitr::kable(head(cz_fw), digits = 4)

## ----reproduction_data, eval=FALSE--------------------------------------------
# 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
# )