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 <liguowei1...@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@r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-devel > [[alternative HTML version deleted]] ______________________________________________ R-devel@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-devel