Displaying 20 results from an estimated 300 matches similar to: "R crash with complex matrix algebra when using EISPACK=TRUE"
2002 Jul 09
1
EISPACK symmetric matrix eigenvalue routines
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
2006 Mar 03
1
NA in eigen()
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
2006 Jun 28
1
calculating the spacial autocorrelation for poisson data
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
1997 May 18
2
R-alpha: Eigenvalue Computation Query
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
2005 May 01
2
eigen() may fail for some symmetric matrices, affects mvrnorm()
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
2007 Mar 18
1
eigen returns NAs from a real matrix
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
2006 Mar 28
2
R crashes during 'eigen'
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,
2012 Apr 11
0
phangorn and calculation of a rate matrix
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,
2001 Jan 13
1
DWIT
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
2010 Mar 19
1
Howto get unnormalized eigenvectors?
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
2003 Apr 18
1
Problem with eigen() and LAPACK
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
2012 Mar 09
1
Eigenvalue calculation of sparse matrices
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
2008 Jul 17
2
Sampling distribution (PDF & CDF) of correlation
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
2004 Jul 27
4
Problems with Lapack's eigen() on 64-bit
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
1997 Aug 25
2
R-alpha: eigen and batch
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
2002 Dec 25
2
inv.logit (package boot) (PR#2394)
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
>
2009 Jul 21
6
[Bug 22867] New: [KMS] black screen on nouveau.ko module load
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
2009 Aug 22
0
expansion of %h in HostName field of ssh_config
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.
2010 Oct 20
0
[Rpy] "lapack routines cannot be loaded" in Ubuntu Linux 9.10
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
2009 Apr 24
1
the puzzle of eigenvector and eigenvalue
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