Nothing
## ----setup, include=FALSE-----------------------------------------------------
knitr::opts_chunk$set(
collapse = TRUE,
comment = "#>"
)
library(dplyr, quietly=TRUE)
library(tidyr, quietly=TRUE)
library(testthat, quietly=TRUE)
## ----read_control_data, eval=FALSE--------------------------------------------
# uk_2010_data <- iotables_download ( source = "uk_2010" )
# save ( uk_2010_data, file = file.path('data-raw', 'uk_2010_data.rda'))
# uk_test_results <- iotables:::uk_2010_results_get ()
# #saved as package data
## ----downloaded_excel_data----------------------------------------------------
library(iotables)
data(uk_2010_data)
uk_siot <- iotable_get (
labelled_io_data = uk_2010_data,
source = "uk_2010_siot" )
uk_published_coeff <- iotable_get (
labelled_io_data = uk_2010_data,
source = "uk_2010_coeff" )
uk_published_inverse <- iotable_get (
labelled_io_data = uk_2010_data,
source = "uk_2010_inverse" )
## ----compare_coeff------------------------------------------------------------
uk_input_coeff <- input_coefficient_matrix_create(data_table = uk_siot)
coeff_comparison <- select (uk_input_coeff, 1 ) %>%
left_join ( uk_published_coeff, by = "prod_na")
test_that("correct data is returned", {
expect_equal(round(uk_input_coeff[,2:8], 8),
round(coeff_comparison [,2:8], 8)) })
## ----compare_inverse----------------------------------------------------------
uk_calculated_inverse <- leontief_inverse_create(uk_input_coeff)
inverse_comparison <- select (uk_calculated_inverse, 1 ) %>%
left_join (uk_calculated_inverse, by = "prod_na")
test_that("correct data is returned", {
expect_equal(round(uk_calculated_inverse[,2:8], 8),
round(inverse_comparison[,2:8], 8)) })
## ----compare_effects----------------------------------------------------------
employment_effect_results <- uk_test_results %>%
select ( uk_row_label, `Employment cost effects`)
primary_inputs_uk <- coefficient_matrix_create(
data_table = uk_siot,
total = 'output',
return_part = 'primary_inputs')
employment_input <- filter (primary_inputs_uk , prod_na == "D1")
employment_effects <- direct_effects_create( employment_input, uk_calculated_inverse ) %>%
gather ( prod, values, !!2:ncol(.)) %>%
mutate ( prod_na = prod ) %>%
select ( -prod ) %>%
left_join ( select ( metadata_uk_2010, prod_na, uk_row_label ),
by = 'prod_na') %>%
left_join ( employment_effect_results, by = 'uk_row_label' ) %>%
filter ( !is.na(uk_row_label )) %>%
select ( prod_na, uk_row_label, values, `Employment cost effects`)
iotables:::create_knitr_table (
data_table = employment_effects[1:10,],
digits = 4,
caption = "Comparison of Calculated And Published Employment Cost Effects",
col.names = c("industry code", "row label", "calculated", "published"),
col_width = c(2,11,3,3))
## ----gva_effects--------------------------------------------------------------
uk_siot2 <- uk_siot %>%
filter ( prod_na %in% c("B2A3G", "D1", "D29X39") ) %>%
summarize_if ( is.numeric, sum, na.rm = TRUE ) %>%
cbind ( data.frame ( prod_na = "GVA"), . ) %>%
rbind ( uk_siot, .)
gva_effect_results <- uk_test_results %>%
select ( uk_row_label, `GVA effects`)
gva_input <- coefficient_matrix_create(
data_table = uk_siot2,
total = 'output',
return_part = 'primary_inputs') %>%
filter ( prod_na == "GVA" )
gva_effects <- direct_effects_create( gva_input,
uk_calculated_inverse ) %>%
gather ( prod, values, !!2:ncol(.)) %>%
mutate ( prod_na = prod ) %>%
select ( -prod ) %>%
left_join ( select ( metadata_uk_2010, prod_na, uk_row_label ),
by = 'prod_na') %>%
left_join ( gva_effect_results, by = 'uk_row_label' ) %>%
filter ( !is.na(uk_row_label )) %>%
select ( prod_na, uk_row_label, values, `GVA effects`)
iotables:::create_knitr_table (
data_table = gva_effects[1:10,],
digits = 4,
caption = "Comparison of Calculated And Published GVA Effects",
col.names = c("industry code", "row label", "calculated", "published"),
col_width = c(2,11,3,3))
## ----compare_emp_multipliers--------------------------------------------------
empc_multiplier_results <- uk_test_results %>%
select ( uk_row_label, `Employment cost multiplier`)
empc_indicator_uk <- coefficient_matrix_create(
data_table = uk_siot,
total = 'output',
return_part = 'primary_inputs') %>%
filter ( prod_na == 'D1')
empc_multipliers <- input_multipliers_create(
input_requirements = empc_indicator_uk,
uk_calculated_inverse) %>%
gather ( prod, values, !!2:ncol(.)) %>%
mutate ( prod_na = prod ) %>%
select ( -prod ) %>%
left_join ( select ( metadata_uk_2010, prod_na, uk_row_label ),
by = 'prod_na') %>%
left_join ( empc_multiplier_results, by = 'uk_row_label' ) %>%
filter ( !is.na(uk_row_label )) %>%
select ( prod_na, uk_row_label, values, `Employment cost multiplier`)
iotables:::create_knitr_table (
data_table = empc_multipliers [1:10,], digits = 4,
caption = "Comparison of Calculated And Published Employment Cost Multipliers",
col.names = c("industry code", "row label", "calculated", "published"),
col_width = c(2,11,3,3))
## ----gva_comp, eval=FALSE-----------------------------------------------------
# gva_multipliers <- input_multipliers_create(
# input_requirements = gva_input,
# uk_calculated_inverse) %>%
# gather ( prod, values, !!2:ncol(.)) %>%
# mutate ( prod_na = prod ) %>%
# select ( -prod ) %>%
# left_join ( select ( metadata_uk_2010, prod_na, uk_row_label ),
# by = 'prod_na') %>%
# left_join ( gva_multiplier_results, by = 'uk_row_label' ) %>%
# filter ( !is.na(uk_row_label )) %>%
# select ( prod_na, uk_row_label, values, `GVA multiplier`)
#
# iotables:::create_knitr_table (
# data_table = gva_multipliers [1:10,], digits = 4,
# caption = "Comparison of Calculated And Published GVA Multipliers",
# col.names = c("industry code", "row label", "calculated", "published"),
# col_width = c(2,11,3,3))
#
# iotables:::create_knitr_table (
# data_table = gva_multipliers [1:10,],
# digits = 4,
# caption = "Comparison of Calculated And Published GVA Multipliers",
# col.names = c("industry code", "row label",
# "calculated", "published"),
# col_width = c(2,11,3,3))
## ----compare_output_multipliers-----------------------------------------------
output_multiplier_results <- uk_test_results %>%
select ( uk_row_label, `Output multiplier`)
uk_output_multipliers <- output_multiplier_create(uk_input_coeff) %>%
gather ( prod, values, !!2:ncol(.)) %>%
mutate ( prod_na = prod ) %>%
select ( -prod ) %>%
left_join ( select ( metadata_uk_2010, prod_na, uk_row_label ),
by = 'prod_na') %>%
left_join ( output_multiplier_results,
by = 'uk_row_label' ) %>%
filter ( !is.na(uk_row_label) ) %>%
select ( prod_na, uk_row_label, values, `Output multiplier`)
iotables:::create_knitr_table (
data_table = uk_output_multipliers [1:10,],
digits = 4,
caption = "Comparison of Calculated And Published Output Multipliers",
col.names = c("industry code", "row label",
"calculated", "published"),
col_width = c(2,11,3,3))
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.