# NOT RUN {
# estimate a model using the houses dataset
fit <- diseq_deterministic_adjustment(
HS | RM | ID | TREND ~
RM + TREND + W + CSHS + L1RM + L2RM + MONTH |
RM + TREND + W + L1RM + MA6DSF + MA3DHF + MONTH,
fair_houses(), correlated_shocks = FALSE,
estimation_options = list(control = list(maxit = 1e+5)))
# get estimated normalized shortages
head(normalized_shortages(fit))
# get estimated relative shortages
head(relative_shortages(fit))
# get the estimated shortage probabilities
head(shortage_probabilities(fit))
# get the estimated shortage indicators
head(shortage_indicators(fit))
# get the estimated shortages
head(shortages(fit))
# get the estimated shortage variance
shortage_standard_deviation(fit)
# }
Run the code above in your browser using DataLab