Displaying 20 results from an estimated 5000 matches similar to: "Q: Problems with eigen() vs. svd()"
2001 Sep 06
1
svd and eigen
Hello List,
i need help for eigen and svd functions. I have a non-symmetric
square matrix. These matrix is not positive (some eigenvalues are
negative). I want to diagonalise these matrix. So, I use svd and
eigen and i compare the results. eigen give me the "good" eigenvalues
(positive and negative). I compare with another software and the
results are the same. BUT, when i use svd,
2008 May 16
1
Dimensions of svd V matrix
Hi,
I'm trying to do PCA on a n by p wide matrix (n < p), and I'd like to
get more principal components than there are rows. However, svd() only
returns a V matrix of with n columns (instead of p) unless the argument
nv=p is set (prcomp calls svd without setting it). Moreover, the
eigenvalues returned are always min(n, p) instead of p, even if nv is set:
> x <-
2010 Sep 22
3
eigen and svd
Dear R-helpers,
could anybody explain me briefly what is the difference between
eigenvectors returned by 'eigen' and 'svd' functions and how they are
related?
Thanks in advance
Ondrej Mikula
2008 May 23
1
SVD on a matix
Hi All,
I performed an svd on a matrix X and saved the first three column of the
left singular matrix U. ( I assume that they correspond to the projection of
the matrix on the first three eigen vectors that corresponds to the first
three largest eigenvalues). I would like to know how much variance is
explained by the first eigenvectors? how can I find that.
Thanks for your help
--
View this
2001 May 19
1
COMPUTING DETERMINANT FROM SVD
Dear R-users,
I computed determinant of a square matrix "var.r" using the SVD output:
detr _ 1
d _ svd(var.r)$d
for (i in 1:length(d)) {
detr _ detr*d[i]
}
print(detr)
30.20886
BUT when I tried :
det(var.r)
I got :
-30.20886
Is this because SVD output will only give absolute of the eigenvalues ?, If
this is the case
how can I get the original eigenvalues?
Thanks,
Agus
2010 May 21
2
Data reconstruction following PCA using Eigen function
Hi all,
As a molecular biologist by training, I'm fairly new to R (and statistics!),
and was hoping for some advice. First of all, I'd like to apologise if my
question is more methodological rather than relating to a specific R
function. I've done my best to search both in the forum and elsewhere but
can't seem to find an answer which works in practice.
I am carrying out
2013 Jun 18
1
eigen(symmetric=TRUE) for complex matrices
R-3.0.1 rev 62743, binary downloaded from CRAN just now; macosx 10.8.3
Hello,
eigen(symmetric=TRUE) behaves strangely when given complex matrices.
The following two lines define 'A', a 100x100 (real) symmetric matrix
which theoretical considerations [Bochner's theorem] show to be positive
definite:
jj <- matrix(0,100,100)
A <- exp(-0.1*(row(jj)-col(jj))^2)
A's being
2003 Feb 06
6
Confused by SVD and Eigenvector Decomposition in PCA
Hey, All
In principal component analysis (PCA), we want to know how many percentage
the first principal component explain the total variances among the data.
Assume the data matrix X is zero-meaned, and
I used the following procedures:
C = covriance(X) %% calculate the covariance matrix;
[EVector,EValues]=eig(C) %%
L = diag(EValues) %%L is a column vector with eigenvalues as the elements
percent
2003 Jul 03
2
SVD and spectral decompositions of a hermitian matrix
Hi:
I create a hermitian matrix and then perform its singular value
decomposition. But when I put it back, I don't get the original
hermitian matrix. I am having the same problem with spectral value
decomposition as well.
I am using R 1.7.0 on Windows. Here is my code:
X <- matrix(rnorm(16)+1i*rnorm(16),4)
X <- X + t(X)
X[upper.tri(X)] <- Conj(X[upper.tri(X)])
Y <-
2009 Nov 25
1
which to trust...princomp() or prcomp() or neither?
According to R help:
princomp() uses eigenvalues of covariance data.
prcomp() uses the SVD method.
yet when I run the (eg., USArrests) data example and compare with my own
"hand-written" versions of PCA I get what looks like the opposite.
Example:
comparing the variances I see:
Using prcomp(USArrests)
-------------------------------------
Standard deviations:
[1] 83.732400 14.212402
2000 Jul 05
0
svd() (Linpack) problems/bug for ill-conditioned matrices (PR#594)
After fixing princomp(), recently,
{tiny negative eigen-values are possible for non-negative
definite matrices}
Fritz Leisch drew my attention to the fact the not only eigen() can be
funny, but also svd().
Adrian Trappleti found that the singular values returned
can be "-0" instead of "0". This will be a problem in something like
sd <- svd(Mat) $ d
2005 Mar 14
1
r: eviews and r // eigen analysis
hi all
i have a question that about the eigen analysis found in R and in
eviews.
i used the same data set in the two packages and found different
answers. which is incorrect?
the data is:
aa ( a correlation matrix)
1 0.9801 0.9801 0.9801 0.9801
0.9801 1 0.9801 0.9801 0.9801
0.9801 0.9801 1 0.9801 0.9801
0.9801 0.9801 0.9801 1 0.9801
0.9801 0.9801 0.9801 0.9801 1
now
> svd(aa)
$d
[1] 4.9204
2007 Oct 17
3
Observations on SVD linpack errors, and a workaround
Lately I'm getting this error quite a bit:
Error in La.svd(x, nu, nv) : error code 1 from Lapack routine 'dgesdd'
I'm running R 2.5.0 on a 64 bit Intel machine running Fedora (8 I think).
Maybe the 64 bit platform is more fragile about declaring convergence.
I'm seeing way more of these errors than I ever have before.
From R-Help I see that this issue comes up from time to
2005 Apr 25
1
The eigen function
I'm using R version 2.0.1 on a Windows 2000 operating system. Here is some
actual code I executed:
> test
[,1] [,2]
[1,] 1000 500
[2,] 500 250
> eigen(test, symmetric=T)$values
[1] 1.250000e+03 -3.153033e-15
> eigen(test, symmetric=T)$values[2] >= 0
[1] FALSE
> eigen(test, symmetric=T, only.values=T)$values
[1] 1250 0
> eigen(test, symmetric=T,
2001 Mar 23
1
eigen segfaults on 0-diml matrices (PR#882)
>From one of my students' simulations:
m <- matrix(1, 0, 0) # 1 to force numeric not logical
eigen(m)
and segfault in TRED2 in src/appl/eigen.f
Easy to fix, but I wonder what else might have been overlooked?
(svd is protected).
--please do not edit the information below--
Version:
platform = sparc-sun-solaris2.7
arch = sparc
os = solaris2.7
system = sparc, solaris2.7
status =
2001 Feb 05
1
SVD of complex matrices
Is there a way to determine the SVD of a complex matrix using R?
(I'm using v1.0.1 and svd() won't do the trick).
I know LAPACK has a function to do this.
Thanks
--
Ben Stapley
Biomolecular Modelling Lab
Imperial Cancer Research
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send
2003 Feb 14
1
eigen() error: R Version 1.6.1 on Mac OS X (PR#2550)
Consider this matrix:
> sg
X1 X2 X3 X4 X5
1 3.240 2.592 2.592 2.592 2.592
2 2.592 3.240 2.592 2.592 2.592
3 2.592 2.592 3.240 2.592 2.592
4 2.592 2.592 2.592 3.240 2.592
5 2.592 2.592 2.592 2.592 3.240
If I compute the eigenvalues of the 'sg' matrix using R Version 1.5.0
(2002-04-29) under Linux (or using Version 1.4.0 (2001-12-19) under
Solaris), I obtain:
>
2006 Jan 18
1
function 'eigen' (PR#8503)
Full_Name: Pierre Legendre
Version: 2.1.1
OS: Mac OSX 10.4.3
Submission from: (NULL) (132.204.120.81)
I am reporting the mis-behaviour of the function 'eigen' in 'base', for the
following input matrix:
A <- matrix(c(2,3,4,-1,3,1,1,-2,0),3,3)
eigen(A)
I obtain the following results, which are incorrect for eigenvalues and
eigenvectors 2 and 3 (incorrect imaginary portions):
2003 Dec 22
1
La.eigen hangs R when NaN is present (PR#6003)
Full_Name: Sundar Dorai-Raj
Version: 1.8.1
OS: Windows 2000 Professional
Submission from: (NULL) (12.64.199.173)
I discovered this problem when trying to use princomp in package:mva when a
column in my matrix was all zeros and I set cor = TRUE (thus division by 0).
Doing so hangs R, never to return. I have to shut down Rterm in the Task Manager
and lose all work from the current image. I tracked
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