similar to: grep

I received several answer to my querry re reading zipped files, R.Gentleman, Liaw Andy, BD Ripley,J Holtman, A.Matt, and others, sorry if I missed your names. a- use pipe to run external unzip to a file and read from there b- use gzfile() to open directly a gzipped file c- use zip.file.extract() solution b works like a charm, unfortunately it open gzipped files and not zipped files, would it be

Hello, I have a problem. I have to read a HUGE file which has to be line processed, so I would like to use scan like p <- scan(pfile,what=c(0,0,0,"",0,""),nlines=1) but it continues to read the first line, so I have to add a skip and increment the skip after each read. It takes forever to read a test file of 1 Mb, the real one is in 60's Mb Is there a way to read line

Hi, I got a code from S that uses 'nlminb' to minimize a function with constraints. Is there a similar function in R? Thanks. R. Heberto Ghezzo Ph.D. Meakins-Christie Labs McGill University Montreal - Canada heberto at meakins.lan.mcgill.ca -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello, I wrote the following function to compute multiple comparisons in a one way anova and randomized blocks anova. aov1 <- function(y,g,s=NULL,comp="mca",meth="Sidak") { # fun <- function(x) c(mean(x,na.rm=T),sd(x,na.rm=T),length(x[!is.na(x)])) # li <- length(unique(g)) cat(" Analysis of Variance with Multiple comparisons\n\n") cat("

Hello, I think this qualify as a bug > x<-c(1,2,3,4,NA,6,7) > mean(x) [1] NA > mean(x,na.rm=T) [1] 3.833333 > sd(x) Error in var(as.vector(x)) : missing observations in cov/cor > sd(x,na.rm=T) Error in sd(x, na.rm = T) : unused argument(s) (na.rm ...) > var(x) Error in var(x) : missing observations in cov/cor > var(x,na.rm=T) [1] 5.366667 > why sd() does not

Hello everybody I am trying to extract some observations from a data frame, ie the subjects that belong to a given group, and although t-test etc work if I try to obtain the number of subjects in the subgroup i get some funny numbers. All the subjects with NA for group are included in the subgroups 'group==1' etc. Is this a bug? for example in Windoze R : Copyright 2001, The R

Hi I am using R65.1 in Windows 95 I have a CSV file from Excell > a<-read.table("c:/heberto/mgc/tst.csv",header=T,sep=",") > attach(a) > a manolo fvcpp fevpp fvvcpp tlcpp rvpp rvtlpp plmaxpp 1 1 99.28 97.67 98.38 91.14 102.9 111.25 117.64 2 1 86.97 68.56 78.89 94.60 112.34 118.53 159.20 3 1 81.12 71.76 88.37 89.16

Hi, Does somebody knows of the existence of a library of functions for solving differential equations, I need a solver for a system of non linear first order differential equations. I think there was something in IMSL but that was a long time ago and probably there are some new algorithms in the new now. Thanks for any pointers. . R. Heberto Ghezzo Ph.D. Meakins-Christie Labs McGill

Hello, another problem. > x<-rep(1,10) > y<-rep(c(1,2),c(5,5)) > z<-seq(1:10) > ab<-data.frame(x,y,z) # now I want to do some work by the value of 'y' > by(ab,y,mean) y: 1 x y z 1 1 3 ------------------------------------------------------------ y: 2 x y z 1 2 8 # I do not want all the means, only the mean of 'z' > by(ab,y,function(x)

I would like to obtain a crosstabulation of means(var, quantiles...) i.e. I have a data frame with Var-i, Var-j, Var-k, Var-X, var-Y I like to have the mean of Var-X for each combination of Var-i,Var-j. One solution would be: by(var-i,Var-j,mean(Var-x)) but I would like it better formatted and with mean,S.Dev,n for each cell? Does anybody have some function to do this or some ideas how to go

Hi, I am relatively new to R programming. I have Guido's version on a Win95 Pentium with 32 Mb RAM. I tried to compile the Package MIXREG of Rolf Turner in S to R almost everything works OK but he has a function po = partial outer product as follows: po <- function(a,b) { array(apply(b,2,"*",a),dim=c(dim(a),dim(b)[2]))) } in S-Plus given x = ( 1 1 1 2 1 3 1 4

Hi R experts, Can somebody suggest what can be wrong? PIII with 128 Mb RAM Using R 1.1 in Win98, the latest version downloaded yesterday. I run update.packages so 'survival', ' mass', and 'nlme' are the latest versions in CRAN > library(survival) Error in dyn.load(x, as.logical(local), as.logical(now)) : unable to load shared library

hello. I am trying to read a big file with different sections each with a different format, actually it is a map in format *.e00 of MapView. I can read the whole file with no problems from R using scan() but it takes too long, some files are 50 meg and have to be read line by line to catch the section codes. So if I do the reading from C would it be faster? if it is, how do I manage with the

Hello, I am trying to modify some functions in an existing package. I can get a function by just typing the name in the console, but how can I get all the functions in a package? the ./R/ subdirectory has files xx.rdb and xx.rdx which I can not decompress with 7-zip or with unzip, unrar etc. The help files are compressed with gzip in the ./man/ folder and can be unzipped easily. Thanks for your

Hi R friends, I havent asked a silly question in a long time so here it is: Reading the last issue of Stat Can J. there is an article in a single degree of freedom test for non aditivity of interactions in anova tables with 1 obs per cell. The authors claim it is more powerfull than Tukey but for the case of multiplicative interaction, which is the alternative studied by Tukey. I tried to

Hello, I think the following qualify as a bug given: x<-c(1,2,3,4,2,4,2,2,4,2,2,2,4,3,2,1) z<-embed(x,4) z [,1] [,2] [,3] [,4] [1,] 4 3 2 1 [2,] 2 4 3 2 [3,] 4 2 4 3 [4,] 2 4 2 4 [5,] 2 2 4 2 [6,] 4 2 2 4 [7,] 2 4 2 2 [8,] 2 2 4 2 [9,] 2 2 2 4 [10,] 4 2 2 2

Hello, I think the following qualify as a bug given: x<-c(1,2,3,4,2,4,2,2,4,2,2,2,4,3,2,1) z<-embed(x,4) z [,1] [,2] [,3] [,4] [1,] 4 3 2 1 [2,] 2 4 3 2 [3,] 4 2 4 3 [4,] 2 4 2 4 [5,] 2 2 4 2 [6,] 4 2 2 4 [7,] 2 4 2 2 [8,] 2 2 4 2 [9,] 2 2 2 4 [10,] 4 2 2 2

I am using R 2.0.2 on a WinXP I am trying to get the code of the Kruskal-Wallis test but > kruskal.test function (x, ...) UseMethod("kruskal.test") <environment: namespace:stats> > ls(3) [1] "acf" "acf2AR" "add.scope" .............. [181] "kruskal.test" "ks.test"

Hello, Using R 1051 in win98 When I try to build a library Addreg with DESCRIPTION, Index and Title in R and subdirectories Man and R. From a Command window I do: c:\> c:\rw1051\bin\Rcmd INSTALL c:\r_src\src\library\addreg --library=c:\r\src\library make: Entering directory '/cygdrive/c/rw1051/src/gnuwin32' make DLLNM= \ -C c:/r_src/src/library/addreg PKG=addreg RHOME=c:/rw1051

hello, I am using R1.2.2 on win98 with 64 meg of RAM I try to read a big csv file and get the error : a <- read.csv("c:/all2.csv",header=T) Error: vector memory exhausted (limit reached?) In addition: Warning message: Reached total allocation of 16Mb: see help(memory.size) Lost warning messages from there on, nothing works File-Exit, ^F4, q() , top right cross, the same message