| general.calculator | R Documentation |
Generate the generalized effect formulae for an autoregressive distributed lag (ADL) model, given pulse effects and shock/treatment history
general.calculator(d.x, d.y, h, limit, pulses)
d.x |
the order of differencing of the x variable in the ADL model. (Generally, this is the same x variable used in |
d.y |
the order of differencing of the y variable in the ADL model. (Generally, this is the same y variable used in |
h |
an integer for the shock/treatment history. |
limit |
an integer for the number of periods (s) to determine the generalized effect (beginning at 0) |
pulses |
a list of pulse effect formulae used to construct the generalized effect formulae. We expect this will be provided by |
general.calculator does no calculation. It generates a list of mpoly formulae that contain variable names that represent the generalized effect in each period. The expectation is that these will be evaluated using coefficients from an object containing an ADL model with corresponding variables. Note: mpoly does not allow variable names with a .; variables passed to general.calculator should not include this character
a list of limit + 1 mpoly formulae containing the generalized effect formula in each period
Soren Jordan, Garrett N. Vande Kamp, and Reshi Rajan
# ADL(1,1)
x.lags <- c("x" = 0, "l_1_x" = 1) # lags of x
y.lags <- c("l_1_y" = 1)
s <- 5
pulse.effects <- pulse.calculator(x.vrbl = x.lags, y.vrbl = y.lags, limit = s)
# Assume that both x and y are in levels and we want a pulse treatment
general.pulse.effects <- general.calculator(d.x = 0, d.y = 0,
h = -1, limit = s, pulses = pulse.effects)
general.pulse.effects
# Apply a step treatment
general.step.effects <- general.calculator(d.x = 0, d.y = 0,
h = 0, limit = s, pulses = pulse.effects)
general.step.effects
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.