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

Adding value added

Browse Source
This commit is contained in:
alex 2026-08-14 16:33:12 -06:00
parent 7b95f4af1a
commit a7f30832d1
1 changed files with 28 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
Leontief.r
View File

@ -18,9 +18,37 @@ 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)
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
###########Value Added
VALUE_ADD_MAT <- colSums(I_MAT[c('EmpComp','PropInc','OthPropInc','TOPI','SUB'),])/colSums(I_MAT)
VALUE_ADD_MAT
colSums(VALUE_ADD_MAT[colnames(TYPE2_A)]*MULTIPLIER_OUTPUT_TYPE2 )
colSums(VALUE_ADD_MAT[colnames(TYPE1_A)]*MULTIPLIER_OUTPUT_TYPE1 )
colSums(MULTIPLIER_OUTPUT_TYPE1*VALUE_ADD_MAT[colnames(TYPE1_A)] )
VA_TYPE2 <- colSums(solve(diag(nrow(TYPE2_A))-TYPE2_A*VALUE_ADD_MAT[colnames(TYPE2_A)]))
VA_TYPE1 <- colSums(solve(diag(nrow(TYPE1_A))-TYPE1_A*VALUE_ADD_MAT[colnames(TYPE1_A)]))
EmpComp_MAT <- I_MAT[c('EmpComp'),]/colSums(I_MAT)
EmpComp_TYPE2 <- colSums(solve(diag(nrow(TYPE2_A))-TYPE2_A*EmpComp_MAT[colnames(TYPE2_A)]))
EmpComp_TYPE2
EmpComp_TYPE1 <- colSums(solve(diag(nrow(TYPE1_A))-TYPE1_A*EmpComp_MAT[colnames(TYPE1_A)]))
min(EmpComp_TYPE1 )
PropInc_MAT <- I_MAT[c('PropInc'),]/colSums(I_MAT)
OthPropInc_MAT <- I_MAT[c('OthPropInc'),]/colSums(I_MAT)
TOPI_MAT <- I_MAT[c('TOPI'),]/colSums(I_MAT)
SUB_MAT <- I_MAT[c('SUB'),]/colSums(I_MAT)