Replace initval blocks by steady_state_model blocks in fs2000 unit tests
Speeds up estimationtime-shift
parent
e07607011c
commit
0104f91f4b
|
@ -52,21 +52,34 @@ P2 = P(+1);
|
||||||
c2 = c(+1);
|
c2 = c(+1);
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
|
Y_obs = 1;
|
||||||
|
P_obs = 1;
|
||||||
|
P2 = P;
|
||||||
|
c2 = c;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -52,23 +52,37 @@ P2 = P(+1);
|
||||||
c2 = c(+1);
|
c2 = c(+1);
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
|
Y_obs = 1;
|
||||||
|
P_obs = 1;
|
||||||
|
P2 = P;
|
||||||
|
c2 = c;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
var e_a; stderr 0.014;
|
var e_a; stderr 0.014;
|
||||||
var e_m; stderr 0.005;
|
var e_m; stderr 0.005;
|
||||||
|
|
|
@ -1,65 +0,0 @@
|
||||||
% computes the steady state of fs2000 analyticaly
|
|
||||||
% largely inspired by the program of F. Schorfheide
|
|
||||||
function [ys,check] = fs2000k_steadystate(ys,exe)
|
|
||||||
global M_
|
|
||||||
|
|
||||||
alp = M_.params(1);
|
|
||||||
bet = M_.params(2);
|
|
||||||
gam = M_.params(3);
|
|
||||||
mst = M_.params(4);
|
|
||||||
rho = M_.params(5);
|
|
||||||
psi = M_.params(6);
|
|
||||||
del = M_.params(7);
|
|
||||||
|
|
||||||
check = 0;
|
|
||||||
|
|
||||||
dA = exp(gam);
|
|
||||||
gst = 1/dA;
|
|
||||||
m = mst;
|
|
||||||
|
|
||||||
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
|
||||||
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
|
||||||
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
|
||||||
n = xist/(nust+xist);
|
|
||||||
P = xist + nust;
|
|
||||||
k = khst*n;
|
|
||||||
|
|
||||||
l = psi*mst*n/( (1-psi)*(1-n) );
|
|
||||||
c = mst/P;
|
|
||||||
d = l - mst + 1;
|
|
||||||
y = k^alp*n^(1-alp)*gst^alp;
|
|
||||||
R = mst/bet;
|
|
||||||
W = l/n;
|
|
||||||
ist = y-c;
|
|
||||||
q = 1 - d;
|
|
||||||
|
|
||||||
e = 1;
|
|
||||||
|
|
||||||
gp_obs = m/dA;
|
|
||||||
gy_obs = dA;
|
|
||||||
|
|
||||||
P_obs = 1;
|
|
||||||
Y_obs = 1;
|
|
||||||
|
|
||||||
P2=P;
|
|
||||||
c2=c;
|
|
||||||
|
|
||||||
ys =[
|
|
||||||
m
|
|
||||||
P
|
|
||||||
c
|
|
||||||
e
|
|
||||||
W
|
|
||||||
R
|
|
||||||
k
|
|
||||||
d
|
|
||||||
n
|
|
||||||
l
|
|
||||||
gy_obs
|
|
||||||
gp_obs
|
|
||||||
Y_obs
|
|
||||||
P_obs
|
|
||||||
y
|
|
||||||
dA
|
|
||||||
P2
|
|
||||||
c2 ];
|
|
|
@ -1,65 +0,0 @@
|
||||||
% computes the steady state of fs2000 analyticaly
|
|
||||||
% largely inspired by the program of F. Schorfheide
|
|
||||||
function [ys,check] = fs2000_b1L1L_steadystate(ys,exe)
|
|
||||||
global M_
|
|
||||||
|
|
||||||
alp = M_.params(1);
|
|
||||||
bet = M_.params(2);
|
|
||||||
gam = M_.params(3);
|
|
||||||
mst = M_.params(4);
|
|
||||||
rho = M_.params(5);
|
|
||||||
psi = M_.params(6);
|
|
||||||
del = M_.params(7);
|
|
||||||
|
|
||||||
check = 0;
|
|
||||||
|
|
||||||
dA = exp(gam);
|
|
||||||
gst = 1/dA;
|
|
||||||
m = mst;
|
|
||||||
|
|
||||||
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
|
||||||
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
|
||||||
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
|
||||||
n = xist/(nust+xist);
|
|
||||||
P = xist + nust;
|
|
||||||
k = khst*n;
|
|
||||||
|
|
||||||
l = psi*mst*n/( (1-psi)*(1-n) );
|
|
||||||
c = mst/P;
|
|
||||||
d = l - mst + 1;
|
|
||||||
y = k^alp*n^(1-alp)*gst^alp;
|
|
||||||
R = mst/bet;
|
|
||||||
W = l/n;
|
|
||||||
ist = y-c;
|
|
||||||
q = 1 - d;
|
|
||||||
|
|
||||||
e = 1;
|
|
||||||
|
|
||||||
gp_obs = m/dA;
|
|
||||||
gy_obs = dA;
|
|
||||||
|
|
||||||
P_obs = 1;
|
|
||||||
Y_obs = 1;
|
|
||||||
|
|
||||||
P2=P;
|
|
||||||
c2=c;
|
|
||||||
|
|
||||||
ys =[
|
|
||||||
m
|
|
||||||
P
|
|
||||||
c
|
|
||||||
e
|
|
||||||
W
|
|
||||||
R
|
|
||||||
k
|
|
||||||
d
|
|
||||||
n
|
|
||||||
l
|
|
||||||
gy_obs
|
|
||||||
gp_obs
|
|
||||||
Y_obs
|
|
||||||
P_obs
|
|
||||||
y
|
|
||||||
dA
|
|
||||||
P2
|
|
||||||
c2 ];
|
|
|
@ -353,8 +353,6 @@ EXTRA_DIST = \
|
||||||
printMakeCheckOctaveErrMsg.m \
|
printMakeCheckOctaveErrMsg.m \
|
||||||
fataltest.m \
|
fataltest.m \
|
||||||
AIM/data_ca1.m \
|
AIM/data_ca1.m \
|
||||||
AIM/fs2000_b1L1L_AIM_steadystate.m \
|
|
||||||
AIM/fs2000_b1L1L_steadystate.m \
|
|
||||||
AIM/fsdat.m \
|
AIM/fsdat.m \
|
||||||
block_bytecode/run_ls2003.m \
|
block_bytecode/run_ls2003.m \
|
||||||
bvar_a_la_sims/bvar_sample.m \
|
bvar_a_la_sims/bvar_sample.m \
|
||||||
|
@ -365,12 +363,10 @@ EXTRA_DIST = \
|
||||||
steady_state/walsh1_old_ss_steadystate.m \
|
steady_state/walsh1_old_ss_steadystate.m \
|
||||||
data/test.xls \
|
data/test.xls \
|
||||||
analytic_derivatives/fsdat_simul.m \
|
analytic_derivatives/fsdat_simul.m \
|
||||||
fs2000/fs2000a_steadystate.m \
|
|
||||||
fs2000/fsdat_simul.m \
|
fs2000/fsdat_simul.m \
|
||||||
k_order_perturbation/run_fs2000kplusplus.m \
|
k_order_perturbation/run_fs2000kplusplus.m \
|
||||||
ls2003/data_ca1.m \
|
ls2003/data_ca1.m \
|
||||||
measurement_errors/data_ca1.m \
|
measurement_errors/data_ca1.m \
|
||||||
measurement_errors/fs2000_corr_me_ml_mcmc/fs2000_corr_ME_steadystate.m \
|
|
||||||
measurement_errors/fs2000_corr_me_ml_mcmc/fsdat_simul.m \
|
measurement_errors/fs2000_corr_me_ml_mcmc/fsdat_simul.m \
|
||||||
missing/simulate_data_with_missing_observations.m \
|
missing/simulate_data_with_missing_observations.m \
|
||||||
objectives/sgu_ex1.mat \
|
objectives/sgu_ex1.mat \
|
||||||
|
@ -393,18 +389,14 @@ EXTRA_DIST = \
|
||||||
ms-sbvar/archive-files/specification_2v2c.dat \
|
ms-sbvar/archive-files/specification_2v2c.dat \
|
||||||
recursive/data_ca1.m \
|
recursive/data_ca1.m \
|
||||||
kalman_filter_smoother/fsdat_simul.m \
|
kalman_filter_smoother/fsdat_simul.m \
|
||||||
kalman_filter_smoother/fs2000a_steadystate.m \
|
|
||||||
identification/kim/kim2_steadystate.m \
|
identification/kim/kim2_steadystate.m \
|
||||||
|
identification/as2007/as2007_steadystate.m \
|
||||||
estimation/fsdat_simul.m \
|
estimation/fsdat_simul.m \
|
||||||
ep/mean_preserving_spread.m \
|
ep/mean_preserving_spread.m \
|
||||||
third_order/comparison_policy_functions_dynare_mathematica.m \
|
third_order/comparison_policy_functions_dynare_mathematica.m \
|
||||||
third_order/policyfunctions.mat \
|
third_order/policyfunctions.mat \
|
||||||
shock_decomposition/example1_calib_shock_decomp_data.mat \
|
shock_decomposition/example1_calib_shock_decomp_data.mat \
|
||||||
shock_decomposition/fsdat_simul.m \
|
shock_decomposition/fsdat_simul.m \
|
||||||
estimation/fs2000_MCMC_jumping_covariance_steadystate.m \
|
|
||||||
estimation/fs2000_initialize_from_calib_steadystate.m \
|
|
||||||
filter_step_ahead/fs2000_filter_step_ahead_bayesian_steadystate.m \
|
|
||||||
filter_step_ahead/fs2000_filter_step_ahead_ML_steadystate.m \
|
|
||||||
loglinear/results_exp.mat \
|
loglinear/results_exp.mat \
|
||||||
smoother2histval/fsdat_simul.m \
|
smoother2histval/fsdat_simul.m \
|
||||||
optimal_policy/Ramsey/find_c.m \
|
optimal_policy/Ramsey/find_c.m \
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -32,21 +32,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,23 +30,33 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
var e_a; stderr 0.014;
|
var e_a; stderr 0.014;
|
||||||
var e_m; stderr 0.005;
|
var e_m; stderr 0.005;
|
||||||
|
@ -72,4 +82,4 @@ varobs gp_obs gy_obs;
|
||||||
|
|
||||||
options_.solve_tolf = 1e-12;
|
options_.solve_tolf = 1e-12;
|
||||||
|
|
||||||
estimation(order=1,datafile=fsdat_simul,nobs=192,loglinear,mh_replic=2000,mh_nblocks=2,mh_jscale=0.8);
|
estimation(order=1,datafile=fsdat_simul,nobs=192,loglinear,mh_replic=3000,mh_nblocks=2,mh_jscale=0.8,moments_varendo,selected_variables_only) y m;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -1,73 +0,0 @@
|
||||||
% computes the steady state of fs2000 analyticaly
|
|
||||||
% largely inspired by the program of F. Schorfheide
|
|
||||||
|
|
||||||
% Copyright (C) 2004-2010 Dynare Team
|
|
||||||
%
|
|
||||||
% This file is part of Dynare.
|
|
||||||
%
|
|
||||||
% Dynare is free software: you can redistribute it and/or modify
|
|
||||||
% it under the terms of the GNU General Public License as published by
|
|
||||||
% the Free Software Foundation, either version 3 of the License, or
|
|
||||||
% (at your option) any later version.
|
|
||||||
%
|
|
||||||
% Dynare is distributed in the hope that it will be useful,
|
|
||||||
% but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
||||||
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
||||||
% GNU General Public License for more details.
|
|
||||||
%
|
|
||||||
% You should have received a copy of the GNU General Public License
|
|
||||||
% along with Dynare. If not, see <http://www.gnu.org/licenses/>.
|
|
||||||
|
|
||||||
function [ys,check] = fs2000_MCMC_jumping_covariance_steadystate(ys,exe)
|
|
||||||
global M_
|
|
||||||
|
|
||||||
alp = M_.params(1);
|
|
||||||
bet = M_.params(2);
|
|
||||||
gam = M_.params(3);
|
|
||||||
mst = M_.params(4);
|
|
||||||
rho = M_.params(5);
|
|
||||||
psi = M_.params(6);
|
|
||||||
del = M_.params(7);
|
|
||||||
|
|
||||||
check = 0;
|
|
||||||
|
|
||||||
dA = exp(gam);
|
|
||||||
gst = 1/dA;
|
|
||||||
m = mst;
|
|
||||||
|
|
||||||
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
|
||||||
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
|
||||||
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
|
||||||
n = xist/(nust+xist);
|
|
||||||
P = xist + nust;
|
|
||||||
k = khst*n;
|
|
||||||
|
|
||||||
l = psi*mst*n/( (1-psi)*(1-n) );
|
|
||||||
c = mst/P;
|
|
||||||
d = l - mst + 1;
|
|
||||||
y = k^alp*n^(1-alp)*gst^alp;
|
|
||||||
R = mst/bet;
|
|
||||||
W = l/n;
|
|
||||||
ist = y-c;
|
|
||||||
q = 1 - d;
|
|
||||||
|
|
||||||
e = 1;
|
|
||||||
|
|
||||||
gp_obs = m/dA;
|
|
||||||
gy_obs = dA;
|
|
||||||
|
|
||||||
ys =[
|
|
||||||
m
|
|
||||||
P
|
|
||||||
c
|
|
||||||
e
|
|
||||||
W
|
|
||||||
R
|
|
||||||
k
|
|
||||||
d
|
|
||||||
n
|
|
||||||
l
|
|
||||||
gy_obs
|
|
||||||
gp_obs
|
|
||||||
y
|
|
||||||
dA ];
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
varobs gp_obs gy_obs;
|
varobs gp_obs gy_obs;
|
||||||
|
|
|
@ -29,23 +29,33 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
var e_a; stderr 0.014;
|
var e_a; stderr 0.014;
|
||||||
var e_m; stderr 0.005;
|
var e_m; stderr 0.005;
|
||||||
|
|
|
@ -1,73 +0,0 @@
|
||||||
% computes the steady state of fs2000 analyticaly
|
|
||||||
% largely inspired by the program of F. Schorfheide
|
|
||||||
|
|
||||||
% Copyright (C) 2004-2010 Dynare Team
|
|
||||||
%
|
|
||||||
% This file is part of Dynare.
|
|
||||||
%
|
|
||||||
% Dynare is free software: you can redistribute it and/or modify
|
|
||||||
% it under the terms of the GNU General Public License as published by
|
|
||||||
% the Free Software Foundation, either version 3 of the License, or
|
|
||||||
% (at your option) any later version.
|
|
||||||
%
|
|
||||||
% Dynare is distributed in the hope that it will be useful,
|
|
||||||
% but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
||||||
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
||||||
% GNU General Public License for more details.
|
|
||||||
%
|
|
||||||
% You should have received a copy of the GNU General Public License
|
|
||||||
% along with Dynare. If not, see <http://www.gnu.org/licenses/>.
|
|
||||||
|
|
||||||
function [ys,check] = fs2000_steadystate(ys,exe)
|
|
||||||
global M_
|
|
||||||
|
|
||||||
alp = M_.params(1);
|
|
||||||
bet = M_.params(2);
|
|
||||||
gam = M_.params(3);
|
|
||||||
mst = M_.params(4);
|
|
||||||
rho = M_.params(5);
|
|
||||||
psi = M_.params(6);
|
|
||||||
del = M_.params(7);
|
|
||||||
|
|
||||||
check = 0;
|
|
||||||
|
|
||||||
dA = exp(gam);
|
|
||||||
gst = 1/dA;
|
|
||||||
m = mst;
|
|
||||||
|
|
||||||
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
|
||||||
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
|
||||||
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
|
||||||
n = xist/(nust+xist);
|
|
||||||
P = xist + nust;
|
|
||||||
k = khst*n;
|
|
||||||
|
|
||||||
l = psi*mst*n/( (1-psi)*(1-n) );
|
|
||||||
c = mst/P;
|
|
||||||
d = l - mst + 1;
|
|
||||||
y = k^alp*n^(1-alp)*gst^alp;
|
|
||||||
R = mst/bet;
|
|
||||||
W = l/n;
|
|
||||||
ist = y-c;
|
|
||||||
q = 1 - d;
|
|
||||||
|
|
||||||
e = 1;
|
|
||||||
|
|
||||||
gp_obs = m/dA;
|
|
||||||
gy_obs = dA;
|
|
||||||
|
|
||||||
ys =[
|
|
||||||
m
|
|
||||||
P
|
|
||||||
c
|
|
||||||
e
|
|
||||||
W
|
|
||||||
R
|
|
||||||
k
|
|
||||||
d
|
|
||||||
n
|
|
||||||
l
|
|
||||||
gy_obs
|
|
||||||
gp_obs
|
|
||||||
y
|
|
||||||
dA ];
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -65,21 +65,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -1,73 +0,0 @@
|
||||||
% computes the steady state of fs2000 analyticaly
|
|
||||||
% largely inspired by the program of F. Schorfheide
|
|
||||||
|
|
||||||
% Copyright (C) 2004-2010 Dynare Team
|
|
||||||
%
|
|
||||||
% This file is part of Dynare.
|
|
||||||
%
|
|
||||||
% Dynare is free software: you can redistribute it and/or modify
|
|
||||||
% it under the terms of the GNU General Public License as published by
|
|
||||||
% the Free Software Foundation, either version 3 of the License, or
|
|
||||||
% (at your option) any later version.
|
|
||||||
%
|
|
||||||
% Dynare is distributed in the hope that it will be useful,
|
|
||||||
% but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
||||||
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
||||||
% GNU General Public License for more details.
|
|
||||||
%
|
|
||||||
% You should have received a copy of the GNU General Public License
|
|
||||||
% along with Dynare. If not, see <http://www.gnu.org/licenses/>.
|
|
||||||
|
|
||||||
function [ys,check] = fs2000_steadystate(ys,exe)
|
|
||||||
global M_
|
|
||||||
|
|
||||||
alp = M_.params(1);
|
|
||||||
bet = M_.params(2);
|
|
||||||
gam = M_.params(3);
|
|
||||||
mst = M_.params(4);
|
|
||||||
rho = M_.params(5);
|
|
||||||
psi = M_.params(6);
|
|
||||||
del = M_.params(7);
|
|
||||||
|
|
||||||
check = 0;
|
|
||||||
|
|
||||||
dA = exp(gam);
|
|
||||||
gst = 1/dA;
|
|
||||||
m = mst;
|
|
||||||
|
|
||||||
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
|
||||||
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
|
||||||
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
|
||||||
n = xist/(nust+xist);
|
|
||||||
P = xist + nust;
|
|
||||||
k = khst*n;
|
|
||||||
|
|
||||||
l = psi*mst*n/( (1-psi)*(1-n) );
|
|
||||||
c = mst/P;
|
|
||||||
d = l - mst + 1;
|
|
||||||
y = k^alp*n^(1-alp)*gst^alp;
|
|
||||||
R = mst/bet;
|
|
||||||
W = l/n;
|
|
||||||
ist = y-c;
|
|
||||||
q = 1 - d;
|
|
||||||
|
|
||||||
e = 1;
|
|
||||||
|
|
||||||
gp_obs = m/dA;
|
|
||||||
gy_obs = dA;
|
|
||||||
|
|
||||||
ys =[
|
|
||||||
m
|
|
||||||
P
|
|
||||||
c
|
|
||||||
e
|
|
||||||
W
|
|
||||||
R
|
|
||||||
k
|
|
||||||
d
|
|
||||||
n
|
|
||||||
l
|
|
||||||
gy_obs
|
|
||||||
gp_obs
|
|
||||||
y
|
|
||||||
dA ];
|
|
|
@ -65,21 +65,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -1,73 +0,0 @@
|
||||||
% computes the steady state of fs2000 analyticaly
|
|
||||||
% largely inspired by the program of F. Schorfheide
|
|
||||||
|
|
||||||
% Copyright (C) 2004-2010 Dynare Team
|
|
||||||
%
|
|
||||||
% This file is part of Dynare.
|
|
||||||
%
|
|
||||||
% Dynare is free software: you can redistribute it and/or modify
|
|
||||||
% it under the terms of the GNU General Public License as published by
|
|
||||||
% the Free Software Foundation, either version 3 of the License, or
|
|
||||||
% (at your option) any later version.
|
|
||||||
%
|
|
||||||
% Dynare is distributed in the hope that it will be useful,
|
|
||||||
% but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
||||||
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
||||||
% GNU General Public License for more details.
|
|
||||||
%
|
|
||||||
% You should have received a copy of the GNU General Public License
|
|
||||||
% along with Dynare. If not, see <http://www.gnu.org/licenses/>.
|
|
||||||
|
|
||||||
function [ys,check] = fs2000_steadystate(ys,exe)
|
|
||||||
global M_
|
|
||||||
|
|
||||||
alp = M_.params(1);
|
|
||||||
bet = M_.params(2);
|
|
||||||
gam = M_.params(3);
|
|
||||||
mst = M_.params(4);
|
|
||||||
rho = M_.params(5);
|
|
||||||
psi = M_.params(6);
|
|
||||||
del = M_.params(7);
|
|
||||||
|
|
||||||
check = 0;
|
|
||||||
|
|
||||||
dA = exp(gam);
|
|
||||||
gst = 1/dA;
|
|
||||||
m = mst;
|
|
||||||
|
|
||||||
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
|
||||||
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
|
||||||
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
|
||||||
n = xist/(nust+xist);
|
|
||||||
P = xist + nust;
|
|
||||||
k = khst*n;
|
|
||||||
|
|
||||||
l = psi*mst*n/( (1-psi)*(1-n) );
|
|
||||||
c = mst/P;
|
|
||||||
d = l - mst + 1;
|
|
||||||
y = k^alp*n^(1-alp)*gst^alp;
|
|
||||||
R = mst/bet;
|
|
||||||
W = l/n;
|
|
||||||
ist = y-c;
|
|
||||||
q = 1 - d;
|
|
||||||
|
|
||||||
e = 1;
|
|
||||||
|
|
||||||
gp_obs = m/dA;
|
|
||||||
gy_obs = dA;
|
|
||||||
|
|
||||||
ys =[
|
|
||||||
m
|
|
||||||
P
|
|
||||||
c
|
|
||||||
e
|
|
||||||
W
|
|
||||||
R
|
|
||||||
k
|
|
||||||
d
|
|
||||||
n
|
|
||||||
l
|
|
||||||
gy_obs
|
|
||||||
gp_obs
|
|
||||||
y
|
|
||||||
dA ];
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,23 +30,33 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
var e_a; stderr 0.014;
|
var e_a; stderr 0.014;
|
||||||
var e_m; stderr 0.005;
|
var e_m; stderr 0.005;
|
||||||
|
|
|
@ -30,23 +30,33 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
var e_a; stderr 0.014;
|
var e_a; stderr 0.014;
|
||||||
var e_m; stderr 0.005;
|
var e_m; stderr 0.005;
|
||||||
|
|
|
@ -30,23 +30,33 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
var e_a; stderr 0.014;
|
var e_a; stderr 0.014;
|
||||||
var e_m; stderr 0.005;
|
var e_m; stderr 0.005;
|
||||||
|
|
|
@ -30,23 +30,33 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
var e_a; stderr 0.014;
|
var e_a; stderr 0.014;
|
||||||
var e_m; stderr 0.005;
|
var e_m; stderr 0.005;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -33,21 +33,32 @@ Y_obs/Y_obs(-1) = gy_obs;
|
||||||
P_obs/P_obs(-1) = gp_obs;
|
P_obs/P_obs(-1) = gp_obs;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gy_obs = exp(gam);
|
||||||
|
gp_obs = exp(-gam);
|
||||||
|
Y_obs=gy_obs;
|
||||||
|
P_obs=gp_obs;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -1,60 +0,0 @@
|
||||||
% computes the steady state of fs2000 analyticaly
|
|
||||||
% largely inspired by the program of F. Schorfheide
|
|
||||||
function [ys,check] = fs2000a_steadystate(ys,exe)
|
|
||||||
global M_
|
|
||||||
|
|
||||||
alp = M_.params(1);
|
|
||||||
bet = M_.params(2);
|
|
||||||
gam = M_.params(3);
|
|
||||||
mst = M_.params(4);
|
|
||||||
rho = M_.params(5);
|
|
||||||
psi = M_.params(6);
|
|
||||||
del = M_.params(7);
|
|
||||||
|
|
||||||
check = 0;
|
|
||||||
|
|
||||||
dA = exp(gam);
|
|
||||||
gst = 1/dA;
|
|
||||||
m = mst;
|
|
||||||
|
|
||||||
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
|
||||||
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
|
||||||
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
|
||||||
n = xist/(nust+xist);
|
|
||||||
P = xist + nust;
|
|
||||||
k = khst*n;
|
|
||||||
|
|
||||||
l = psi*mst*n/( (1-psi)*(1-n) );
|
|
||||||
c = mst/P;
|
|
||||||
d = l - mst + 1;
|
|
||||||
y = k^alp*n^(1-alp)*gst^alp;
|
|
||||||
R = mst/bet;
|
|
||||||
W = l/n;
|
|
||||||
ist = y-c;
|
|
||||||
q = 1 - d;
|
|
||||||
|
|
||||||
e = 1;
|
|
||||||
|
|
||||||
gp_obs = m/dA;
|
|
||||||
gy_obs = dA;
|
|
||||||
|
|
||||||
P_obs = 1;
|
|
||||||
Y_obs = 1;
|
|
||||||
|
|
||||||
ys =[
|
|
||||||
m
|
|
||||||
P
|
|
||||||
c
|
|
||||||
e
|
|
||||||
W
|
|
||||||
R
|
|
||||||
k
|
|
||||||
d
|
|
||||||
n
|
|
||||||
l
|
|
||||||
gy_obs
|
|
||||||
gp_obs
|
|
||||||
Y_obs
|
|
||||||
P_obs
|
|
||||||
y
|
|
||||||
dA ];
|
|
|
@ -35,18 +35,27 @@ M-M(-1)+d = l;
|
||||||
y = k(-1)^alp*(A*n)^(1-alp);
|
y = k(-1)^alp*(A*n)^(1-alp);
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
gA = exp(gam);
|
||||||
gM = mst;
|
gst = 1/gA;
|
||||||
P = 2.25;
|
gM = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
W = 4;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
R = 1.02;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
d = 0.85;
|
n = xist/(nust+xist);
|
||||||
n = 0.19;
|
P = xist + nust;
|
||||||
l = 0.86;
|
k = khst*n;
|
||||||
y = 0.6;
|
|
||||||
gA = exp(gam);
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
|
c = mst/P;
|
||||||
|
d = l - mst + 1;
|
||||||
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -59,18 +59,27 @@ M-M(-1)+d = l;
|
||||||
y = k(-1)^alp*(A*n)^(1-alp);
|
y = k(-1)^alp*(A*n)^(1-alp);
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
gA = exp(gam);
|
||||||
gM = mst;
|
gst = 1/gA;
|
||||||
P = 2.25;
|
gM = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
W = 4;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
R = 1.02;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
d = 0.85;
|
n = xist/(nust+xist);
|
||||||
n = 0.19;
|
P = xist + nust;
|
||||||
l = 0.86;
|
k = khst*n;
|
||||||
y = 0.6;
|
|
||||||
gA = exp(gam);
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
|
c = mst/P;
|
||||||
|
d = l - mst + 1;
|
||||||
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -28,21 +28,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -28,21 +28,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -28,21 +28,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -28,21 +28,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -28,21 +28,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -1,6 +1,6 @@
|
||||||
function [ys,check1]=as2007_steadystate(junk, ys)
|
function [ys,check1]=as2007_steadystate(ys,exo)
|
||||||
|
|
||||||
global M_ options_
|
global M_
|
||||||
|
|
||||||
for j=1:size(M_.param_names,1)
|
for j=1:size(M_.param_names,1)
|
||||||
eval([deblank(M_.param_names(j,:)),' = M_.params(j);'])
|
eval([deblank(M_.param_names(j,:)),' = M_.params(j);'])
|
||||||
|
@ -21,12 +21,14 @@ YGR=gam_steady;
|
||||||
INFL = pi_steady;
|
INFL = pi_steady;
|
||||||
INT = pi_steady+rr_steady+4*gam_steady;
|
INT = pi_steady+rr_steady+4*gam_steady;
|
||||||
|
|
||||||
ys=[
|
%% end own model equations
|
||||||
pie
|
|
||||||
y
|
for iter = 1:length(M_.params) %update parameters set in the file
|
||||||
R
|
eval([ 'M_.params(' num2str(iter) ') = ' M_.param_names(iter,:) ';' ])
|
||||||
g
|
end
|
||||||
z
|
|
||||||
YGR
|
NumberOfEndogenousVariables = M_.orig_endo_nbr; %auxiliary variables are set automatically
|
||||||
INFL
|
for ii = 1:NumberOfEndogenousVariables
|
||||||
INT];
|
varname = deblank(M_.endo_names(ii,:));
|
||||||
|
eval(['ys(' int2str(ii) ') = ' varname ';']);
|
||||||
|
end
|
||||||
|
|
|
@ -1,6 +1,6 @@
|
||||||
function [ys,check1]=kim2_steadystate(junk, ys)
|
function [ys,check1]=kim2_steadystate(ys,exo)
|
||||||
|
|
||||||
global M_ options_
|
global M_
|
||||||
|
|
||||||
for j=1:size(M_.param_names,1)
|
for j=1:size(M_.param_names,1)
|
||||||
eval([deblank(M_.param_names(j,:)),' = M_.params(j);'])
|
eval([deblank(M_.param_names(j,:)),' = M_.params(j);'])
|
||||||
|
@ -19,10 +19,14 @@ i=delta*k;
|
||||||
c=(((a*k^alph)^(1+theta)-s*(i/s)^(1+theta))/(1-s))^(1/(1+theta))*(1-s);
|
c=(((a*k^alph)^(1+theta)-s*(i/s)^(1+theta))/(1-s))^(1/(1+theta))*(1-s);
|
||||||
lam = (1-s)^theta/c^(1+theta)/(1+theta);
|
lam = (1-s)^theta/c^(1+theta)/(1+theta);
|
||||||
|
|
||||||
ys=[
|
%% end own model equations
|
||||||
c
|
|
||||||
k
|
for iter = 1:length(M_.params) %update parameters set in the file
|
||||||
i
|
eval([ 'M_.params(' num2str(iter) ') = ' M_.param_names(iter,:) ';' ])
|
||||||
a
|
end
|
||||||
lam
|
|
||||||
];
|
NumberOfEndogenousVariables = M_.orig_endo_nbr; %auxiliary variables are set automatically
|
||||||
|
for ii = 1:NumberOfEndogenousVariables
|
||||||
|
varname = deblank(M_.endo_names(ii,:));
|
||||||
|
eval(['ys(' int2str(ii) ') = ' varname ';']);
|
||||||
|
end
|
||||||
|
|
|
@ -31,21 +31,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
m = mst;
|
dA = exp(gam);
|
||||||
P = 2.25;
|
gst = 1/dA;
|
||||||
c = 0.45;
|
m = mst;
|
||||||
e = 1;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
W = 4;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
R = 1.02;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
k = 6;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -31,21 +31,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
m = mst;
|
dA = exp(gam);
|
||||||
P = 2.25;
|
gst = 1/dA;
|
||||||
c = 0.45;
|
m = mst;
|
||||||
e = 1;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
W = 4;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
R = 1.02;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
k = 6;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -31,21 +31,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
m = mst;
|
dA = exp(gam);
|
||||||
P = 2.25;
|
gst = 1/dA;
|
||||||
c = 0.45;
|
m = mst;
|
||||||
e = 1;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
W = 4;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
R = 1.02;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
k = 6;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
m = mst;
|
dA = exp(gam);
|
||||||
P = 2.25;
|
gst = 1/dA;
|
||||||
c = 0.45;
|
m = mst;
|
||||||
e = 1;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
W = 4;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
R = 1.02;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
k = 6;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
m = mst;
|
dA = exp(gam);
|
||||||
P = 2.25;
|
gst = 1/dA;
|
||||||
c = 0.45;
|
m = mst;
|
||||||
e = 1;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
W = 4;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
R = 1.02;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
k = 6;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
m = mst;
|
dA = exp(gam);
|
||||||
P = 2.25;
|
gst = 1/dA;
|
||||||
c = 0.45;
|
m = mst;
|
||||||
e = 1;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
W = 4;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
R = 1.02;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
k = 6;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -34,26 +34,36 @@ P_1 = P;
|
||||||
AUXv = bet*P*(alp*exp(-alp*(gam+log(e)))*k(-1)^(alp-1)*n^(1-alp)+(1-del)*exp(-(gam+log(e))))/(c(+1)*P(+1)*m);
|
AUXv = bet*P*(alp*exp(-alp*(gam+log(e)))*k(-1)^(alp-1)*n^(1-alp)+(1-del)*exp(-(gam+log(e))))/(c(+1)*P(+1)*m);
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
m = mst;
|
dA = exp(gam);
|
||||||
m_1=mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
P_1 = 2.25;
|
m_1=mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
k = 6;
|
P = xist + nust;
|
||||||
d = 0.85;
|
P_1 = P;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
AUXv = 1;
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
|
AUXv = 1/(c*m);
|
||||||
end;
|
end;
|
||||||
|
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
var e_a; stderr 0.014;
|
var e_a; stderr 0.014;
|
||||||
var e_m; stderr 0.005;
|
var e_m; stderr 0.005;
|
||||||
|
|
|
@ -34,24 +34,33 @@ P_1 = P;
|
||||||
AUXv = bet*P*(alp*exp(-alp*(gam+log(e)))*k(-1)^(alp-1)*n^(1-alp)+(1-del)*exp(-(gam+log(e))))/(c(+1)*P(+1)*m);
|
AUXv = bet*P*(alp*exp(-alp*(gam+log(e)))*k(-1)^(alp-1)*n^(1-alp)+(1-del)*exp(-(gam+log(e))))/(c(+1)*P(+1)*m);
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
m = mst;
|
dA = exp(gam);
|
||||||
m_1=mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
P_1 = 2.25;
|
m_1=mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
k = 6;
|
P = xist + nust;
|
||||||
d = 0.85;
|
P_1 = P;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
AUXv = 1;
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
|
AUXv = 1/(c*m);
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -30,21 +30,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -31,21 +31,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -31,21 +31,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -65,21 +65,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -33,21 +33,32 @@ Y_obs/Y_obs(-1) = gy_obs;
|
||||||
P_obs/P_obs(-1) = gp_obs;
|
P_obs/P_obs(-1) = gp_obs;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
|
Y_obs = gy_obs;
|
||||||
|
P_obs = gp_obs;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -1,60 +0,0 @@
|
||||||
% computes the steady state of fs2000 analyticaly
|
|
||||||
% largely inspired by the program of F. Schorfheide
|
|
||||||
function [ys,check] = fs2000a_steadystate(ys,exe)
|
|
||||||
global M_
|
|
||||||
|
|
||||||
alp = M_.params(1);
|
|
||||||
bet = M_.params(2);
|
|
||||||
gam = M_.params(3);
|
|
||||||
mst = M_.params(4);
|
|
||||||
rho = M_.params(5);
|
|
||||||
psi = M_.params(6);
|
|
||||||
del = M_.params(7);
|
|
||||||
|
|
||||||
check = 0;
|
|
||||||
|
|
||||||
dA = exp(gam);
|
|
||||||
gst = 1/dA;
|
|
||||||
m = mst;
|
|
||||||
|
|
||||||
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
|
||||||
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
|
||||||
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
|
||||||
n = xist/(nust+xist);
|
|
||||||
P = xist + nust;
|
|
||||||
k = khst*n;
|
|
||||||
|
|
||||||
l = psi*mst*n/( (1-psi)*(1-n) );
|
|
||||||
c = mst/P;
|
|
||||||
d = l - mst + 1;
|
|
||||||
y = k^alp*n^(1-alp)*gst^alp;
|
|
||||||
R = mst/bet;
|
|
||||||
W = l/n;
|
|
||||||
ist = y-c;
|
|
||||||
q = 1 - d;
|
|
||||||
|
|
||||||
e = 1;
|
|
||||||
|
|
||||||
gp_obs = m/dA;
|
|
||||||
gy_obs = dA;
|
|
||||||
|
|
||||||
P_obs = 1;
|
|
||||||
Y_obs = 1;
|
|
||||||
|
|
||||||
ys =[
|
|
||||||
m
|
|
||||||
P
|
|
||||||
c
|
|
||||||
e
|
|
||||||
W
|
|
||||||
R
|
|
||||||
k
|
|
||||||
d
|
|
||||||
n
|
|
||||||
l
|
|
||||||
gy_obs
|
|
||||||
gp_obs
|
|
||||||
Y_obs
|
|
||||||
P_obs
|
|
||||||
y
|
|
||||||
dA ];
|
|
|
@ -63,21 +63,30 @@ gy_obs = dA*y/y(-1);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
varobs gp_obs gy_obs;
|
varobs gp_obs gy_obs;
|
||||||
|
|
|
@ -1,73 +0,0 @@
|
||||||
% computes the steady state of fs2000 analyticaly
|
|
||||||
% largely inspired by the program of F. Schorfheide
|
|
||||||
|
|
||||||
% Copyright (C) 2004-2013 Dynare Team
|
|
||||||
%
|
|
||||||
% This file is part of Dynare.
|
|
||||||
%
|
|
||||||
% Dynare is free software: you can redistribute it and/or modify
|
|
||||||
% it under the terms of the GNU General Public License as published by
|
|
||||||
% the Free Software Foundation, either version 3 of the License, or
|
|
||||||
% (at your option) any later version.
|
|
||||||
%
|
|
||||||
% Dynare is distributed in the hope that it will be useful,
|
|
||||||
% but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
||||||
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
||||||
% GNU General Public License for more details.
|
|
||||||
%
|
|
||||||
% You should have received a copy of the GNU General Public License
|
|
||||||
% along with Dynare. If not, see <http://www.gnu.org/licenses/>.
|
|
||||||
|
|
||||||
function [ys,check] = fs2000_corr_ME_steadystate(ys,exe)
|
|
||||||
global M_
|
|
||||||
|
|
||||||
alp = M_.params(1);
|
|
||||||
bet = M_.params(2);
|
|
||||||
gam = M_.params(3);
|
|
||||||
mst = M_.params(4);
|
|
||||||
rho = M_.params(5);
|
|
||||||
psi = M_.params(6);
|
|
||||||
del = M_.params(7);
|
|
||||||
|
|
||||||
check = 0;
|
|
||||||
|
|
||||||
dA = exp(gam);
|
|
||||||
gst = 1/dA;
|
|
||||||
m = mst;
|
|
||||||
|
|
||||||
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
|
||||||
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
|
||||||
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
|
||||||
n = xist/(nust+xist);
|
|
||||||
P = xist + nust;
|
|
||||||
k = khst*n;
|
|
||||||
|
|
||||||
l = psi*mst*n/( (1-psi)*(1-n) );
|
|
||||||
c = mst/P;
|
|
||||||
d = l - mst + 1;
|
|
||||||
y = k^alp*n^(1-alp)*gst^alp;
|
|
||||||
R = mst/bet;
|
|
||||||
W = l/n;
|
|
||||||
ist = y-c;
|
|
||||||
q = 1 - d;
|
|
||||||
|
|
||||||
e = 1;
|
|
||||||
|
|
||||||
gp_obs = m/dA;
|
|
||||||
gy_obs = dA;
|
|
||||||
|
|
||||||
ys =[
|
|
||||||
m
|
|
||||||
P
|
|
||||||
c
|
|
||||||
e
|
|
||||||
W
|
|
||||||
R
|
|
||||||
k
|
|
||||||
d
|
|
||||||
n
|
|
||||||
l
|
|
||||||
gy_obs
|
|
||||||
gp_obs
|
|
||||||
y
|
|
||||||
dA ];
|
|
|
@ -31,21 +31,30 @@ gy_obs = dA*y/y(-2);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -31,21 +31,30 @@ gy_obs = dA*y/y(-2);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -31,21 +31,30 @@ gy_obs = dA*y/y(-2);
|
||||||
gp_obs = (P/P(-1))*m(-1)/dA;
|
gp_obs = (P/P(-1))*m(-1)/dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
initval;
|
steady_state_model;
|
||||||
k = 6;
|
dA = exp(gam);
|
||||||
m = mst;
|
gst = 1/dA;
|
||||||
P = 2.25;
|
m = mst;
|
||||||
c = 0.45;
|
khst = ( (1-gst*bet*(1-del)) / (alp*gst^alp*bet) )^(1/(alp-1));
|
||||||
e = 1;
|
xist = ( ((khst*gst)^alp - (1-gst*(1-del))*khst)/mst )^(-1);
|
||||||
W = 4;
|
nust = psi*mst^2/( (1-alp)*(1-psi)*bet*gst^alp*khst^alp );
|
||||||
R = 1.02;
|
n = xist/(nust+xist);
|
||||||
d = 0.85;
|
P = xist + nust;
|
||||||
n = 0.19;
|
k = khst*n;
|
||||||
l = 0.86;
|
|
||||||
y = 0.6;
|
l = psi*mst*n/( (1-psi)*(1-n) );
|
||||||
gy_obs = exp(gam);
|
c = mst/P;
|
||||||
gp_obs = exp(-gam);
|
d = l - mst + 1;
|
||||||
dA = exp(gam);
|
y = k^alp*n^(1-alp)*gst^alp;
|
||||||
|
R = mst/bet;
|
||||||
|
W = l/n;
|
||||||
|
ist = y-c;
|
||||||
|
q = 1 - d;
|
||||||
|
|
||||||
|
e = 1;
|
||||||
|
|
||||||
|
gp_obs = m/dA;
|
||||||
|
gy_obs = dA;
|
||||||
end;
|
end;
|
||||||
|
|
||||||
shocks;
|
shocks;
|
||||||
|
|
|
@ -2,12 +2,12 @@ var y c k m n R pi z u;
|
||||||
varexo e sigma;
|
varexo e sigma;
|
||||||
// sigma stands for phi in the eq 2.37 p.69
|
// sigma stands for phi in the eq 2.37 p.69
|
||||||
|
|
||||||
parameters alpha beta delta gamm phi1 eta a b rho phi2 Psi thetass;
|
parameters alphha betta delta gamm phi1 eta a b rho phi2 Psi thetass;
|
||||||
//phi1 stands for capital phi in eq.2.68 and 2.69
|
//phi1 stands for capital phi in eq.2.68 and 2.69
|
||||||
//phi2 stands for lowercase phi in eq. 2.66
|
//phi2 stands for lowercase phi in eq. 2.66
|
||||||
|
|
||||||
alpha = 0.36;
|
alphha = 0.36;
|
||||||
beta = 0.989;
|
betta = 0.989;
|
||||||
gamm = 0.5;
|
gamm = 0.5;
|
||||||
delta = 0.019;
|
delta = 0.019;
|
||||||
phi1 = 2;
|
phi1 = 2;
|
||||||
|
@ -21,17 +21,17 @@ thetass = 1.0125;
|
||||||
|
|
||||||
model;
|
model;
|
||||||
|
|
||||||
(a*exp(c)^(1-b)+(1-a)*exp(m)^(1-b))^((b-phi1)/(1-b))*a*exp(c)^(-b) = (a*exp(c)^(1-b)+(1-a)*exp(m)^(1-b))^((b-phi1)/(1-b))*(1-a)*exp(m)^(-b)+beta*(a*exp(c(+1))^(1-b)+(1-a)*exp(m(+1))^(1-b))^((b-phi1)/(1-b))*a*exp(c(+1))^(-b)/(1+pi(+1));
|
(a*exp(c)^(1-b)+(1-a)*exp(m)^(1-b))^((b-phi1)/(1-b))*a*exp(c)^(-b) = (a*exp(c)^(1-b)+(1-a)*exp(m)^(1-b))^((b-phi1)/(1-b))*(1-a)*exp(m)^(-b)+betta*(a*exp(c(+1))^(1-b)+(1-a)*exp(m(+1))^(1-b))^((b-phi1)/(1-b))*a*exp(c(+1))^(-b)/(1+pi(+1));
|
||||||
|
|
||||||
Psi*(1-exp(n))^(-eta)/(a*exp(c)^(-b)*(a*exp(c)^(1-b) + (1-a)*exp(m)^(1-b))^((b-phi1)/(1-b))) = (1-alpha)*exp(y)/exp(n);
|
Psi*(1-exp(n))^(-eta)/(a*exp(c)^(-b)*(a*exp(c)^(1-b) + (1-a)*exp(m)^(1-b))^((b-phi1)/(1-b))) = (1-alphha)*exp(y)/exp(n);
|
||||||
|
|
||||||
(a*exp(c)^(1-b)+(1-a)*exp(m)^(1-b))^((b-phi1)/(1-b))*a*exp(c)^(-b) = beta*exp(R(+1))*(a*exp(c(+1))^(1-b)+(1-a)*exp(m(+1))^(1-b))^((b-phi1)/(1-b))*a*exp(c(+1))^(-b);
|
(a*exp(c)^(1-b)+(1-a)*exp(m)^(1-b))^((b-phi1)/(1-b))*a*exp(c)^(-b) = betta*exp(R(+1))*(a*exp(c(+1))^(1-b)+(1-a)*exp(m(+1))^(1-b))^((b-phi1)/(1-b))*a*exp(c(+1))^(-b);
|
||||||
|
|
||||||
exp(R) = alpha*exp(y)/exp(k(-1)) + 1-delta;
|
exp(R) = alphha*exp(y)/exp(k(-1)) + 1-delta;
|
||||||
|
|
||||||
exp(k) = (1-delta)*exp(k(-1))+exp(y)-exp(c);
|
exp(k) = (1-delta)*exp(k(-1))+exp(y)-exp(c);
|
||||||
|
|
||||||
exp(y) = exp(z)*exp(k(-1))^alpha*exp(n)^(1-alpha);
|
exp(y) = exp(z)*exp(k(-1))^alphha*exp(n)^(1-alphha);
|
||||||
|
|
||||||
exp(m) = exp(m(-1))*(u+thetass)/(1+pi);
|
exp(m) = exp(m(-1))*(u+thetass)/(1+pi);
|
||||||
|
|
||||||
|
|
|
@ -1,32 +1,28 @@
|
||||||
function [ys,check] = walsh1_old_ss_steadystate(ys0,exo)
|
function [ys,check] = walsh1_old_ss_steadystate(ys,exo)
|
||||||
global M_
|
global M_
|
||||||
|
|
||||||
check = 0;
|
% read out parameters to access them with their name
|
||||||
|
NumberOfParameters = M_.param_nbr;
|
||||||
params = M_.params;
|
for ii = 1:NumberOfParameters
|
||||||
alpha = params(1);
|
paramname = deblank(M_.param_names(ii,:));
|
||||||
beta = params(2);
|
eval([ paramname ' = M_.params(' int2str(ii) ');']);
|
||||||
delta = params(3);
|
end
|
||||||
gamm = params(4);
|
% initialize indicator
|
||||||
phi1 = params(5);
|
check = 0;
|
||||||
eta = params(6);
|
|
||||||
a = params(7);
|
|
||||||
b = params(8);
|
%% Enter model equations here
|
||||||
rho = params(9);
|
|
||||||
phi2 = params(10);
|
|
||||||
Psi = params(11);
|
|
||||||
thetass = params(12);
|
|
||||||
|
|
||||||
pi = thetass-1;
|
pi = thetass-1;
|
||||||
en = 1/3;
|
en = 1/3;
|
||||||
eR = 1/beta;
|
eR = 1/betta;
|
||||||
y_k = (1/alpha)*(1/beta-1+delta);
|
y_k = (1/alphha)*(1/betta-1+delta);
|
||||||
ek = en*y_k^(-1/(1-alpha));
|
ek = en*y_k^(-1/(1-alphha));
|
||||||
ec = ek*(y_k-delta);
|
ec = ek*(y_k-delta);
|
||||||
em = ec*(a/(1-a))^(-1/b)*((thetass-beta)/thetass)^(-1/b);
|
em = ec*(a/(1-a))^(-1/b)*((thetass-betta)/thetass)^(-1/b);
|
||||||
ey = ek*y_k;
|
ey = ek*y_k;
|
||||||
Xss = a*ec^(1-b)*(1+(a/(1-a))^(-1/b)*((thetass-beta)/thetass)^((b-1)/b));
|
Xss = a*ec^(1-b)*(1+(a/(1-a))^(-1/b)*((thetass-betta)/thetass)^((b-1)/b));
|
||||||
Psi = (1-alpha)*(ey/en)*Xss^((b-phi1)/(1-b))*a*ec^(-b)*(1-en)^eta;
|
Psi = (1-alphha)*(ey/en)*Xss^((b-phi1)/(1-b))*a*ec^(-b)*(1-en)^eta;
|
||||||
n = log(en);
|
n = log(en);
|
||||||
k = log(ek);
|
k = log(ek);
|
||||||
m = log(em);
|
m = log(em);
|
||||||
|
@ -36,6 +32,14 @@ function [ys,check] = walsh1_old_ss_steadystate(ys0,exo)
|
||||||
z = 0;
|
z = 0;
|
||||||
u = 0;
|
u = 0;
|
||||||
|
|
||||||
ys = [y c k m n R pi z u]';
|
%% end own model equations
|
||||||
M_.params(11) = Psi;
|
|
||||||
|
for iter = 1:length(M_.params) %update parameters set in the file
|
||||||
|
eval([ 'M_.params(' num2str(iter) ') = ' M_.param_names(iter,:) ';' ])
|
||||||
|
end
|
||||||
|
|
||||||
|
NumberOfEndogenousVariables = M_.orig_endo_nbr; %auxiliary variables are set automatically
|
||||||
|
for ii = 1:NumberOfEndogenousVariables
|
||||||
|
varname = deblank(M_.endo_names(ii,:));
|
||||||
|
eval(['ys(' int2str(ii) ') = ' varname ';']);
|
||||||
|
end
|
||||||
|
|
Loading…
Reference in New Issue