From f698c5e51ae349dfc3d775a4b3ff2b4930d175af Mon Sep 17 00:00:00 2001 From: Alex Gebben Work Date: Mon, 17 Aug 2026 17:07:19 -0600 Subject: [PATCH] Changing flow of code --- Leontief.r | 26 ++++++++++++++----- .../.~lock.Wyoming_ixi_39_naics2_2024.csv# | 1 - .../.~lock.Wyoming_leontief_naics2_2024.xlsx# | 2 +- 3 files changed, 20 insertions(+), 9 deletions(-) delete mode 100644 Wyoming_2024_SAM_v1.7/.~lock.Wyoming_ixi_39_naics2_2024.csv# diff --git a/Leontief.r b/Leontief.r index 67c3465..543d90f 100644 --- a/Leontief.r +++ b/Leontief.r @@ -11,27 +11,39 @@ TOTALS <- ifelse(TOTALS==0,0.000000,TOTALS) %>% as.numeric #Column division A_MAT <- t(t(I_MAT)/TOTALS) ACCOUNT_RANGE <- (LEN-13):LEN -TYPE1_A <- A_MAT[-ACCOUNT_RANGE,-ACCOUNT_RANGE] -MULTIPLIER_OUTPUT_TYPE1 <- solve(diag(LEN-14)-TYPE1_A) +KEEP_TYPE1 <- !(colnames(A_MAT) %in% c('EmpComp','PropInc','OthPropInc','TOPI','SUB','HH','FICA','K','FED','SG','LG','INV','RUSA','ROW')) +A_MAT1 <- A_MAT[KEEP_TYPE1,KEEP_TYPE1] + +A_MAT1[1,1] +MULTIPLIER_OUTPUT_TYPE1 <- solve(diag(nrow(A_MAT1))-A_MAT1) TYPE1_MULTIPLIERS_TOTAL <- colSums(MULTIPLIER_OUTPUT_TYPE1 ) -DEMAND_RANGE <- ACCOUNT_RANGE[-1:-6] -TYPE2_A <- A_MAT[-DEMAND_RANGE,-DEMAND_RANGE] -MULTIPLIER_OUTPUT_TYPE2 <- solve(diag(LEN-8)-TYPE2_A) +KEEP_TYPE2 <- !(colnames(A_MAT) %in% c('FICA','K','FED','SG','LG','INV','RUSA','ROW')) +A_MAT2 <- A_MAT[KEEP_TYPE2,KEEP_TYPE2] + + +MULTIPLIER_OUTPUT_TYPE2 <- solve(diag(nrow(A_MAT2))-A_MAT2) MULTIPLIER_OUTPUT_TYPE2 HH_VALUE_ACCOUNTS <- (nrow(MULTIPLIER_OUTPUT_TYPE2 )-5):nrow(MULTIPLIER_OUTPUT_TYPE2 ) TYPE2_MULTIPLIERS_TOTAL <- colSums(MULTIPLIER_OUTPUT_TYPE2[-HH_VALUE_ACCOUNTS ,]) TYPE2_MULTIPLIERS_TOTAL - +names(TYPE2_MULTIPLIERS_TOTAL ) +INDUCED <- TYPE2_MULTIPLIERS_TOTAL[!(names(TYPE2_MULTIPLIERS_TOTAL) %in% c('EmpComp','PropInc','OthPropInc','TOPI','SUB','HH'))]-TYPE1_MULTIPLIERS_TOTAL +INDIRECT <- TYPE1_MULTIPLIERS_TOTAL-1 +INDIRECT +INDUCED +TYPE1_MULTIPLIERS_TOTAL ###########Value Added VALUE_ADD_MAT <- colSums(I_MAT[c('EmpComp','PropInc','OthPropInc','TOPI','SUB'),])/colSums(I_MAT) +VALUE_ADD_MAT (I_MAT[c('EmpComp'),])/colSums(I_MAT) (I_MAT[c('PropInc'),])/colSums(I_MAT) (I_MAT[c('OthPropInc'),])/colSums(I_MAT) (I_MAT[c('TOPI'),])/colSums(I_MAT) (I_MAT[c('SUB'),])/colSums(I_MAT) - +nrow(MULTIPLIER_OUTPUT_TYPE2) +MULTIPLIER_OUTPUT_TYPE2 length( TYPE2_MULTIPLIERS_TOTAL) length( I_MAT[c('EmpComp'),]) /colSums(I_MAT))) diff --git a/Wyoming_2024_SAM_v1.7/.~lock.Wyoming_ixi_39_naics2_2024.csv# b/Wyoming_2024_SAM_v1.7/.~lock.Wyoming_ixi_39_naics2_2024.csv# deleted file mode 100644 index f685ce0..0000000 --- a/Wyoming_2024_SAM_v1.7/.~lock.Wyoming_ixi_39_naics2_2024.csv# +++ /dev/null @@ -1 +0,0 @@ -,alex,pop-os,14.08.2026 10:07,file:///home/alex/.config/libreoffice/4; \ No newline at end of file diff --git a/Wyoming_2024_SAM_v1.7/.~lock.Wyoming_leontief_naics2_2024.xlsx# b/Wyoming_2024_SAM_v1.7/.~lock.Wyoming_leontief_naics2_2024.xlsx# index f997bbf..c39a321 100644 --- a/Wyoming_2024_SAM_v1.7/.~lock.Wyoming_leontief_naics2_2024.xlsx# +++ b/Wyoming_2024_SAM_v1.7/.~lock.Wyoming_leontief_naics2_2024.xlsx# @@ -1 +1 @@ -,alex,pop-os,14.08.2026 14:19,file:///home/alex/.config/libreoffice/4; \ No newline at end of file +,alex,pop-os,17.08.2026 17:02,file:///home/alex/.config/libreoffice/4; \ No newline at end of file