From wdunlap at tibco.com  Mon Feb  2 17:07:27 2015
From: wdunlap at tibco.com (William Dunlap)
Date: Mon, 2 Feb 2015 08:07:27 -0800
Subject: [Rd] error code 1 from Lapack routine 'dsyevr'
In-Reply-To: <1422846495154-4702639.post@n4.nabble.com>
References: <1422713746884-4702571.post@n4.nabble.com>
	<54CDDB89.5090406@stats.ox.ac.uk>
	<1422846495154-4702639.post@n4.nabble.com>
Message-ID: <CAF8bMcYVoS2d=ybry9imR0omRj7Q1iH3F84cmOKmX-CMZPkrSQ@mail.gmail.com>

You can start diagnosing the problem by capturing the matrix that caused
eigen() to stop.   You can do this in a variety of ways; here is one
    trace(eigen, quote(lastEigenX <<- x))
After setting the trace, make your offending function call and after the
error the global variable 'lastEigenX' will contain the bad matrix.  Take
the trace off of eigen
   untrace(eigen)
If all goes well then
   eigen(lastEigenX)
will cause the same error and you or someone on this list can see what
is odd about that matrix (e.g., by looking at its singular values).




Bill Dunlap
TIBCO Software
wdunlap tibco.com

On Sun, Feb 1, 2015 at 7:08 PM, eigen <liguowei1991 at gmail.com> wrote:

> Thank you for your reply. Do you have any idea of how to get rid of the
> errors? I tried Null function to calculate eigenvectors and nearPD to get
> approximate positive definite matrix first but they also had errors.
>
>
>
> --
> View this message in context:
> http://r.789695.n4.nabble.com/error-code-1-from-Lapack-routine-dsyevr-tp4702571p4702639.html
> Sent from the R devel mailing list archive at Nabble.com.
>
> ______________________________________________
> R-devel at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>

	[[alternative HTML version deleted]]



[Rd] error code 1 from Lapack routine 'dsyevr'
William Dunlap