How can I estimate a Vector Autoregressive (VAR) Model by OLS?

12 次查看（过去 30 天）
Valentina
Valentina 2013-10-17
评论： Spyridon Kariofylis 2020-7-10
How can I estimate a VAR Model with the vgxvarx function by Ordinary Least Squares ( OLS )?
By default vgxvarx uses the Maximum Likelihood and I can't find how to change it.
I tried with the default MLE:
Spec = vgxset('n',5,'nAR',1,'Constant',false);
EstSpec = vgxvarx(Spec, Y, [], Y0);
but I get this error:
Error using mvregress (line 425)
Covariance is not positive-definite.
Error in vgxvarx (line 521)
[x,Q,~,xvar] = mvregress(D,R, 'covtype',covartype, 'varformat',varformat, ...
Here:
http://www.mathworks.co.uk/help/econ/vgxvarx.html
it mentions to set 'MaxIter' to 1 for OLS, but it's not clear how to use the OLS approach.
  2 个评论
Shashank Prasanna
Shashank Prasanna 2013-10-17
Would you be willing to share why you don't want to use the MLE approach in the Econometrics toolbox?
Valentina
Valentina 2013-10-17
编辑：Valentina 2013-10-18
Because with MLE I get the above-mentioned error, while I was able to estimate the model with the same data, using a function in the R software that uses OLS.

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Hang Qian
Hang Qian 2014-3-30
编辑：Hang Qian 2014-3-30
Yes, estimation of a VAR(p) model by OLS is possible using the vgxvarx functionality. The vgxvarx uses maximum likelihood for rigorous treatment of missing values and presample values.
If the data are complete and presample values are specified (using the first p values of the data), vgxvarx will produce an estimator identical to the OLS estimator. For example, consider a VAR(2) model with 3 variables,
Y = rand(100,3);
Spec = vgxset('n',3,'nAR',2);
EstSpec = vgxvarx(Spec,Y(3:100,:),[],Y(1:2,:));
OLS1 = [EstSpec.AR{1},EstSpec.AR{2}]'
OLS2 = [Y(2:end-1,:),Y(1:end-2,:)] \ Y(3:end,:)
The second estimator is the raw OLS estimator.
norm(OLS1-OLS2) suggests that vgxvarx reproduces the raw OLS estimator.
Hang Qian

类别

在 Help Center 和 File Exchange 中查找有关 Conditional Mean Models 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by