VARselect | R Documentation |
The function returns infomation criteria and final prediction error for sequential increasing the lag order up to a VAR(p)-proccess. which are based on the same sample size.
VARselect(y, lag.max = 10, type = c("const", "trend", "both", "none"),
season = NULL, exogen = NULL)
y |
Data item containing the endogenous variables |
lag.max |
Integer for the highest lag order (default is
|
type |
Type of deterministic regressors to include. |
season |
Inlusion of centered seasonal dummy variables (integer value of frequency). |
exogen |
Inlusion of exogenous variables. |
Estimates a VAR by OLS per equation. The model is of the following form:
\bold{y}_t = A_1 \bold{y}_{t-1} + \ldots + A_p \bold{y}_{t-p} + CD_t
+ \bold{u}_t
where \bold{y}_t
is a K \times 1
vector of endogenous
variables and u_t
assigns a spherical disturbance term of the
same dimension. The coefficient matrices A_1, \ldots, A_p
are of
dimension K \times K
. In addition, either a constant and/or a
trend can be included as deterministic regressors as well as centered
seasonal dummy variables and/or exogenous variables (term CD_T
, by
setting the type
argument to the corresponding value and/or
setting season
to the desired frequency (integer) and/or providing a
matrix object for exogen
, respectively. The default for type
is
const
and for season
and exogen
the default is
set to NULL
.
Based on the same sample size the following information criteria and
the final prediction error are computed:
AIC(n) = \ln \det(\tilde{\Sigma}_u(n)) + \frac{2}{T}n K^2 \quad,
HQ(n) = \ln \det(\tilde{\Sigma}_u(n)) + \frac{2 \ln(\ln(T))}{T}n K^2 \quad,
SC(n) = \ln \det(\tilde{\Sigma}_u(n)) + \frac{\ln(T)}{T}n K^2 \quad,
FPE(n) = \left ( \frac{T + n^*}{T - n^*} \right )^K
\det(\tilde{\Sigma}_u(n)) \quad ,
with \tilde{\Sigma}_u (n) = T^{-1} \sum_{t=1}^T \bold{\hat{u}}_t
\bold{\hat{u}}_t'
and n^*
is the total number of the
parameters in each equation and n
assigns the lag order.
A list with the following elements:
selection |
Vector with the optimal lag number according to each criterium. |
criteria |
A matrix containing the values of the criteria up to
|
Bernhard Pfaff
Akaike, H. (1969), Fitting autoregressive models for prediction, Annals of the Institute of Statistical Mathematics, 21: 243-247.
Akaike, H. (1971), Autoregressive model fitting for control, Annals of the Institute of Statistical Mathematics, 23: 163-180.
Akaike, H. (1973), Information theory and an extension of the maximum likelihood principle, in B. N. Petrov and F. Csáki (eds.), 2nd International Symposium on Information Theory, Académia Kiadó, Budapest, pp. 267-281.
Akaike, H. (1974), A new look at the statistical model identification, IEEE Transactions on Automatic Control, AC-19: 716-723.
Hamilton, J. (1994), Time Series Analysis, Princeton University Press, Princeton.
Hannan, E. J. and B. G. Quinn (1979), The determination of the order of an autoregression, Journal of the Royal Statistical Society, B41: 190-195.
Lütkepohl, H. (2006), New Introduction to Multiple Time Series Analysis, Springer, New York.
Quinn, B. (1980), Order determination for a multivariate autoregression, Journal of the Royal Statistical Society, B42: 182-185.
Schwarz, G. (1978), Estimating the dimension of a model, Annals of Statistics, 6: 461-464.
VAR
data(Canada)
VARselect(Canada, lag.max = 5, type="const")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.