estimation — computes estimation.
estimation
[(OPTIONS)] ;
datafile =
FILENAME
: the datafile (a .m file or a .mat file)
nobs = INTEGER
: the number of observations to be used (default: all observations in the file)
first_obs = INTEGER
: the number of the first observation to be used (default = 1)
prefilter = 1: the estimation procedure demeans the data (default=0, no prefiltering)
presample = INTEGER
: the number of observations to be skipped before evaluating the likelihood (default = 1)
loglinear: computes a log--linear approximation of the model instead of a linear (default) approximation. The data must correspond to the definition of the variables used in the modelx.
nograph: no graphs should be plotted
lik_init: INTEGER
: type of initialization of Kalman filter.
1 (default): for stationary models, the initial matrix of variance of the error of forecast is set equal to the unconditional variance of the state variables.
2: for nonstationary models: a wide prior is used with an initial matrix of variance of the error of forecast diagonal with 10 on the diagonal.
conf_sig = {INTEGER | DOUBLE}
: the level for the confidence intervals reported in the results (default = 0.90)
mh_replic = INTEGER
: number of replication for Metropolis Hasting algorithm. For the time being, mh_replic should be larger than 1200 (default = 20000.)
mh_nblocks = INTEGER
: number of paralletl chains for Metropolis Hasting algorithm (default = 2).
mh_drop = DOUBLE
: the fraction of initially generated parameter vectors to be dropped before using posterior simulations (default = 0.5)
mh_jscale = DOUBLE
: the scale to be used for the jumping distribution in MH algorithm. The default value is rarely satisfactory. This option must be tune to obtain, ideally, an accpetation rate of 25% in the Metropolis-Hastings algorithm (default = 0.2).
mh_init_scale=DOUBLE
: the scale to be used for drawing the initial value of the Metropolis-Hastings chain (default=2*mh_scale).
mode_file=FILENAME
: name of the file containing previous value for the mode. When computing the mode, Dynare stores the mode (xparam1
) and the hessian (hh
) in a file called
.MODEL NAME
_mode
mode_compute=INTEGER
: specifies the optimizer for the mode computation.
0: the mode isn't computed. mode_file must be specified
1: uses Matlab fmincon.
2: uses Lester Ingber's Adaptive Simulated Annealing.
3: uses Matlab fminunc.
4 (default): uses Chris Sim's csminwel.
mode_check: when mode_check is set, Dynare plots the posterior density for values around the computed mode for each estimated parameter in turn. This is helpful to diagnose problems with the optimizer.
prior_trunc=DOUBLE
: probability of extreme values of the prior density that is ignored when computing bounds for the parameters (default=1e-32).
load_mh_file: when load_mh_file is declared, Dynare adds to previous Metropolis-Hastings simulations instead of starting from scratch.
optim=(fmincon options
): can be used to set options for fmincon, the optimizing function of Matlab Optimizaiton toolbox. Use Matlab syntax for these options
(default: ('display','iter','LargeScale','off','MaxFunEvals',100000,'TolFun',1e-8,'TolX',1e-6))
nodiagnostic: doesn't compute the convergence diagnostics for Metropolis (default: diagnostics are computed and displayed).
bayesian_irf triggers the computation of the posterior distribution of IRFs. The length of the IRFs are controlled by the irf option
moments_varendo triggers the computation of the posterior distribution of the theoretical moments of the endogenous variablesas well as variance decomposition
filtered_vars triggers the computation of the posterior distribution of filtered endogenous variables and shocks
smoother triggers the computation of the posterior distribution of smoothered endogenous variables and shocks
forecast = INTEGER
computes the posterior distribution of a forecast on INTEGER
periods after the end of the sample used in estimation
tex requests the printing of results and graphs in TeX tables and graphics that can be later directly included in Latex files (not yet implemented)
All options for stoch_simul
If no jscale parameter is used in estimated_params, the procedure uses mh_jscale for all parameters. If mh_jscale option isn't set, the procedure uses 0.2 for all parameters.
results from posterior optimization (also for maximum likelihood)
marginal log density
mean and shortest confidence interval from posterior simulation
Metropolis-Hastings convergence graphs that still need to be documented
graphs with prior, posterior and mode
graphs of smoothed shocks, smoothed observation errors, smoothed and historical variables
Many estimation results are made avaiable in Matlab structure oo_
. It is a structure of structures and the field names are selfexplanatory (to be expanded).
If you know how to compute the steady state for your model, you can provide a Matlab function doing the computation instead of using steady. The function should be called with the name of the .mod
file followed by _steadystate
. See fs2000a_steadystate.m
in examples/fs2000
directory.