# NOT RUN {
IT17 <- matrix(c(
1.47, 6.47, 0.57, 2.99, 0.12 * 0.60 / (0.60 + 12.10 + 2.23 + 1.45),
0.12 * 12.10 / (0.60 + 12.10 + 2.23 + 1.45),
0.12 * 2.23 / (0.60 + 12.10 + 2.23 + 1.45),
0.12 * 1.45 / (0.60 + 12.10 + 2.23 + 1.45),
2.18, 76.32, 12.83, 43, 13.30 * 0.60 / (0.60 + 12.10 + 2.23 + 1.45),
13.30 * 12.10 / (0.60 + 12.10 + 2.23 + 1.45),
13.30 * 2.23 / (0.60 + 12.10 + 2.23 + 1.45),
13.30 * 1.45 / (0.60 + 12.10 + 2.23 + 1.45),
0.82, 19.47, 23.33, 34.88, 2.96 * 0.60 / (0.60 + 12.10 + 2.23 + 1.45),
2.96 * 12.10 / (0.60 + 12.10 + 2.23 + 1.45),
2.96 * 2.23 / (0.60 + 12.10 + 2.23 + 1.45),
2.96 * 1.45 / (0.60 + 12.10 + 2.23 + 1.45),
6.53, 13.92, 21.88, 0, 0, 0, 0, 0,
0.23, 4.05, 6.76, 0, 0, 0, 0, 0,
0, 6.43, 3.40, 0, 0, 0, 0, 0,
0.13, 8.87, 10.46, 0, 0, 0, 0, 0,
0, 0, 0, 1.45, 0, 0, 0, 0
), 8, 8, TRUE)
OT17 <- matrix(c(
11.02, 0, 0, 0, 0.60, 0, 0, 0,
0, 135.53, 0, 0, 0, 12.10, 0, 0,
0, 0, 79.23, 0, 0, 0, 2.23, 0,
0, 0, 0, 42.33, 0, 0, 0, 0,
0, 0, 0, 11.04, 0, 0, 0, 0,
0.34, 0, 0, 9.49, 0, 0, 0, 0,
0, 0, 0, 19.46, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 1.45
), 8, 8, TRUE)
rownames(IT17) <- rownames(OT17) <-
c("agri", "manu", "serv", "lab", "cap", "tax", "dividend", "bond.ROW")
colnames(IT17) <- colnames(OT17) <- c(
"sector.agri", "sector.manu", "sector.serv", "sector.hh",
"sector.FT.agri", "sector.FT.manu", "sector.FT.serv", "sector.FT.bond.ROW"
)
ge <- gemInputOutputTable_8_8(
IT = IT17,
OT = OT17
)
#### technology progress
IT.TP <- IT17
IT.TP ["lab", "sector.manu"] <-
IT.TP ["lab", "sector.manu"] * 0.9
geTP <- gemInputOutputTable_8_8(
IT = IT.TP,
OT = OT17
)
geTP$z / ge$z
geTP$p / ge$p
geTP$value.added
prop.table(geTP$value.added) - prop.table(ge$value.added)
#### capital accumulation
OT.CA <- OT17
OT.CA["cap", "sector.hh"] <- OT.CA["cap", "sector.hh"] * 1.1
geCA <- gemInputOutputTable_8_8(
IT = IT17,
OT = OT.CA
)
geCA$z / ge$z
geCA$p / ge$p
geCA$p
geCA$value.added
prop.table(geCA$value.added) - prop.table(ge$value.added)
#### tax change
OT.TC <- OT17
OT.TC["tax", "sector.agri"] <- OT.TC["tax", "sector.agri"] * 2
geTC <- gemInputOutputTable_8_8(
IT = IT17,
OT = OT.TC
)
geTC$z / ge$z
geTC$p / ge$p
##
IT.TC2 <- IT17
IT.TC2["tax", "sector.manu"] <- IT.TC2["tax", "sector.manu"] * 0.8
geTC2 <- gemInputOutputTable_8_8(
IT = IT.TC2,
OT = OT17
)
geTC2$z / ge$z
geTC2$p / ge$p
# }
Run the code above in your browser using DataLab