# \donttest{
fit <- diseq_basic(
HS | RM | ID | TREND ~
RM + TREND + W + CSHS + L1RM + L2RM + MONTH |
RM + TREND + W + L1RM + MA6DSF + MA3DHF + MONTH,
fair_houses(),
correlated_shocks = FALSE
)
# get estimated aggregate demand
aggregate_demand(fit)
# simulate the deterministic adjustment model
model <- simulate_model(
"diseq_deterministic_adjustment", list(
# observed entities, observed time points
nobs = 500, tobs = 3,
# demand coefficients
alpha_d = -0.6, beta_d0 = 9.8, beta_d = c(0.3, -0.2), eta_d = c(0.6, -0.1),
# supply coefficients
alpha_s = 0.2, beta_s0 = 4.1, beta_s = c(0.9), eta_s = c(-0.5, 0.2),
# price equation coefficients
gamma = 0.9
),
seed = 1356
)
# estimate the model object
fit <- estimate(model)
# get estimated aggregate demand
aggregate_demand(fit)
# get estimated aggregate demand
aggregate_supply(fit)
# }
Run the code above in your browser using DataLab