Sometimes the software complains about some matrix not being positive definite. Is it the singularity that's causing issues (because the iterative algorithm involves inverting matrices)? Also, what does this say about the data or model we're trying to fit?