similar to: R crash with complex matrix algebra when using EISPACK=TRUE

Can someone confirm that the EISPACK routines for eigenvalues of symmetric matrix are in base R. They seem to be, but I can't seem to locate where they are in the src tree. Thanks. Chong Gu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or

Hi, I am using eigen to get an eigen decomposition of a square, symmetric matrix. For some reason, I am getting a column in my eigen vectors (the 52nd column out of 601) that is a column of all NAs. I am using the option, symmetric=T for eigen. I just discovered that I do not get this behavior when I use the option EISPACK=T. With EISPACK=T, the 52nd eigenvector is (up to rounding error) a

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20060628/2ff7714d/attachment.pl

I have been looking at the "eigen" function and have reintroduced the ability to compute (right) eigenvalues and vectors for non-symmetric matrices. I've also made "eigen" complex capable. The code is based on the eispack entry points RS, RG, CH, CG (which is what S appears to use too). The problem with both the S and R implementations is that they consume huge amounts

Hi all, Recently our statistics students noticed that their Gibbs samplers were crashing due to some NaNs in some parameters. The NaNs came from mvrnorm (Ripley & Venables' MASS package multivariate normal sampling function) and with some more investigation it turned out that they were generated by function eigen, the eigenvalue computing function. The problem did not seem to happen

Hi, All: Attached please find a symmetric, indefinite matrix for which 'eigen(...)$vectors' included NAs: > load("eigenBug.Rdata") > sum(is.na(eigen(eigenBug)$vectors)) [1] 5670 > sessioninfo() Error: could not find function "sessioninfo" > sessionInfo() R version 2.4.1 (2006-12-18) i386-pc-mingw32 locale: LC_COLLATE=English_United

Hi all, Hi, When I want to compute the eigenvalues & eigenvectors of a specific matrix, R crashes (i.e. it stops responding to any input). I've tried it with different versions of R (2.1.1, 2.2.0, 2.2.1) - all with crashing as result. What I did before the crash was: M <- as.matrix(read.table("thematrix",header=T)) eigen(M) If, instead of eigen(M), I use eigen(M,

Hi, I'm trying to calculate a ratematrix for a RNA aligment (U instead of T) in order to use it as a ratematrix in Profidst (a phylogenetic program which takes into account both the primary sequence and the secondary structure of the RNA, in my case rRNA). The sequence-structure aligment has been made in 4SALE (a java app) and saved as one-letter encoded (using a 12 letters alphabet,

Hi, I'm interested in helping out with the Vorbis and Tarkin projects (in my already limited time!). I don't have an extensive background in all of the mathematics involved, but it looks like my mathematics from undergrad Physics is good enough to get a grasp on it. In the past few days I've started doing some reading on the mathematical basics. As part of my reading, I found

Hi, I try to calculate the angle between two first eigenvectors of different covariance matrices of biological phenotypic traits for different populations. My issue here is, that all possibilities to do so seem to normalize the eigenvectors to length 1. Although the helpfile of eigen() states, that using eigen(, symmetric = FALSE, EISPACK =TRUE) skips normalization this is (I guess) not applicable

Hi all, when testing the new improvements in the new 1.7.0-version I stumbled over the following: >eigen(matrix(c(0,.3,2,.9),2,2)) Error in eigen(matrix(c(0,.3,2,.9),2,2)) : LAPACK routine DGEEV gave error code -13 >eigen(matrix(c(0,.3,2,.9),2,2),EISPACK=TRUE) $values [1] 1.3458236 -0.4458236 $vectors [,1] [,2] [1,] -1.1436890 -0.9760443 [2,] -0.7696018

Dear all, I am currently working on the calculation of eigenvalues (and -vectors) of large matrices. Since these are mostly sparse matrices and I remember some specific functionalities in MATLAB for sparse matrices, I started a research how to optimize the calculation of eigenvalues of a sparse matrix. The function eigen itself works with the LAPACK library which has no special handling for

Hi all, I'm looking for an analytic method to obtain the PDF & CDF of the sampling distribution of a given correlation (rho) at a given sample size (N). I've attached code describing a monte carlo method of achieving this, and while it is relatively fast, an analytic solution would obviously be optimal. get.cors <- function(i, x, y, N){ end=i*N

I'm only now realizing that we have severe problems with R on our AMD 'Opteron' and 'Athlon64' clients running Redhat Enterprise with all 64-bit libraries (AFAICS). The Lapack problem happens for R-patched and R-devel both on the Opteron and the Athlon64. Here are platform details: o "gcc -v" and "g77 -v" both end with the line gcc version 3.2.3

Batch: Putting q(save=F) at the end of my file does not work in my context because I can no longer source the file without quitting. I have that quit statement in my .First so that I always quit that way interactively. The problem is that it is ignored in batch. eigen: The crash occurs on my 586 running Red Hat Linux 2.0.27 but not on my son's 486 running SLackware Linix 2.0.29. We both

Full_Name: Ole Christensen Version: 1.6.1 OS: linux-gnu Submission from: (NULL) (130.225.18.176) In package boot : > inv.logit(800) [1] NaN where it should have been 1. The problem is caused by exp(x) returning Inf when x is large. One way of fixing the problem [there may be better ways] would be to include the line out[x > 709] <- 1 in inv.logit() Cheers Ole >

http://bugs.freedesktop.org/show_bug.cgi?id=22867 Summary: [KMS] black screen on nouveau.ko module load Product: xorg Version: git Platform: x86-64 (AMD64) OS/Version: Linux (All) Status: NEW Severity: normal Priority: medium Component: Driver/nouveau AssignedTo: nouveau at lists.freedesktop.org

Hi It seems like a nice idea to expand a %h in the HostName field in ssh_config to the host given as argument to ssh. often one would have a entry in their ssh_config like: camel* User ... It's okay if the system knows that the host camel01 fx points to the right host. But what if the actual host is camel01.daimi.au.dk this is something you wouldn't like to write in your terminal.

On 20 October 2010 at 09:57, Aman Thakral wrote: | Thanks for your reply. I just created a link to | /usr/local/lib/R/lib/libRlapack.so in /usr/local/lib. It works now. | I'm just curious, is this the "right" way to do it? Is there a better way? I am coming to this a little late but ... it is a band aid to a self-inflicted wound. You only have libRlapack because you built R

Dear all I am so glad the R can provide the efficient calculate about eigenvector and eigenvalue. However, i have some puzzle about the procedure of eigen. Fristly, what kind of procedue does the R utilize such that the eigen are obtained? For example, A=matrix(c(1,2,4,3),2,2) we can define the eigenvalue lamda, such as det | 1-lamda 4 | =0 | 2 3-lamda | then