On Thu, Sep 08, 2011 at 09:49:12AM -0500, Greg Sterijevski wrote:
> Cholesky, in my opinion, is not robust as you have discovered. When it
> encounters a non-psd matrix it gives up. Maybe that is the correct course of
> action, but I still think that when you are using the getCovariance to
> estimate the curvature in the neighborhood of a point it would be okay to
> take the generalized inverse and not worry too much about why your matrix is
> bordering on non-PSDness... If your optimization stops on that point, that
> is another story and should be flagged.
> 

In fact, in my current case (which is a unit test), the point is supposed to
be the minimum, by construction.  I'm trying to figure out where the
problem comes from (namely, whether the Jacobian matrix is correct)...

Gilles

> 
> > A possibly more robust option here is to use Cholesky decomposition,
> > > which is known to be stable for symmetric positive definite
> > > matrices, which the covariance matrix being inverted here should
> > > be.  The exceptions thrown will be different; but they will give
> > > more specific information about what is wrong with the covariance
> > > matrix.
> >
> > I've tried it with my problem, and it also throws an exception.
> > However, I would like to obtain the covariance matrix anyway, because I've
> > no other clue as to what might be wrong.
> > So I think that, at least, users should be able to set the positive
> > definiteness threshold in order to avoid raising an exception.
> >
> >
> >

---------------------------------------------------------------------
To unsubscribe, e-mail: dev-unsubscr...@commons.apache.org
For additional commands, e-mail: dev-h...@commons.apache.org

Reply via email to