Bayesian estimation: Log data density [Laplace approximation] is NaN

[jpfeifer]

You are not handling parameter dependence correctly. In e.g.

vy_obs =log(vy)-log(vy_ss);

the value of vy_ss is fixed before estimation and will not be updated when the parameters change.
See e.g.

Replication Adolfson 2007

You cannot estimate the model with that mod-file. You are neglecting parameter dependence. After running estimation, Dynare will update the estimated parameters but not take any dependence into account that you did not specify. Rather, your calibration only once updates the other parameters depending on the estimated one. That’s why you should use model-local variables (the ones with the pound operator). See Remark 4 (Parameter dependence and the use of model-local variables) in Pfeifer(2013): “…