dynardl.all.plots {dynamac}

Combine all of the potential plots of a simulated response in a dynardl model

Description

Combine all of the potential plots of a simulated response in a dynardl model

Usage

dynardl.all.plots(
  x,
  type = "area",
  bw = FALSE,
  last.period = NULL,
  start.period = 1,
  tol = (abs(x$model$ymean) * 0.01),
  abs.errors = "none",
  ylim = NULL,
  xlab = NULL,
  ylab = NULL,
  ...
)

Arguments

x

a dynardl model with a simulation to be plotted. Since all plots includes absolute cumulative differences, fullsims must be TRUE in the dynardl simulation

type

whether the plot should be an area plot (area) or a spike plot (spike)

bw

should the colors be in black and white (for publication)? The default is FALSE

last.period

when deciding when to stop calculating the absolute value of the shocks to the dependent variable, you can specify a specific period in which to stop calculating absolute cumulative differences. Specify a tol or a last.period. If both are specified, last.period overrides tol

start.period

which period of the simulation to begin the plot with. You can view the equilibriating behavior of the dependent variable, or you can skip forward in time (maybe to just before the shock). The default is 1 (the first period of the simulation)

tol

when deciding when to stop calculating the absolute value of the shocks to the dependent variable, you can specify the minimum amount of movement required to qualify as a non-noise change over time periods (for calculating absolute cumulative differences). The default is 0.1 percent of the mean of the dependent variable. Specify a tol or a last.period. If both are specified, last.period overrides tol

abs.errors

when calculating confidence for the absolute cumulative effect, should differences accumulate in each time time period (cumulate, which could be explosive if the error in the model is large), should differences be observed at each time (within.period, which will have smaller values in equilibrium than when changing), or should only the values be plotted (none)

ylim

a user-defined y-limit to be used instead of the default (for instance, for shared axes. Use caution, as it will be passed to all plots)

xlab

a user-defined x-label to be used instead of the default (use caution, as it will be passed to all plots)

ylab

a user-defined y-label to be used instead of the default (use caution, as it will be passed to all plots)

...

other arguments to be passed to the call to plot. Use caution, as they will be passed to all plots

Details

When running dynardl, simulate must be TRUE so that there is a simulation to plot. Also, fullsims must be TRUE as the plot will contain absolute cumulative differences. See dynardl.simulation.plot for arguments to the individual plotting types

Value

a 2 x 3 grid of the plots of the simulated dynardl model effects plots

Author(s)

Soren Jordan and Andrew Q. Philips

Examples

# Using the ineq data in dynamac
# Shocking Income Top 10
# Not run: simulations are time-intensive to estimate as an example
## Not run: 
set.seed(1)
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", fullsims = TRUE)

# Shows all of the potential responses
dynardl.all.plots(ardl.model)	
# Same plot, but with spikeplot
dynardl.all.plots(ardl.model, type = "spike")  
# Grayscale plots
dynardl.all.plots(ardl.model, bw = TRUE)	 

## End(Not run)