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

Bug fix: thanks to Johannes Pfeifer

time-shift
Marco Ratto 2011-06-06 15:21:59 +02:00
parent 4b85444b18
commit 2c977073d2
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
matlab/identification_analysis.m
View File

@ -177,7 +177,7 @@ if info(1)==0,
cmm = siJ(:,ind1)*((AHess(ind1,ind1))\siJ(:,ind1)');
chh = siH(:,ind1)*((AHess(ind1,ind1))\siH(:,ind1)');
ind1=ind1(ind1>offset);
clre = siLRE(:,ind1)*((AHess(ind1,ind1))\siLRE(:,ind1)');
clre = siLRE(:,ind1-offset)*((AHess(ind1,ind1))\siLRE(:,ind1-offset)');
rhoM=sqrt(1./diag(inv(tildaM(indok,indok))));
% deltaM = deltaM.*abs(params');
flag_score=1;
@ -201,7 +201,7 @@ if info(1)==0,
ind1=find(ide_hess.ind0);
chh = siH(:,ind1)*((MIM(ind1,ind1))\siH(:,ind1)');
ind1=ind1(ind1>offset);
clre = siLRE(:,ind1)*((MIM(ind1,ind1))\siLRE(:,ind1)');
clre = siLRE(:,ind1-offset)*((MIM(ind1,ind1))\siLRE(:,ind1-offset)');
rhoM(indok)=sqrt(1./diag(inv(tildaM(indok,indok))));
normaliz(indok) = (sqrt(diag(inv(tildaM(indok,indok))))./deltaM(indok))'; %sqrt(diag(inv(MIM(indok,indok))))';
% deltaM = deltaM.*abs(params');