Re: [R] reading very large files

2007-02-04 Thread juli g. pausas
i386, mingw32 status major 2 minor 4.1 year 2006 month 12 day18 svn rev40228 language R version.string R version 2.4.1 (2006-12-18) Juli On 03/02/07, Marc Schwartz [EMAIL PROTECTED] wrote: On Sat, 2007-02-03 at 19:06 +0100, juli g

Re: [R] reading very large files

2007-02-04 Thread juli g. pausas
that was only 1/4th done after 25 h! Many thanks Juli On 02/02/07, juli g. pausas [EMAIL PROTECTED] wrote: Hi all, I have a large file (1.8 GB) with 900,000 lines that I would like to read. Each line is a string characters. Specifically I would like to randomly select 3000 lines. For smaller

[R] reading very large files

2007-02-02 Thread juli g. pausas
Hi all, I have a large file (1.8 GB) with 900,000 lines that I would like to read. Each line is a string characters. Specifically I would like to randomly select 3000 lines. For smaller files, what I'm doing is: trs - scan(myfile, what= character(), sep = \n) trs- trs[sample(length(trs), 3000)]

[R] repeated measures with random effects

2005-10-04 Thread juli g. pausas
: aov(H ~ B * A %in% B * Time + Error(id) ) where id is the factor coded for the repeated subjects. But my question is how to include the random factor C ? Any help would be appreciated. Thanks Juli -- Juli G. Pausas CEAM UA http://www.ceam.es/lass/pausas.htm

[R] APE: compar.gee( )

2004-02-08 Thread juli g. pausas
Dear all, I don't understand the following behaviour: Running compar.gee (in library ape ) with and without the option 'data', it give me different results Example: Start R load(eiber.RData) ls() [1] gee.na mydata mytree library(ape) # runnig with the option data= mydata

[R] subset select within a function

2004-01-21 Thread juli g. pausas
Dear all, I'd like to subset a df within a function, and use select for choosing the variable. Something like (simplified example): mydf - data.frame(a= 0:9, b= 10:19) ttt - function(vv) { tmpdf - subset(mydf, select= vv) mean(tmpdf$vv) } ttt(mydf$b) But this is not the correct way. Any

[R] barplot plot together

2003-12-18 Thread juli g. pausas
for setting the 2 plots at the same scale? I have unsuccessfully tried: plot(t, data, type=b, xlim=c(0,7)) plot(t, data, type=b, xlim=c(min(t),max(t))) (R 1.8.1, for Windows) Thanks Juli -- Juli G. Pausas Centro de Estudios Ambientales del Mediterraneo (CEAM) C/ Charles R. Darwin 14, Parc

Re: [R] barplot plot together

2003-12-18 Thread juli g. pausas
), axes=F, xlab=, ylab=) axis(4); mtext(Number of fires, side=4, line=3, col=1) legend(2, 155, Number, lty= 1, pch= 19, bty=n, cex=0.8) legend(2, 163, Area, fill= 2, bty=n, cex=0.8) Cheers Juli Martin Maechler wrote: juli == juli g pausas [EMAIL PROTECTED] on Thu, 18 Dec 2003 10:54:08 +0100

[R] NA, deleting rows

2003-12-18 Thread juli g. pausas
Dear colleges, I do not understand the following behaviour: aa - data.frame(a1= 1:10, a2= c(rep(NA, 5), 1:5) ) aa a1 a2 1 1 NA 2 2 NA 3 3 NA 4 4 NA 5 5 NA 6 6 1 7 7 2 8 8 3 9 9 4 10 10 5 aa[!aa$a2==1, ] # removing rows with a2==1 a1 a2 NA NA NA NA.1 NA NA NA.2 NA NA

[R] the name of a variable in a function

2003-09-17 Thread juli g. pausas
Dear collegues, How can I get the name of a variable (and not the variable) within a function ? For instance, in the following function, I'd like to create a variable in the dataframe df with the same name to the variable passed in var: prova - function( var ) { df -

[R] plots to metafile and x/y ratio

2003-03-15 Thread juli g. pausas
Dear all, I have a problem which I'm not sure if it is due to my machine or if I'm doing something wrong. I'm ploting a little map using plot(val, type=l, asp=1) where val is the object with the latitude/longitude data. In the screen the figure is perfect, and I can copy to the clipboard and

[R] any european maps

2003-03-14 Thread juli g. pausas
Dear all, I've seen that there are some maps, or at least costlines in R. The oz package is the map of Autralia and in the fields package there is a US map and a world map. This world map allows to select smaller sections such Europe: library(fields) world(xlim=c(-10,18),ylim=c(36,60))

[R] variables in a correlation matrix

2003-03-12 Thread juli g. pausas
Hi, Given a correlation matrix, how can I know which variables have certain correlation values? aa - matrix( rnorm(25, 1,1), 5,5) colnames(aa) - c(v1, v2, v3, v4, v5) aacor - cor(aa) For instance, which variables are negatively correlated? aacor[aacor0] # provide the r values, but how I get

[R] hist() basic question

2003-03-08 Thread juli g. pausas
Hi, This is a very basic question, but I would like to undestand hist(). I thought that the hist( , freq=FALSE) should provide the relative frequencies (probabilities), and so they should sum 1, however: set.seed(2) ah - hist(rnorm(100), freq=F) sum(ah$intensities) [1] 2 set.seed(2) bh -