Nothing
## ----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)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.