Fixing bug in Ghosh B

USEPA · Jan 30, 2025 · e66e69b · e66e69b
1 parent d3d3543
commit e66e69b
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/R/CalculationFunctions.R b/R/CalculationFunctions.R
@@ -691,9 +691,9 @@ calculateGhoshB <- function(model,use_domestic_requirements=FALSE) {
   x <- model$x # total industry output
   A <- model$A
   if(model$specs$CommodityorIndustryType == "Commodity") {
-    B <- solve(diag(x)) %*% A %*% diag(x)
-  } else if(model$specs$CommodityorIndustryType == "Industry") {
     B <- solve(diag(q)) %*% A %*% diag(q)
+  } else if(model$specs$CommodityorIndustryType == "Industry") {
+    B <- solve(diag(x)) %*% A %*% diag(x)
   }
 
   row.names(B) <- model$Commodities$Code_Loc