I added init="adapt_diag" and the sampler running smoothly. Thanks!
And it looks like can catch the volatile of r
The volatility of v looks better too.
I am not sure I finished the correct model. Because the more smooth curve generated by Waskom’s Code. ( \hat{i} is the inverse of v_cap in my model).


