Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

derive an M matrix for IEF models #317

Merged
merged 9 commits into from
Oct 30, 2024
Merged

derive an M matrix for IEF models #317

Changes from 1 commit
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
a35b2c7
derive an M matrix for IEF models #311
bl-young Sep 4, 2024
119ea9c
use domestic production demand and import consumption demand to creat…
bl-young Sep 10, 2024
04b7a8e
revise approach for deriving M to use q and y_m
bl-young Sep 10, 2024
7c7f855
derive m from use table
bl-young Sep 10, 2024
adb2daf
fix error in deriving M matrix for 2R models
bl-young Sep 10, 2024
b80c38b
ignore location when removing imports
bl-young Sep 10, 2024
c9762f4
validate that M is between M_m and M_d
bl-young Sep 11, 2024
e232ca1
move M matrix test to validation and print any failures
bl-young Sep 13, 2024
0910392
Merge branch 'develop' into issue_311
bl-young Oct 30, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 2 additions & 2 deletions R/ExternalImportFactors.R
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -135,14 +135,14 @@ buildModelwithImportFactors <- function(model, configpaths = NULL) {
#' @param model, An EEIO model object with model specs and crosswalk table loaded
#' @return An M matrix of flows x sector
deriveMMatrix <- function(model) {
# logging::loginfo("Deriving M matrix (total emissions and resource use per dollar) ...")
logging::loginfo("Deriving M matrix (total emissions and resource use per dollar) ...")
q <- model$q
loc <- grepl(model$specs$ModelRegionAcronyms[1], model$FinalDemandMeta$Code_Loc)
bl-young marked this conversation as resolved.
Show resolved Hide resolved
import_code <- model$FinalDemandMeta[model$FinalDemandMeta$Group=="Import" & loc, "Code_Loc"]
# derive total imports (m) from the Use table
U_m <- model$U - model$U_d
# Exclude imports col when calculating total imports
m <- head(rowSums(U_m[, !(colnames(U_m) %in% import_code)]), -3) # drop VA
m <- rowSums(U_m[model$Commodities$Code_Loc, !(colnames(U_m) %in% import_code)]) # drop VA

dr <- q / (q + m)
mr <- 1 - dr
Expand Down