<p>Hi,<br>
I am just playing around with hmm_train and experience some problems with the number of gaussians. I don't have a strategy on how to decide the number of gaussians, yet. So I just tried some. <br>
With values around five everything seems to be fine. But with 10 and above I get <br>
error: Mat::col(): index out of bounds<br>
terminate called after throwing an instance of 'std::logic_error'<br>
I run it like this:<br>
./hmm_train -i observationPKM.csv -n 13 -t gmm -g 10 -o hmm.xml -l labels.csv <br>
Also the Debug informations look like this:<br>
[DEBUG] Covariance matrix is not positive definite.  Adding perturbation.<br>
[DEBUG] Covariance matrix is not positive definite.  Adding perturbation.<br>
[DEBUG] Covariance matrix is not positive definite.  Adding perturbation.<br>
[DEBUG] Covariance matrix is not positive definite.  Adding perturbation.<br>
[DEBUG] Covariance matrix is not positive definite.  Adding perturbation.<br>
[DEBUG] Covariance matrix is not positive definite.  Adding perturbation.<br>
[DEBUG] Covariance matrix is not positive definite.  Adding perturbation.<br>
[DEBUG] Point 4 assigned to empty cluster 0.<br>
[DEBUG] Point 5 assigned to empty cluster 2.<br>
[DEBUG] Point 6 assigned to empty cluster 3.<br>
...<br>
Is this normal? And is there an upper limit to the number of gaussian? <br>
Greetings<br>
Davud</p>

<p style="font-size:small;-webkit-text-size-adjust:none;color:#666;">&mdash;<br>Reply to this email directly or <a href="https://github.com/mlpack/mlpack/issues/479">view it on GitHub</a>.<img alt="" height="1" src="https://github.com/notifications/beacon/AJ4bFHUDZKSrwNeVHSYxpoKgjleklcXhks5pHzISgaJpZM4GmccO.gif" width="1" /></p>
<div itemscope itemtype="http://schema.org/EmailMessage">
<div itemprop="action" itemscope itemtype="http://schema.org/ViewAction">
  <link itemprop="url" href="https://github.com/mlpack/mlpack/issues/479"></link>
  <meta itemprop="name" content="View Issue"></meta>
</div>
<meta itemprop="description" content="View this Issue on GitHub"></meta>
</div>