Hey all, I’m seeing a weird issue where my Stan model fails to initialize with a very simple model: parameters { cholesky_factor_corr[100] ccm; } model { ccm ~ lkj_corr_cholesky(1); } because lkj_corr_cholesky(1); on any of the constrained cholesky_factor_corr matrices returns -nan, I think b…

For Cholesky factors that big, it is possible that the diagonal elements result in the square root of a slightly negative number when symbolically they should be slightly positive. If you search for the “onion method” on the Stan-users Google group, there is some code for a more numerically stable w…

-nan results with lkj_corr_cholesky

Developers

bgoodri July 22, 2020, 1:35am 3

onion.stan (1.1 KB) is the 2020 version.

4 Likes

Using the onion method in occupancy hierarchical model

Topic		Replies	Views
Cholesky_factor_corr initializing with invalid values Developers	8	2849	July 14, 2021
Cholesky_decompose error Modeling fitting-issues	5	1141	May 24, 2019
Initialization between (-2, 2) failed after 100 attempts Modeling specification	4	397	May 22, 2023
Bad LKJ Initialization in model from "Fast hierarchical Gaussian processes" Modeling specification	2	1455	March 22, 2018
Newly-Specified BHM Returns "Initialization Failed" (Edited) Modeling specification	1	263	September 27, 2023

-nan results with lkj_corr_cholesky

Related Topics