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	2860	July 14, 2021
Question on cholesky_factor_corr and the resulting correlation matrix Modeling	4	1056	January 8, 2020
Covariance issues with Dynamic Factor Model with large N Modeling covariance , time-series	7	547	January 17, 2024
Initialization Issues with Cholesky Factorization of Covariance Matrix Modeling performance , covariance , prior-choice	1	985	April 25, 2020
Cholesky_decompose error Modeling fitting-issues	5	1158	May 24, 2019

-nan results with lkj_corr_cholesky

Related topics