| title | date | output | editor_options | ||||||
|---|---|---|---|---|---|---|---|---|---|
Validating USEEIOv2.0 Model |
2021-12-30 |
|
|
This document presents validation results of USEEIOv2.0 model.
model <- buildModel(modelname)
#> 2021-12-30 09:29:51 INFO::Begin model initialization...
#> 2021-12-30 09:29:52 INFO::Initializing IO tables...
#> 2021-12-30 09:29:52 INFO::Initializing Gross Output tables...
#> 2021-12-30 09:29:54 INFO::Initializing Chain Price Index tables...
#> 2021-12-30 09:29:55 INFO::Initializing model satellite tables...
#> 2021-12-30 09:29:55 INFO::Loading Water withdrawals flows from DataCommons...
#> 2021-12-30 09:29:55 INFO::Loading Criteria and Hazardous Air Emissions flows from DataCommons...
#> 2021-12-30 09:30:06 INFO::Loading Point source industrial releases to ground flows from DataCommons...
#> 2021-12-30 09:30:06 INFO::Loading Point source releases to water flows from DataCommons...
#> 2021-12-30 09:30:07 INFO::Loading Greenhouse Gases flows from DataCommons...
#> 2021-12-30 09:30:07 INFO::Loading Land use flows from DataCommons...
#> 2021-12-30 09:30:07 INFO::Loading Mineral extraction flows from DataCommons...
#> 2021-12-30 09:30:07 INFO::Loading Energy extraction flows from DataCommons...
#> 2021-12-30 09:30:08 INFO::Loading Nitrogen and Phosphorus Releases from Agriculture flows from DataCommons...
#> 2021-12-30 09:30:08 INFO::Loading Pesticide releases flows from DataCommons...
#> 2021-12-30 09:30:08 INFO::Loading Commercial non-hazardous waste excluding construction activities flows from DataCommons...
#> 2021-12-30 09:30:09 INFO::Loading Commercial non-hazardous waste from construction activities flows from DataCommons...
#> 2021-12-30 09:30:09 INFO::Loading Commercial RCRA-defined hazardous waste flows from DataCommons...
#> 2021-12-30 09:30:11 INFO::Loading Employment flows from DataCommons...
#> 2021-12-30 09:30:12 INFO::Generating Value Added flows...
#> 2021-12-30 09:30:12 INFO::Initializing model indicators...
#> 2021-12-30 09:30:12 INFO::Getting Greenhouse Gases indicator from DataCommons...
#> 2021-12-30 09:30:12 INFO::Getting Acidification Potential indicator from DataCommons...
#> 2021-12-30 09:30:12 INFO::Getting Eutrophication Potential indicator from DataCommons...
#> 2021-12-30 09:30:12 INFO::Getting Freshwater Ecotoxicity Potential indicator from DataCommons...
#> 2021-12-30 09:30:18 INFO::Getting Human Health - Cancer indicator from DataCommons...
#> 2021-12-30 09:30:20 INFO::Getting Human Health - Noncancer indicator from DataCommons...
#> 2021-12-30 09:30:21 INFO::Getting Human Health Toxicity indicator from DataCommons...
#> 2021-12-30 09:30:23 INFO::Getting Human Health - Respiratory Effects indicator from DataCommons...
#> 2021-12-30 09:30:24 INFO::Getting Ozone Depletion indicator from DataCommons...
#> 2021-12-30 09:30:24 INFO::Getting Smog Formation Potential indicator from DataCommons...
#> 2021-12-30 09:30:25 INFO::Getting Freshwater withdrawals indicator from DataCommons...
#> 2021-12-30 09:30:25 INFO::Getting Land use indicator from DataCommons...
#> 2021-12-30 09:30:25 INFO::Getting Hazardous Air Pollutants indicator from DataCommons...
#> 2021-12-30 09:30:26 INFO::Getting Pesticides indicator from DataCommons...
#> 2021-12-30 09:30:26 INFO::Getting Nonrenewable Energy Use indicator from DataCommons...
#> 2021-12-30 09:30:26 INFO::Getting Renewable Energy Use indicator from DataCommons...
#> 2021-12-30 09:30:26 INFO::Getting Energy Use indicator from DataCommons...
#> 2021-12-30 09:30:27 INFO::Getting Minerals and Metals Use indicator from DataCommons...
#> 2021-12-30 09:30:27 INFO::Getting Value Added indicator from useeior...
#> 2021-12-30 09:30:27 INFO::Getting Jobs Supported indicator from useeior...
#> 2021-12-30 09:30:27 INFO::Getting Commercial RCRA Hazardous Waste indicator from useeior...
#> 2021-12-30 09:30:27 INFO::Getting Commercial Municipal Solid Waste indicator from useeior...
#> 2021-12-30 09:30:28 INFO::Getting Commercial Construction and Demolition Debris indicator from useeior...
#> 2021-12-30 09:30:28 INFO::Loading demand vectors ...
#> 2021-12-30 09:30:28 INFO::Loading CompleteProduction demand vector...
#> 2021-12-30 09:30:28 INFO::Loading DomesticProduction demand vector...
#> 2021-12-30 09:30:28 INFO::Loading CompleteConsumption demand vector...
#> 2021-12-30 09:30:28 INFO::Loading DomesticConsumption demand vector...
#> 2021-12-30 09:30:33 INFO::Building commodity-by-commodity A matrix (direct requirements)...
#> 2021-12-30 09:30:33 INFO::Building commodity-by-commodity A_d matrix (domestic direct requirements)...
#> 2021-12-30 09:30:34 INFO::Calculating L matrix (total requirements)...
#> 2021-12-30 09:30:34 INFO::Calculating L_d matrix (domestic total requirements)...
#> 2021-12-30 09:30:34 INFO::Building B matrix (direct emissions and resource use per dollar)...
#> 2021-12-30 09:30:37 INFO::Building C matrix (characterization factors for model indicators)...
#> 2021-12-30 09:30:42 INFO::Calculating D matrix (direct environmental impacts per dollar)...
#> 2021-12-30 09:30:42 INFO::Calculating M matrix (total emissions and resource use per dollar)...
#> 2021-12-30 09:30:42 INFO::Calculating M_d matrix (total emissions and resource use per dollar from domestic activity)...
#> 2021-12-30 09:30:42 INFO::Calculating N matrix (total environmental impacts per dollar)...
#> 2021-12-30 09:30:42 INFO::Calculating N_d matrix (total environmental impacts per dollar from domestic activity)...
#> 2021-12-30 09:30:42 INFO::Calculating Rho matrix (price year ratio)...
#> 2021-12-30 09:30:42 INFO::Calculating Phi matrix (producer over purchaser price ratio)...
#> 2021-12-30 09:30:42 INFO::Model build complete.Validate that commodity output can be recalculated (within 1%) with the model total requirements matrix (L) and demand vector (y) for US production
econval <- compareOutputandLeontiefXDemand(model, tolerance = 0.01)
cat(paste("Number of sectors passing:",econval$N_Pass))Number of sectors passing: 403
cat(paste("Number of sectors failing:",econval$N_Fail))Number of sectors failing: 2
cat(paste("Sectors failing:", paste(econval$Failure$rownames, collapse = ", ")))Sectors failing: S00402/US, S00300/US
Validate that commodity output can be recalculated (within 1%) with model total domestic requirements matrix (L_d) and model demand (y) for US production
econval <- compareOutputandLeontiefXDemand(model,use_domestic=TRUE, tolerance = 0.01)
cat(paste("Number of sectors passing:",econval$N_Pass))Number of sectors passing: 403
cat(paste("Number of sectors failing:",econval$N_Fail))Number of sectors failing: 2
cat(paste("Sectors failing:", paste(econval$Failure$rownames, collapse = ", ")))Sectors failing: S00402/US, S00300/US
Validate that flow totals by commodity (E_c) can be recalculated (within 1%) using the model satellite matrix (B), market shares matrix (V_n), total requirements matrix (L), and demand vector (y) for US production
modelval <- compareEandLCIResult(model, tolerance = 0.01)
cat(paste("Number of flow totals by commodity passing:",modelval$N_Pass))Number of flow totals by commodity passing: 1080540
cat(paste("Number of flow totals by commodity failing:",modelval$N_Fail))Number of flow totals by commodity failing: 0
#cat(paste("Sectors failing:", paste(modelval$Failure$variable, collapse = ", ")))Validate that flow totals by commodity (E_c) can be recalculated (within 1%) using the model satellite matrix (B), market shares matrix (V_n), total domestic requirements matrix (L_d), and demand vector (y) for US production
dom_val <- compareEandLCIResult(model,use_domestic=TRUE, tolerance = 0.01)
cat(paste("Number of flow totals by commodity passing:",dom_val$N_Pass))Number of flow totals by commodity passing: 1080540
cat(paste("Number of flow totals by commodity failing:",dom_val$N_Fail))Number of flow totals by commodity failing: 0
cat(paste("Sectors with flow totals failing:", paste(dom_val$Failure$variable, collapse = ", ")))Sectors with flow totals failing:
q_x_val <- compareCommodityOutputXMarketShareandIndustryOutputwithCPITransformation(model, tolerance = 0.01)
cat(paste("Number of flow totals by commodity passing:",q_x_val$N_Pass))Number of flow totals by commodity passing: 403
cat(paste("Number of flow totals by commodity failing:",q_x_val$N_Fail))Number of flow totals by commodity failing: 2
cat(paste("Sectors with flow totals failing:", paste(q_x_val$Failure$rownames, collapse = ", ")))Sectors with flow totals failing: S00402/US, S00300/US
q_val <- compareCommodityOutputandDomesticUseplusProductionDemand(model, tolerance = 0.01)
cat(paste("Number of flow totals by commodity passing:",q_val$N_Pass))Number of flow totals by commodity passing: 404
cat(paste("Number of flow totals by commodity failing:",q_val$N_Fail))Number of flow totals by commodity failing: 1
cat(paste("Sectors with flow totals failing:", paste(q_val$Failure$rownames, collapse = ", ")))Sectors with flow totals failing: S00300/US