In USEPA/useeior: USEEIO R modeling software

This document presents validation results of r modelname model.

Build and Calculate Model

model <- buildModel(modelname)

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))
cat(paste("Number of sectors failing:",econval$N_Fail))
cat(paste("Sectors failing:", paste(econval$Failure$rownames, collapse = ", ")))

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))
cat(paste("Number of sectors failing:",econval$N_Fail))
cat(paste("Sectors failing:", paste(econval$Failure$rownames, 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 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))
cat(paste("Number of flow totals by commodity failing:",modelval$N_Fail))
#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))
cat(paste("Number of flow totals by commodity failing:",dom_val$N_Fail))
cat(paste("Sectors with flow totals failing:", paste(dom_val$Failure$variable, collapse = ", ")))

Validate that commodity output are properly transformed to industry output via MarketShare

q_x_val <- compareCommodityOutputXMarketShareandIndustryOutputwithCPITransformation(model, tolerance = 0.01)
cat(paste("Number of flow totals by commodity passing:",q_x_val$N_Pass))
cat(paste("Number of flow totals by commodity failing:",q_x_val$N_Fail))
cat(paste("Sectors with flow totals failing:", paste(q_x_val$Failure$rownames, collapse = ", ")))

Validate that commodity output equals to domestic use plus production demand

q_val <- compareCommodityOutputandDomesticUseplusProductionDemand(model, tolerance = 0.01)
cat(paste("Number of flow totals by commodity passing:",q_val$N_Pass))
cat(paste("Number of flow totals by commodity failing:",q_val$N_Fail))
cat(paste("Sectors with flow totals failing:", paste(q_val$Failure$rownames, collapse = ", ")))

USEPA/useeior documentation built on April 12, 2024, 1:36 p.m.