VARMAX Procedure

Vector Autoregressive Model with Exogenous Variables

A VAR process can be affected by other observable variables that are determined outside the system of interest. Such variables are called exogenous (independent) variables. Exogenous variables can be stochastic or nonstochastic. The process can also be affected by the lags of exogenous variables. A model used to describe this process is called a VARX(p,s) model.

The VARX(p,s) model is written as

StartLayout 1st Row  bold y Subscript t Baseline equals bold-italic delta plus sigma-summation Underscript i equals 1 Overscript p Endscripts normal upper Phi Subscript i Baseline bold y Subscript t minus i Baseline plus sigma-summation Underscript i equals 0 Overscript s Endscripts normal upper Theta Subscript i Superscript asterisk Baseline bold x Subscript t minus i Baseline plus bold-italic epsilon Subscript t EndLayout

where bold x Subscript t Baseline equals left-parenthesis x Subscript 1 t Baseline comma ellipsis comma x Subscript r t Baseline right-parenthesis prime is an r-dimensional time series vector and normal upper Theta Subscript i Superscript asterisk is a k times r matrix.

For example, a VARX(1,0) model is

StartLayout 1st Row  bold y Subscript t Baseline equals bold-italic delta plus normal upper Phi 1 bold y Subscript t minus 1 Baseline plus normal upper Theta 0 Superscript asterisk Baseline bold x Subscript t Baseline plus bold-italic epsilon Subscript t EndLayout

where bold y Subscript t Baseline equals left-parenthesis y Subscript 1 t Baseline comma y Subscript 2 t Baseline comma y Subscript 3 t Baseline right-parenthesis prime and bold x Subscript t Baseline equals left-parenthesis x Subscript 1 t Baseline comma x Subscript 2 t Baseline right-parenthesis prime.

The following statements fit the VARX(1,0) model to the given data:

data grunfeld;
   input year y1 y2 y3 x1 x2 x3;
   label y1='Gross Investment GE'
         y2='Capital Stock Lagged GE'
         y3='Value of Outstanding Shares GE Lagged'
         x1='Gross Investment W'
         x2='Capital Stock Lagged W'
         x3='Value of Outstanding Shares Lagged W';
datalines;
1935  33.1 1170.6  97.8 12.93  191.5   1.8
1936  45.0 2015.8 104.4 25.90  516.0    .8
1937  77.2 2803.3 118.0 35.05  729.0   7.4
1938  44.6 2039.7 156.2 22.89  560.4  18.1

   ... more lines ...   

/*--- Vector Autoregressive Process with Exogenous Variables ---*/

proc varmax data=grunfeld;
   model y1-y3 = x1 x2 / p=1 lagmax=5
                         printform=univariate
                         print=(impulsx=(all) estimates);
run;

The VARMAX procedure output is shown in Figure 23 through Figure 25.

Figure 23 shows the descriptive statistics for the dependent (endogenous) and independent (exogenous) variables with labels.

Figure 23: Descriptive Statistics for the VARX(1, 0) Model

The VARMAX Procedure

Number of Observations 20
Number of Pairwise Missing 0

Simple Summary Statistics
Variable Type N Mean Standard
Deviation
Min Max Label
y1 Dependent 20 102.29000 48.58450 33.10000 189.60000 Gross Investment GE
y2 Dependent 20 1941.32500 413.84329 1170.60000 2803.30000 Capital Stock Lagged GE
y3 Dependent 20 400.16000 250.61885 97.80000 888.90000 Value of Outstanding Shares GE Lagged
x1 Independent 20 42.89150 19.11019 12.93000 90.08000 Gross Investment W
x2 Independent 20 670.91000 222.39193 191.50000 1193.50000 Capital Stock Lagged W


Figure 24 shows the parameter estimates for the constant, the lag zero coefficients of exogenous variables, and the lag one AR coefficients. From the schematic representation of parameter estimates, the significance of the parameter estimates can be easily verified. The symbol "C" means the constant and "XL0" means the lag zero coefficients of exogenous variables.

Figure 24: Parameter Estimates for the VARX(1, 0) Model

The VARMAX Procedure

Type of Model VARX(1,0)
Estimation Method Least Squares Estimation

Constant
Variable Constant
y1 -12.01279
y2 702.08673
y3 -22.42110

XLag
Lag Variable x1 x2
0 y1 1.69281 -0.00859
  y2 -6.09850 2.57980
  y3 -0.02317 -0.01274

AR
Lag Variable y1 y2 y3
1 y1 0.23699 0.00763 0.02941
  y2 -2.46656 0.16379 -0.84090
  y3 0.95116 0.00224 0.93801

Schematic Representation
Variable/Lag C XL0 AR1
y1 . +. ...
y2 + .+ ...
y3 - .. +.+
+ is > 2*std error,  - is < -2*std error,  . is between,  * is N/A


Figure 25 shows the parameter estimates and their significance.

Figure 25: Parameter Estimates for the VARX(1, 0) Model Continued

Model Parameter Estimates
Equation Parameter Estimate Standard
Error
t Value Pr > |t| Variable
y1 CONST1 -12.01279 27.47108 -0.44 0.6691 1
  XL0_1_1 1.69281 0.54395 3.11 0.0083 x1(t)
  XL0_1_2 -0.00859 0.05361 -0.16 0.8752 x2(t)
  AR1_1_1 0.23699 0.20668 1.15 0.2722 y1(t-1)
  AR1_1_2 0.00763 0.01627 0.47 0.6470 y2(t-1)
  AR1_1_3 0.02941 0.04852 0.61 0.5548 y3(t-1)
y2 CONST2 702.08673 256.48046 2.74 0.0169 1
  XL0_2_1 -6.09850 5.07849 -1.20 0.2512 x1(t)
  XL0_2_2 2.57980 0.50056 5.15 0.0002 x2(t)
  AR1_2_1 -2.46656 1.92967 -1.28 0.2235 y1(t-1)
  AR1_2_2 0.16379 0.15193 1.08 0.3006 y2(t-1)
  AR1_2_3 -0.84090 0.45304 -1.86 0.0862 y3(t-1)
y3 CONST3 -22.42110 10.31166 -2.17 0.0487 1
  XL0_3_1 -0.02317 0.20418 -0.11 0.9114 x1(t)
  XL0_3_2 -0.01274 0.02012 -0.63 0.5377 x2(t)
  AR1_3_1 0.95116 0.07758 12.26 0.0001 y1(t-1)
  AR1_3_2 0.00224 0.00611 0.37 0.7201 y2(t-1)
  AR1_3_3 0.93801 0.01821 51.50 0.0001 y3(t-1)


The fitted model is given as

StartLayout 1st Row 1st Column Start 3 By 1 Matrix 1st Row  y Subscript 1 t 2nd Row  y Subscript 2 t 3rd Row  y Subscript 3 t EndMatrix 2nd Column equals 3rd Column Start 6 By 1 Matrix 1st Row  negative 12.013 2nd Row  left-parenthesis 27.471 right-parenthesis 3rd Row  702.086 4th Row  left-parenthesis 256.480 right-parenthesis 5th Row  negative 22.421 6th Row  left-parenthesis 10.312 right-parenthesis EndMatrix plus Start 6 By 2 Matrix 1st Row 1st Column 1.693 2nd Column negative 0.009 2nd Row 1st Column left-parenthesis 0.544 right-parenthesis 2nd Column left-parenthesis 0.054 right-parenthesis 3rd Row 1st Column negative 6.099 2nd Column 2.580 4th Row 1st Column left-parenthesis 5.078 right-parenthesis 2nd Column left-parenthesis 0.501 right-parenthesis 5th Row 1st Column negative 0.023 2nd Column negative 0.013 6th Row 1st Column left-parenthesis 0.204 right-parenthesis 2nd Column left-parenthesis 0.020 right-parenthesis EndMatrix StartBinomialOrMatrix x Subscript 1 t Choose x Subscript 2 t EndBinomialOrMatrix 2nd Row 1st Column Blank 2nd Column plus 3rd Column Start 6 By 3 Matrix 1st Row 1st Column 0.237 2nd Column 0.008 3rd Column 0.029 2nd Row 1st Column left-parenthesis 0.207 right-parenthesis 2nd Column left-parenthesis 0.016 right-parenthesis 3rd Column left-parenthesis 0.049 right-parenthesis 3rd Row 1st Column negative 2.467 2nd Column 0.164 3rd Column negative 0.841 4th Row 1st Column left-parenthesis 1.930 right-parenthesis 2nd Column left-parenthesis 0.152 right-parenthesis 3rd Column left-parenthesis 0.453 right-parenthesis 5th Row 1st Column 0.951 2nd Column 0.002 3rd Column 0.938 6th Row 1st Column left-parenthesis 0.078 right-parenthesis 2nd Column left-parenthesis 0.006 right-parenthesis 3rd Column left-parenthesis 0.018 right-parenthesis EndMatrix Start 3 By 1 Matrix 1st Row  y Subscript 1 comma t minus 1 2nd Row  y Subscript 2 comma t minus 1 3rd Row  y Subscript 3 comma t minus 1 EndMatrix plus Start 3 By 1 Matrix 1st Row  epsilon Subscript 1 t 2nd Row  epsilon Subscript 2 t 3rd Row  epsilon Subscript 3 t EndMatrix EndLayout
Last updated: June 19, 2025