[R] Generalized SVD
Hi,
I need to computed the GSVD of two matrices. For doing so I used in my C
code the lapack routine dggsvd. But when I source my file gsvd.R I have
the following error:
Error in eval.with.vis(expr, envir, enclos) :
BLAS/LAPACK routine 'DGGSVD' gave error code -1
Is there a problem with the parameters passed through the gsvd.R script?
Is there a problem within the C script?
I will really appreciate your help.
Thanks,
Philippe
--
Philippe Hupé
Institut Curie, CNRS UMR 144 / Service Bioinformatique
Laboratoire de Transfert (4ème étage)
26 rue d'Ulm
75005 Paris - France
Email : [EMAIL PROTECTED]
Tél :+33 (0)1 44 32 42 75
Fax :+33 (0)1 42 34 65 28
website : http://bioinfo.curie.fr
vA - c(1.0, 2.0, 3.0,
3.0, 2.0, 1.0,
4.0, 5.0, 6.0,
7.0, 8.0, 8.0)
A - matrix(vA,4,3, byrow=TRUE)
vB - c(-2.0, -3.0, 3.0,
4.0, 6.0, 5.0)
B - matrix(vB,2,3, byrow=TRUE)
if(ncol(A)!=ncol(B))
{
stop(A and B matrix does not have the same column number)
}
gsvd - .C(GSVD,
JOBU = as.character('U'),
JOBV = as.character('V'),
JOBQ = as.character('Q'),
M = as.integer(nrow(A)),
N = as.integer(ncol(A)),
P = as.integer(nrow(B)),
K = integer(1),
L = integer(1),
A = as.double(A),
LDA = as.integer(nrow(A)),
B = as.double(B),
LDB = as.integer(nrow(B)),
ALPHA = double(ncol(A)),
BETA = double(ncol(A)),
U = double(nrow(A)*nrow(A)),
LDU = as.integer(nrow(A)),
V = double(nrow(B)*nrow(B)),
LDV = as.integer(nrow(B)),
Q = double(ncol(A)*ncol(A)),
LDQ = as.integer(ncol(A)),
WORK = double(max(c(3*ncol(A),nrow(A),nrow(B)))+ncol(A)),
IWORK = integer(ncol(A)),
INFO = integer(1),
PACKAGE=GSVD)
GSVD_1.0.tar.gz
Description: PostScript document
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] * creating vignettes ... ERROR
Hi, I have the following error when building a package * creating vignettes ... ERROR Error: chunk 3 Error in pmatch(x, table, duplicates.ok) : argument is not of mode character Error in buildVignettes(dir = .) : Error: chunk 3 Error in pmatch(x, table, duplicates.ok) : argument is not of mode character Execution halted I use R-2.0.1 under debian linux 2.6.9 Thanks for any idea Philippe -- Philippe Hupé UMR 144 - Service Bioinformatique Institut Curie Laboratoire de Transfert (4ème étage) 26 rue d'Ulm 75005 Paris - France Email : [EMAIL PROTECTED] Tél :+33 (0)1 44 32 42 75 Fax :+33 (0)1 42 34 65 28 __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] aggregate and median
I am trying to use the function aggregate with the median function but I get the following error: Error in FUN(X[[1]], ...) : Argument INDEX When I replace median by mean, it works perfectly Can someone tell me where the problem comes from? Thx I am running R 2.0.0 on SunOS 5.9 -- Philippe Hupé UMR 144 - Service Bioinformatique Institut Curie Laboratoire de Transfert (4ème étage) 26 rue d'Ulm 75005 Paris - France Email : [EMAIL PROTECTED] Tél :+33 (0)1 44 32 42 75 Fax :+33 (0)1 42 34 65 28 __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] Colored label on dendrogram
Hi, Can someone tell me how to put labels of differente colors on a tree dendrogram. Thanks -- Philippe Hupé UMR 144 - Service Bioinformatique Institut Curie Laboratoire de Transfert (4ème étage) 26 rue d'Ulm 75005 Paris - France Email : [EMAIL PROTECTED] Tél :+33 (0)1 44 32 42 75 Fax :+33 (0)1 42 34 65 28 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] problem with aggregate
Hi, I am trying to do the following aggregation : data - data.frame(a=rep(2,10),b=rep(a,10)) aggregate(data$a, by=list(a1=data$a, b1=data$b), NROW) but I have the following error message Error in names-.default(`*tmp*`, value = c(names(by), names(x))) : names attribute [3] must be the same length as the vector [2] In addition: Warning message: row names were found from a short variable and have been discarded in: data.frame(w, lapply(y, unlist, use.names = FALSE)) It works if there are more than 2 modalities in data$b. I can test if there is only one modality of course but there is no reason for the function not working in this particular case. I use R 1.9.0 under Solaris -- Philippe Hupé UMR 144 - Service Bioinformatique Institut Curie Laboratoire de Transfert (4ème étage) 26 rue d'Ulm 75005 Paris - France Email : [EMAIL PROTECTED] Tél :+33 (0)1 44 32 42 75 Fax :+33 (0)1 42 34 65 28 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
[R] how to insert a double quote with the paste function?
Hello, I would like to write in a variable the following string : Hello World If I do cat(Hello ,\World\), I have the good result on the screen but it can not be affected to a variable. If I do the same thing with paste paste(Hello ,\World\) it does not work since \ seems to be not recognized. So what is the trick to do that. Thanks, Philippe -- -- Philippe Hupé Institut Curie - Equipe Bioinformatique 26, rue d'Ulm - 75005 PARIS France +33 (0)1 44 32 42 75 [EMAIL PROTECTED] mailto:[EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] How to speed the shut down of the graphic device
Andrew C. Ward a écrit :
Dear Philippe,
Perhaps you could try a different graphics device (maybe
postscript). On my machine, the time differences were all 1
second rather than the 3 you reported.If 300s is really
too long for you, you could get a new computer or run your
script on a faster one.
I have a computer with 3GHz processor (Pentium IV) and 2 Go of RAM so I
don't think this is a matter of computer performance :).
Andrew C. Ward
CAPE Centre
Department of Chemical Engineering
The University of Queensland
Brisbane Qld 4072 Australia
[EMAIL PROTECTED]
Quoting Philippe Hupé [EMAIL PROTECTED]:
hello,
I use R1.7.1 under winXP and I am running the following
script example :
for (i in 1:10)
{
x - rnorm(100)
png( paste(D:/essai,i,.png,sep=))
plot(x)
t1 - Sys.time()
dev.off()
t2 - Sys.time()
print(t2-t1)
}
at each step, it takes about 3 seconds to shut down the
graphic device.
I want to generate about one hundred of image and of
course it takes too
much time. Is there any trick ?
Philippe
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
.
--
--
Philippe Hupé
Institut Curie - Equipe Bioinformatique
26, rue d'Ulm - 75005 PARIS France
+33 (0)1 42 34 65 29
[EMAIL PROTECTED] mailto:[EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] How to speed the shut down of the graphic device
Prof Brian Ripley a écrit :
On Mon, 4 Aug 2003, Philippe Hupé wrote:
Andrew C. Ward a écrit :
Dear Philippe,
Perhaps you could try a different graphics device (maybe
postscript). On my machine, the time differences were all 1
second rather than the 3 you reported.If 300s is really
too long for you, you could get a new computer or run your
script on a faster one.
I have a computer with 3GHz processor (Pentium IV) and 2 Go of RAM so I
don't think this is a matter of computer performance :).
My lowly 1.4GHz P4M laptop does it in 2.5 secs, so there does seem to be a
performance problem with your computer.
BTW, the way to time a command is to use system.time()
Andrew C. Ward
CAPE Centre
Department of Chemical Engineering
The University of Queensland
Brisbane Qld 4072 Australia
[EMAIL PROTECTED]
Quoting Philippe Hupé [EMAIL PROTECTED]:
hello,
I use R1.7.1 under winXP and I am running the following
script example :
for (i in 1:10)
{
x - rnorm(100)
png( paste(D:/essai,i,.png,sep=))
plot(x)
t1 - Sys.time()
dev.off()
t2 - Sys.time()
print(t2-t1)
}
at each step, it takes about 3 seconds to shut down the
graphic device.
I want to generate about one hundred of image and of
course it takes too
much time. Is there any trick ?
Philippe
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
.
I have done the same thing under Debian/sid wiht R1.7.1 and it takes
almost 0 second !! Then, is there any problem with winXP ?
--
--
Philippe Hupé
Institut Curie - Equipe Bioinformatique
26, rue d'Ulm - 75005 PARIS France
+33 (0)1 42 34 65 29
[EMAIL PROTECTED] mailto:[EMAIL PROTECTED]
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] how to set the variable name in a loop
Hello,
I would like to have variables whose name are var1, var2, ... in a loop :
for (i in 1:10)
{
var(i) - i # where var(i) is var1,
}
Thanks in advance
Philippe
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] How to speed the shut down of the graphic device
hello,
I use R1.7.1 under winXP and I am running the following script example :
for (i in 1:10)
{
x - rnorm(100)
png( paste(D:/essai,i,.png,sep=))
plot(x)
t1 - Sys.time()
dev.off()
t2 - Sys.time()
print(t2-t1)
}
at each step, it takes about 3 seconds to shut down the graphic device.
I want to generate about one hundred of image and of course it takes too
much time. Is there any trick ?
Philippe
__
[EMAIL PROTECTED] mailing list
https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] symbol size on a plot
Hi, I would like to get from a plot the size of the symbols plotted. Imagine I have the following plot function : plot(1:2,1:2, pch=15, cex=4) I would like the get the values SIZE1 and SIZE2 so that if I plot the following rectangle : rect(1.5,1.5, 1.5+SIZE1, 1.5+SIZE2) then the size of this square is exactely the same as the one of the symbols that have been plotted. Thanks for any idea. -- -- Philippe Hupé Institut Curie - Equipe Bioinformatique 26, rue d'Ulm - 75005 PARIS France +33 (0)1 42 34 65 29 [EMAIL PROTECTED] mailto:[EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] How to get pixel position of a plot
Hi, I would like to plot a graph on the jpeg device for example and to write a table associated to this plot where I have the pixel coordonates of each plotted point so that I could include the jpeg image onto a html page and get all the information about each point when the mouse button is pressed. The indentify() can do it but on the window device... Thanks for any idea. Philippe -- -- Philippe Hupé Institut Curie - Equipe Bioinformatique 26, rue d'Ulm - 75005 PARIS France +33 (0)1 42 34 65 29 [EMAIL PROTECTED] mailto:[EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Markov chain simulation
Hi, Does anybody know a function to simulate a Markov chain given a probability transition matrix and an initial state ? Thanks. Philippe -- -- Philippe Hupé Institut Curie - Equipe Bioinformatique 26, rue d'Ulm - 75005 PARIS France +33 (0)1 42 34 65 29 [EMAIL PROTECTED] mailto:[EMAIL PROTECTED] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Rcmd check problem
I have the following message error when I try to buid my package : C:\Documents and Settings\PhilippeRcmd check mypackage * checking for working latex ...Error: environment variable TMPDIR not set (or set to unusable value) and no default available. at C:\PROGRA~1\R\rw1071\share\perl/R/Utils.pm line 165 Does someone have any idea ? Thanks OS Windows XP miktex R 1.7.1 ActivePerle 5.8.0 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Fortran compilation problem
Hello, I have problem when i am trying to install a package containing fortran source code. I give below the result of installation for one package (I have always the same message : /usr/bin/ld: can not find -lreadline): g77 -mieee-fp -fPIC -g -O2 -c akima.new.f -o akima.new.o g77 -mieee-fp -fPIC -g -O2 -c idbvip.f -o idbvip.o g77 -mieee-fp -fPIC -g -O2 -c idcldp.f -o idcldp.o g77 -mieee-fp -fPIC -g -O2 -c idgrid.f -o idgrid.o g77 -mieee-fp -fPIC -g -O2 -c idlctn.f -o idlctn.o g77 -mieee-fp -fPIC -g -O2 -c idpdrv.f -o idpdrv.o g77 -mieee-fp -fPIC -g -O2 -c idptip.f -o idptip.o g77 -mieee-fp -fPIC -g -O2 -c idptli.f -o idptli.o g77 -mieee-fp -fPIC -g -O2 -c idsfft.f -o idsfft.o g77 -mieee-fp -fPIC -g -O2 -c idtang.f -o idtang.o g77 -mieee-fp -fPIC -g -O2 -c idxchg.f -o idxchg.o g77 -mieee-fp -fPIC -g -O2 -c tripack.f -o tripack.o g77 -mieee-fp -fPIC -g -O2 -c ttidbs.f -o ttidbs.o gcc -shared -o akima.so akima.new.o idbvip.o idcldp.o idgrid.o idlctn.o idpdrv.o idptip.o idptli.o idsfft.o idtang.o idxchg.o tripack.o ttidbs.o -L/usr/lib/gcc-lib/i386-linux/2.95.4 -lreadline -ldl -lncurses -lg2c-pic -lm -L/usr/lib/R/bin -lR /usr/bin/ld: can not find -lreadline collect2: ld returned 1 exit status make: *** [akima.so] Erreur 1 ERROR: compilation failed for package 'akima' Thanks in advance. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
Re: [R] Fortran compilation problem
I have the 1.6.2 version installed on Debian Linux woody (deb packages have been downloading from CRAN) Le Dimanche 30 Mars 2003 20:04, [EMAIL PROTECTED] a écrit : You have not given us even minimal information (R version, platform, etc). My guess is that you installed R from an RPM, and you do not have (something like) readline-devel installed. However, this is easily fixed by editing R_HOME/etc/Makeconf and removing '-lreadline -ldl -lncurses' from the FLIBS macro. This has come up many times, so a search of the archives should be revealing. On Sun, 30 Mar 2003, Philippe Hupé wrote: Hello, I have problem when i am trying to install a package containing fortran source code. I give below the result of installation for one package (I have always the same message : /usr/bin/ld: can not find -lreadline): g77 -mieee-fp -fPIC -g -O2 -c akima.new.f -o akima.new.o g77 -mieee-fp -fPIC -g -O2 -c idbvip.f -o idbvip.o g77 -mieee-fp -fPIC -g -O2 -c idcldp.f -o idcldp.o g77 -mieee-fp -fPIC -g -O2 -c idgrid.f -o idgrid.o g77 -mieee-fp -fPIC -g -O2 -c idlctn.f -o idlctn.o g77 -mieee-fp -fPIC -g -O2 -c idpdrv.f -o idpdrv.o g77 -mieee-fp -fPIC -g -O2 -c idptip.f -o idptip.o g77 -mieee-fp -fPIC -g -O2 -c idptli.f -o idptli.o g77 -mieee-fp -fPIC -g -O2 -c idsfft.f -o idsfft.o g77 -mieee-fp -fPIC -g -O2 -c idtang.f -o idtang.o g77 -mieee-fp -fPIC -g -O2 -c idxchg.f -o idxchg.o g77 -mieee-fp -fPIC -g -O2 -c tripack.f -o tripack.o g77 -mieee-fp -fPIC -g -O2 -c ttidbs.f -o ttidbs.o gcc -shared -o akima.so akima.new.o idbvip.o idcldp.o idgrid.o idlctn.o idpdrv.o idptip.o idptli.o idsfft.o idtang.o idxchg.o tripack.o ttidbs.o -L/usr/lib/gcc-lib/i386-linux/2.95.4 -lreadline -ldl -lncurses -lg2c-pic -lm -L/usr/lib/R/bin -lR /usr/bin/ld: can not find -lreadline collect2: ld returned 1 exit status make: *** [akima.so] Erreur 1 ERROR: compilation failed for package 'akima' __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Spatial Clustering
Hello, Does anyone know any package which can perform spatial clustering ? Thanks. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] iterative reweighting in LOESS function
Hello, Is there an iterative reweighting option in LOESS function from modreg package to perform robust fitting in the presence of outliers in the data ? Thanks in advance. Philippe. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Median Regression Model
Hello, Is there any function to perform median regression model ? Thanks in advance. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
[R] Least-squares means
Is there any function to compute Least-squares means from a linear model as lsmeans does in SAS. Thanks. Philippe. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help
