similar to: Overlaying graphs from different datasets with ggplot

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)

For a given random variable rv, for instance, rv = rnorm(1000), I plot its density curve and calculate some quantiles: plot(density(rv)) P10P50P90 = = quantile(rv,probs = c(10,50,90)/100) I would like to color the area between P10 and P90 and under the curve and mark the P50 on the curve. > rv = rnorm(1000) > plot(density(rv)) > P10P50P90 = = quantile(rv,probs = c(10,50,90)/100) Could

Hi, I have a dataset called 'fish'. fish$Species returns extract 1. When I use fish$Species != c("CRA","PHC"), i.e. I want all species except "CRA" and "PHC", I get extract 2 which is blatantly wrong. Can anyone see what I'm doing wrong? Regards, James EXTRACT 1 > fish$Species [1] ALB ALB ALB ALB ALB ALB ALB ALB ALB ALB ALB ALB ALB

I have written a few different summary functions. I want to calculate the statistics by groups and I am having trouble getting the output as a dataframe. I have attached one example with a small dataset that calculates summary stats and percentiles, I have others that calculate upper confidence limits etc. I would like the output to be converted to a dataframe with one of the columns as the

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

We're using Samba 1.9.18p10 on a Linux 2.0.35 box as a fileserver chiefly for software development under VC++. Builds are done from the command line, not the IDE. I'm seeing rather slow build performance under NT4SP3. Build performance is fine under Win95. Build performance is also good when building from Samba 1.9.18p8 running on a Sun UltraSparc. To give some concrete figures, building

Hi Kurt, I've read the new document. It's a very fine piece of technical documentation. However, it doesn't help me in the present predicament. Going now by this book I'm only getting more frustrated. I'm stuck! For example adding driver with rcpclient invariably ends in "Segmentation fault" (Signal 11) immediately after pop_sec_ctx, whose push counterpart was

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! I have a problem with samba. I run the freebsd samba-des port (samba-des-1.9.17.4) in a pretty much vanilla setting; home dir, a few shared dirs, no printing (yet). The problem is performance. While Samba is not terribly slow it's still too slow. Copying large files takes about half a minute/meg on an Ethernet LAN with no other activity. Trying to read a file using netscape

Some day I may figure out how ggplot2 works. I am trying to plot 5 columns of data on a graph (similar to a simple matplot) =========================================================================== library(ggplot2) bmi <- structure(list(pct = 2:21, P10 = c(14.6, 14.5, 14.2, 13.9, 13.7, 13.7, 13.9, 14.2, 14.5, 14.8, 15.3, 15.9, 16.6, 17.2, 17.8, 18.1, 18.3, 18.4, 18.5, 18.6), P25 =

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,

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

NICs with Linux. > > I think I have tracked down the problem though. even > though it may not > > look like it, I think that the 3com NIC is bad. I put a 10 > baseT only > > Realtek card in a cliant, and it works fine.. Is there any > known problems > > with 3COM 905b's running at 10 Mbits and SAMBA. anyway, I > think I will > > buy Intel NICs

Hi, I use CentOS as a firewall/proxy/webserver/fileserver in my small network. As the small-spec machine with CentOS is heavily loaded/used I can't afford downtime. 20 GB, pentium II with only 128 MB RAM. However I want to know the news on 4.5 is it due soon? Can I gain more by running CentOS 3.x range on such an old machine like mine? Will CentOS 5.0 mean you need a minimum of 512 RAM?

Hi Kurt, I've read the new document. It's a very fine piece of technical documentation. However, it doesn't help me in the present predicament. Going now by this book I'm only getting more frustrated. What are your versions of Samba, CUPS(?), OS, and Printer again? 2.2.8a, 1.1.18, SuSE 8.2 w/2.4.20-62, Lexmark C720 I'm stuck! I would

Hi Kurt, I tried to explain my problem to you at Göttingen but it looks like I'm the only one who gets this error. I've just reinstalled my box with SuSE 8.2, kernel 2.4.20-62, samba 2.2.8a, and cups 1.1.18, moved all config files, databases, scripts and jobs to it and it looks like samba has become even better and faster. No one noticed anything wrong except I just can't come to

hello all, when running ./configure in my [new] R-0.90 directory, the process stop at the following step : (...) checking whether g77 and gcc agree on int and double... can anybody tell me what's the problem (I believe it is due to g77 but what can I do?) thanks, Mathieu I'm using Redhat 6.0 on a Intel P90 -- ----------------------------------------------------------------------

We used syslinux package on IPCop with success to install from floppy, from cdrom, from usb key and pxe. I upgrade the distribution nine month ago from 2.11 to more recent versions (3.11, later 3.20). We add one consistent report of a laptop machine that boot with 2.11 but fail with 3.11 and 3.20 Symptom is that the machine don't display anything after loading the kernel and initrd. with

> ... > I used the latest and greatest SuSE 9.0 Professional... > I then installed all the latest patches via YaST. That > gives me a kernel of 2.4.21 (-144 in SuSE speak) and > Samba 2.2.8a > > I had the configuration backed up on another box, so I > used that as the base for Samba 2.2.8a. I have tried > chmod, chown of various directories, making profile >

Hello everyone, I am working with plm package and I have problem with random and within models, which are giving errors which says "empty model". However, the model is not empty. In the source code for plm.fit, where the error originates it says something like (writing from the top of my head...) X <- model.matrix(formula,data, lhs=1,...) if (ncol(X) == 0) stop("empty