[R] what are the limits on R array sizes?

2006-02-08 Thread Mike Miller
I have some computers with a massive amount of memory, and I have some jobs that could use very large matrix sizes. Can R handle matrices of larger than 2GB? If I were to create a matrix of 1,000,000 x 1,000, it would use about 8GB. Can R work with an array of that size if I have compiled R

Re: [R] open source and R

2005-11-13 Thread Mike Miller
On Sun, 13 Nov 2005, Roger Bivand wrote: On Sun, 13 Nov 2005, Robert wrote: If I do not know C or FORTRAN, how can I fully understand the package or possibly improve it? By learning enough to see whether that makes a difference for your purposes. Life is hard, but that's what makes life

Re: [R] problems with for: warnings and segfault

2005-11-11 Thread Mike Miller
On Fri, 11 Nov 2005, Ronaldo Reis-Jr. wrote: Segmentation fault This is a R bug or an error in my for function? All seg faults are bugs. Mike __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] How to find statistics like that.

2005-11-10 Thread Mike Miller
On Thu, 10 Nov 2005, Ruben Roa wrote: A statistic is any real-valued or vector-valued function whose domain includes the sample space of a random sample. The p-value is a real-valued function and its domain includes the sample space of a random sample. The p-value has a sampling

Re: [R] Is there no LU-Decomposition in R?

2005-11-09 Thread Mike Miller
On Wed, 9 Nov 2005, Christian Hinz wrote: I need the LU_Decomposition for my conversion of the DIRK/SDIRK (singly diagonally Runge Kutta implicit methods) algorithm. I think you need the Matrix package: http://cran.r-project.org/src/contrib/Descriptions/Matrix.html Mike

Re: [R] writing R shell scripts?

2005-11-09 Thread Mike Miller
On Wed, 9 Nov 2005, Henrik Bengtsson wrote: Put everything in curly brackets as my above example show. Ah ha! This leads to the latest version of my little one-liner (called doR): #!/bin/sh echo output - { $1 }; write.table(file=stdout(), .Last.value, row.names=FALSE, col.names=FALSE); q() |

[R] R-and-octave (was Element-by-element multiplication operator?)

2005-11-09 Thread Mike Miller
On Wed, 9 Nov 2005, Gabor Grothendieck wrote: See: http://cran.r-project.org/doc/contrib/R-and-octave-2.txt That is a really beautiful page. If the author is willing to update it to include some of the new Octave functionality -- specifically N-dimensional arrays -- I'm sure I can help

Re: [R] How to find statistics like that.

2005-11-09 Thread Mike Miller
On Wed, 9 Nov 2005, Gao Fay wrote: Hi there, Suppose mu is constant, and error is normally distributed with mean 0 and fixed variance s. I need to find a statistics that: Y_i = mu + beta1* I1_i beta2*I2_i + beta3*I1_i*I2_i + +error, where I_i is 1 Y_i is from group A, and 0 if Y_i is from

Re: [R] writing R shell scripts?

2005-11-08 Thread Mike Miller
Many thanks for the suggestions. Here is a related question: When I do things like this... echo matrix(rnorm(25*2),c(25,2)) | R --slave --no-save ...I get the desired result except that I would like to suppress the [,1] row and column labels so that only the values go to stdout. What is the

Re: [R] writing R shell scripts?

2005-11-08 Thread Mike Miller
On Wed, 9 Nov 2005, Henrik Bengtsson wrote: What you really want to do might be solved by write.table(), e.g. x - matrix(rnorm(25*2),c(25,2)); write.table(file=stdout(), x, row.names=FALSE, col.names=FALSE); Thanks. That does what I want. There is one remaining problem for my echo method.

Re: [R] writing R shell scripts?

2005-11-08 Thread Mike Miller
On Wed, 9 Nov 2005, Henrik Bengtsson wrote: A note of concern: When writing batch scripts like this, be explicit and use the print() statement. A counter example to compare echo 1; 2 | R --slave --no-save and echo print(1); print(2) | R --slave --no-save I guess you are saying

[R] writing R shell scripts?

2005-11-07 Thread Mike Miller
I'm new to the list. I've used R and S-PLUS a bit for about 15 years but am now working to make R my main program for all numerical and statistical computing. I also use Octave for this kind of work and I recommend it (it is also under the GPL). Here's my question: In Octave I can write