Proposed specification for estimation functions
This includes factorization of common code and extensions for (unexpected) structural change
Estimated parameters description
xparam1 is a column n by 1 vector holding the current values of the estimated parameters (n is the number of parameters).
parameters_description is a (global?) strucure of arrays describing the estimated parameters with the following fields:
type is a n by one vector of integers (1,2,3,4,5) where:
1 -> standard deviation of a structural shock.
2 -> standard deviation of a measurement error.
3 -> correlation betwwen two structural shocks.
4 -> correlation between two measurement errors.
5 -> "deep" parameter.
ID1 is a n by 1 vector of indices targeting to the (parameter, exogenous variable or endogenous variable) names (in M_.param_names, M_.exo_names and M_.endo_names).
ID2 is a n by 1 vector of NaNs or indices (second exogenous or endogenous variable index for the estimated correlations).
StartPeriod n by 1 (beginning of sub-period)
EndPeriod n by 1 (end of sub-period)
Prior is a sub-structure describing the prior belief on a parameter. The fiels are:
shape is an n by 1 vector of integers (1,2,3,4,5,6) where:
1 -> Beta density.
2 -> Gamma density.
3 -> Gaussian density.
4 -> Inverse gamma (type 1) density.
5 -> Uniform density.
6 -> Inverse gamma (type 2) density.
mean is an n by 1 vector of doubles (mean of the prior density).
mode is an n by 1 vector of doubles (mode of the prior density).
standard is an n by 1 vector of doubles (standard deviation of the prior density).
lower_bound is an n by 1 vector of doubles (and Inf/-Inf) (prior density domain).
upper_bound is an n by 1 vector of doubles (and Inf/-Inf) (prior density domain).
fhp is an n by 1 vector of doubles (first hypermarameter
for the beta and gamma priors, 
for the inverse gamma prior, lower bound for the uniform prior). shp is an n by 1 vector of doubles (second hypermarameter
for the beta and gamma priors, 
for the inverse gamma prior, upper bound for the uniform prior).
lower_bound n ny 1 vector of doubles (effective lower bound for the estimated parameter).
upper_bound n by 1 vector of doubles (effective upper bound for the estimated parameter).
nvx integer scalar (number of estimated structural shock standard deviations).
nvx_id nvx by 1 vector of integers.
nvn integer scalar (number of estimated measurement error standard deviations).
nvn_id nvn by 1 vector of integers.
ncx integer scalar (number of estimated structural shock correlations).
ncx_id ncx by 1 vector of integers.
ncn integer scalar (number of estimated measurement error correlations).
ncn_id ncn by 1 vector of integers.
np integer scalar (number of estimated "deep" parameters).
np_id np by 1 vector of integers.
Remark 1. The considered subsamples for different parameters may be different.
Remark 2. StartPeriod and EndPeriod specify the closed sub sample used to estimate the parameter. If these two fields are empty the parameter is estimated using the whole sample.
Remark 3. The field Prior is empty if the model is estimated by maximum lilkelihood, simulated moments or indirect inference.
Remark 4. xxx_id vectors select subset of xparam1 (for instance, xparam_1(np_id) selects all the lines of xparam1 corresponding to the estimated "deep" parameters).
Functions
- Likelihood/posterior evaluation (separate two functions ?)
- INPUTS: xparam1, data, parameters_specification, ys, coeff_trends, M_.params, M_.Sigma_e, M_.H, list_of_observed_variables
- OUTPUTS: likelihood/posterior value, info
- Parameter updates (for period 1):
- INPUTS: xparam1, ys, coeff_trends,M_.params, M_.Sigma_e, M_.H
- OUTUTS: coeff_trends, params, Sigma_e, H
- reduced form (for period 1):
- INPUTS: params
- OUTPUTS: T R ys
- Computes steady state
- Computes linear solution for state variables + observed variables
- data_filtering (for period 1):
- INPUTS: data, ys, coeff_trends, StartPeriod, EndPeriod
- OUTPUTS: centered_data for subperiod
- Kalman filter initialization
- INPUTS: Sigma_e T R list_of_observed_variables
- OUTPUTS: a0 P0 Z
- Remark: Z is computed only for nonstationary models. mfs is computed once before these functions
- Diffuse Kalman filter recursions takes place here ?
- Kalman filter
- INPUTS: a0 Po Z T R centered_data Sigma_e H - OUTPUTS: log-likelihood, a(t|t-1), P(t|t-1)
if period < NP:
- parameter updates
- reduced form
- data_filtering
- priors evaluation