FCVARboot | R Documentation |
FCVARboot
generates a distribution of a likelihood ratio
test statistic using a wild bootstrap, following the method of
Boswijk, Cavaliere, Rahbek, and Taylor (2016). It takes two sets
of options as inputs to estimate the model under the null and the
unrestricted model.
FCVARboot(x, k, r, optRES, optUNR, B)
x |
A matrix of variables to be included in the system. |
k |
The number of lags in the system. |
r |
The cointegrating rank. |
optRES |
An S3 object of class |
optUNR |
An S3 object of class |
B |
The number of bootstrap samples. |
A list FCVARboot_stats
containing the estimation results,
including the following parameters:
LRbs
A B x 1 vector of simulated likelihood ratio statistics
pv
An approximate p-value for the likelihood ratio statistic based on the bootstrap distribution.
H
A list containing the likelihood ratio test results.
It is identical to the output from FCVARhypoTest
, with one addition,
namely H$pvBS
which is the bootstrap p-value
mBS
The model estimates under the null hypothesis.
mUNR
The model estimates under the alternative hypothesis.
Boswijk, Cavaliere, Rahbek, and Taylor (2016) "Inference on co-integration parameters in heteroskedastic vector autoregressions," Journal of Econometrics 192, 64-85.
FCVARoptions
to set default estimation options.
FCVARestn
is called to estimate the models under the null and alternative hypotheses.
Other FCVAR postestimation functions:
FCVARhypoTest()
,
GetCharPolyRoots()
,
MVWNtest()
,
plot.FCVAR_roots()
,
summary.FCVAR_roots()
,
summary.MVWN_stats()
opt <- FCVARoptions() opt$gridSearch <- 0 # Disable grid search in optimization. opt$dbMin <- c(0.01, 0.01) # Set lower bound for d,b. opt$dbMax <- c(2.00, 2.00) # Set upper bound for d,b. opt$constrained <- 0 # Impose restriction dbMax >= d >= b >= dbMin ? 1 <- yes, 0 <- no. x <- votingJNP2014[, c("lib", "ir_can", "un_can")] opt$plotRoots <- 0 optUNR <- opt optRES <- opt optRES$R_Beta <- matrix(c(1, 0, 0), nrow = 1, ncol = 3) set.seed(42) FCVARboot_stats <- FCVARboot(x, k = 2, r = 1, optRES, optUNR, B = 2) # In practice, set the number of bootstraps so that (B+1)*alpha is an integer, # where alpha is the chosen level of significance. # For example, set B = 999 (but it takes a long time to compute).
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.