stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

GSA: cleanup of doc and inclusion in build system

time-shift
Sébastien Villemot 2011-03-17 15:48:31 +01:00
parent ebe28543fe
commit 05ebb51102
15 changed files with 21 additions and 9086 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
configure.ac
View File

@ -174,6 +174,7 @@ AC_CONFIG_FILES([Makefile
doc/userguide/Makefile
doc/parallel/Makefile
doc/internals/Makefile
doc/gsa/Makefile
tests/Makefile
matlab/dynare_version.m
windows/dynare-version.nsi

2
doc/Makefile.am
View File

@ -1,4 +1,4 @@
SUBDIRS = preprocessor macroprocessor userguide parallel internals
SUBDIRS = preprocessor macroprocessor userguide parallel internals gsa
PDF_TARGETS =

18
doc/gsa/Makefile.am Normal file
View File

@ -0,0 +1,18 @@
if HAVE_PDFLATEX
if HAVE_BIBTEX
pdf-local: gsa.pdf
endif
endif
SRC = gsa.tex marco.bib
EXTRA_DIST = $(SRC)
gsa.pdf: $(SRC)
$(PDFLATEX) gsa
$(BIBTEX) gsa
$(PDFLATEX) gsa
$(PDFLATEX) gsa
clean-local:
rm -f *.pdf *.log *.aux *.toc *.lof *.blg *.bbl *.out *~

24
doc/gsa/dyn_tbx.aux
View File

@ -1,24 +0,0 @@
\relax
\citation{Ratto_CompEcon_2008}
\@writefile{toc}{\contentsline {section}{\numberline {1}Introduction}{2}}
\newlabel{s:intro}{{1}{2}}
\@writefile{toc}{\contentsline {section}{\numberline {2}Use of the Toolbox}{2}}
\@writefile{toc}{\contentsline {section}{\numberline {3}List of options}{3}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.1}Sampling options}{3}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.2}Stability mapping}{3}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.3}Reduced form mapping}{4}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.4}Mapping the fit}{5}}
\citation{Iskrev2009}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.5}Screening analysis}{6}}
\@writefile{toc}{\contentsline {subsection}{\numberline {3.6}Identification analysis (under development)}{6}}
\citation{Iskrev2009}
\@writefile{toc}{\contentsline {section}{\numberline {4}Directory structure}{8}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.1}Binary data files}{8}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.2}Stability analysis}{8}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.3}RMSE analysis}{9}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.4}Reduced form mapping}{10}}
\bibstyle{plainnat}
\bibdata{marco}
\bibcite{Iskrev2009}{{1}{2009}{{Iskrev}}{{}}}
\@writefile{toc}{\contentsline {subsection}{\numberline {4.5}Screening analysis}{11}}
\bibcite{Ratto_CompEcon_2008}{{2}{2008}{{Ratto}}{{}}}

18
doc/gsa/dyn_tbx.bbl
View File

@ -1,18 +0,0 @@
\begin{thebibliography}{2}
\providecommand{\natexlab}[1]{#1}
\providecommand{\url}[1]{\texttt{#1}}
\expandafter\ifx\csname urlstyle\endcsname\relax
\providecommand{\doi}[1]{doi: #1}\else
\providecommand{\doi}{doi: \begingroup \urlstyle{rm}\Url}\fi
\bibitem[Iskrev(2009)]{Iskrev2009}
Nikolay Iskrev.
\newblock Local identification in \textsc{DSGE} models.
\newblock mimeo, 2009.
\bibitem[Ratto(2008)]{Ratto_CompEcon_2008}
M.~Ratto.
\newblock Analysing dsge models with global sensitivity analysis.
\newblock \emph{Computational Economics}, 31:\penalty0 115--139, 2008.
\end{thebibliography}

30
doc/gsa/dyn_tbx.blg
View File

@ -1,30 +0,0 @@
This is BibTeX, Version 0.99cThe top-level auxiliary file: dyn_tbx.aux
The style file: plainnat.bst
Database file #1: marco.bib
I was expecting a `,' or a `}'---line 921 of file marco.bib
:
: pages = {},
(Error may have been on previous line)
I'm skipping whatever remains of this entry
I was expecting a `,' or a `}'---line 1177 of file marco.bib
:
: organisation = {Los Alamos National Laboratories}
(Error may have been on previous line)
I'm skipping whatever remains of this entry
I was expecting a `,' or a `}'---line 1364 of file marco.bib
:
: organisation = {}
(Error may have been on previous line)
I'm skipping whatever remains of this entry
You're missing a field name---line 2210 of file marco.bib
:
: %@ARTICLE{kn:NIH,
(Error may have been on previous line)
I'm skipping whatever remains of this entry
You're missing a field name---line 2211 of file marco.bib
:
: % author = {NIH},
(Error may have been on previous line)
I'm skipping whatever remains of this entry
Warning--empty institution in Iskrev2009
(There were 5 error messages)

BIN
doc/gsa/dyn_tbx.dvi
View File

Binary file not shown.

284
doc/gsa/dyn_tbx.log
View File

@ -1,284 +0,0 @@
This is pdfeTeX, Version 3.141592-1.30.6-2.2 (MiKTeX 2.5) (preloaded format=latex 2006.9.25) 2 OCT 2009 16:32
entering extended mode
**C:/SVN/dyn_util/latex/dyn_tbx.tex
(C:\SVN\dyn_util\latex\dyn_tbx.tex
LaTeX2e <2005/12/01>
Babel <v3.8g> and hyphenation patterns for english, dumylang, nohyphenation, ge
rman, ngerman, french, loaded.
("C:\Program Files\MiKTeX 2.5\tex\latex\base\article.cls"
Document Class: article 2005/09/16 v1.4f Standard LaTeX document class
("C:\Program Files\MiKTeX 2.5\tex\latex\base\size12.clo"
File: size12.clo 2005/09/16 v1.4f Standard LaTeX file (size option)
)
\c@part=\count79
\c@section=\count80
\c@subsection=\count81
\c@subsubsection=\count82
\c@paragraph=\count83
\c@subparagraph=\count84
\c@figure=\count85
\c@table=\count86
\abovecaptionskip=\skip41
\belowcaptionskip=\skip42
\bibindent=\dimen102
)
("C:\Program Files\MiKTeX 2.5\tex\latex\amsfonts\amssymb.sty"
Package: amssymb 2002/01/22 v2.2d
("C:\Program Files\MiKTeX 2.5\tex\latex\amsfonts\amsfonts.sty"
Package: amsfonts 2001/10/25 v2.2f
\@emptytoks=\toks14
\symAMSa=\mathgroup4
\symAMSb=\mathgroup5
LaTeX Font Info: Overwriting math alphabet `\mathfrak' in version `bold'
(Font) U/euf/m/n --> U/euf/b/n on input line 132.
))
("C:\Program Files\MiKTeX 2.5\tex\latex\amsmath\amsmath.sty"
Package: amsmath 2000/07/18 v2.13 AMS math features
\@mathmargin=\skip43
For additional information on amsmath, use the `?' option.
("C:\Program Files\MiKTeX 2.5\tex\latex\amsmath\amstext.sty"
Package: amstext 2000/06/29 v2.01
("C:\Program Files\MiKTeX 2.5\tex\latex\amsmath\amsgen.sty"
File: amsgen.sty 1999/11/30 v2.0
\@emptytoks=\toks15
\ex@=\dimen103
))
("C:\Program Files\MiKTeX 2.5\tex\latex\amsmath\amsbsy.sty"
Package: amsbsy 1999/11/29 v1.2d
\pmbraise@=\dimen104
)
("C:\Program Files\MiKTeX 2.5\tex\latex\amsmath\amsopn.sty"
Package: amsopn 1999/12/14 v2.01 operator names
)
\inf@bad=\count87
LaTeX Info: Redefining \frac on input line 211.
\uproot@=\count88
\leftroot@=\count89
LaTeX Info: Redefining \overline on input line 307.
\classnum@=\count90
\DOTSCASE@=\count91
LaTeX Info: Redefining \ldots on input line 379.
LaTeX Info: Redefining \dots on input line 382.
LaTeX Info: Redefining \cdots on input line 467.
\Mathstrutbox@=\box26
\strutbox@=\box27
\big@size=\dimen105
LaTeX Font Info: Redeclaring font encoding OML on input line 567.
LaTeX Font Info: Redeclaring font encoding OMS on input line 568.
\macc@depth=\count92
\c@MaxMatrixCols=\count93
\dotsspace@=\muskip10
\c@parentequation=\count94
\dspbrk@lvl=\count95
\tag@help=\toks16
\row@=\count96
\column@=\count97
\maxfields@=\count98
\andhelp@=\toks17
\eqnshift@=\dimen106
\alignsep@=\dimen107
\tagshift@=\dimen108
\tagwidth@=\dimen109
\totwidth@=\dimen110
\lineht@=\dimen111
\@envbody=\toks18
\multlinegap=\skip44
\multlinetaggap=\skip45
\mathdisplay@stack=\toks19
LaTeX Info: Redefining \[ on input line 2666.
LaTeX Info: Redefining \] on input line 2667.
)
("C:\Program Files\MiKTeX 2.5\tex\latex\graphics\graphicx.sty"
Package: graphicx 1999/02/16 v1.0f Enhanced LaTeX Graphics (DPC,SPQR)
("C:\Program Files\MiKTeX 2.5\tex\latex\graphics\keyval.sty"
Package: keyval 1999/03/16 v1.13 key=value parser (DPC)
\KV@toks@=\toks20
)
("C:\Program Files\MiKTeX 2.5\tex\latex\graphics\graphics.sty"
Package: graphics 2006/02/20 v1.0o Standard LaTeX Graphics (DPC,SPQR)
("C:\Program Files\MiKTeX 2.5\tex\latex\graphics\trig.sty"
Package: trig 1999/03/16 v1.09 sin cos tan (DPC)
)
("C:\Program Files\MiKTeX 2.5\tex\latex\00miktex\graphics.cfg"
File: graphics.cfg 2005/12/29 v1.2 MiKTeX 'graphics' configuration
)
Package graphics Info: Driver file: dvips.def on input line 90.
("C:\Program Files\MiKTeX 2.5\tex\latex\graphics\dvips.def"
File: dvips.def 1999/02/16 v3.0i Driver-dependant file (DPC,SPQR)
))
\Gin@req@height=\dimen112
\Gin@req@width=\dimen113
)
("C:\Program Files\MiKTeX 2.5\tex\latex\natbib\natbib.sty"
Package: natbib 2006/08/18 7.4 (PWD)
\bibhang=\skip46
\bibsep=\skip47
LaTeX Info: Redefining \cite on input line 541.
\c@NAT@ctr=\count99
)
("C:\Program Files\MiKTeX 2.5\tex\latex\psfrag\psfrag.sty"
Package: psfrag 1998/04/11 v3.04 PSfrag (MCG)
\pfg@pcount=\count100
\pfg@scount=\count101
\pfg@temp=\write3
)
("C:\Program Files\MiKTeX 2.5\tex\latex\setspace\setspace.sty"
Package: setspace 2000/12/01 6.7 Contributed and Supported LaTeX2e package
Package: `setspace' 6.7 <2000/12/01>
) ("C:\Program Files\MiKTeX 2.5\tex\latex\rotating\rotating.sty"
Package: rotating 1997/09/26 v2.13 Rotation package
("C:\Program Files\MiKTeX 2.5\tex\latex\base\ifthen.sty"
Package: ifthen 2001/05/26 v1.1c Standard LaTeX ifthen package (DPC)
)
\c@r@tfl@t=\count102
\rot@float@box=\box28
) (dyn_tbx.aux)
LaTeX Font Info: Checking defaults for OML/cmm/m/it on input line 20.
LaTeX Font Info: ... okay on input line 20.
LaTeX Font Info: Checking defaults for T1/cmr/m/n on input line 20.
LaTeX Font Info: ... okay on input line 20.
LaTeX Font Info: Checking defaults for OT1/cmr/m/n on input line 20.
LaTeX Font Info: ... okay on input line 20.
LaTeX Font Info: Checking defaults for OMS/cmsy/m/n on input line 20.
LaTeX Font Info: ... okay on input line 20.
LaTeX Font Info: Checking defaults for OMX/cmex/m/n on input line 20.
LaTeX Font Info: ... okay on input line 20.
LaTeX Font Info: Checking defaults for U/cmr/m/n on input line 20.
LaTeX Font Info: ... okay on input line 20.
LaTeX Font Info: Try loading font information for U+msa on input line 41.
("C:\Program Files\MiKTeX 2.5\tex\latex\amsfonts\umsa.fd"
File: umsa.fd 2002/01/19 v2.2g AMS font definitions
)
LaTeX Font Info: Try loading font information for U+msb on input line 41.
("C:\Program Files\MiKTeX 2.5\tex\latex\amsfonts\umsb.fd"
File: umsb.fd 2002/01/19 v2.2g AMS font definitions
) [1
]
Overfull \hbox (4.3pt too wide) in paragraph at lines 80--81
[][][]\OT1/cmtt/m/n/12 dynare_sensitivity(option1=<opt1_val>,option2=<opt2_val>
,...)
[]
[2]
Overfull \hbox (7.07413pt too wide) in paragraph at lines 121--133
[][]
[]
[3]
Underfull \hbox (badness 10000) in paragraph at lines 136--141
[]
Overfull \hbox (26.50142pt too wide) in paragraph at lines 147--177
[][][]
[]
Underfull \hbox (badness 10000) in paragraph at lines 147--177
[]
Underfull \hbox (badness 10000) in paragraph at lines 178--180
[]
[4]
LaTeX Font Info: Try loading font information for OMS+cmr on input line 204.
("C:\Program Files\MiKTeX 2.5\tex\latex\base\omscmr.fd"
File: omscmr.fd 1999/05/25 v2.5h Standard LaTeX font definitions
)
LaTeX Font Info: Font shape `OMS/cmr/m/n' in size <12> not available
(Font) Font shape `OMS/cmsy/m/n' tried instead on input line 204.
Overfull \hbox (26.51085pt too wide) in paragraph at lines 216--236
[][][]
[]
[5]
Overfull \hbox (44.59268pt too wide) in paragraph at lines 255--270
[][][]
[]
[6] [7]
Underfull \hbox (badness 10000) in paragraph at lines 291--299
[]
Underfull \hbox (badness 10000) in paragraph at lines 291--299
[]
[8]
Overfull \hbox (12.08382pt too wide) in paragraph at lines 363--365
[][][]\OT1/cmtt/m/n/12 <DYNARE_file>_rmse_post*.fig\OT1/cmr/m/n/12 : save re-su
lts for the anal-y-sis us-ing Metropo-
[]
Overfull \hbox (8.24295pt too wide) in paragraph at lines 375--378
[][][]\OT1/cmtt/m/n/12 <DYNARE_file>_rmse_prior_<name of observedseries>_corr_*
.fig\OT1/cmr/m/n/12 : for
[]
[9]
Overfull \hbox (5.72563pt too wide) in paragraph at lines 384--389
\OT1/cmr/m/it/12 in red \OT1/cmr/m/n/12 the cdf of the log-posterior cor-re-spo
nd-ing to the best 10% RMSE's,
[]
Overfull \hbox (0.88243pt too wide) in paragraph at lines 394--397
[][][]\OT1/cmtt/m/n/12 <DYNARE_file>_rmse_prior_lik_SA_*.fig\OT1/cmr/m/n/12 : w
hen []\OT1/cmtt/m/n/12 lik_only=1\OT1/cmr/m/n/12 , this shows
[]
Overfull \hbox (7.05743pt too wide) in paragraph at lines 397--400
[][][]\OT1/cmtt/m/n/12 <DYNARE_file>_rmse_prior_post_SA_*.fig\OT1/cmr/m/n/12 :
when []\OT1/cmtt/m/n/12 lik_only=1\OT1/cmr/m/n/12 , this shows
[]
Overfull \hbox (9.38109pt too wide) in paragraph at lines 413--419
[][][]\OT1/cmtt/m/n/12 <DYNARE_file>_redform_<endo name>_vs_shocks_*.fig\OT1/cm
r/m/n/12 : shows bar charts
[]
Overfull \hbox (26.57321pt too wide) in paragraph at lines 425--430
[][]\OT1/cmr/m/n/12 Detailed re-sults of the anal-y-ses are shown in the sub-fo
lder []\OT1/cmtt/m/n/12 \GSA\redform_stab\OT1/cmr/m/n/12 ,
[]
[10]
Overfull \hbox (9.38109pt too wide) in paragraph at lines 463--468
[][][]\OT1/cmtt/m/n/12 <DYNARE_file>_redform_<endo name>_vs_shocks_*.fig\OT1/cm
r/m/n/12 : shows bar charts
[]
(dyn_tbx.bbl [11]) [12] (dyn_tbx.aux) )
Here is how much of TeX's memory you used:
1989 strings out of 95420
21421 string characters out of 1187324
74420 words of memory out of 1073000
5051 multiletter control sequences out of 60000
14651 words of font info for 57 fonts, out of 1000000 for 2000
14 hyphenation exceptions out of 8191
29i,13n,24p,335b,294s stack positions out of 5000i,500n,10000p,200000b,32768s
Output written on dyn_tbx.dvi (12 pages, 31640 bytes).

BIN
doc/gsa/dyn_tbx.pdf
View File

Binary file not shown.

3707
doc/gsa/dyn_tbx.ps
View File

File diff suppressed because it is too large Load Diff

486
doc/gsa/dyn_tbx.tex.bak
View File

@ -1,486 +0,0 @@
% ----------------------------------------------------------------
% AMS-LaTeX Paper ************************************************
% **** -----------------------------------------------------------
\documentclass[12pt,a4paper]{article}
\usepackage{amssymb,amsmath}
\usepackage[dvips]{graphicx}
\usepackage{natbib}
\usepackage{psfrag}
\usepackage{setspace}
\usepackage{rotating}
%\singlespacing (interlinea singola)
%\onehalfspacing (interlinea 1,5)
%\doublespacing (interlinea doppia)
%\bibpunct{(}{)}{;}{a}{,}{,}
\bibpunct[, ]{(}{)}{;}{a}{,}{,}
%\pagestyle{headings}
% ----------------------------------------------------------------
\begin{document}
% ----------------------------------------------------------------
\title{Sensitivity Analysis Toolbox for DYNARE}%
\author{Marco Ratto\\
European Commission, Joint Research Centre \\
TP361, IPSC, via E. Fermi 1\\21020 Ispra
(VA) Italy\\
\texttt{marco.ratto@jrc.it}
\thanks{The author gratefully thanks Christophe Planas, Kenneth Judd, Michel Juillard,
Alessandro Rossi, Frank Schorfheide and the participants to the
Courses on Global Sensitivity Analysis for Macroeconomic
Models (Ispra, 2006-2007-2008) for interesting discussions and
helpful suggestions.}}
%%% To have the current date inserted, use \date{\today}:
%%% To insert a footnote, add thanks in the date/title/author fields:
\date{\today}
%\date{\today \thanks{Authors gratefully acknowledge the
%contribution by ... for ...}}
\maketitle %\tableofcontents
%\doublespacing
%-----------------------------------------------------------------------
\begin{abstract}
\noindent The Sensitivity Analysis Toolbox for DYNARE is a set of
MATLAB routines for the analysis of DSGE models with global
sensitivity analysis. The routines are thought to be used within
the DYNARE v4 environment.
\begin{description}
\item \textbf{Keywords}: Stability Mapping , Reduced form solution, DSGE models,
Monte Carlo filtering, Global Sensitivity Analysis.
\end{description}
\end{abstract}
\newpage
% ----------------------------------------------------------------
\section{Introduction} \label{s:intro}
The Sensitivity Analysis Toolbox for DYNARE is a collection of
MATLAB routines implemented to answer the following questions: (i)
Which is the domain of structural coefficients assuring the
stability and determinacy of a DSGE model? (ii) Which parameters
mostly drive the fit of, e.g., GDP and which the fit of inflation?
Is there any conflict between the optimal fit of one observed
series versus another one? (iii) How to represent in a direct,
albeit approximated, form the relationship between structural
parameters and the reduced form of a rational expectations model?
The discussion of the methodologies and their application is
described in \cite{Ratto_CompEcon_2008}.
\section{Use of the Toolbox}
The DYNARE parser now recognizes sensitivity analysis commands.
The syntax is based on a single command:
\vspace{0.5cm}
\verb"dynare_sensitivity(option1=<opt1_val>,option2=<opt2_val>,...)"
\vspace{0.5cm} \noindent with a list of options described in the
next section.
With respect to the previous version of the toolbox, in order to
work properly, the sensitivity analysis Toolbox \emph{no longer}
needs that the DYNARE estimation environment is set-up.
Therefore, \verb"dynare_sensitivity" is the only command to run to
make a sensitivity analysis on a DSGE model\footnote{Of course,
when the user needs to perform the mapping of the fit with a
posterior sample, a Bayesian estimation has to be performed
beforehand}.
\section{List of options}
\subsection{Sampling options}
\begin{tabular}{r|l|l}
% after \\ : \hline or \cline{col1-col2} \cline{col3-col4} ...
option name & default & description \\ \hline
\verb"Nsam"& 2048& Size of MC sample \\
\verb"ilptau"& 1& 1 = use $LP_\tau$ quasi-Monte Carlo \\
& & 0 = use LHS Monte Carlo \\
\verb"pprior"& 1& 1 = sample from prior distributions\\
& & 0 = sample from multivariate
normal \\
& & \hspace{0.5 cm} $N(\hat{\theta},\Sigma)$, $\hat{\theta}$ is posterior mode \\
& & \hspace{0.5 cm} $\Sigma = H^{-1}$, $H$ is Hessian at the mode\\
\verb"prior_range"& 1& 1 = sample \textit{uniformly} from prior ranges\\
& & 0 = sample from prior distributions: \\
& & this requires MATLAB Statistics Toolbox\\
\verb"morris"& 0& 0 = no Morris sampling for screening \\
& & 1 = Morris sampling for screening \\
\verb"morris_nliv"& 6& number of levels in Morris design\\
\verb"morris_ntra"& 20& number of trajectories in Morris design\\
\verb"ppost"& 0& 0 = don't use Metropolis posterior sample\\
& & 1 = use Metropolis posterior sample: this \\
& & \hspace{0.5 cm} overrides any other sampling option! \\ \hline
\end{tabular}
\subsection{Stability mapping}
\begin{tabular}{r|l|l}
option name & default & description \\ \hline
\verb"stab"& 1& 1 = perform stability mapping \\
& & 0 = no stability mapping is performed\\
\verb"load_stab"& 0& 0 = generate a new sample\\
& & 1 = load a previously created sample \\
\verb"alpha2_stab"& 0.4& critical value for correlations $\rho$ in filtered samples:\\
& & plot couples of parameters with \\
& & $|\rho|>$\verb"alpha2_stab"\\
\verb"ksstat"& 0.1& critical value for Smirnov statistics $d$: \\
& & plot parameters with $d>$\verb"ksstat"\\ \hline
\end{tabular}
\newpage
\subsection{Reduced form mapping}% and identification}
The mapping of the reduced form solution forces the use of samples
from prior ranges or prior distributions, i.e.:
\\
\verb"options_.opt_gsa.pprior=1;"\\
\verb"options_.opt_gsa.ppost=0;"\\
It uses 250 samples to optimize smoothing parameters and 1000
samples to compute the fit. The rest of the sample is used for
out-of-sample validation. \vspace{0.5cm}
\begin{tabular}{r|l|l}
option name & default & description \\ \hline
\verb"redform"& 0& 0 = don't prepare MC sample of \\
& & reduced form matrices \\
& & 1 = prepare MC sample of \\
& & reduced form matrices \\
\verb"load_redform"& 0& 0 = estimate the mapping of \\
& & reduced form model\\
& & 1 = load previously estimated mapping\\
\verb"logtrans_redform"& 0& 0 = use raw entries\\
& & 1 = use log-transformed entries \\
\verb"threshold_redform"& []& [] = don't filter MC entries \\
& & of reduced form coefficients\\
& & [\verb"max" \verb"max"] = analyse filtered \\
& & entries within the range [\verb"max" \verb"max"] \\
\verb"ksstat_redform"& 0.1& critical value for Smirnov statistics $d$ \\
& & when reduced form entries are filtered\\
\verb"alpha2_redform"& 0.3& critical value for correlation $\rho$ \\
& & when reduced form entries are filtered\\
\verb"namendo"& () & list of endogenous variables \\
& : & jolly character to indicate ALL endogenous \\
\verb"namlagendo"& () & list of lagged endogenous variables:\\
& & analyse entries [\verb"namendo"$\times$\verb"namlagendo"]\\
& : & jolly character to indicate ALL endogenous \\
\verb"namexo"& ()& list of exogenous variables:\\
& & analyse entries
[\verb"namendo"$\times$\verb"namexo"]\\
& : & jolly character to indicate ALL exogenous \\\hline
\end{tabular}
\vspace{0.5cm} \\
One can also load a previously estimated mapping
with a new MC sample, to look at the forecast for the new MC sample.\\
\subsection{Mapping the fit}
The RMSE analysis can be performed with different types of
sampling options:
\begin{enumerate}
\item when \verb"pprior=1" and \verb"ppost=0", the Toolbox analyses the RMSE's for
the MC sample obtained by sampling parameters from their prior
distributions (or prior ranges): this analysis provides some hints
about what parameter drives the fit of which observed series,
prior to the full estimation;
\item when \verb"pprior=0" and \verb"ppost=0", the Toolbox analyses the RMSE's for
a multivariate normal MC sample, with covariance matrix based on
the inverse Hessian at the optimum: this analysis is useful when
ML estimation is done (i.e. no Bayesian estimation);
\item when \verb"ppost=1" the Toolbox analyses
the RMSE's for the posterior sample obtained by DYNARE's
Metropolis procedure.
\end{enumerate}
The use of cases 2. and 3. requires an estimation step beforehand!
To facilitate the sensitivity analysis after estimation, the
\verb"dynare_sensitivity" command also allows to indicate some
options of \verb"dynare_estimation". These are:
\begin{itemize}
\item \verb"datafile"
\item \verb"mode_file"
\item \verb"first_obs"
\item \verb"nobs"
\item \verb"prefilter"
\item \verb"presample"
\item \verb"loglinear"
\end{itemize}
\vspace{1cm}
\begin{tabular}{r|l|l}
option name & default & description \\ \hline
\verb"rmse"& 0& 0 = no RMSE analysis\\
& & 1 = do RMSE analysis \\
\verb"load_rmse"& 0& 0 = make a new RMSE analysis\\
& & 1 = load previous RMSE analysis \\
\verb"lik_only"& 0& 0 = compute RMSE's for all observed series\\
& & 1 = compute only likelihood and posterior \\
\verb"var_rmse"& varobs& list of observed series to be considered\\
\verb"pfilt_rmse"& 0.1& filtering threshold for RMSE's: default it to\\
& & filter the best 10\% for each observed series\\
\verb"istart_rmse"& 1& start computing RMSE's from \verb"istart_rmse":\\
& & use 2 to avoid big initial error \\
\verb"alpha_rmse"& 0.002& critical value for Smirnov statistics $d$:\\
& & plot parameters with $d>$\verb"alpha_rmse"\\
\verb"alpha2_rmse"& 1& critical value for correlation $\rho$\\
& & plot couples of parameters with
$|\rho|>$\verb"alpha2_rmse"\\
\verb"glue"& 0& prepare for GLUE graphical interface\\\hline
\end{tabular}
\subsection{Screening analysis}
The screening analysis does not require any additional options
with respect to those listed in the `Sampling options':
\verb"morris", \verb"morris_nliv", \verb"morris_ntra". The Toolbox
performs all the analyses required and displays results.
\subsection{Identification analysis (under development)}
Setting the option \verb"identification=1", an identification
analysis based on theoretical moments is performed. Sensitivity plots are provided that
allow to infer which parameters are most likely to be less
identifiable.
Prerequisite for properly running all the identification routines, is the keyword
\verb"identification;" in the DYNARE model file. This keyword triggers the computation of analytic derivatives of the model with respect to estimated parameters and shocks. This is required for options \verb"morris=2" below, which implement \cite{Iskrev2009} identification analysis.
\vspace{1cm}
\begin{tabular}{r|l|l}
option name & default & description \\ \hline
\verb"identification"& 0 & 0 = no identification analysis \\
& & 1 = performs identification analysis:\\
& & this forces \verb"redform"=0 and default\verb"morris"=1\\
\verb"morris"& 1 & 1 = Screening analysis (Type II error)\\
& & 0 = ANOVA Mapping (Type I error)\\
& & 2 = Analytic derivatives (similar to Type II error)\\
\verb"morris_nliv"& 6 & number of levels in Morris design\\
\verb"morris_ntra"& 20& number of trajectories in Morris design\\
\verb"trans_ident"& 0 & data transformation for ANOVA Mapping\\
& & 0 = no transformation\\
& & 1 = log-transformation\\
& & 2 = rank-transformation\\\hline
\end{tabular}
\vspace{1cm}
\noindent For example, the following commands in the DYNARE model file
\vspace{1cm}
\noindent\verb"identification;"\\
\verb"dynare_sensitivity(identification=1, morris=2);"
\vspace{1cm}
\noindent trigger the identification analysis using analytic derivatives \citep{Iskrev2009}, jointly with the mapping of the acceptable region.
The identification analysis with derivatives can also be triggered by the commands
\vspace{1cm}
\noindent\verb"identification;"\\
\verb"dynare_identification;"
\vspace{1cm}
This does not do the mapping of acceptable regions for the model and uses the standard random sampler of DYNARE. It completely offsets any use of the sensitivity analysis toolbox.
\newpage
\section{Directory structure}
Sensitivity analysis results are saved on the hard-disk of the
computer. The Toolbox uses a dedicated folder called \verb"GSA",
located in \\
\\
\verb"<DYNARE_file>\GSA", \\
\\
where \verb"<DYNARE_file>.mod" is the name of the DYNARE model
file.
\subsection{Binary data files}
A set of binary data files is saved in the \verb"GSA" folder:
\begin{description}
\item[]\verb"<DYNARE_file>_prior.mat": this file stores
information about the analyses performed sampling from the prior
ranges, i.e. \verb"pprior=1" and \verb"ppost=0";
\item[]\verb"<DYNARE_file>_mc.mat": this file stores
information about the analyses performed sampling from
multivariate normal, i.e. \verb"pprior=0" and \verb"ppost=0";
\item[]\verb"<DYNARE_file>_post.mat": this file stores information
about analyses performed using the Metropolis posterior sample,
i.e. \verb"ppost=1".
\end{description}
\begin{description}
\item[]\verb"<DYNARE_file>_prior_*.mat": these files store
the filtered and smoothed variables for the prior MC sample,
generated when doing RMSE analysis (\verb"pprior=1" and
\verb"ppost=0");
\item[]\verb"<DYNARE_file>_mc_*.mat": these files store
the filtered and smoothed variables for the multivariate normal MC
sample, generated when doing RMSE analysis (\verb"pprior=0" and
\verb"ppost=0").
\end{description}
\subsection{Stability analysis}
Figure files \verb"<DYNARE_file>_prior_*.fig" store results for
the stability mapping from prior MC samples:
\begin{description}
\item[]\verb"<DYNARE_file>_prior_stab_SA_*.fig": plots of the Smirnov
test analyses confronting the cdf of the sample fulfilling
Blanchard-Kahn conditions with the cdf of the rest of the sample;
\item[]\verb"<DYNARE_file>_prior_stab_indet_SA_*.fig": plots of the Smirnov
test analyses confronting the cdf of the sample producing
indeterminacy with the cdf of the original prior sample;
\item[]\verb"<DYNARE_file>_prior_stab_unst_SA_*.fig": plots of the Smirnov
test analyses confronting the cdf of the sample producing unstable
(explosive roots) behaviour with the cdf of the original prior
sample;
\item[]\verb"<DYNARE_file>_prior_stable_corr_*.fig": plots of
bivariate projections of the sample fulfilling Blanchard-Kahn
conditions;
\item[]\verb"<DYNARE_file>_prior_indeterm_corr_*.fig": plots of
bivariate projections of the sample producing indeterminacy;
\item[]\verb"<DYNARE_file>_prior_unstable_corr_*.fig": plots of
bivariate projections of the sample producing instability;
\item[]\verb"<DYNARE_file>_prior_unacceptable_corr_*.fig": plots of
bivariate projections of the sample producing unacceptable
solutions, i.e. either instability or indeterminacy or the
solution could not be found (e.g. the steady state solution could
not be found by the solver).
\end{description}
Similar conventions apply for \verb"<DYNARE_file>_mc_*.fig" files,
obtained when samples from multivariate normal are used.
\subsection{RMSE analysis}
Figure files \verb"<DYNARE_file>_rmse_*.fig" store results for the
RMSE analysis.
\begin{description}
\item[]\verb"<DYNARE_file>_rmse_prior*.fig": save results for
the analysis using prior MC samples;
\item[]\verb"<DYNARE_file>_rmse_mc*.fig": save results for
the analysis using multivariate normal MC samples;
\item[]\verb"<DYNARE_file>_rmse_post*.fig": save results for
the analysis using Metropolis posterior samples.
\end{description}
The following types of figures are saved (we show prior sample to
fix ideas, but the same conventions are used for multivariate
normal and posterior):
\begin{description}
\item[]\verb"<DYNARE_file>_rmse_prior_*.fig": for each parameter, plots the cdf's
corresponding to the best 10\% RMES's of each observed series;
\item[]\verb"<DYNARE_file>_rmse_prior_dens_*.fig": for each parameter, plots the pdf's
corresponding to the best 10\% RMES's of each observed series;
\item[]\verb"<DYNARE_file>_rmse_prior_<name of observedseries>_corr_*.fig": for each observed series plots the
bi-dimensional projections of samples with the best 10\% RMSE's,
when the correlation is significant;
\item[]\verb"<DYNARE_file>_rmse_prior_lnlik*.fig": for each observed
series, plots \emph{in red} the cdf of the log-likelihood
corresponding to the best 10\% RMSE's, \emph{in green} the cdf of
the rest of the sample and \emph{in blue }the cdf of the full
sample; this allows to see the presence of some idiosyncratic
behaviour;
\item[]\verb"<DYNARE_file>_rmse_prior_lnpost*.fig": for each observed
series, plots \emph{in red} the cdf of the log-posterior
corresponding to the best 10\% RMSE's, \emph{in green} the cdf of
the rest of the sample and \emph{in blue }the cdf of the full
sample; this allows to see idiosyncratic behaviour;
\item[]\verb"<DYNARE_file>_rmse_prior_lnprior*.fig": for each observed
series, plots \emph{in red} the cdf of the log-prior corresponding
to the best 10\% RMSE's, \emph{in green} the cdf of the rest of
the sample and \emph{in blue }the cdf of the full sample; this
allows to see idiosyncratic behaviour;
\item[]\verb"<DYNARE_file>_rmse_prior_lik_SA_*.fig": when
\verb"lik_only=1", this shows the Smirnov tests for the filtering
of the best 10\% log-likelihood values;
\item[]\verb"<DYNARE_file>_rmse_prior_post_SA_*.fig": when
\verb"lik_only=1", this shows the Smirnov test for the filtering
of the best 10\% log-posterior values.
\end{description}
\subsection{Reduced form mapping}
In the case of the mapping of the reduced form solution, synthetic
figures are saved in the \verb"\GSA" folder:
\begin{description}
\item[]\verb"<DYNARE_file>_redform_<endo name>_vs_lags_*.fig":
shows bar charts of the sensitivity indices for the \emph{ten most
important} parameters driving the reduced form coefficients of the
selected endogenous variables (\verb"namendo") versus lagged
endogenous variables (\verb"namlagendo"); suffix \verb"log"
indicates the results for log-transformed entries;
\item[]\verb"<DYNARE_file>_redform_<endo name>_vs_shocks_*.fig":
shows bar charts of the sensitivity indices for the \emph{ten most
important} parameters driving the reduced form coefficients of the
selected endogenous variables (\verb"namendo") versus exogenous
variables (\verb"namexo"); suffix \verb"log" indicates the results
for log-transformed entries;
\item[]\verb"<DYNARE_file>_redform_GSA(_log).fig": shows bar chart of
all sensitivity indices for each parameter: this allows to notice
parameters that have a minor effect for \emph{any} of the reduced
form coefficients,
\end{description}
Detailed results of the analyses are shown in the subfolder
\verb"\GSA\redform_stab", where the detailed results of the
estimation of the single functional relationships between
parameters $\theta$ and reduced form coefficient are stored in
separate directories named as:
\begin{description}
\item[]\verb"<namendo>_vs_<namlagendo>": for the entries of the
transition matrix;
\item[]\verb"<namendo>_vs_<namexo>": for entries of the matrix of
the shocks.
\end{description}
Moreover, analyses for log-transformed entries are denoted with
the following suffixes ($y$ denotes the generic reduced form
coefficient):
\begin{description}
\item[]\verb"log": $y^*=\log(y)$;
\item[]\verb"minuslog": $y^*=\log(-y)$;
\item[]\verb"logsquared": $y^*=\log(y^2)$ for symmetric fat tails;
\item[]\verb"logskew": $y^*=\log(|y+\lambda|)$ for asymmetric fat tails.
\end{description}
The optimal type of transformation is automatically selected
without the need of any user's intervention.
\subsection{Screening analysis}
The results of the screening analysis with Morris sampling design
are stored in the subfolder \verb"\GSA\SCREEN". The data file
\verb"<DYNARE_file>_prior" stores all the information of the
analysis (Morris sample, reduced form coefficients, etc.).
Screening analysis merely concerns reduced form coefficients.
Similar synthetic bar charts as for the reduced form analysis with
MC samples are saved:
\begin{description}
\item[]\verb"<DYNARE_file>_redform_<endo name>_vs_lags_*.fig":
shows bar charts of the elementary effect tests for the \emph{ten
most important} parameters driving the reduced form coefficients
of the selected endogenous variables (\verb"namendo") versus
lagged endogenous variables (\verb"namlagendo");
\item[]\verb"<DYNARE_file>_redform_<endo name>_vs_shocks_*.fig":
shows bar charts of the elementary effect tests for the \emph{ten
most important} parameters driving the reduced form coefficients
of the selected endogenous variables (\verb"namendo") versus
exogenous variables (\verb"namexo");
\item[]\verb"<DYNARE_file>_redform_screen.fig": shows bar chart of
all elementary effect tests for each parameter: this allows to
identify parameters that have a minor effect for \emph{any} of the
reduced form coefficients.
\end{description}
% ----------------------------------------------------------------
\bibliographystyle{plainnat}
%\bibliographystyle{amsplain}
%\bibliographystyle{alpha}
\bibliography{marco}
\newpage
\end{document}
% ----------------------------------------------------------------

0
doc/gsa/dyn_tbx.tex → doc/gsa/gsa.tex
View File

2260
doc/gsa/marco.bib
View File

File diff suppressed because it is too large Load Diff

2275
doc/gsa/marco.bib.bak
View File

File diff suppressed because it is too large Load Diff

2
windows/dynare.nsi
View File

@ -120,7 +120,7 @@ SectionEnd
Section "Documentation and examples (Dynare and Dynare++)"
SetOutPath $INSTDIR\doc
File ..\doc\manual.pdf ..\doc\guide.pdf ..\doc\userguide\UserGuide.pdf ..\doc\bvar-a-la-sims.pdf ..\doc\dr.pdf ..\doc\macroprocessor\macroprocessor.pdf ..\doc\preprocessor\preprocessor.pdf ..\doc\parallel\parallel.pdf
File ..\doc\manual.pdf ..\doc\guide.pdf ..\doc\userguide\UserGuide.pdf ..\doc\bvar-a-la-sims.pdf ..\doc\dr.pdf ..\doc\macroprocessor\macroprocessor.pdf ..\doc\preprocessor\preprocessor.pdf ..\doc\parallel\parallel.pdf ..\doc\gsa\gsa.pdf
SetOutPath $INSTDIR\doc\manual-html
File ..\doc\manual-html\*.html