similar to: unknown compression

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

Hello, I have a data frame with the following variables: ID, X1,X2,X3,X4,X5,Y1,Y2,Y3,Y4,Y5 and some other that do not matter, some of the X and Y can be missing (NA). I want to compute the slope of the linear regression Y ~ X for each subject, so using apply(DF,1,FUN,ra.rm=TRUE) now How do I define FUN? The X are different for each subject. Thanks for any help R.Heberto Ghezzo Ph.D. Montreal -

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

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 have a big data frame where consecutive time dates and corresponding observed values for each subject (ID) are on a line. I want to compute the linear slope for each subject. I would like to use apply but I do not know how to express the corresponding function. An example using a loop follows # # create dummy data set There are missing values a <- c(1,2,3,4, 1,1,1,1, 2,2,3,3,

hello, i tried to run playwith but : > library(playwith) Loading required package: lattice Loading required package: cairoDevice Loading required package: gWidgetsRGtk2 Loading required package: gWidgets Error in inDL(x, as.logical(local), as.logical(now), ...) : unable to load shared object 'H:/R/cran/RGtk2/libs/i386/RGtk2.dll': LoadLibrary failure: The specified procedure could

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,

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 , I do not know anything abount Ubunto, but I found a Portable Ubunto for Windows and since so many people prefer Linux to Windows I decided to give it a try. It runs very nicely, so I tried to load R, following Instructions in CRAN I added the line deb http://probability.ca/cran/bin/linux/ubuntu hardy/ to /etc/apt/sources.list and then from a console I did sudo apt-get update sudo apt-get

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

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)

Hello, I never used Sweave before but now I try and got: > > rnwfile <- system.file("Sweave", "example-1.Rnw", package = "utils") > Sweave(rnwfile) Writing to file example-1.tex Processing code chunks with options ... 1 : echo term verbatim 2 : term verbatim pdf You can now run (pdf)latex on 'example-1.tex' >

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, if I do: x <- c(0.5,1,3,6,10,20,40) y <- 10-log(x)+rnorm(7,0,0.05) r1 <- lm(y ~ log(x)) plot(log(x),y) abline(r1) # I get a nice plot with the regression line almost over the points. but: plot(x,y,log="x") abline(r1) gives me exactly the same plot for the points but the regression line is completely off ! I would like the plot with the real values of X on

Hello, Working with Windows 7 in a HP laptop with R-2.10.1 I download and installed ImageMagick-6.3.7.7-Q16-Windows-dll.exe and GTK 2.12.9-win32-2, then downloaded and installed from local file EBImage_3.2.0.zip and I got: > library(EBImage) Loading required package: abind Error in inDL(x, as.logical(local), as.logical(now), ...) : unable to load shared library