# NOT RUN {
# initialize the basic model using the houses dataset
model <- new(
"diseq_basic", # model type
c("ID", "TREND"), "HS", "RM", # keys, quantity, and price variables
"RM + TREND + W + CSHS + L1RM + L2RM + MONTH", # demand specification
"RM + TREND + W + L1RM + MA6DSF + MA3DHF + MONTH", # supply specification
fair_houses(), # data
correlated_shocks = FALSE # allow shocks to be correlated
)
# estimate the model object (BFGS is used by default)
est <- estimate(model)
# get estimated aggregate demand
aggregate_demand(model, est@coef)
# 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
est <- estimate(model)
# get estimated aggregate demand
aggregate_demand(model, est@coef)
# get estimated aggregate demand
aggregate_supply(model, est@coef)
# }
Run the code above in your browser using DataLab