# NOT RUN {
es.agri <- 0.2 # the elasticity of substitution
es.manu <- 0.3
es.serv <- 0.1
es.VA.agri <- 0.25
es.VA.manu <- 0.5
es.VA.serv <- 0.8
d.agri <- c(260, 345, 400, 200, 160)
d.manu <- c(320, 390, 365, 250, 400)
d.serv <- c(150, 390, 320, 400, 210)
d.hh <- c(635, 600, 385, 0, 0)
# d.hh <- c(635, 600, 100, 0, 0)
D <- cbind(d.agri, d.manu, d.serv, d.hh)
dst.agri <- Node$new("sector.agri",
type = "SCES", sigma = 1 - 1 / es.agri,
alpha = 1,
beta = prop.table(c(
sum(d.agri[1:3]),
sum(d.agri[4:5])
))
)
dst.agri$AddChild("cc1.agri",
type = "Leontief",
a = prop.table(d.agri[1:3])
)$
AddChild("agri")$AddSibling("manu")$AddSibling("serv")$
parent$
AddSibling("cc2.agri",
type = "SCES", sigma = 1 - 1 / es.VA.agri,
alpha = 1,
beta = prop.table(d.agri[4:5])
)$
AddChild("lab")$AddSibling("cap")
##
dst.manu <- Node$new("sector.manu",
type = "SCES", sigma = 1 - 1 / es.manu,
alpha = 1,
beta = prop.table(c(
sum(d.manu[1:3]),
sum(d.manu[4:5])
))
)
dst.manu$AddChild("cc1.manu",
type = "Leontief",
a = prop.table(d.manu[1:3])
)$
AddChild("agri")$AddSibling("manu")$AddSibling("serv")$
parent$
AddSibling("cc2.manu",
type = "SCES", sigma = 1 - 1 / es.VA.manu,
alpha = 1,
beta = prop.table(d.manu[4:5])
)$
AddChild("lab")$AddSibling("cap")
##
dst.serv <- Node$new("sector.serv",
type = "SCES", sigma = 1 - 1 / es.serv,
alpha = 1,
beta = prop.table(c(
sum(d.serv[1:3]),
sum(d.serv[4:5])
))
)
dst.serv$AddChild("cc1.serv",
type = "Leontief",
a = prop.table(d.serv[1:3])
)$
AddChild("agri")$AddSibling("manu")$AddSibling("serv")$
parent$
AddSibling("cc2.serv",
type = "SCES", sigma = 1 - 1 / es.VA.serv,
alpha = 1,
beta = prop.table(d.serv[4:5])
)$
AddChild("lab")$AddSibling("cap")
##
dst.hh <- Node$new("sector.hh",
type = "SCES", sigma = -1,
alpha = 1,
beta = prop.table(d.hh[1:3])
)
dst.hh$AddChild("agri")$AddSibling("manu")$AddSibling("serv")
dstl <- list(
dst.agri,
dst.manu,
dst.serv,
dst.hh
)
ge <- gemInputOutputTable_5_4(dstl)
#### labor supply increase
geLSI <- gemInputOutputTable_5_4(dstl, supply.labor = 850 * 1.08)
geLSI$p
geLSI$z / ge$z
# }
Run the code above in your browser using DataLab