84 lines
2.8 KiB
R
84 lines
2.8 KiB
R
library(tidyverse)
|
|
I_MAT <- read.csv("Wyoming_2024_SAM_v1.7/Wyoming_ixi_39_naics2_2024.csv")
|
|
NAMES <- I_MAT[,1]
|
|
LEN <- nrow(I_MAT)
|
|
I_MAT <- I_MAT[,-1]
|
|
I_MAT <- as.matrix(I_MAT,nrow=LEN)
|
|
colnames(I_MAT) <- NAMES
|
|
rownames(I_MAT) <- NAMES
|
|
TOTALS <- colSums(I_MAT )
|
|
TOTALS <- ifelse(TOTALS==0,0.000000,TOTALS) %>% as.numeric
|
|
#Column division
|
|
A_MAT <- t(t(I_MAT)/TOTALS)
|
|
ACCOUNT_RANGE <- (LEN-13):LEN
|
|
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 )
|
|
|
|
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)))
|
|
|
|
|
|
TYPE1_MULTIPLIERS_TOTAL
|
|
|
|
|
|
colSums(I_MAT[(colnames(I_MAT) %in% 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)
|
|
|
|
|
|
|
|
|