stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Quicker with partial correlations.

time-shift
Marco Ratto 2010-03-08 16:17:00 +01:00
parent 8a7bdefe2a
commit e66b7b5ee2
1 changed files with 11 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
matlab/identification_checks.m
View File

@ -144,21 +144,27 @@ PcoH = NaN(npar,npar);
PcoJ = NaN(npar,npar);
PcoGP = NaN(npar0,npar0);
for ii = 1:size(H1,2);
for jj = 1:size(H1,2);
PcoH(ind1(ii),ind1(ii)) = 1;
for jj = ii+1:size(H1,2);
PcoH(ind1(ii),ind1(jj)) = [cosn([H1(:,ii),H1(:,jj)])];
PcoH(ind1(jj),ind1(ii)) = PcoH(ind1(ii),ind1(jj));
end
end
for ii = 1:size(JJ1,2);
for jj = 1:size(JJ1,2);
PcoJ(ind2(ii),ind2(ii)) = 1;
for jj = ii+1:size(JJ1,2);
PcoJ(ind2(ii),ind2(jj)) = [cosn([JJ1(:,ii),JJ1(:,jj)])];
PcoJ(ind2(jj),ind2(ii)) = PcoJ(ind2(ii),ind2(jj));
end
end
for ii = 1:size(gp1,2);
for jj = 1:size(gp1,2);
PcoGP(ind3(ii),ind3(jj)) = [cosn([gp1(:,ii),gp1(:,jj)])];
end
PcoGP(ind3(ii),ind3(ii)) = 1;
for jj = ii+1:size(gp1,2);
PcoGP(ind3(ii),ind3(jj)) = [cosn([gp1(:,ii),gp1(:,jj)])];
PcoGP(ind3(jj),ind3(ii)) = PcoGP(ind3(ii),ind3(jj));
end
end