1. Introduction
In order to analyze the structure and/or forecast realizations of a stochastic process, or an observation within the time series, the models currently developed by econometricians have a certain degree of complexity. Thus, econometric estimation and the analysis of stochastic processes to explain economic phenomena through models are increasingly relevant, mainly in the context of time series, taking into consideration that various processes focus on understanding the dynamic structure of the series and on the possibility of forecasting its dynamic pattern of temporal behavior or the extrapolation of a stochastic process [
1,
2,
3,
4] where the lags of the variables involved play a key role in terms of the autoregressive models, as is the case in the estimation of autoregressive vector models (VARs), the central theme in this work.
These models are expressed through differential equations, since each variable is explained by the lags of both itself and the remaining variables.
It should be noted that each of the variables involved must meet the assumption of stationarity as a particular state of statistical equilibrium, where their probability distributions remain stable over time [
5,
6]. This implies that once the system is interrupted by some type of shocks, it will adjust back to equilibrium [
7] or the shocks gradually disappear.
In estimating these models, it is accepted and often required that, if the estimators meet the tests, then they are the best linearly unbiased estimators (BLUEs). But what happens if the tests applied to the model are not fulfilled? Are the estimators not valid for the analysis? Does the model have to be scrapped?
2. Empirical Obtaining of Estimators in a VAR Model
It is important to point out that when there are large samples, the assumptions of normality, homoscedasticity, and the absence of autocorrelation in the errors are hardly fulfilled. This occurs regularly when using short duration data, for example monthly, quarterly, as well as long periods in the analysis. This could mean a limitation of the model that leads to strong criticism in this regard; however, as Wooldridge mentioned in his modern approach, given the law of large numbers, an asymptotic normality is assumed, due to the size of the sample and in terms of the homoscedasticity and absence of autocorrelation in the errors, the results obtained allow themselves to be the best linearly unbiased estimators, pointed out by Guarati and Porter [
8].
In this regard, it has agreed in common that for the forecasting purposes, the VAR models are required to fulfill the assumptions of the estimation: normality, homoscedasticity, and the absence of autocorrelation in the errors. However, if the estimated model is used only to analyze the structural changes of the economic variables, the requirement could be relaxed to the stability of the solutions; that is, the convergence in terms of the dynamic analysis, which can be determined by estimating the inverse roots of the characteristic polynomial of the autoregressive vector.
It should be noted that, since these are unrestricted models, the main advantage is that there will be no specification errors in the empirical estimation, in addition to the fact that the long-term cointegration solution is exempt from the problem of spuriousness or meaningless regressions, as it is defined by Granger and Newbold [
9], with the initial idea owed to Yule [
10].
Therefore, the VAR models not only provide a better estimate of forecasts compared to static ones, but also could be analyzed in a dynamic and structural manner where the importance of a shock of one variable on the others is revealed, with the relaxation of the related assumptions.
To apply the structural analysis of systems of simultaneous differential equations derived from the VAR approach, the series corresponding to Mexican exports to the US market and imports of Chinese origin are used with an annual periodicity from 2001 to 2020 (
Figure 1).
It should be noted that the study Mexican trade with its two most important trading partners during the last 20 years is of great importance given the two facts. Firstly, China’s entry into the World Trade Organization, an unprecedented event that contributed to the expansion of the Chinese products around the world markets, and eventually has become the second largest trading partner for Mexico since 2003. Secondly, at the beginning of 2020 COVID-19 pandemic has generated an important structural change in the trading among the three countries.
It is evident that both series are non-stationary since they have a trend. In order to estimate a VAR model, it is necessary for the series to be integrated in the same order to be transformed into stationariy; in this particular case, they are I(1).
Figure 2 shows the stationarity of the two series after the respective first difference, which was confirmed by augmented Dickey–Fuller [
11] and Phillips–Perron [
12] tests.
Subsequently, the VAR models were estimated with one and two lags respectively, estimations that accredited the stability tests within the unit circle and the residuals model’s presented normality, homoscedasticity and an absence of autocorrelation.
2.1. Solution of the VAR System of Equations with One Lag
The equation derived from the normalized equation estimated in the Eviews 12 Student Version Lite of S&P Global, New York, NY, USA with a model that shows stability within the unit circle (inverse roots of the autoregressive characteristic polynomial) is as follows:
To simplify the notation, suppose that
and
, so that
Particular solutions that can be supposed for the two variables are
and
, which implies that
and
; then,
Reducing yields are the following:
Solving the system yields
and
, while normalizing yields the following:
Note that 1.8 is very close to the solution provided by the software (1.76).
To obtain the complementary solutions, the following can been carried out:
We assume that and where γ and δ are constants; therefore, and .
Substituting into the difference equations that are now homogeneous,
Multiply everything by
, arriving at the following:
Solving the equation, the values of are and .
The complementary solutions are as follows:
Consequently, the general solutions are as follows:
2.2. Solution of the VAR System of Equations with Two Lags
Regrouping yields the following:
Solving the system results in and .
Normalizing then results in the following:
In this case, the solution given by the software is 1.67; therefore, it is consistent.
Following this reasoning, we could generalize as follows:
Let be a VAR model of
variables with
lags:
We express the VAR model in difference equations:
Supposing that , and , then , , and .
Analogously, if
,
and
then,
Factoring the constants yields the following:
In matrix form, this is represented as
where
,
.
This results in a system of simultaneous linear equations with the number of variables equal to the number of equations, i.e., it will always be a square matrix.
The solution of the system constitutes the vector k of the long-term cointegration equation.
To achieve the above, that is, the general solutions converging to their particular solutions, the complementary solutions would have to be dynamically stable, or they would converge to nullity.
represent the number of variables.
Since all the solutions of the characteristic equation are real numbers, .
To fulfill the above, it is required that
, and the particular solutions will be as follows:
On the other hand, to achieve , no doubt the necessary and sufficient condition is that the series involved in the system of difference equations has to be stationary, which is consistent with a single time series.