alex/Tapestry_analysis
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Changing flow of code

Browse Source
This commit is contained in:
Alex Gebben Work 2026-08-17 17:07:19 -06:00
parent 4edb3480be
commit f698c5e51a
3 changed files with 20 additions and 9 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
Leontief.r
View File

@ -11,27 +11,39 @@ TOTALS <- ifelse(TOTALS==0,0.000000,TOTALS) %>% as.numeric
#Column division #Column division
A_MAT <- t(t(I_MAT)/TOTALS) A_MAT <- t(t(I_MAT)/TOTALS)
ACCOUNT_RANGE <- (LEN-13):LEN ACCOUNT_RANGE <- (LEN-13):LEN
TYPE1_A <- A_MAT[-ACCOUNT_RANGE,-ACCOUNT_RANGE] KEEP_TYPE1 <- !(colnames(A_MAT) %in% c('EmpComp','PropInc','OthPropInc','TOPI','SUB','HH','FICA','K','FED','SG','LG','INV','RUSA','ROW'))
MULTIPLIER_OUTPUT_TYPE1 <- solve(diag(LEN-14)-TYPE1_A) 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 ) TYPE1_MULTIPLIERS_TOTAL <- colSums(MULTIPLIER_OUTPUT_TYPE1 )
DEMAND_RANGE <- ACCOUNT_RANGE[-1:-6] KEEP_TYPE2 <- !(colnames(A_MAT) %in% c('FICA','K','FED','SG','LG','INV','RUSA','ROW'))
TYPE2_A <- A_MAT[-DEMAND_RANGE,-DEMAND_RANGE] A_MAT2 <- A_MAT[KEEP_TYPE2,KEEP_TYPE2]
MULTIPLIER_OUTPUT_TYPE2 <- solve(diag(LEN-8)-TYPE2_A)
MULTIPLIER_OUTPUT_TYPE2 <- solve(diag(nrow(A_MAT2))-A_MAT2)
MULTIPLIER_OUTPUT_TYPE2 MULTIPLIER_OUTPUT_TYPE2
HH_VALUE_ACCOUNTS <- (nrow(MULTIPLIER_OUTPUT_TYPE2 )-5):nrow(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 <- colSums(MULTIPLIER_OUTPUT_TYPE2[-HH_VALUE_ACCOUNTS ,])
TYPE2_MULTIPLIERS_TOTAL 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 Added
VALUE_ADD_MAT <- colSums(I_MAT[c('EmpComp','PropInc','OthPropInc','TOPI','SUB'),])/colSums(I_MAT) 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('EmpComp'),])/colSums(I_MAT)
(I_MAT[c('PropInc'),])/colSums(I_MAT) (I_MAT[c('PropInc'),])/colSums(I_MAT)
(I_MAT[c('OthPropInc'),])/colSums(I_MAT) (I_MAT[c('OthPropInc'),])/colSums(I_MAT)
(I_MAT[c('TOPI'),])/colSums(I_MAT) (I_MAT[c('TOPI'),])/colSums(I_MAT)
(I_MAT[c('SUB'),])/colSums(I_MAT) (I_MAT[c('SUB'),])/colSums(I_MAT)
nrow(MULTIPLIER_OUTPUT_TYPE2)
MULTIPLIER_OUTPUT_TYPE2
length( TYPE2_MULTIPLIERS_TOTAL) length( TYPE2_MULTIPLIERS_TOTAL)
length( I_MAT[c('EmpComp'),]) length( I_MAT[c('EmpComp'),])
/colSums(I_MAT))) /colSums(I_MAT)))

1
Wyoming_2024_SAM_v1.7/.~lock.Wyoming_ixi_39_naics2_2024.csv#
View File

@ -1 +0,0 @@
,alex,pop-os,14.08.2026 10:07,file:///home/alex/.config/libreoffice/4;

2
Wyoming_2024_SAM_v1.7/.~lock.Wyoming_leontief_naics2_2024.xlsx#
View File

@ -1 +1 @@
,alex,pop-os,14.08.2026 14:19,file:///home/alex/.config/libreoffice/4; ,alex,pop-os,17.08.2026 17:02,file:///home/alex/.config/libreoffice/4;