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.
On 02 Feb 2015, at 04:08 , 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. >How could we? All we know is that you are having trouble running some unspecified code on some unspecified data on an unspecified platform. It is a bit like someone calling from Finland and telling you that they are lost looking for the University and are standing at a roadsign pointing to "Keskus".> -- > 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-- Peter Dalgaard, Professor, Center for Statistics, Copenhagen Business School Solbjerg Plads 3, 2000 Frederiksberg, Denmark Phone: (+45)38153501 Email: pd.mes at cbs.dk Priv: PDalgd at gmail.com
fortune candidate On Mon, Feb 2, 2015 at 8:54 AM, peter dalgaard <pdalgd at gmail.com> wrote:> > On 02 Feb 2015, at 04:08 , 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. > > > > How could we? All we know is that you are having trouble running some > unspecified code on some unspecified data on an unspecified platform. > > It is a bit like someone calling from Finland and telling you that they > are lost looking for the University and are standing at a roadsign pointing > to "Keskus". > > > > -- > > 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 > > -- > Peter Dalgaard, Professor, > Center for Statistics, Copenhagen Business School > Solbjerg Plads 3, 2000 Frederiksberg, Denmark > Phone: (+45)38153501 > Email: pd.mes at cbs.dk Priv: PDalgd at gmail.com > > ______________________________________________ > R-devel at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-devel >[[alternative HTML version deleted]]
Sorry for being unclear. I am using R x64 3.1.2 for simulation. Basically I got a matrix of genetic data and calculated the kernel matrix. Now I have trouble applying eigen function on the kernel matrix for a few simulations but not most of them. I tried other methods but none of them worked. I guess my kernel matrix may not be very good. Thank you. -- View this message in context: http://r.789695.n4.nabble.com/error-code-1-from-Lapack-routine-dsyevr-tp4702571p4702676.html Sent from the R devel mailing list archive at Nabble.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]]
On 2 February 2015 at 10:07, William Dunlap <wdunlap at tibco.com> wrote: <snip>> > 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).Preferably *not* this list as this doesn't really seem to be about developing R or with/for it. R-Help would have been a more appropriate venue for a question such as this. G> > 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]] > > ______________________________________________ > R-devel at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-devel >-- Gavin Simpson, PhD [[alternative HTML version deleted]]