我用ML估计正确，但贝叶斯估计出错了，Mod 文件和出错信息如下，烦请老师指点：

var x y;
varexo e,u;

parameters rho beta;
beta=0.9;
rho=0.9;
model(linear);
y=beta*y(+1)+x+u;
x=rho*x(-1)+e;
end;
shocks;
var e; stderr 0.1;
var u; stderr 0.1;
end;


stoch_simul(periods=900,order=1,irf=40,drop=400,aim_solver) y x;

save Simul_data y x;
varobs y x;
estimated_params;
rho, normal_pdf, 0.7, 0.2;
end;

//estimation(datafile=Simul_data, mode_check, prefilter=1, order=1, mode_compute=1) y x;
estimation(datafile=Simul_data,prefilter=1,mh_replic=500,mh_nblocks=2,mh_jscale=0.4,mode_compute=1,order=1,bayesian_irf,irf=30,mode_check) y x;

出错信息：

                                Max     Line search  Directional  First-order
Iter F-count        f(x)   constraint   steplength   derivative   optimality Procedure
    0      2      15050.8       -1.272                                         
Error in AIM: aimcode=3 : Aim: too many big roots.
??? Error using ==> mtimes
Inner matrix dimensions must agree.

Error in ==> DsgeLikelihood at 172
    Pstar =
    lyapunov_symm(T,R*Q*R',options_.qz_criterium,options_.lyapunov_complex_threshold);
   
Error in ==> nlconst at 814
                f = feval(funfcn{3},x,varargin{:});

Error in ==> fmincon at 720
    [X,FVAL,LAMBDA,EXITFLAG,OUTPUT,GRAD,HESSIAN]=...

Error in ==> dynare_estimation_1 at 180
            [xparam1,fval,exitflag,output,lamdba,grad,hessian_fmincon] = ...

Error in ==> dynare_estimation at 62
    dynare_estimation_1(var_list,varargin{:});

Error in ==> simplest at 113
dynare_estimation(var_list_);

Error in ==> dynare at 120
evalin('base',fname) ;

把下面的命令复制到.mod，然后跑一遍看。
var x y;
varexo e,u;

parameters rho beta;
beta=0.9;
rho=0.9;
model(linear);
y=beta*y(+1)+x+u;
x=rho*x(-1)+e;
end;
shocks;
var e; stderr 0.1;
var u; stderr 0.1;
end;


stoch_simul(periods=900,order=1,irf=40,drop=400, nograph) y x;

save Simul_data y x;
varobs y x;
estimated_params;
rho, normal_pdf, 0.7, 0.2;
end;

//estimation(datafile=Simul_data, mode_check, prefilter=1, order=1, mode_compute=1) y x;
estimation(datafile=Simul_data,prefilter=1,mh_replic=500,mh_nblocks=2,mh_jscale=0.4,mode_compute=1,order=1,bayesian_irf,irf=30,mode_check) y x;

我去掉了"aim_solver"。你可能会得到下面这张图：

这个图横轴上的红点，都是违反BK condition的地方。

谢谢老师，对比了一下文件，出错在这里：
stoch_simul(periods=900,order=1,irf=40,drop=400,nograph) y x;

我原来的文件里是：
stoch_simul(periods=900,order=1,irf=40,drop=400,aim_solver) y x;

在最后一个参数的差别，不解，看了reference manual，这里这样讲的：
Use the Anderson-Moore Algorithm (AIM) to compute the decision rules, instead
of using Dynare’s default method based on a generalized Schur decomposition. This
option is only valid for first order approximation. See AIM website for more details
on the algorithm.
难道说这个算法有问题？

这个算法没有问题。我感觉问题出在模型某个设置上面，模型的某个地方和这个算法有冲突，但是模型能够被这个算法解出saddle-path，但这个算法在参数估计的时候和模型发生冲突，导致某个地方的矩阵维度坍塌下来了，第一个红色错误信息，
??? Error using ==> mtimes
Inner matrix dimensions must agree.
这个技术错误是来自Matlab的，不是Dynare。我需要想一下，如果不是程序的bug，那就一定是某个地方我模型设置和AIM冲突了。

非常感谢，感觉有点入门了。dynare感觉还不很成熟，我在dynare中某个地方块注释忘了删去，它只在运行时报告。如
上面的程序，我在红色的地方漏了删除一个块注释符：
var x y;
varexo e,u;

parameters rho beta;
beta=0.9;
rho=0.9;
model(linear);
y=beta*y(+1)+x+u;
x=rho*x(-1)+e;
end;
shocks;
var e; stderr 0.1;
var u; stderr 0.1;
end;

*/
stoch_simul(periods=900,order=1,irf=40,drop=400,nograph) y x;

save Simul_data y x;
varobs y x;
estimated_params;
rho, normal_pdf, 0.7, 0.2;
end;

estimation(datafile=Simul_data, mode_check, prefilter=1, order=1, mode_compute=1) y x;
//estimation(datafile=Simul_data,prefilter=1,mh_replic=500,mh_nblocks=2,mh_jscale=0.4,mode_compute=1,order=1,bayesian_irf,irf=30,mode_check) y x;
，它只在运行时报告说：
??? Error: File: simplest.m Line: 83 Column: 1
Unexpected MATLAB operator.

Error in ==> dynare at 120
evalin('base',fname) ;
有些功能应该放在dynare中完成为好，如dynare的注释和matlab是不同的，检查这方面的错误应由dynare完成；在参数估计时如果未增加冲击变量，可能会出现kalman滤波时的条件不满足，这在预处理时似乎也能由dynare提醒用户，但它没这么做或没做到。

我很高兴你能从中学到些东西。

我后面还打算写3个关于DSGE的帖子，准备把一个新凯恩斯的benchmark模型解构了放在这里，你后面可以看到会碰到的问题远远比我们能想象到的多。光是一个calibration就会带来无数的问题，后面的参数估计更是复杂。相比来看前面算optimality condition，算saddle-path难度都算低的。

后面你会看到更多的关于dynare应用的技巧，比如和matlab之间的高度融合做一些循环运算和编程设计。

感谢楼主，下载好好学习。

来捧场了

谢谢楼主分享 虽然不懂 但还是先收藏了

你好 ，你照着我写的这个note做一遍，你自然就看到Metropolis-Hasting算法了。

DSGE模型挺好的

留着以后看

自己编了一个mod，出现以下问题
There are 5 eigenvalue(s) larger than 1 in modulus
for 4 forward-looking variable(s)
Blanchard Kahn conditions are not satisfied: no stable equilibrium
请问各位高手如何解决啊，多谢

又见最新note，支持

这个帖子太好啦

虽然不是很懂，但是感觉很厉害的样子

学习了！

请问可以申请读你的博士吗？

看来高手还挺多。

你能传mod上来么？

出现以下错误信息：

There are 5 eigenvalue(s) larger than 1 in modulus
for 4 forward-looking variable(s)

The rank conditions ISN'T verified!

??? Error using ==> print_info at 40
Blanchard Kahn conditions are not satisfied: no stable equilibrium
Error in ==> stoch_simul at 71
    print_info(info, options_.noprint);
Error in ==> doctor at 264
info = stoch_simul(var_list_);
Error in ==> dynare at 120
evalin('base',fname) ;

运行出来的完整信息如下：
Configuring Dynare ...
[mex] Generalized QZ.
[mex] Sylvester equation solution.
[mex] Kronecker products.
[mex] Sparse kronecker products.
[mex] Bytecode evaluation.
[mex] k-order perturbation solver.
[mex] k-order solution simulation.

Starting Dynare (version 4.2.5).
Starting preprocessing of the model file ...
WARNING: doctor.mod:78.13-21: dr_algo option is now deprecated, and may be removed in a future version of Dynare
Found 18 equation(s).
Evaluating expressions...done
Computing static model derivatives:
- order 1
Computing dynamic model derivatives:
- order 1
- order 2
Processing outputs ...done
Preprocessing completed.
Starting MATLAB/Octave computing.


STEADY-STATE RESULTS:

c1                  0
c2                  0
y                   0
l1                  0
w1                  0
l2                  0
w2                  0
x                   0
b1                  0
b2                  0
h1                  0
h2                  0
pi                  0
r                   0
q                   0
A                   0
d                   0
er                  0

EIGENVALUES:
         Modulus             Real        Imaginary

      1.666e-016      -1.666e-016                0
             0.9              0.9                0
            0.95             0.95                0
            0.98             0.98                0
           1.007            1.007                0
           1.045            1.045                0
           1.202             1.18           0.2289
           1.202             1.18          -0.2289
           2.003            2.003                0


There are 5 eigenvalue(s) larger than 1 in modulus
for 4 forward-looking variable(s)

The rank conditions ISN'T verified!

??? Error using ==> print_info at 40
Blanchard Kahn conditions are not satisfied: no stable equilibrium

Error in ==> stoch_simul at 71
    print_info(info, options_.noprint);

Error in ==> doctor at 264
info = stoch_simul(var_list_);

Error in ==> dynare at 120
evalin('base',fname) ;
谢谢！

你把下面的程序到.mod里，你自己跑一遍看。
var  c1 c2 y l1  w1  l2 w2 x b1 b2 h1 h2 pi r q A d ;
varexo espor ed ea er;
%predetermined_parameter q;

%%
parameters beta1 alpha xx rr m2 dd gmae hh1 eta beta2 theta c1toy c2toy qh1toy b2toy k b1toy qh2toy deta ruod ruoa ruor R;
alpha=0.64;
beta1=0.99;
dd=0.08;
xx=1.05;
rr=0.3;
hh1=0.1;
gmae=0.99;
eta=0.54;
beta2=0.98;
theta=0.75;
deta=0.7;
ruod=0.98;
ruoa=0.95;
ruor=0.9;
R=1/beta1;
m2=0.9;
b1toy=0;

c2toy=(1-alpha)*xx/(1+(rr-1)*beta1*m2*dd/(1-gmae));
qh2toy=c2toy*(dd/(1-gmae));
b2toy=beta1*m2*qh2toy;
c1toy=1-c2toy;
qh1toy=c1toy*dd/hh1;
k=(1-theta)/theta*(1-beta1*theta);

%%
initval;
c1=0;
c2=0;
y=0;
l1=0;
w1=0;
  l2=0;
w2=0;
x=0;
b1=0;
b2=0;
h1=0;
h2=0;
pi=0;
r=0;
q=0;
A=0;
d=0;
%er=0;
end;

%%
model(linear);
y=c1toy*c1+c2toy*c2;
c1=-(eta-1)*l1+w1;
c2=-(eta-1)*l2+w2;
y=A+alpha*l1+(1-alpha)*l2;
x=y-l1-w1;
l2+w2=y-x;
b1toy*b1=b2toy*b2;
h1*(-qh1toy)=h2*qh2toy;
c1=c1(+1)+pi(+1)-r;
q=(1-beta1)*(d-h1)+beta1*(q(+1)-c1(+1))+c1;
b2toy*b2=c2toy*c2+qh2toy*deta*h2+rr*b2toy*(b2(-1)+r(-1)-pi)-(1-alpha)*(y-x)/xx;
q=(beta2+m2*(beta1-beta2))*q(+1)+(1-(beta2+m2*(beta1-beta2)))*(d-h2)+(1-m2*beta1)*(c2-(beta2-m2*beta2)/(1-m2*beta1)*c2(+1))-m2*beta1*(r-pi(+1));
b2=q(+1)+h2+pi(+1)-r;
%pi=beta1*pi(+1)-k*x;
pi=beta1*pi(+1)-k*x(+1);         %what is 'x'? what does theory says about it?
r=0.95*pi+0.95*y+er;
d=ruod*d(-1)+ed;
%A=ruoa*A(-1)+ea;
A=0.9*A(-1)+ea;
%er=ruor*er(-1)+espor;
end;


%%
steady;
check;
resid;

shocks;
%var espor;
%stderr 0.8;
var ed;
stderr 0.08;
var ea;
stderr 2.5;
var er;
stderr 2.5;
end;

%%
stoch_simul(irf=20,nograph,order=1, periods=500) c1 c2 y l1  w1  l2 w2 x b1 b2 h1 h2 pi r q A d;
%stoch_simul(dr_algo=0,periods=1000); datatomfile('simudata',[]);

我仅仅只是把它改得能跑完模拟，但是问题很多。有几个变量都没法模拟出数值来。你如果参数设置完全是按照谋篇论文来的，那么可能参数设置问题不大，但如果你模型有任何修改的地方，都可以导致系统不稳定。我不知道你的.mod根据什么论文来做的，但是我认为你需要找错的话，先从模型逻辑上出发。

我没法帮你找出具体的错误，因为这需要对你模型的逻辑非常熟悉。需要你自己再自己读一下论文，同时看有那些变量根据理论是需要定为predetermined或者是nonpredetermined。模型的推导如果有错误也会导致模型解不出来。