- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 19470 个
- 通用积分
- 34.0385
- 学术水平
- 20 点
- 热心指数
- 41 点
- 信用等级
- 26 点
- 经验
- 15873 点
- 帖子
- 394
- 精华
- 1
- 在线时间
- 924 小时
- 注册时间
- 2008-7-5
- 最后登录
- 2024-4-9
副教授
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
- 19470 个
- 通用积分
- 34.0385
- 学术水平
- 20 点
- 热心指数
- 41 点
- 信用等级
- 26 点
- 经验
- 15873 点
- 帖子
- 394
- 精华
- 1
- 在线时间
- 924 小时
- 注册时间
- 2008-7-5
- 最后登录
- 2024-4-9
| 开心 2017-3-10 22:48:30 |
---|
签到天数: 2 天 连续签到: 1 天 [LV.1]初来乍到
|
20论坛币
我的代码如下,可是我不知道把冲击ed设置为多少算是一个单位的标准差的冲击,请前辈赐教,就是标红那行设定值为多少。
var y pi i w zd; %共5个
varexo ed ;
parameters eta fai delta_m delta_p sigma kappa zeta gamma omega_pi omega_i;
eta = 633.6624;
fai = 14.61883;
delta_m = 0.4777149;
delta_p = 0.4955435;
sigma = 0.0003956;
kappa = 1924.567;
zeta = 0.1521255;
gamma = -100.5403;
omega_pi = 1;
omega_i = 1.906987;
model;
y = delta_m * y(-1) + delta_p * y(+1) - sigma * ( i - pi(+1) ) + ed;% 需求正向冲击
i = i(-1) + sigma * y / omega_i + omega_pi * ( kappa + gamma * w ) * sigma * pi / omega_i;
pi = kappa * y + zeta * w + gamma * w * y;%可能正可能负,可能恒定
w = -kappa * y / ( fai * ( gamma * y + zeta ) ) ; % eta * y - fai * zd;%eta * y + fai * zd;
zd = w/fai - eta * y / fai;
end;
initval;%把稳态值赋予给内生变量
y = 0;
pi = 0;
i = 0;
w = 0;
zd = 0;
end;
shocks;
%var ea =0.1^2;%技术冲击
%var ez =0.1^2;%冲击
var ed = 0.001^2;%冲击
end;
resid(1);%用于判断给定的初始值是不是稳态值
steady;
check;
%when used in loops, noprint is a good option.
stoch_simul(irf=30,periods=5000,replic=1000);
|
|