# NOT RUN {
# initialize the 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 demanded and supplied quantities
head(cbind(demanded_quantities(model, est@coef),
supplied_quantities(model, est@coef)))
# }
Run the code above in your browser using DataLab