dynare 4.6.1 C++ used and complie LATEX file use MikTex

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

Do it best!
Renmin University of China. Economy Department.
Author: Daniel tulips liu

copyright
@ This implementation was written by Johannes Pfeifer. In case you spot mistakes,
email me at jpfeifer@gmx.de

DSGE dynare RBC model and latex file

this .mod file is complied with dynare 4.6.1 ,for old dynare if you use C++ complier than you need install a C++ complier.

but ,when someone use dynare 4.6.1, the dynare contain a gcc C++ complier ,you can set dynary like this ,

declare *** model(use_dll)***
when use this use_dll ，dynare will complie C++ file .

if use old dynare version ,for example ,dynare 4.5.7,.
Microsoft Visual studio ,

dynare RBC.mod mscv % this argment mscv tell dynare the C++ complier is Microsoft Visual studio C++ .

dynare RBC.mod mingw % mingw is tell dynare that the C++ complier is MinGW C++ complier.

Dynare file

***RBC_baseline_first_diff_bayesian.mod ***

/*
 * This file estimates a baseline RBC model with TFP and government spending shocks on simulated data. 
 * The model setup is described in Handout_RBC_model.pdf and resembles the one in King/Rebelo (1999): 
 * Resuscitating Real Business Cycles, Handbook of Macroeconomics, Volume 1 and
 * Romer (2012), Advanced macroeconomics, 4th edition
 * The driving processes are estimated as AR(1)-processes on linearly detrended data.
 *
 * This implementation was written by Johannes Pfeifer. In case you spot mistakes,
 * email me at jpfeifer@gmx.de
 *
 * Please note that the following copyright notice only applies to this Dynare 
 * implementation of the model.
 */

/*
 * Copyright (C) 2014-2016 Johannes Pfeifer
 *
 * This 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.
 *
 * It 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.
 *
 * For a copy of the GNU General Public License,
 * see <http://www.gnu.org/licenses/>.
 */

// define variables

var y           ${y}$ (long_name='output')
    c           ${c}$ (long_name='consumption')
    k           ${k}$ (long_name='capital')
    l           ${l}$ (long_name='hours')
    z           ${z}$ (long_name='TFP')
    ghat        ${\hat g}$ (long_name='government spending')
    r           ${r}$ (long_name='annualized interest rate')
    w           ${w}$ (long_name='real wage')
    invest      ${i}$ (long_name='investment') 
    log_y       ${\log(y)}$ (long_name='log output')
    log_k       ${\log(k)}$ (long_name='log capital stock')
    log_c       ${\log(c)}$ (long_name='log consumption')
    log_l       ${\log(l)}$ (long_name='log labor')
    log_w       ${\log(w)}$ (long_name='log real wage')
    log_invest  ${\log(i)}$ (long_name='log investment')
    c_obs       ${\Delta c}$ (long_name='consumption growth rate')
    g_obs       ${\Delta g}$ (long_name='government spending growth rate')
    y_obs       ${\Delta y}$ (long_name='output growth rate')
    ;

varexo eps_z ${\varepsilon_z}$ (long_name='TFP shock')
       eps_g ${\varepsilon_g}$ (long_name='government spending shock')
    ;
    
parameters 
    beta    ${\beta}$   (long_name='discount factor')
    psi     ${\psi}$    (long_name='labor disutility parameter')
    sigma   ${\sigma}$  (long_name='risk aversion')
    delta   ${\delta}$  (long_name='depreciation rate')
    alpha   ${\alpha}$  (long_name='capital share')
    rhoz    ${\rho_z}$  (long_name='persistence TFP shock')
    rhog    ${\rho_g}$  (long_name='persistence G shock')
    gammax  ${\gamma_x}$ (long_name='composite growth rate')
    gshare  ${\frac{G}{Y}}$ (long_name='government spending share')
    n       ${n}$       (long_name='population growth')
    x       ${x}$       (long_name='technology growth (per capita output growth)')
    i_y     ${\frac{I}{Y}}$ (long_name='investment-output ratio')
    k_y     ${\frac{K}{Y}}$ (long_name='capital-output ratio')
    g_ss    ${\bar G}$ (long_name='government spending in steady state')
    ;

%----------------------------------------------------------------
% set parameter values 
%----------------------------------------------------------------
sigma=1;
alpha= 0.33;
i_y=0.25;
k_y=10.4;
x=0.0055;
n=0.0027;
rhoz=0.97;
rhog=0.989;
gshare=0.2038;
/* end of set parameters values */

%----------------------------------------------------------------
% enter model equations
%----------------------------------------------------------------

model(use_dll); // this part set use_dll will genarate C++ file and complie it
[name='Euler equation']
c^(-sigma)=beta/gammax*c(+1)^(-sigma)*
    (alpha*exp(z(+1))*(k/l(+1))^(alpha-1)+(1-delta));
[name='Labor FOC']
psi*c^sigma*1/(1-l)=w;
[name='Law of motion capital'] 
gammax*k=(1-delta)*k(-1)+invest;
[name='resource constraint']
y=invest+c+g_ss*exp(ghat);
[name='production function']
y=exp(z)*k(-1)^alpha*l^(1-alpha);
[name='real wage/firm FOC labor']
w=(1-alpha)*y/l;
[name='annualized real interest rate/firm FOC capital']
r=4*alpha*y/k(-1);
[name='exogenous TFP process']
z=rhoz*z(-1)+eps_z;
[name='government spending process']
ghat=rhog*ghat(-1)+eps_g;
[name='Definition log output']
log_y = log(y);
[name='Definition log capital']
log_k = log(k);
[name='Definition log consumption']
log_c = log(c);
[name='Definition log hours']
log_l = log(l);
[name='Definition log wage']
log_w = log(w);
[name='Definition log investment']
log_invest = log(invest);
[name='Definition consumption growth']
c_obs=log_c-log_c(-1);
[name='Definition output growth']
y_obs=log_y-log_y(-1);
[name='Definition G growth']
g_obs=ghat-ghat(-1);
end;  
/*  end of model defination */

steady_state_model;
    //Do Calibration
    //calibrate the model to steady state labor of 0.33,i.e. compute the corresponding steady state values
    // and the labor disutility parameter by hand;
    gammax=(1+n)*(1+x);
    delta=i_y/k_y-x-n-n*x;
    beta=(1+x)*(1+n)/(alpha/k_y+(1-delta));
    l=0.33;
    k = ((1/beta*(1+n)*(1+x)-(1-delta))/alpha)^(1/(alpha-1))*l; 
    invest = (x+n+delta+n*x)*k;
    y=k^alpha*l^(1-alpha);
    g=gshare*y;
    g_ss=g;
    c = (1-gshare)*k^(alpha)*l^(1-alpha)-invest;
    psi=(1-alpha)*(k/l)^alpha*(1-l)/c^sigma;
    w = (1-alpha)*y/l;
    r = 4*alpha*y/k;
    log_y = log(y);
    log_k = log(k);
    log_c = log(c);
    log_l = log(l);
    log_w = log(w);
    log_invest = log(invest);
    z = 0; 
    ghat =0;
    c_obs=0;
    g_obs=0;
    y_obs=0;
end;

//set shock variances
shocks;
    var eps_z=0.0068^2;// in data
    var eps_g=0.0105^2;// in data
end;

//check the starting values for the steady state
resid;

// compute steady state given the starting values
steady;
// check Blanchard-Kahn-conditions
check;

//simulate data for estimation
stoch_simul(order=1,periods=200)  g_obs c_obs;
datatomfile('first_diff_data',[]);


estimated_params;
rhog, beta_pdf,0.7,0.1;
rhoz, beta_pdf,0.7,0.1;
stderr eps_z, inv_gamma_pdf, 0.01, 0.1;
stderr eps_g, inv_gamma_pdf, 0.01, 0.1;
end;

options_.TeX=1;

varobs g_obs c_obs;

estimation(datafile=first_diff_data,mh_jscale=1.5,mode_check,
mh_replic=40000,mh_nblocks=1,bayesian_irf,smoother,
conditional_variance_decomposition=[1 2 4]) log_y log_k log_c log_l log_w r z ghat g_obs c_obs;


%----------------------------------------------------------------
% generate LaTeX output
%----------------------------------------------------------------

write_latex_dynamic_model;
write_latex_parameter_table;
write_latex_definitions;
write_latex_prior_table;
generate_trace_plots(1);
collect_latex_files;
if system(['pdflatex -halt-on-error -interaction=batchmode ' M_.fname '_TeX_binder.tex'])
    error('TeX-File did not compile.')
end

run dynare 4.6.1 do not need set C++ complier

matlab console

matlab: dynare RBC_baseline_first_diff_bayesian.mod debug

Latex file

\documentclass[12pt]{article}
\usepackage[margin=2cm]{geometry}
\usepackage{psfrag}
\usepackage{graphicx}
\usepackage{epstopdf}
\usepackage{longtable,booktabs}
\usepackage{amsmath,amsfonts}
\usepackage{breqn}
\usepackage{float,morefloats,caption}
\begin{document}
\include{RBC_baseline_first_diff_bayesian/latex/RBC_baseline_first_diff_bayesian_latex_definitions}
\include{RBC_baseline_first_diff_bayesian/latex/RBC_baseline_first_diff_bayesian_latex_parameters}
\include{RBC_baseline_first_diff_bayesian/latex/RBC_baseline_first_diff_bayesian_priors_table}
\include{RBC_baseline_first_diff_bayesian/latex/dynamic_content}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_BayesianIRF}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_Constant}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_MCMC_inefficiency_factors}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_Posterior_Mean_1}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_Posterior_Mean_2}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_Posterior_Mode_1}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_Posterior_Mode_2}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_PriorsAndPosteriors}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_SmoothedShocks}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_SmoothedVariables}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_UpdatedVariables}
\include{RBC_baseline_first_diff_bayesian/Output/RBC_baseline_first_diff_bayesian_geweke}
\include{RBC_baseline_first_diff_bayesian/graphs/RBC_baseline_first_diff_bayesian_CheckPlots}
\include{RBC_baseline_first_diff_bayesian/graphs/RBC_baseline_first_diff_bayesian_Priors}
\include{RBC_baseline_first_diff_bayesian/graphs/RBC_baseline_first_diff_bayesian_TracePlot_Posterior_blck_1}
\include{RBC_baseline_first_diff_bayesian/graphs/RBC_baseline_first_diff_bayesian_TracePlot_SE_eps_g_blck_1}
\include{RBC_baseline_first_diff_bayesian/graphs/RBC_baseline_first_diff_bayesian_TracePlot_SE_eps_z_blck_1}
\include{RBC_baseline_first_diff_bayesian/graphs/RBC_baseline_first_diff_bayesian_TracePlot_rhog_blck_1}
\include{RBC_baseline_first_diff_bayesian/graphs/RBC_baseline_first_diff_bayesian_TracePlot_rhoz_blck_1}
\end{document}

complie latex file in matlab IDE

install MikTex on your compuet

use command in matlab consose

!pdflatex RBC_baseline_first_diff_bayesian_TeX_binder.tex % use MikTex complie .tex file

Enjoy

yours. Daniel tulips liu . Renmin University of China.

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：MikTex dynare LaTeX comp file

parameters
    beta    ${\beta}$   (long_name='discount factor')
    psi     ${\psi}$    (long_name='labor disutility parameter')
    sigma   ${\sigma}$  (long_name='risk aversion')
    delta   ${\delta}$  (long_name='depreciation rate')
    alpha   ${\alpha}$  (long_name='capital share')
    rhoz    ${\rho_z}$  (long_name='persistence TFP shock')
    rhog    ${\rho_g}$  (long_name='persistence G shock')
    gammax  ${\gamma_x}$ (long_name='composite growth rate')
    gshare  ${\frac{G}{Y}}$ (long_name='government spending share')
    n       ${n}$       (long_name='population growth')
    x       ${x}$       (long_name='technology growth (per capita output growth)')
    i_y     ${\frac{I}{Y}}$ (long_name='investment-output ratio')
    k_y     ${\frac{K}{Y}}$ (long_name='capital-output ratio')
    g_ss    ${\bar G}$ (long_name='government spending in steady state')
    ;

There are something wrong when running to the code ‘write_latex_dynamic_model’ since we can't find the write_latex_dynamic_model.m file in the package of dynare 4.6.4 version.

Computer;

Console ,  控制台。