Closed rossfadely closed 12 years ago
oh I understand why that is! Remember me saying to Dilip "why is this positive definite?" It relies on the psis being not-insane -- same with the lambdas. In general the safe move is to initialize everything to a small variance -- maybe even set by the smallest eigenvalue of the empirical data variance tensor....
see initialization thread
The culprit seems to be the psi initializations - too large of variance causes the initial cov to have negative logdet sometimes.
Setting the initialization of the psis to 0.5 * datavar seems to be ok for 100 seeds..