similar to: extract numerical variables from a data frame

Hi, I just encountered a problem in R that may easily be fixed: If one uses attach for a data.frame e.g. 10000 times and forgets detach, then R gets incredibly slow (less then 10% of the original speed). My system: platform powerpc-apple-darwin6.0 arch powerpc os darwin6.0 system powerpc, darwin6.0 status major 1 minor 6.1 year 2002

I tried to define a function like: fnx <- function(x, by.vars=Month) print(by(x, by.vars, summary)) But this doesn't work (does not find x$Month; unlike other functions, such as subset(), the INDICES argument to "by" does not look for variables in dataset x. Is fully documented, but I forget every time). So I tried using "with": fnxx <- function(x, by.vars=Month)

Hello, First, sorry for my poor english, I will try to be understood. It's the first time I try this "r-help mailing list" and I hope it will be a success. I am working on heteroscedasticity analysis. I would like to get the "Box-Ljung" and the "Lagrange multipliers" test. I found the first one in the library "ts", but I can't find the second one.

Hello Anyone that now of any function in R that can calculate copulas? Or if anyone have any code avaible I would be more than interested. Thank you in advance /Thomas ______________________________________________ R-help at stat.math.ethz.ch mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Greetings everyone, I am new at R. My questions is rather basic. Looking R manuals looks like there should be a way to read MS excel files into R. Could somebody tell me which library should I use for that? Thanks in advance Victor H. Marín Laboratorio de Modelación Ecológica Depto. de Ciencias Ecológicas, Fac. de Ciencias, Universidad de Chile Casilla 653 Santiago, Chile

R version 1.8.1, OS Windows 98 Dear colleagues, if I import vegetation data (first row with column labels and first column with row labels) like 7MYRGERM;7AGRGIGA;7DRYOCTO;5MYRGERM;7SALELEA;7CHOCHON;7SALNIG?;....... t401;5;2;2;3;4;2;2;2;1;2;1;2;2;1;2;2;2;1;2;1;0;0;...... t403;3;0;0;6;4;0;3;0;0;3;0;0;0;0;3;0;0;0;2;0;2;0;..... with read.table("data.file", header=TRUE,

Hello, I am planning a study with the main point to evaluate the interaction of two treatments, but for ethical reasons one cell is empty, that with patients receaving no treatment at all Treatment B

Dear all, I have found the following (for me) incomprehensible behaviour of ISOdate (POSIXct): > ISOdate(1900,6,16) [1] "1900-06-15 14:00:00 Westeurop?ische Sommerzeit" > ISOdate(1950,6,16) [1] "1950-06-16 14:00:00 Westeurop?ische Sommerzeit" Note that in the first case I get the 15th of June back, not the 16th as I would have expected! This happened under R-1.7.1 on

bug.report() tells me to email to r-bugs@r-project.org, whereas the Web site http://www.r-project.org/ points me to r-bugs@biostat.ku.dk. Which should I believe? Simon Fear Senior Statistician Syne qua non Ltd Tel: +44 (0) 1379 644449 Fax: +44 (0) 1379 644445 email: Simon.Fear@synequanon.com web: http://www.synequanon.com Number of attachments included with this message: 0 This

Dear R users, I want to retrieve "..." argument values within a function. Here is a small exmaple: myfunc <- function(x, ...) { if (hasArg(ylim)) a <- ylim plot(x, ...) } x <- rnorm(100) myfunc(x, ylim=c(-0.5, 0.5)) Error in myfunc(x, ylim = c(-0.5, 0.5)) : Object "ylim" not found > I need to retrieve values of "ylim" (if it is defined

I've searched on CRAN for axes, axis, and other terms I've already forgotten, without (re)discovering the reason for S using "non-joining" axes by default, instead of box("l"). MASS points me towards Cleveland (1993) but I don't have ready access to this any more. Could someone give me a one-liner to justify this choice to a sceptic? It's something to do

Hello, Pass by reference appears to be a topic which comes up from time to time, but I wasn't able to find something in the R-help archives which tells how to accomplish it. I have a problem that you may have seen before -- R runs out of memory when processing large matrices. Part of the problem for me is that I am using some large matrices as function arguments, and these are modified,

Is there a good way in R to impute values which exist, but are less than the detection level for an assay? Thanks, Jonathan Williams OPTIMA Radcliffe Infirmary Woodstock Road OXFORD OX2 6HE Tel +1865 (2)24356

Dear all can anybady tell me how to seperate R's input and output screen?I mean just like SAS or some others, the commands does not mixed with results. thank you

Hi I have a huge array with series of data. For each cell in the array I fit a linear model, either using lm() or gls() with lm() there is no problem, but with gls() I get an error: Error in glsEstimate(glsSt, control = glsEstControl) : computed gls fit is singular, rank 2 as soon as there are data like this: > y1 <- c(0,0,0,0) > x1 <- c(0,1,1.3,0) > gls(y1~x1)

Hello, I've a data frame with 15 colums and 6000 rows, and I need the data in a single vector of size 90000 for ttest. Is there such a conversion function in R, or would I have to write my own loop over the colums? thanks for your help + kind regards Arne

Hi everyone. look at this: x <- 1:4 x[c(1,2,1)] <- (1:3) print(x[1]) I get 3, but isn't NA more appropriate? [1 would be as sensible]. FWIW, the equivalent Fortran 95 statement is flagged as an error. R-intro, section 2.7, says that in such cases the assignment is carried out "in order" which might support getting 3. To my way of thinking, the concept of "in

Bonjour, Je suis ?tudiant stagiaire ? Paris et je rencontre quelques difficult?s en programmation R. J'ai une data frame compos?e de 4 colonnes et 250 lignes et dont chaque ligne est une famille. J'ai fait un tirage al?atoire avec remise des familles 250 fois ce qui m'am?ne ? une nouvelle dataframe. A cette nouvelle dataframe, j'applique un programme qui calcule 2 param?tre X1

what does this mean in R-1.8.1 release notes? o median() no longer `works' for odd-length factor variables.

does anyone know how to use the apply function on a two-variable t-test? i've tried everything. apply(data[1:4],1,t.test) # works, but only a one-variable test apply((data[1:4],data[5:9]),1,t.test) # returns a syntax error i've also tried using the for loop: for (i in 1:1000) t.test(data[i,1:4],data[i,5:9]) doesn't print anything! it processes for awhile, so i know it's doing