I’m working with a model with around 50 relevant exogenous variables (i.e. excluding AR processes, “made-up” variables from other endos., etc.), and about 35 parameters, I mention this so you know approximately the size of the model.

Most of the parameters are chosen with prior assumptions or past literature, what I’m left with is about 17 parameter values that I want to calibrate with moment-matching. I already have the calibration for 3 of then through an analytic ss (exogenizing ss target for variables and endogenizing parameter values such that those targets are met).

With the other 14 parameters I want that my model matches second-moments such as standard deviations, correlations and autocorrelations. I already have a routine for doing this, but so far I have not had luck matching those data moments, tried different initial parametrizations and the furthest I got is a norm of the errors vector of around 0.03, but with the biggest relative deviation of theoretical moment from data moment of 100%, which is huge.

Thus, I would be very grateful if you may give me some advice on how to choose those initial guesses for the numeric solver, and also if it would be advisable to try something else, like dropping some of the moments to match or change of strategy. Thanks!

PD: I’m using Matlab’s `fsolve()`

to get a solution of the vector equation `F=0`

by choosing parameters to calibrate, where the elements of `F`

are the difference between model and data moments.