Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Exogenous regressors #4

Merged
merged 12 commits into from
Aug 4, 2017
Prev Previous commit
Next Next commit
"Changed the notation for the dimension of the exogenous regressors c…
…ontained in X

(m instead of w)"
  • Loading branch information
jm4474 committed Jul 29, 2017
commit 00bee7ceea4bafac2a3eee2efec8f52d04bd4da8
6 changes: 3 additions & 3 deletions functions/RForm/CovAhat_Sigmahat_Gamma.m
Original file line number Diff line number Diff line change
Expand Up @@ -22,7 +22,7 @@

k = size(Z,2);

w = size(X,2) - (n*p);
m = size(X,2) - (n*p);

XSVARp = X;

Expand All @@ -41,7 +41,7 @@

%Each 2-D page contains [eta_t W_t', eta_t X_t', eta_t eta_t'-Sigma, eta_tZ_t'-Gamma];

vecAss1 = reshape(auxeta,[(n*w)+(p*(n^2))+n^2+(n*k),1,T1aux]);
vecAss1 = reshape(auxeta,[(n*m)+(p*(n^2))+n^2+(n*k),1,T1aux]);

%Each 2-D page contains [vec(eta_tW_t'); vec(eta_tX_t') ; vec(eta_t*eta_t'-Sigma) ; vec(eta_tZ_t'-Gamma)]

Expand Down Expand Up @@ -73,7 +73,7 @@

Q2 = Z'*XSVARp/T1aux;

Shat = [kron([zeros(n*p,w),eye(n*p)]*(Q1^(-1)),eye(n)), zeros((n^2)*p,n^2+(k*n)); ...
Shat = [kron([zeros(n*p,m),eye(n*p)]*(Q1^(-1)),eye(n)), zeros((n^2)*p,n^2+(k*n)); ...
zeros(n*(n+1)/2,((n^2)*p)+n), V, zeros(n*(n+1)/2,k*n);...
-kron(Q2*(Q1^(-1)),eye(n)), zeros(k*n,n^2),eye(k*n) ];

Expand Down