Hello Stanimaniacs,
Once I asked how to code inverse transformation by Stan. And then, I learned Newton method.
After that, I noticed my transformation failed when skewness parameter is larger or lesser. So I tried coding inverse transformation by bisection method.
But even though the transformation goes well, checking the expose_stan_function by rstan, the optimization failed by the model. When I saw the error in console, LS failed
was printed.
Anyone knows the problem in my code?
The operated codes were located in Google Colab.
The plots of transformation result by rstan were also attached. First one is by bisection’s method, and second one is by Netwon’s method. This transformation uses nu=0.9, and it failed by Netwon’s method as the gradient become steep. So I implemented by bisection method (and it failed!!).