Element of Newton isn't finite!

[amanda]

Hello everyone,
I try to reproduce a DSGE model with trend inflation (non linear). I have 19 equations for 19 endogenous variables with 3 shocks in the model. but I cant understand how to solve this message error…some help plz.

??? Error using ==> lnsrch1 at 53
Some element of Newton direction isn’t finite. Jacobian maybe singular or there is a
problem with initial values

Error in ==> solve1 at 129
[x,f,fvec,check]=lnsrch1(xold,fold,g,p,stpmax,func,j1,j2,varargin{:});

Error in ==> dynare_solve at 130
[x,info]=solve1(func,x,j1(r(i):r(i+1)-1),j2(r(i):r(i+1)-1),jacobian_flag,
bad_cond_flag, varargin{:});

[jpfeifer]

Please post the mod-file.

[amanda]

Thx for your reply,

here is the mod-file…
trendinf.mod (4.95 KB)

[jpfeifer]

You have

[code]b=(b(-1)^(rhob))*(Bs^(1-rhob))*exp(eb); %eq13

a=(As^(1-rhoa))*(a(-1)^(rhoa))*exp(ea); %eq14[/code]

with starting values 0. Thus, you have the problem that the derivative with respect to b(-1) (and a(-1)) will be b(-1)^(0.8-1)=1/b(-1)^0.2. But b has steady state 0. Thus 1/b(-1)^0.2=1/0=Inf. This yields the error.

[amanda]

Hi,

So, I’ve corrected the issue with initial values, but still, I do not get the solution, again the same problems…

[jpfeifer]

Please try the newest Dynare unstable. It should provide you with more explicit error messages.

[difficultModels2013]

Hi,

I am trying to run stochastic simulation for Canada TOTEM model but I receive a similar error message:

??? Error using ==> lnsrch1 at 53
Some element of Newton direction isn’t finite. Jacobian maybe singular or there is a problem with initial values

Error in ==> solve1 at 127
[x,f,fvec,check]=lnsrch1(xold,fold,g,p,stpmax,func,j1,j2,varargin{:});

Error in ==> dynare_solve at 126
[x,info]=solve1(func,x,j1(r(i):r(i+1)-1),j2(r(i):r(i+1)-1),jacobian_flag, bad_cond_flag, varargin{:});

Error in ==> resol at 102
[steady_state,check1] = dynare_solve(fh,steady_state,options_.jacobian_flag,…

Error in ==> stoch_simul at 66
[oo_.dr, info] = resol(oo_.steady_state,0);

Error in ==> canada at 2188
info = stoch_simul(var_list_);

Error in ==> dynare at 120
evalin(‘base’,fname) ;

canada.mod (50.3 KB)

[jpfeifer]

Put resid(1); before steady and follow the equations with residuals.

[difficultModels2013]

Thank you,

I’ve tried this and noticed that the only equation with residual significantly greater than 0 is:

lpxrow_prow = 0.85*lpxrow_prow(-1)+lpxrow_prow_shk.

lpxrow_prow_shk is a parameter that’s calibrated to be zero in the mod file and initial value for the lpxrow_prow is 2.36808551e-01.
If I change initial value for lpxrow_prow than it creates a chain reaction and other equations start giving non zero residuals. So, I’ve decided to
change the calibrated lpxrow_prow_shk value from 0 to 3.5521383e-02. Then stoch_simul(periods=2100) works. However, it took about 3.5 hours and in the end
all of the moments and correlation of simulated variables were NaNs. I’ve attached the adjusted mod file. What am I missing?

AUTOCORRELATION OF SIMULATED VARIABLES

VARIABLE 1 2 3 4 5
cbal NaN NaN NaN NaN NaN
cbal_cor NaN NaN NaN NaN NaN
defn_yn NaN NaN NaN NaN NaN
defn_yn_ss NaN NaN NaN NaN NaN
den_w NaN NaN NaN NaN NaN
deprk NaN NaN NaN NaN NaN
dr1rown NaN NaN NaN NaN NaN
gbn_cor NaN NaN NaN NaN NaN
gbn_yn NaN NaN NaN NaN NaN
gbn_yn_ss NaN NaN NaN NaN NaN
gn_yn_cor NaN NaN NaN NaN NaN
gn_yn_ss NaN NaN NaN NaN NaN
inff NaN NaN NaN NaN NaN
infprow NaN NaN NaN NaN NaN
infq NaN NaN NaN NaN NaN
infq2 NaN NaN NaN NaN NaN
infq_pg NaN NaN NaN NaN NaN
infq_row NaN NaN NaN NaN NaN
infrow NaN NaN NaN NaN NaN
infrow2 NaN NaN NaN NaN NaN
infrow_ss NaN NaN NaN NaN NaN
la NaN NaN NaN NaN NaN
la0_cor NaN NaN NaN NaN NaN
la5_cor NaN NaN NaN NaN NaN
la9_cor NaN NaN NaN NaN NaN
ladotrow NaN NaN NaN NaN NaN
la_inv NaN NaN NaN NaN NaN
la_tot NaN NaN NaN NaN NaN
lc NaN NaN NaN NaN NaN
lc0_cor NaN NaN NaN NaN NaN
lc995_cor NaN NaN NaN NaN NaN
lcfl NaN NaN NaN NaN NaN
lchm NaN NaN NaN NaN NaN
lckl NaN NaN NaN NaN NaN
lcklc NaN NaN NaN NaN NaN
lcklcm NaN NaN NaN NaN NaN
lcm NaN NaN NaN NaN NaN
lcom NaN NaN NaN NaN NaN
lcomc NaN NaN NaN NaN NaN
lcomd NaN NaN NaN NaN NaN
lcomdc NaN NaN NaN NaN NaN
lcomdg NaN NaN NaN NaN NaN
lcomdinv NaN NaN NaN NaN NaN
lcomdx NaN NaN NaN NaN NaN
lcomx NaN NaN NaN NaN NaN
lcomy NaN NaN NaN NaN NaN
lcshare NaN NaN NaN NaN NaN
lcshare_add NaN NaN NaN NaN NaN
lcshare_sreq NaN NaN NaN NaN NaN
lcx8 NaN NaN NaN NaN NaN
lc_cor NaN NaN NaN NaN NaN
lc_cor_tot NaN NaN NaN NaN NaN
leec NaN NaN NaN NaN NaN
leeg NaN NaN NaN NaN NaN
leeinv NaN NaN NaN NaN NaN
leex NaN NaN NaN NaN NaN
lforex NaN NaN NaN NaN NaN
lforexn0_cor NaN NaN NaN NaN NaN
lforexn_cor NaN NaN NaN NaN NaN
lforexn_cor_tot NaN NaN NaN NaN NaN
lg NaN NaN NaN NaN NaN
lgbn_r NaN NaN NaN NaN NaN
lgkl NaN NaN NaN NaN NaN
lgklc NaN NaN NaN NaN NaN
lgklcm NaN NaN NaN NaN NaN
lgm NaN NaN NaN NaN NaN
lhaw NaN NaN NaN NaN NaN
lhawc NaN NaN NaN NaN NaN
lhawcom NaN NaN NaN NaN NaN
lhawcom_sreq NaN NaN NaN NaN NaN
lhawc_sreq NaN NaN NaN NaN NaN
lhawg NaN NaN NaN NaN NaN
lhawg_sreq NaN NaN NaN NaN NaN
lhawinv NaN NaN NaN NaN NaN
lhawinv_sreq NaN NaN NaN NaN NaN
lhawx NaN NaN NaN NaN NaN
lhawx_sreq NaN NaN NaN NaN NaN
lhaw_sreq NaN NaN NaN NaN NaN
lindex_w NaN NaN NaN NaN NaN
linv NaN NaN NaN NaN NaN
linvc NaN NaN NaN NaN NaN
linvcom NaN NaN NaN NaN NaN
linvg NaN NaN NaN NaN NaN
linvinv NaN NaN NaN NaN NaN
linvkl NaN NaN NaN NaN NaN
linvklc NaN NaN NaN NaN NaN
linvklcm NaN NaN NaN NaN NaN
linvm NaN NaN NaN NaN NaN
linvx NaN NaN NaN NaN NaN
linv_add NaN NaN NaN NaN NaN
linv_cor NaN NaN NaN NaN NaN
linv_cor1 NaN NaN NaN NaN NaN
lk NaN NaN NaN NaN NaN
lkcom NaN NaN NaN NaN NaN
lkg NaN NaN NaN NaN NaN
lkinv NaN NaN NaN NaN NaN
lktotal NaN NaN NaN NaN NaN
lkx NaN NaN NaN NaN NaN
lk_cor NaN NaN NaN NaN NaN
ll NaN NaN NaN NaN NaN
llabshare NaN NaN NaN NaN NaN
llabshare_add NaN NaN NaN NaN NaN
llabshare_sreq NaN NaN NaN NaN NaN
llagrlpc NaN NaN NaN NaN NaN
llc NaN NaN NaN NaN NaN
llcom NaN NaN NaN NaN NaN
llg NaN NaN NaN NaN NaN
llinv NaN NaN NaN NaN NaN
llx NaN NaN NaN NaN NaN
lm NaN NaN NaN NaN NaN
lm0_cor NaN NaN NaN NaN NaN
lmanx NaN NaN NaN NaN NaN
lmcg_r NaN NaN NaN NaN NaN
lmcinv_r NaN NaN NaN NaN NaN
lmcm_r NaN NaN NaN NaN NaN
lmcsum_r NaN NaN NaN NaN NaN
lmcx_r NaN NaN NaN NaN NaN
lmc_r NaN NaN NaN NaN NaN
lmpk NaN NaN NaN NaN NaN
lmpkcom NaN NaN NaN NaN NaN
lmpkg NaN NaN NaN NaN NaN
lmpkinv NaN NaN NaN NaN NaN
lmpkx NaN NaN NaN NaN NaN
lmpl NaN NaN NaN NaN NaN
lm_cor NaN NaN NaN NaN NaN
lm_cor_tot NaN NaN NaN NaN NaN
lpcomdsum_r NaN NaN NaN NaN NaN
lpcomd_r NaN NaN NaN NaN NaN
lpcomrow_prow NaN NaN NaN NaN NaN
lpcom_r NaN NaN NaN NaN NaN
lpcrow_r NaN NaN NaN NaN NaN
lpc_cor NaN NaN NaN NaN NaN
lpc_r NaN NaN NaN NaN NaN
lpfx_r NaN NaN NaN NaN NaN
lpg_r NaN NaN NaN NaN NaN
lpinvsum_r NaN NaN NaN NaN NaN
lpinv_r NaN NaN NaN NaN NaN
lpmanx_r NaN NaN NaN NaN NaN
lpm_fitsum_r NaN NaN NaN NaN NaN
lpm_fit_r NaN NaN NaN NaN NaN
lpm_r NaN NaN NaN NaN NaN
lprow2_cor NaN NaN NaN NaN NaN
lprow2_r NaN NaN NaN NaN NaN
lprow_cor NaN NaN NaN NaN NaN
lpxrow_prow NaN NaN NaN NaN NaN
lpxrow_r NaN NaN NaN NaN NaN
lpx_r NaN NaN NaN NaN NaN
lp_r NaN NaN NaN NaN NaN
lqcom_r NaN NaN NaN NaN NaN
lqg_r NaN NaN NaN NaN NaN
lqinv_r NaN NaN NaN NaN NaN
lqx_r NaN NaN NaN NaN NaN
lq_r NaN NaN NaN NaN NaN
ltdnr_r NaN NaN NaN NaN NaN
ltincr_r NaN NaN NaN NaN NaN
ltransf_r NaN NaN NaN NaN NaN
lw NaN NaN NaN NaN NaN
lwn_r NaN NaN NaN NaN NaN
lwsum_r NaN NaN NaN NaN NaN
lx NaN NaN NaN NaN NaN
lxdc0_cor NaN NaN NaN NaN NaN
lxdc_cor NaN NaN NaN NaN NaN
lxdc_cor_tot NaN NaN NaN NaN NaN
lxdg_cor NaN NaN NaN NaN NaN
lxdinv_cor NaN NaN NaN NaN NaN
lxdx_cor NaN NaN NaN NaN NaN
lxf_cor NaN NaN NaN NaN NaN
lxkl NaN NaN NaN NaN NaN
lxklc NaN NaN NaN NaN NaN
lxklcm NaN NaN NaN NaN NaN
lxm NaN NaN NaN NaN NaN
lxw0_cor NaN NaN NaN NaN NaN
lxw_cor NaN NaN NaN NaN NaN
lxw_cor_tot NaN NaN NaN NaN NaN
lxw_r NaN NaN NaN NaN NaN
lx_cor NaN NaN NaN NaN NaN
ly NaN NaN NaN NaN NaN
lykl NaN NaN NaN NaN NaN
lyrow NaN NaN NaN NaN NaN
lyrow_cor NaN NaN NaN NaN NaN
lyrow_gap NaN NaN NaN NaN NaN
lyrow_sreq NaN NaN NaN NaN NaN
ly_gap NaN NaN NaN NaN NaN
nfa NaN NaN NaN NaN NaN
nfa_ss NaN NaN NaN NaN NaN
num_w NaN NaN NaN NaN NaN
pcomrow_prow_ss NaN NaN NaN NaN NaN
pertarget NaN NaN NaN NaN NaN
pertran NaN NaN NaN NaN NaN
profit NaN NaN NaN NaN NaN
r1 NaN NaN NaN NaN NaN
r1n NaN NaN NaN NaN NaN
r1n_ss NaN NaN NaN NaN NaN
r1n_tran NaN NaN NaN NaN NaN
r1row NaN NaN NaN NaN NaN
r1rown NaN NaN NaN NaN NaN
r1rown_des NaN NaN NaN NaN NaN
r1row_ss NaN NaN NaN NaN NaN
rfn NaN NaN NaN NaN NaN
rf_prem NaN NaN NaN NaN NaN
rgb NaN NaN NaN NaN NaN
rgbn NaN NaN NaN NaN NaN
rgbn_cor NaN NaN NaN NaN NaN
rgbn_ss NaN NaN NaN NaN NaN
rgb_ss NaN NaN NaN NaN NaN
rhn NaN NaN NaN NaN NaN
rh_prem NaN NaN NaN NaN NaN
risk NaN NaN NaN NaN NaN
risk_cor NaN NaN NaN NaN NaN
risk_dyn NaN NaN NaN NaN NaN
risk_ss NaN NaN NaN NaN NaN
rkn NaN NaN NaN NaN NaN
target NaN NaN NaN NaN NaN
tdn NaN NaN NaN NaN NaN
tdn_cor NaN NaN NaN NaN NaN
tdn_ss NaN NaN NaN NaN NaN
tinc NaN NaN NaN NaN NaN
tincr_yn_ss NaN NaN NaN NaN NaN
tinc_cor NaN NaN NaN NaN NaN
tinc_ss NaN NaN NaN NaN NaN
transf_r_cor NaN NaN NaN NaN NaN
transf_yn NaN NaN NaN NaN NaN
transf_yn_ss NaN NaN NaN NaN NaN
u NaN NaN NaN NaN NaN
ucom NaN NaN NaN NaN NaN
ug NaN NaN NaN NaN NaN
uinv NaN NaN NaN NaN NaN
usact NaN NaN NaN NaN NaN
usum NaN NaN NaN NaN NaN
utotal NaN NaN NaN NaN NaN
ux NaN NaN NaN NaN NaN
w_inf NaN NaN NaN NaN NaN
interest NaN NaN NaN NaN NaN
inflation NaN NaN NaN NaN NaN
inflationq NaN NaN NaN NaN NaN
outputgap NaN NaN NaN NaN NaN
fispol NaN NaN NaN NaN NaN
canada2.mod (50.7 KB)

[jpfeifer]

Your model is simply huge. And I guess there is still a mistake somewhere. Dynare won’t be able to help you here. Relying on pure numerical force is unlikely to solve the issue.

I also don’t get your point. If lpxrow_prow_shk is really 0, then the equation

directly implies that lpxrow_prow is also 0 in steady state. How come you have anything different?