# NOT RUN {
# Using the ineq data in dynamac
# Shocking Income Top 10
ardl.model <- dynardl(concern ~ incshare10 + urate, data = ineq,
lags = list("concern" = 1, "incshare10" = 1),
diffs = c("incshare10", "urate"),
lagdiffs = list("concern" = 1),
ec = TRUE, simulate = TRUE, range = 30,
shockvar = "incshare10", graph = FALSE)
area.simulation.graph(ardl.model) # Shows absolute levels
area.simulation.graph(ardl.model, changes = TRUE) # Shows changes from mean level
area.simulation.graph(ardl.model, bw = TRUE) # Grayscale plots
# }
Run the code above in your browser using DataLab