similar to: Observations on SVD linpack errors, and a workaround

Hi, I am looking for C/C++ codes for computing generalized inverse of a matrix. Can anyone help me in this regard? Thanks, Mahbub.

Dear All, Under Debian GNU/Linux La.svd (with method = "dgesdd") sometimes gives the error "Error in La.svd(data, nu = 0, nv = min(nrow, ncol), method = "dgesdd") : LAPACK routine DGESDD gave error code -12" It seems not to depend on the data per se, but on the relationship between numbers of rows and columns. For example, if the number of columns is 100,

Hi, I wonder where I can find an example of using a function in LAPACK library in a user's own C code. I wrote a C program which will be compiled and linked to produce a DLL file and then loaded into R. I hope to use a function from LAPACK library, for example, dgesdd, in the program. Following R manual, I call the function by F77_CALL(dgesdd) in the program. The program can be compiled

Dear R helpers, I am working with R 2.4.1 GUI 1.18 (4038) for MacOSX. I have a matrix of 10 000 genes and try to run the following commands: > model.mix<-makeModel (data=data, formula=~Dye+Array+Sample+Time, random=~Array+Sample) > anova.mix<-fitmaanova (data, model.mix) > test.mix<-matest (data, model=model.mix, term="Time", n.perm=100, test.method=c(1,0,1,1))

Error in La.svd(x, nu, nv) : error code 1 from Lapack routine ‘dgesdd’ what resources are there to track down errors like this [[alternative HTML version deleted]]

Full_Name: Hugh Miller Version: 2.8.1 OS: XP Submission from: (NULL) (128.250.24.101) Hi there, I've been looking at approaches that use the projection pursuit regression function fairly (ppr) heavily, and have discovered that it crashes my R system on occasion. It only happens with the inputs are pathological in some way I don't understand. I have pasted such an example below. Any

I am attempting to perform a pca on a data frame of dimension 5000x19, but when I execute pcapres<-prcomp(pres,center=TRUE) the following error message is returned: Error in La.svd(x, nu, nv, method) : error code 18 from Lapack routine dgesdd Where am I going wrong? I am running R-1.8.0 on Debian. Regards, Laura

Hi, I'm just wondering if it is planned to include the Lapack routine DGESDD (and friends) in R-1.4.0? This is faster (supposedly by a factor of ~6 for large matrices) than DGESVD which is currently (R-1.3.1) called by La.svd. And if it is not in the plans yet, is there a chance it could be? I've added it to my local version of R-1.3.1 and so far see a factor of 4 improvement over

Dear all I have a big matrix of standardized values (dimensions 285x5829) and R fails to calculate the principal components using prcomp() with the following error message: pc <- prcomp(my.matrix) Error in La.svd(x, nu, nv, method) : error code 1 from Lapack routine dgesdd Is the matrix too big? I'm using R-1.8.1 under Unix (Solaris8) and Linux(Suse 8.2). I tried to perform a principal

Hello R-users, during one of my analyses that involve a SVD, I get the following error message: Error in La.svd(x, nu, nv, method) : error code 1 from Lapack routine dgesdd With a search on the R web site, I only found references to error codes 17 and 3 for this particular routine. I also found the Lapack web site, but could not find a list of the possible error messages. If somebody knows what

Well, CentOS 5.3, which amounts to the same thing. I recently decided to upgrade my main research machine from Fedora Core 8 -> CentOS 5.3. Basically, I was looking to move to a distro with longer 'term-of-life' than the release schedule for Fedora currently allows. The machine is a multi-Opteron box, so both 32- and 64-bit apps natively supported. Since I do a lot of 'linear

Dear all, Has anyone seen the same problem? I tried compiling R-1.9.0 alpha (2004-03-16) on our Opteron box running SUSE Linux ES8. I ran: ./configure --enable-R-shlib --with-blas=goto and got: Source directory: . Installation directory: /usr/local C compiler: gcc -m64 -O2 -g -msse2 -march=k8 -Wall --pedantic C++ compiler: g++ -m64 -O2 -g

> Problem with this is client needs to listen to the call recordings and my interface will only display .wav or .mp3 so they will moan if they have to wait until the next day for today's recordings If you're up to writing a bit of shell script, and are running on Linux, you could automate the conversion process so that it happens as soon as the recording is completed. Look at the

Dear all, The 'make check' step fails for the pacakge mva on IBM AIX. The tail of the Rout log file looks like: > for(factors in 2:4) print(update(Harman23.FA, factors = factors)) Call: factanal(factors = factors, covmat = Harman23.cor) Uniquenesses: height arm.span forearm lower.leg weight 0.170 0.107 0.166

Dear all, I've been learning biplot (Gabriel, 1971) and some days ago I sent for this list a procedural function with invitation for a collaborative package. Jari Oksanen made some suggestions and I agree with all. So, I reworked the function under the object-oriented programming (OOP/S3). I think it is now a good frame for more resources. Below it is the function and a small script to

Dear Philipp, this is just a tentative answer because debugging is really not possible without a reproducible example (or, at a very bare minimum, the output from traceback()). Anyway, thank you for reporting this interesting numerical issue; I'll try to replicate some similar behaviour on a similarly dimensioned artificial dataset when I have some time (which might not be soon). As for now,

Dear R community, I have a data matrix (531X314), and would like to apply the prcomp. However, I got this error Lapack message. I am using R3.2.2. I googled a bit and found that it might be related to converge issue. ?Just wonder if there is a way to get around it? Thank you very much! Ace On Thursday, December 29, 2016 11:44 AM, Ista Zahn <istazahn at gmail.com> wrote: Use

Commit-ID: f97573c5c1a47595c5a059b84b74f54cb2cd8d6d Gitweb: http://git.kernel.org/?p=libs/klibc/klibc.git;a=commit;h=f97573c5c1a47595c5a059b84b74f54cb2cd8d6d Author: Stephen Kitt <steve at sk2.org> AuthorDate: Mon, 6 Jun 2016 22:45:27 +0800 Committer: Ben Hutchings <ben at decadent.org.uk> CommitDate: Fri, 25 Jan 2019 02:57:21 +0000 [klibc] jobs: Handle string-based job

Commit-ID: 8115b5d06f75ea543d01355384cf5556aa388481 Gitweb: http://git.kernel.org/?p=libs/klibc/klibc.git;a=commit;h=8115b5d06f75ea543d01355384cf5556aa388481 Author: Stephen Kitt <steve at sk2.org> AuthorDate: Mon, 6 Jun 2016 22:45:27 +0800 Committer: Ben Hutchings <ben at decadent.org.uk> CommitDate: Sat, 28 Mar 2020 21:42:54 +0000 [klibc] dash: jobs: Handle string-based

Here''s the script. I would have posted it to a website, but I figure this''ll be better so it is always on the list for people in the future in case they wanted to see it. If you have any questions about any of it, please let me know. Mike and I aren''t the cleanest of programmers. Cron the script to run whenever you need it to. :) #!/usr/bin/perl # # TC Helper Script: