stepan
/
quest-iii-estimation
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Remove switch parameter DDYN.

master
Stéphane Adjemian (Argos) 2023-10-12 18:03:10 +02:00
parent 3ce318f4bb
commit c378472dca
Signed by: stepan
GPG Key ID: A6D44CB9C64CE77B
2 changed files with 5 additions and 19 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Q3.mod
View File

@ -20,7 +20,7 @@ varexo
E_EPS_C E_EPS_ETA E_EPS_ETAM E_EPS_ETAX E_EPS_EX E_EPS_IG E_EPS_INOMW E_EPS_L E_EPS_LOL E_EPS_M E_EPS_PPI E_EPS_PW E_EPS_RPREME E_EPS_RPREMK E_EPS_TR E_EPS_W E_EPS_Y E_EPS_YW E_EPS_G;
parameters
A1E A2E ALPHAX ALPHAE ALPHAGE BETAE BGADJ1 BGADJ2 BGTAR DELTAE DELTAGE DGEX DGIM DGPM DGPX DDYN E_EX_INOMW E_EX_R E_EX_RW G1E GAMI2E GAMIE GAMLE GAMPE GAMPME GAMPXE GAMWE GP0 GPCPI0 GPOP0 GPW0 GSLAG GVECM
A1E A2E ALPHAX ALPHAE ALPHAGE BETAE BGADJ1 BGADJ2 BGTAR DELTAE DELTAGE DGEX DGIM DGPM DGPX E_EX_INOMW E_EX_R E_EX_RW G1E GAMI2E GAMIE GAMLE GAMPE GAMPME GAMPXE GAMWE GP0 GPCPI0 GPOP0 GPW0 GSLAG GVECM
GFLAG GLAGFLAG IGLAGFLAG IGFLAG TRFLAG TWFLAG GEXOFLAG IGEXOFLAG TREXOFLAG GAMIFLAG SLCFLAG GSN GTFP0 GY0 GYW0 HABE HABLE IGSLAG IGVECM ILAGE INFLAGE IG1E IGSN ISN KAPPAE L0 LOL LYWY0 OMEGE RHOCE RHOETA RHOETAM RHOETAX RHOEXE RHOGE RHOIG RHOL0 RHOLE RHOLOL RHOPPI1 RHOPPI2 RHOPPI3 RHOPPI4 RHOPCPM RHOPWPX RHORPE RHORPK RHOUCAP0 RII RIP RIX RPI RPP RPX RXI RXP RXX RXY RPREME RPREMK SE SFPE SFPME SFPXE SFWE SIGC SIGEXE SIGIME SLC SSC TAUE TP THETAE TINFE TR1E TRSN RHOTR TYE1 TYE2 TVAT TW0 TW1 UCAP0 WRLAG ZETE interestq_exog inflationannual_exog;
//estimated parameters (mean posterior distribution)
@ -103,7 +103,6 @@ TW1 = 0.8;
ZETE = 0.4000;
// SWITCHES for various simulations
DDYN = 1;
GFLAG=1;
GLAGFLAG=1;
IGLAGFLAG=1;
@ -121,7 +120,6 @@ A1E = 0.0669; % this is actually a function of the steady state and es
OMEGE = 1.4836; % this is actually a function of the steady state and estimated parameters to get L0=0.65;
GSN = 0.203;
IGSN = 0.025;
//ISN = 0.17; % used only when DDYN=0, and RPREMK would be a function of ISN;
GPCPI0 = 0;
GP0 = 0.005;
GPW0 = GP0;

20
Q3_steadystate.m
View File

@ -6,9 +6,6 @@ for j=1:size(M_.param_names,1)
end
check1=0;
if isempty(DDYN), DDYN=1; end
E_EX_RW = E_EX_INOMW - GPW0;
E_EX_R = 1/BETAE-1;
@ -67,15 +64,10 @@ YWY = exp(LYWY0);
E_LIK = log(DELTAKE+GY0+GPCPI0);
E_LIKG = log(DELTAKGE+GY0+GPCPI0);
IK = (DELTAKE+GY0+GPCPI0);
if DDYN==0,
E_LYKPPI = E_LIK-log(ISN);
YKN = IK/ISN;
else
YKN = (E_EX_R+RPREMK+DELTAE)/(1-TP)/(1-TAUE)/(1-ALPHAE);
E_LYKPPI=log(YKN);
ISN=IK/YKN;
M_.params(strcmp('ISN', M_.param_names)) = ISN;
end
YKN = (E_EX_R+RPREMK+DELTAE)/(1-TP)/(1-TAUE)/(1-ALPHAE);
E_LYKPPI=log(YKN);
ISN=IK/YKN;
M_.params(strcmp('ISN', M_.param_names)) = ISN;
E_L =L0;
@ -206,9 +198,5 @@ A1E = (1-TAUE)*(1-ALPHAE)*exp(E_LYKPPI)/UCAP0;
M_.params(strcmp('A1E', M_.param_names)) = A1E;
M_.params(strcmp('OMEGE', M_.param_names)) = OMEGE;
if DDYN==0
RPREMK = (1-TAUE)*(1-ALPHAE)*(1-TP)*exp(E_LYKPPI) - DELTAE-E_EX_R;
M_.params(strcmp('RPREMK', M_.param_names)) = RPREMK;
end
params = M_.params;