I experience some errors in my estimation (attached the code and data if anyone would like to take a look) :

78.3% of the prior support gives unique saddle-path solution.
0.0% of the prior support gives explosive dynamics.
21.6% of the prior support gives indeterminacy.
Smirnov statistics in driving acceptable behaviour
psi_pi d-stat = 0.524 p-value = 0.000
psi_pis d-stat = 0.524 p-value = 0.000

POSTERIOR KERNEL OPTIMIZATION PROBLEM!
(minus) the hessian matrix at the āmodeā is not positive definite!
=> posterior variance of the estimated parameters are not positive.
You should try to change the initial values of the parameters using
the estimated_params_init block, or use another optimization routine.
Warning: The results below are most likely wrong!

After estimation, 4 of my shocks have 0 std err and t statistic. only one has 0.9879 1.0512 . most constants provide the same pair of values for (std err, t stat) :0 .
also, the smooth shocks graph looks like a white noise for 2 of them .
this is definitely a bad estimation, but i don t understand what is the cause: constants or missmatch between the data/variab ?

If the mle analysis is performed, although dynare notifies me that the estimation is most likely wrong due to the negative hessian, why arenāt confidence intervals computed ?
what else can i use to test the strength of the model ?
iām also not sure what algorhytm to select in the "estimation"function if i donāt perform bayesian estim. so far i have used only 4 and 6.

y corresponds to the o.g. (log deviation of the trend from gdp)
i - interest rate ( diff(1) to obtain a stationary distrib)
pi - inflation rate, obtained as log diff of price index
s - real exchange rate (log diff for the same reason - stationarity) .

i tried to estimate parameters without the shocks block in order to identify the param that have issues.
the identification strength & the sensivity component with moments inf matrix plots show only blue columns (relative to param value) ; relative to prior std which should have been plotted with red are smaller than the blue ones and they don t overlap orā¦ ?

All parameters are identified. That is not the source of your problems. The observation equations are. See below. That something is severely wrong can be seen in the fact that many parameters like psi_y or nui are at their upper bounds at the mode. That always results in a non-positive definite Hessian as you did not find an interior maximum.

Without Bayesian estimation and a specified prior, identification cannot display identification strenght relative to the prior mean.

Confidence intervals for what? You get mean and standard deviation for the parameter. CI for the parameters directly follow. For IRFs, you need to add the bayesian_irf-option

Finding the mode to your model is hard. Use mode_compute=6,9 or 10.

Your observation equations are wrong. The way you do it currently, e.g. pi has mean 0 in the model, but is not demeaned in the data. Similarly, the interest rate in the model is not in first differences, but your data is. You cannot say they are equal. You need to fix that. See Pfeifer(2013): āA Guide to Specifying Observation Equations for the Estimation of DSGE Modelsā sites.google.com/site/pfeiferecon/Pfeifer_2013_Observation_Equations.pdf for details.