similar to: Matrix:::qr.qy and signature(qr = "sparseQR", y = "dgCMatrix")

Dear list, Can anyone tell me how to obtain the rank of a sparse Matrix, for example from package Matrix (class dgCMatrix)? Here is an example of QR decomposition of a sparse matrix (from the sparseQR class help). library(Matrix) data(KNex) mm <- KNex$mm str(mmQR <- qr(mm)) Similarly, using the functions/classes from the relatively new MatrixModels package: library(MatrixModels)

With this example set.seed(123) A <- matrix(runif(40), nrow = 8) y <- 1:nrow(A) A.laqr <- qr(A, LAPACK=TRUE) both qr.qy(A.laqr,y) and qr.qty(A.laqr,y) give the respective error messages Error in qr.qy(A.laqr, y) : 'b' must be a numeric matrix Error in qr.qty(A.laqr, y) : 'b' must be a numeric matrix However when Lapack is not used as in A.liqr <- qr(A,

Hi, I am having problem building R 2.10, and it does not seem anyone else is having a similar problem... My system information: uname -m = x86_64 uname -r = 2.6.31-ARCH uname -s = Linux gcc (GCC) 4.4.1 icc (ICC) 11.0 20081105 The first problem I encounter seems to be with icc and wctype.h during ./configure. export CC="icc -std=c99" export CXX=icpc export OBJC="icc

Hi, I came across a case when the dimensions of matrices returned by qr() operated on a sparse matrix does not coincide with the initial matrix. Here is a spinet code that should produce an example (one of many that I could provide): ?m=205 ?n=199 ?set.seed(7); ?a=matrix(rnorm(m*n), m, n) ?a[sample(seq(m*n), m*(n-4))]=0 ?a=as(a, "Matrix") ?qa=qr(a); ?stopifnot(nrow(qa at R) ==

Hi, I tried applying http://people.freebsd.org/~anholt/dri/files/drm-stable-20030411.diff to stable from today (using device radeondrm and options DRM_LINUX). Kernel reports drm seems to be there: drm0: <ATI Radeon QY RV100 VE> port 0x9000-0x90ff mem 0xed000000-0xed00ffff,0xe$ info: [drm] Initialized radeon 1.8.0 20020828 on minor 0 Now, when I try to start XFree86 everything looks fine

salutations! i am doing some longitudinal modeling with sem and thought calculating some "simple variables" would make my model more readable. this is the smallest subset of my model that illustrates the resulting problem. i have 2 observed exogenous variables (c1, d2) and 4 observed endogenous variables (dc1, dd1, dc2, dd2). c1 is the observed state at time 1, dc1 is the change in c

Hi, One of my pbx's seems to be having some new issues. crackling interference on the zap channels running through the channel bank, and I noticed that these happen when I hit an "irq miss" in zttool: Current Alarms: No alarms. Sync Source: Digium Wildcard T100P T1/PRI C IRQ Misses: 82695 Bipolar Viol: 0 Tx/Rx Levels: 0/8 3 Total/Conf/Act:

Hi, we have installed R-2.8.1 using the current pgi compiler (8.0.2) for AMD64 on a SLES9 system. When I try to install "Matrix" everything is fine until the last step. make[1]: Leaving directory `/tmp/R.INSTALL.TW3399/Matrix/src/AMD' pgCC -L/usr/lib64 -L/usr/X11R6/lib64 -pgf90libs -o Matrix.so CHMfactor.o Csparse.o TMatrix_as.o Tsparse.o init.o Mutils.o chm_common.o

64-bit R-2.12.0 was installed on Sun SPARC Solaris 10. Compiler used is solstudio12.2. Attached is the configure script. I then tried to install a recommended package called Matrix. The compilation failed with the following messages, ....... CC -library=stlport4 -G -L/opt/csw/lib/sparcv9 -L/opt/solstudio12.2/prod/lib/v9 -o Matrix.so CHMfactor.o Csparse.o TMatrix_as.o Tsparse.o init.o Mutils.o

Something probably obivous but I don't see it. I needed to find the first 1 or two digits of some 5 and 6 digit numbers since they identified research sites while the rest of the number was the plot id. I converted the numbers to characters, got the first 1 or 2 characters as appropriate and went to add the new vector to the data.frame. For some reason R is insisting on turning the

Dear Developers, motivated by the new Sun Studio checks I compiled R-devel and several of our packages with Sun Studio 12 on Fedora x64. Everything worked fine and R-devel runs, with the exception of package Matrix where compilation crashes with the following message. The error occurs during building of the recommended packages and also if Matrix is compiled separately: [...] CC -G -lCstd

--- John Kane <jrkrideau at yahoo.ca> wrote: > Date: Thu, 19 Apr 2007 15:50:36 -0400 (EDT) > From: John Kane <jrkrideau at yahoo.ca> > Subject: RE: [R] Character coerced to factor and I > cannot get it back > To: Jorge Cornejo-Donoso <jorgecornejo at uach.cl> > > > --- Jorge Cornejo-Donoso <jorgecornejo at uach.cl> > wrote: > > >

--- John Kane <jrkrideau at yahoo.ca> wrote: > Date: Fri, 20 Apr 2007 10:47:45 -0400 (EDT) > From: John Kane <jrkrideau at yahoo.ca> > Subject: Re: [R] Character coerced to factor and I > cannot get it back > To: jim holtman <jholtman at gmail.com> > > Thanks Jim, > > I can live with it as a factor or I will do as you > suggest. What is bothering

System details: -------------- Sun Blade 1000, Solaris 8, Sun ONE Compiler Suite (rebranded Forte 7), R versions 1.6.0 and 1.6.2 - suspect others are affected, too Problem description: ------------------- Using the recommended env var setting from p.15 of the Admin Guide, the configure script fails the sgemm sunperf test, even though the sunperf libraries are available Not linking with

I need to apply a yearly inflation factor to some wages and supply some simple sums by work category. I have gone at it with a brute force "for" loop approach which seems okay as it is a small dataset. It looks a bit inelegant and given all the warnings in the Intro to R, etc, about using loops I wondered if anyone could suggest something a bit simpler or more efficent? Example:

Dear All, Suppose I have a program as below: Outside is a loop for simulation (with random generated data), inside there are several sapply()'s (10~100) over the data and something else, but these sapply's have to be sequential. And each sapply do not involve very intensive calculation (a few seconds only). So the outside loop takes minutes to finish one iteration. I guess the better way

Buenas tardes, He construido la función “myfun” al objeto de considerar aquellas persones que a partir de una determinada fecha de Apertura tienen como mínimo 65 años. Se tiene su fecha de nacimiento, su fecha de inicio en la institución y su fecha de salida de la misma. Doy vueltas al script y no acabo se saber cómo poder aplicar de un modo eficiente las instrucciones “if” ó bien “ifelse”, y me

Hello, I have a question about IMAP : I want read a IMAP mailbox in my application. I have found 2 pages about this : http://wiki.rubyonrails.org/rails/pages/HowToReceiveEmailsWithActionMailer/versions/203 http://www.codefetch.com/search?qy=login&lang=ruby But, Rails don''t want use ''net/imap''. So I can''t =( Thank you for your help pmt -- Posted via

Good evening list, I'm looking for an R implementation of the "Shuffled Complex Evolution??? (SCE-UA) algorithm after Duan et al. (1993). Does anybody know if there is an extension/ package existing that contains it? Thanks very much for your help! Cheers, Simon Duan QY, Gupta KV, Sorooshian S (1993) Shuffled Complex Evolution Approach for Effective and Efficient Global Minimization. In

-------- Original Message -------- Subject: Sun Solaris 8 compile problem. Date: Thu, 06 Nov 2003 23:03:36 -0700 From: Jack Matthews <jmatthews@t-graphic.com> To: R-help@lists.R-project.org I am new to the list was wondering if anyone can help me find the reason and remedy for compiling R on an ultrasparc 60 solaris 8 with gcc 3.2.2 in 64 bit mode? Fails in checking for a fortran