Displaying 20 results from an estimated 4000 matches similar to: "Query regarding SVD of binary matrix:"
2002 Oct 29
0
patch to mva:prcomp to use La.svd instead of svd (PR#2227)
Per the discussion about the problems with prcomp() when n << p, which
boils down to a problem with svd() when n << p,
here is a patch to prcomp() which substitutes La.svd() instead of svd().
-Greg
(This is really a feature enhancement, but submitted to R-bugs to make sure
it doesn't get lost. )
*** R-1.6.0/src/library/mva/R/prcomp.R Mon Aug 13 17:41:50 2001
---
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
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 <-
2007 Feb 05
0
strange error message get from La.svd(X)
Generator Microsoft Word 11 (filtered medium) Hi,
I'm the mannova package maintainer. We used La.svd(X, method="dgesvd") in maanova package before. After R-2.3.0, the old La.svd() method was deprecated for option method="dgesvd". I changed maanova code correspondingly, which will call method="dgesdd" instead. But after that, we keep getting below error message
2002 Nov 17
1
SVD for reducing dimensions
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
Hi all, this is probably simple and I'm just doing something stupid, sorry
about that :-)
I'm trying to convert words (strings of letters) into a fairly small
dimensional space (say 10, but anything between about 5 and 50 would be ok),
which I will call a feature vector. The the distance between two words
represents the similarity of the
2008 Apr 15
1
SVD of a variance matrix
Hello!
I suppose this is more a matrix theory question than a question on R,
but I will give it a try...
I am using La.svd to compute the singular value decomposition (SVD) of
a variance matrix, i.e., a symmetric nonnegative definite square
matrix. Let S be my variance matrix, and S = U D V' be its SVD. In my
numerical experiments I always got U = V. Is this necessarily the
case? Or I might
2001 Nov 06
1
R-devel & ATLAS generates Dr. Watson on NT (was RE: Look, Wa tson! La.svd & ATLAS)
Prof. Bates & R-devel,
I've done more test with the following results:
I have two versions of ATLAS 3.2.1. One was compiled on my old Thinkpad
600E (PII), the other was compiled on my new Thinkpad T22 (PIIISSE1).
I compiled R-devel dated 10/31, 11/01 and 11/04, linked against either of
the two ATLAS libs. All gave Dr. Watson when given this code:
La.svd(matrix(runif(1e5), 1e3,
2012 Dec 05
1
Understanding svd usage and its necessity in generalized inverse calculation
Dear R-devel:
I could use some advice about matrix calculations and steps that might
make for faster computation of generalized inverses. It appears in
some projects there is a bottleneck at the use of svd in calculation
of generalized inverses.
Here's some Rprof output I need to understand.
> summaryRprof("Amelia.out")
$by.self
self.time self.pct
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,
2002 Apr 01
3
svd, La.svd (PR#1427)
(I tried to send this earlier, but it doesnt seem to have come through,
due to
problems on my system)
Hola:
Both cannot be correct:
> m <- matrix(1:4, 2)
> svd(m)
$d
[1] 5.4649857 0.3659662
$u
[,1] [,2]
[1,] -0.5760484 -0.8174156
[2,] -0.8174156 0.5760484
$v
[,1] [,2]
[1,] -0.4045536 0.9145143
[2,] -0.9145143 -0.4045536
> La.svd(m)
$d
[1]
2011 Sep 13
1
SVD Memory Issue
I am trying to perform Singular Value Decomposition (SVD) on a Term Document
Matrix I created using the 'tm' package. Eventually I want to do a Latent
Semantic Analysis (LSA).
There are 5677 documents with 771 terms (the DTM is 771 x 5677). When I try
to do the SVD, it runs out of memory. I am using a 12GB Dual core Machine
with Windows XP and don't think I can increase the memory
2000 Aug 10
1
svd error (PR#631)
--=====================_24736660==_
Content-Type: text/plain; charset="iso-8859-1"; format=flowed
Content-Transfer-Encoding: quoted-printable
SVD-Error on
R 1.1.0
Windows 98
I get the following error applying svd on a positive definite matrix :
> sk2
[,1] [,2] [,3] [,4] [,5]
[1,] 1.0460139783 0.084356992 -2.810553e-04
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
2012 Dec 06
1
svd(X, LINPACK=TRUE) alters its input
Ordinary functions should not alter their inputs but in R-2.15.2
svd(LINPACK=TRUE,X) does. (It worked in 2.15.0 but not in 2.15.1
or 2.15.2 and became deprecated in 2.15.2.)
> X <- matrix(c(1,2,3, 5,7,11, 13,17,19), 3, 3)
> X
[,1] [,2] [,3]
[1,] 1 5 13
[2,] 2 7 17
[3,] 3 11 19
> svd(X, LINPACK=TRUE)$d
[1] 31.9718214 2.3882717 0.3143114
Warning message:
2013 Jan 14
1
ginv / LAPACK-SVD causes R to segfault on a large matrix.
Dear R-help list members,
I am hoping to get you help in reproducing a problem I am having That is
only reproducible on a large-memory machine. Whenever I run the following
lines, get a segfault listed below:
*** caught segfault ***
address 0x7f092cc46e40, cause 'invalid permissions'
Traceback:
1: La.svd(x, nu, nv)
2: svd(X)
3: ginv(bigmatrix)
Here is the code that I run:
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
2010 Sep 02
0
using R's svd from outside R
Hi,
I have to compute the singular value decomposition of rather large
matrices. My test matrix is 10558 by 4255 and it takes about three
minutes in R to decompose on a 64bit quadruple core linux machine. (R is
running svd in parallel, all four cores are at their maximum load while
doing this.) I tried several blas and lapack libraries as well as the
gnu scientific library in my C++ programm.
2000 May 10
4
Q: Problems with eigen() vs. svd()
At 01:37 PM 5/10/00 +0200, ralle wrote:
>Hi,
>I have a problem understanding what is going on with eigen() for
>nonsymmetric matrices.
>Example:
>h<-rnorm(6)
>> dim(h)<-c(2,3)
>> c<-rnorm(6)
"c" is not a great choice of identifier!
>> dim(c)<-c(3,2)
>> Pi<-h %*% c
>> eigen(Pi)$values
>[1] 1.56216542 0.07147773
These could
2001 Nov 02
1
Look, Watson! La.svd & ATLAS
Dear R-devel,
I had attempted to compile r-devel (dated Oct. 31, 2001) on WinNT with link
to ATLAS, with mostly success. However, when I tried the following, I got a
visit from Dr. Watson:
R : Copyright 2001, The R Development Core Team
Version 1.4.0 Under development (unstable) (2001-10-31)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under
2009 Aug 09
1
Inaccuracy in svd() with R ubuntu package
On two laptops running 32-bit kubuntu, I have found that svd(), invoked
within R 2.9.1 as supplied with the current ubuntu package, returns very
incorrect results when presented with complex-valued input. One of the
laptops is a Dell D620, the other a MacBook Pro. I've also verified the
problem on a 32-bit desktop. On these same systems, R compiled from
source provides apparently