similar to: ggplot2 graphing multiple lines of data

Does something like this get you close: x <- list() keys <- LETTERS[1:6] # create for (i in keys){ x[[i]] <- data.frame(a=1:5, b=1:5, c=1:5) } # output output <- file('tempxx.txt', 'w') for (i in keys){ write.table(i, row.names=FALSE, col.names=FALSE, file=output, quote=FALSE) write.table(x[[i]], file=output, quote=FALSE) } close(output) On Wed, Jul 2,

Dear List, This is just a programming problem which i cannot seem to figure out. I am trying to get a set of power from a test (say, kolmogorov smirnov) out of a distribution (say, G-K distribution) as follows. I am trying to reduce to pain of writing the whole set of data points (p# below) using "for" loop. However, I seem to have some problem in it as the output "M" does not

Hi all, I would like to set up an iterative read & write sequence to avoid reading and writing each file one at a time. Hundreds of data sets to re-calculate.? The code I have works well individually, but would like to set up an iterative read, calculate and write changing the input and output file names each iteration. I? think I have read that there is an R? feature using

Hello, everyone I have a piece of code that looks like this: mrets <- merge(mrets, BMM.SR=apply(mrets, 1, MyFunc, ret="BMM.AV120", stdev="BMM.SD120")) mrets <- merge(mrets, GM1.SR=apply(mrets, 1, MyFunc, ret="GM1.AV120", stdev="GM1.SD120")) mrets <- merge(mrets, IYC.SR=apply(mrets, 1, MyFunc, ret="IYC.AV120",

Hi Bob, are these performance regressions real? They look pretty serious. Ciao, Duncan. On 10/12/11 09:40, llvm-testresults at cs.uiuc.edu wrote: > > bwilson__llvm-gcc_PROD__i386 nightly tester results > > URL http://llvm.org/perf/db_default/simple/nts/332/ > Nickname bwilson__llvm-gcc_PROD__i386:4 > Name curlew.apple.com > > Run ID Order Start Time End Time >

While trying to train randomForest with my dataset, I am ending up with the following error Error in randomForest.default(datatrain, classtrain) : length of response must be the same as predictors My data looks like: A,B,C,D,Class 1,2,1,2,cl1 1,2,1,2,cl1 3,2,1,2,cl2 3,2,1,2,cl2 3,2,1,2,cl2 3,2,1,2,cl2 3,2,1,2,cl2 3,2,1,2,cl2 3,2,1,2,cl2 3,2,12,3,cl2 3,2,1,2,cl2 Actual dataset has around 4000

Yes, they are real. I re-ran the two tests with the biggest execution time regressions, and the results were completely reproducible. On Oct 12, 2011, at 1:24 AM, Duncan Sands wrote: > Hi Bob, are these performance regressions real? They look pretty serious. > > Ciao, Duncan. > > On 10/12/11 09:40, llvm-testresults at cs.uiuc.edu wrote: >> >>

Am I missing something obvious on how to draw multi-line plots in base graphics? In ggplot2, I can do: data(Oxboys, package = "nlme") library(ggplot2) qplot(age, height, data = Oxboys, geom = "line", group = Subject) But in base graphics, the best I can come up with is this: with(Oxboys, plot(age, height, type = "n")) lapply(split(Oxboys[c("age",

Dear all, I've got a problem in including square variables in lme and lmlist functions. I've tried to work on Oxboys data of Pinheiro and Bates'book, which consist of the heights of 26 boys, each mesured on nine different occasions : > Oxboys Grouped Data: height ~ age | Subject Subject age height Occasion 1 1 -1.0000 140.50 1 2 1 -0.7479 143.40

>> If you're looking to share the tape drive so you ................ >You could probably do something clever with a magic >script in the Samba config, but it still wouldn't >work with native Windows backup utilities. (Are >there native Windows backup utilities that can >access a remote tape drive?) You could probably get >it to work with tar under Cygwin if

Dear Martin, Have you try to create a new variable for age squared, say agesq? If you fit the model using this new variable you should get the coefficients. So your new model is something like height~age+agesq I hope this helps, Saghir > -----Original Message----- > From: MARTIN Ludovic [SMTP:martinl@mathinfo.ens.univ-reims.fr] > Sent: Wednesday, 07 May, 2003 2:02 PM > To:

Dear R-users, I recently began using the ggplot2 package and I am still in the process of getting used to it. My goal would be to plot on the same grid a number of curves derived from two distinct datasets. The first dataset (called molten.data) looks like this : Column names : Perc, Week, Weight P10 21 333.3554 P90 21 486.0480 P10 22 452.6347 P90 22 563.8263 P10 23 575.0960 P90

Dear All, I've got a problem with ks.test. I've two realy large vectors, that I'd like to test, but I get an overflow, and the p-value cannot be calculated: > length(genomesv) [1] 390025 > length(scopv) [1] 140002 > ks.test(genomesv, scopv) Two-sample Kolmogorov-Smirnov test data: genomesv and scopv D = 0.2081, p-value = NA alternative hypothesis: two.sided

Full_Name: Vaidotas Zemlys Version: 2.9.0 OS: Ubuntu 8.10 Submission from: (NULL) (213.197.173.50) Steps to reproduce the bug: library(nlme) data(Oxboys) qm=lmList(height~age|Subject,data=Oxboys) grid=with(Oxboys,expand.grid(age=seq(min(age),max(age),length=50),Subject=levels(Subject))) res <- predict(qm,grid,se=TRUE) Erreur dans if (pool) { : l'argument est de longueur nulle res <-

Hello experts I have a thing which I can't vectorize. Can any gurus out there help? This is my truncated dataset; it is a dataframe of five possums and which tree they slept in on three days. R> dataset p14 p88 p82 p90 p98 1 6 7 5 1 3 2 8 9 10 1 11 3 8 7 12 13 14 I can manipulate it with simple functions fine: R> f function(n){n+1000} R> f(dataset)

Hello All, I'd like to create diagrams for dynamic web-pages via a cgi-script that itself uses R. Sadly, I found out that the png and jpeg devices need an X-server :-( , so I cannot use these nice devices. Does any one have any experience using R for on-the-flycreation of graphics for web-pages pages? Any suggestions for workarounds? thanks very much for help, Arne -- Arne Mueller

Hello, I've just started R, and I'm getting a bit mad using it. I've managed to produce a barplot with the labels for the ytick marks placed horizontally (perpendicular to the y-axis) usiing par(las=1). The problem is that most of my labels are in part beyond the plotting area because they are rather long (e.g. "H. sapiens", "D. melanogaster" ...). What is the

I have a TDM400P with some FXO ports, and I wanted to connect the two POTS lines from my Pipeline-75 ISDN router into the FXO interfaces on my Asterisk server. Hooked it up, seemed fine, called in and it answered. The problem is when the call is hung up on, the FXO port never drops. So of course then the P75 just holds the line off hook and you get a busy. So it's good for the first

Dear All, I'm getting confused with the concept R uses to do regression using lm. I'm afmiliar with gnuplot and the build-in fit command, but couldn't get R's lm to work on my data. I know that my data follows a powerlaw or maybe an exponential function, and I'd like to determine the best fitting factors for the formula: a*x^b where b < 0. I've tried thge follwoing:

Hello, I have just installed a fresh copy of FreeBSD 4.0-RELEASE on a P75, and want to compile openssh-2.1.1p2, but the compile fails with: gcc -g -O2 -Wall -I/usr/local/ssl/include -DETCDIR=\"/usr/local/etc\" -DSSH_PROGRAM=\"/usr/local/bin/ssh\" -DSSH_ASKPASS_DEFAULT=\"/usr/local/libexec/ssh/ssh-askpass\" -DHAVE_CONFIG_H -c loginrec.c loginrec.c: In function