关于dynare程序调试

giant 发表于 2013-10-3 22:38
另外问一下，书上稳态的产出大于消费，那么这样的模型还是steady的么？

你说的对，我验证了一下，产出和消费的缺口是投资，太谢谢了！我试试全部对数线性化后看看程序会跑出什么效果。

线性对数化完了以后出了这个报错，大神帮忙看看：
Error using evalin
Undefined function or variable 'line1004'.

Error in dynare (line 120)
evalin('base',fname) ;

程序如下：

var rc rd rl m dp c y w h k z a g;
varexo ze ae ge;
parameters beta delta theta b rhoz rhoa rhog;
beta=0.99;
delta=0.025;
theta=0.36;
b=-2.5805;
rhoz=0.005;
rhoa=0.005;
rhog=0.005;

model;
w=beta*0.0351*rc(-1)+w(-1);
z+w(-1)=c+rd;
0.9004*c=1.6676*m(-1)-0.7672*dp;
1.6676*m+12.4193*k=0.7672*(w(-1)+h)+0.4359*(rc+k(-1))+12.1088*k(-1)+0.7749*(rd+dp);
y=a+theta*k(-1)+(1-theta)*h;
1.2109*y=0.7749*(rl+w(-1)+h)+0.4359*(rc+k(-1));
rc+w(-1)=a+theta*(k(-1)-h);
rc=a+(theta-1)*(k(-1)-h);
0.7672*1.0101*(rl-rd)+1.6676*g=0;
0.7672*dp+1.6676*m=0.7672*(w(-1)+h);
z=rhoz*z(-1)+ze;
a=rhoa*a(-1)+ae;
g=rhog*g(-1)+ge;
end;

initval;
rc=0;
rd=0;
rl=0;
m=0;
dp=0;
c=0;
y=0;
w=0;
k=0;
h=0;
z=0;
a=0;
g=0;
end;
steady;

shocks;
var ze;stderr 0.009;
var ae;stderr 0.009;
var ge;stderr 0.099;
end;
stoch_simul;

giant 发表于 2013-10-4 12:55
线性对数化完了以后出了这个报错，大神帮忙看看：
Error using evalin
Undefined function or variable ' ...

怎么你线性化之后没有forward-looking variables呢？比如c(+1)

rastila 发表于 2013-10-4 15:58
怎么你线性化之后没有forward-looking variables呢？比如c(+1)

我看dynare userguide里说如果遇到这种变量，应该修改timing，把类似k(+1)都表达为k，所以……。
另外，事实上我的模型里有两个方程是关于c的（请看我附的图），有一个是c（+1），但我认为两个式子其实是一个等式，所以就删掉，不知道这样做对不对？

giant 发表于 2013-10-4 16:07
我看dynare userguide里说如果遇到这种变量，应该修改timing，把类似k(+1)都表达为k，所以……。
另外 ...

程序现在显示的是这个：2nd and 3rd order approximation not implemented for purely backward models
是不是就是你说的模型里没有forward-looking变量的原因啊？
Starting Dynare (version 4.3.3).
Starting preprocessing of the model file ...
Found 13 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:

rc                  0
rd                  0
rl                  0
m                   0
dp                  0
c                   0
y                   0
w                   0
h                   0
k                   0
z                   0
a                   0
g                   0
Error using stochastic_solvers (line 160)
2nd and 3rd order approximation not implemented for purely backward models

Error in resol (line 118)
    [dr,info] = stochastic_solvers(dr,check_flag,M,options,oo);

Error in stoch_simul (line 76)
    [oo_.dr,info,M_,options_,oo_] = resol(0,M_,options_,oo_);

Error in basic1004 (line 164)
info = stoch_simul(var_list_);

Error in dynare (line 120)
evalin('base',fname) ;

把model；改成model(linear);

rastila 发表于 2013-10-4 19:12
把model；改成model(linear);

改了，出了这个提示：Blanchard Kahn conditions are not satisfied: no stable equilibrium
我特别怕不符合BK条件，所以对原来的模型都没怎么改。就增加了一个外部变量（总共3个，分别是技术、消费和货币增长率冲击），总共13个变量

哎，有心而无力啊。我用dynare都是手动对数线性化，一来是以前习惯了，速度反而很快，二来是有时候要推loss function。
lz不妨试试手动线性化了再求解。

我先关注下这个帖子

猪人 发表于 2013-10-5 12:06
哎，有心而无力啊。我用dynare都是手动对数线性化，一来是以前习惯了，速度反而很快，二来是有时候要推loss ...

我就是手动线性化的啊，只是在原有模型里增加了一个冲击而已，不知为何BK条件不满足了？loss function是什么我不太清楚，麻烦解释一下，谢谢。