In USEPA/useeior: USEEIO R modeling software

This document presents comparison results of r modelname_pair[1] and r modelname_pair[2] model.

Compare flow totals between two models

mA <- buildModel(modelname_pair[1])
mB <- buildModel(modelname_pair[2])

# Compare flow totals
model_com <- compareFlowTotals(mA, mB)
cat(paste("Number of flow totals by commodity passing:",model_com$N_Pass))
cat(paste("Number of flow totals by commodity failing:",model_com$N_Fail))
#cat(paste("Sectors with flow totals failing:", paste(model_com$Failure$rownames, collapse = ", ")))

if (!is.null(model_com[["FlowDifference"]])) {
  cat("There are flow differences between", mA$specs$Model, "and", mB$specs$Model, "\n\n")
  cat(names(model_com[["FlowDifference"]])[1], "are\n\n")
  print(unlist(model_com[["FlowDifference"]][1], use.names = FALSE))
  cat("\n", names(model_com[["FlowDifference"]])[2], "are\n\n")
  print(unlist(model_com[["FlowDifference"]][2], use.names = FALSE))
}