# \donttest{
zeta <- 1.25
dst.firm1 <- node_new(
"prod2",
type = "CD", alpha = 2,
beta = c(0.5, 0.5),
"lab1", "prod1"
)
dst.firm2.1 <- node_new(
"prod3.1",
type = "CD", alpha = 2,
beta = c(0.5, 0.5),
"lab2.1", "prod2.1"
)
dst.firm2.2 <- node_new(
"prod3.2",
type = "CD", alpha = 1,
beta = c(0.5, 0.5),
"lab2.2", "prod2.2"
)
dst.firm3 <- node_new(
"prod1",
type = "Leontief",
a = c(1, 1)*zeta,
"prod3.1", "prod3.2"
)
dst.consumer <- node_new(
"util",
type = "CD", alpha = 1,
beta = c(1/3, 1/3, 1/3),
"prod1","prod2.1","prod2.2"
)
ge <- sdm2(
A = c(dst.firm1, dst.firm2.1,dst.firm2.2,
dst.firm3,dst.consumer),
B = matrix(c(0, 0, 0, 1,0,
1, 0, 0, 0,0,
1, 0, 0, 0,0,
0, 1, 0, 0,0,
0, 0, 1, 0,0,
0, 0, 0, 0,0,
0, 0, 0, 0,0,
0, 0, 0, 0,0),8,5,TRUE),
S0Exg = matrix(c(NA, NA, NA, NA,NA,
NA, NA, NA, NA,NA,
NA, NA, NA, NA,NA,
NA, NA, NA, NA,NA,
NA, NA, NA, NA,NA,
NA, NA, NA, NA,100,
NA, NA, NA, NA,100,
NA, NA, NA, NA,100),8,5,TRUE),
names.commodity = c("prod1", "prod2.1","prod2.2",
"prod3.1","prod3.2",
"lab1", "lab2.1","lab2.2"),
names.agent = c("firm1", "firm2.1","firm2.2",
"firm3.1",
"consumer"),
numeraire = "lab1",
policy = makePolicyMeanValue(30),
ts = TRUE
)
ge$p
ge$z
ge$D
ge$S
ge$DV
ge$SV
# }
Run the code above in your browser using DataLab