Re: [Rd] R-console vs. bash console (execution halted)

2008-03-11 Thread Roger Bivand
On Tue, 11 Mar 2008, Prof Brian Ripley wrote: > Your first use is non-interactive. I have no idea what > > R CMD source("/data/myscript.R") > > is supposed to do (it does not work for me), but source() when run in RGui > is interactive. > > See ?options and in particular 'error' and the reference

Re: [Rd] R-console vs. bash console (execution halted)

2008-03-11 Thread Prof Brian Ripley
Your first use is non-interactive. I have no idea what R CMD source("/data/myscript.R") is supposed to do (it does not work for me), but source() when run in RGui is interactive. See ?options and in particular 'error' and the reference to ?stop for the different ways errors are handled in int

[Rd] R-console vs. bash console (execution halted)

2008-03-11 Thread Lienert Christophe
Dear Everyone, I am using R 2.6.2 on my Redhat AS4. I installed the RPM offered on the cran website. I would like to use R from my bash console with the following command: /path/to/R --vanilla --slave --args < /data/myscript.R this script invokes a range of functions and tools, among others to

[Rd] Rtools and GCC4 problem

2008-03-11 Thread Joe Byers
>> I am trying to compile rseries from Whit Armstrong and a colleague of mine found a problem with using GCC4 I get the following error when compiling rseries g++-sjlj -Ic:/R/include -O2 -Wall -c Rutils.cpp -o Rutils.o Rutils.cpp: In function 'double* getColPointer(SEXPREC*, int)': Rutils

Re: [Rd] NA warnings for r() {aka "patch for random.c"}

2008-03-11 Thread Berwin A Turlach
G'day Martin, On Tue, 11 Mar 2008 18:07:35 +0100 Martin Maechler <[EMAIL PROTECTED]> wrote: > > "BAT" == Berwin A Turlach <[EMAIL PROTECTED]> > > on Tue, 11 Mar 2008 13:19:46 +0800 writes: [...] > BAT> The first two lines give identical results, as one could > BAT> reasonably

Re: [Rd] NA warnings for r() {aka "patch for random.c"}

2008-03-11 Thread Martin Maechler
> "BAT" == Berwin A Turlach <[EMAIL PROTECTED]> > on Tue, 11 Mar 2008 13:19:46 +0800 writes: BAT> G'day Martin and others, BAT> On Mon, 10 Mar 2008 12:06:01 +0100 BAT> Martin Maechler <[EMAIL PROTECTED]> wrote: >> > "BAT" == Berwin A Turlach <[EMAIL PROTECTED]>

Re: [Rd] write.table with row.names=FALSE unnecessarily slow?

2008-03-11 Thread Martin Maechler
MartinMo> write.table with large data frames takes quite a long time MartinMo> system.time({ MartinMo> + write.table(df, '/tmp/dftest.txt', row.names=FALSE) MartinMo> + }, gcFirst=TRUE) MartinMo> user system elapsed MartinMo> 97.302 1.532 98.837 MartinMo> A r

Re: [Rd] crossprod is slower than t(AA)%*BB

2008-03-11 Thread Prof Brian Ripley
Please, you can't compare calls starting with different memory settings (and the garbage collector does adjust to recent usage). Try them in the other order or in parallel sessions. >From the times you quote it looks to me as if you would benefit from using an optimized BLAS. Here are some res

[Rd] Wishlist: Factor correlations in factanal (PR#10931)

2008-03-11 Thread uhkeller
Full_Name: Ulrich Keller Version: 2.6.2 OS: Ubuntu 7.10 Submission from: (NULL) (158.64.77.190) Most statistical packages report factor correlations for oblique factor rotations. R's factanal() does not. John Fox posted some modifications to R-devel back in 2005 that implement this, but unfortuna

[Rd] crossprod is slower than t(AA)%*BB

2008-03-11 Thread Ole Fredslund Christensen
Dear Rdevelopers The background for this email is that I was helping a PhD student to improve the speed of her R code. I suggested to replace calls like t(AA)%*% BB by crossprod(AA,BB) since I expected this to be faster. The surprising result to me was that this change actually made her code slo

Re: [Rd] write.table with row.names=FALSE unnecessarily slow?

2008-03-11 Thread Prof Brian Ripley
This is a pretty extreme case: why not use write() to write a single column? (It's a bit faster than your patched timing.) In a more realistic test of 10 columns of 1 million rows I see a speedup from 12.2 to 9.7 seconds. So I'll add the patch, but think that significant speedups will be quite