I think I found a little typo in Arellano model, although the result wouldn't change much.
To calculate value when default, defval = u(ae, ydef) + β * (θ * EVc[zero_ind, iy] + (1-θ) * EVd[1, iy]), I think EVc should be EV to make it consistent with the theory in the text and the code in Python.
Hi,
I think I found a little typo in Arellano model, although the result wouldn't change much.
To calculate value when default,
defval = u(ae, ydef) + β * (θ * EVc[zero_ind, iy] + (1-θ) * EVd[1, iy])
, I thinkEVc
should beEV
to make it consistent with the theory in the text and the code in Python.