similar to: scan

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 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("

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

I've got a PhD student starting this year. She'll be working on data mining. She has asked for a reading list while she's still in her home country, which is a really good sign. Her cosupervisor and I don't (think we) have any problem with choosing things that are specifically about data mining, but there are some statistical ideas (sampling, exploratory-vs-confirmatory,

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, I have a problem with 'grep' I read several lines from a file p <- readLines(pfile,n=8) now I have to check that a key word "JABBERWOCKY" is not in that segment a <- grep("JABBERWOCKY",p) which gives a 1 if the key is in the text or 'numeric(0) how can I test it if(a==1) gives Ok in one case and error:missing value where logical needed in the other.

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"

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 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

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)

Colleagues ---------------------------------- System info: R version rw1020 on NT ESS using emacs ver. 20.4 ---------------------------------- To keep my scripts short and tidy, I want to call sub-functions and return variables from these for subsequent use in the calling function. e.g. "test.smooth.opc.time.series" <- function() { load.opc.ascii.files() }

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

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

Windows binaries at CRAN (bin/windows/Windows-NT/base) have been updated to R-1.1.0. See below for a list of Windows-specific changes. We thanks all the people who checked over pre-test versions. guido masarotto (for the R-core team) Windows-specific changes to R ============================= There is now a GUI preferences editor on the Edit menu in Rgui. A data entry editor is now

Windows binaries at CRAN (bin/windows/Windows-NT/base) have been updated to R-1.1.0. See below for a list of Windows-specific changes. We thanks all the people who checked over pre-test versions. guido masarotto (for the R-core team) Windows-specific changes to R ============================= There is now a GUI preferences editor on the Edit menu in Rgui. A data entry editor is now

Hello, I have some problems with RODBC and Excel in Win98 I am using R 1.6.2 and just upgraded RODBC to the last version on CRAN. I have an Excel file with columns Number, Name, Sex, Age, FEV1 on Sheet 1 and Number, Age, FEV1, Name, Sex on Sheet 2. Now I open the channel to the file > chan1 <- odbcConnectExcel("c:/testOdbc.xls") > tables(chan1) and the list appears with the 2

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