similar to: axis annotation

I saw in the archive a post from Mark Wilkinson (Feb 1, 2003), saying that some of his R-generated png plots came out overlapping. I am seeing the same thing (with R 1.6.2 on Linux i686). My input file generated 4 plot files. The first two were fine, but the last 2 featured a weird overlay of the remaining graphs. The problem is not seen with postscript of pdf output. -- Allin Cottrell

I think I'm missing something tonight in the usage of 'legend': plot(0, type="n") legend(locator(1), month.abb[1:5], pch=15, col=1:5) gives me something similar to what I want. But legend(locator(1), month.abb[1:5], fill=T, col=1:5) gives me 5 black boxes. What am I doing wrong? Thank you, > version _ platform i686-pc-linux-gnu arch i686 os

Hi, I posted a while ago about 'overlap' problems using png/jpeg. If what Patrick Connolly suggests is truly happening, I think the following may be related. My new problem is with the following code (the overlap is still there if I use png() instead of postscript(), compounded by the difficulty described below): tmp <- matrix(runif(16000), nrow=16)

When I execute the following code, it works just like I want it to: three pages of nine (or fewer) plots. However, when I execute the code with the first and last lines uncommented, I get three pages (files), but the 2nd & 3rd pages have overlapping plots. It's like a new page wasn't created. I'm pretty sure I've either misplaced or left out a crucial call to some

Hi, I apologize if this has previously been posted. I've just subscribed to the R-help digest. I'm writing a plotting function that uses layout() to plot several different plots on the same device. This function uses plot(), image(), and a custom function that uses text(). Each cell of the layout needs different par() parameters, so what I'd like to do is pass them as lists:

Hello, it wanted to know how I can extract of a dates frame the values peaks according to an interval that I establish. For example if dates are: 1 23 2 4 3 56 4 7 5 99 6 33 extract the date i wanted to divide into intervals of 2 an d to take alone the numbers 23, 56 and 99 of those 3 intervals. Thanks Ruben

Hi, I want to compute the lower triangle of a square matrix (optionally, sans diagonal). With for() loops I can do something like this: ## 5 by 5 matrix rtn for (j in 1:5) { for (k in 1:j) { if (j != k) { ## optional rtn[j, k] <- my.func(j, k) } } } I'd like to do this with apply(). Is there some way I can do this kind of 'short-circuit'? Thanks, Mark Wilkinson

I have two CSV files (exported from Excel), say file1 and file2. The have the same number of rows, and each has several columns, with names on the first line; and some of the columns in file1 are repeated in file2. Using the "foreign" package, I can read these in separately to dataframes say d1 and d2 with > d1<-read.csv("file1") >

I have been using R and Perl. When I read in a text file using the read.table option, and I try to mathematically manipulate the individual elements in the table, I keep getting an "object is not subsettable" error. If I try to use a different method, it works, but takes too much time(basically, I then need to read in values individually into R instead of as a 2D array, so the number of

Hi, I have been using the prcomp function to perform PCA on my example microarray data, (stored in metric text files) which looks like this: 1a 1b 1c 1d 1e 1f ...................................................4r 4s 4t g1 1.2705 1.2766 ...........................................................2.0298 g2 0.1631

Dear R experts, I am trying to save three plots using tiff graphics devices; however the following code only produces two files (Rplot002.tif and Rplot003.tif) showing figures 1 and 3. Here is a simplified ex code tiff(filename ="Rplot%03d.tif",width=24,height=20,units="cm",res=300, pointsize=10, compression = "lzw") plot(1) mtext("Fig

Hello. I run samba 3.0.21c on RHEL AS 4 EMT_64 U3 and I've experienced a strange problem where smbd will go (and stay) in the D state. While waiting for IO (possibly), these smbd processes do not respond to client requests, thus samba appears to be down. I could not kill these processes via kill -9 , so I rebooted which worked. The only error I could see out of the ordinary would be

Hi, I've just started using R and RSPerl. I have some code as follows: &R::initR("--no-save"); &R::call("t.test", (\@array1, \@array2)); where @array1 and @array2 are both 1-dimensional arrays in Perl having 54675 elements each. On execution the output is as follows: Calling R function name `t.test', # arguments: 3 1) Arg type 3 Got a reference to a

R-Devel-ers: I have an S4 method that simply formats an object: setGeneric("formatMe", function(x) standardGeneric("formatMe")) setMethod("formatMe", "ANY", function(x) format(x)) If I issue the above in an R session, then define an S4 class with its own format method, I get the desired result: > setClass("A",contains="numeric")

Hi, I am using function mars() in package mda to find knots in a whole bunch of predictor variables. I hope to be able to replicate all or some of the basis functions that the MARS software from Salford Systems creates. When I ran mars() on a small dataset, I was able to get the knots. However, when I tried running mars() on a larger dataset (145 predictor variables), for a different

Hi All, I am running quantile regression in a "for loop" starting with 1 variable and adding a variable at a time reaching a maximum of 20 variables. I get the following warning messages after my "for" loop runs. Should I be concerned about these messages? I am building predictive models and am not interested in inference. Warning messages: 1: In

Hi, I am trying to read a SAS version 9.1.3 SAS dataset into R (to preserve the SAS labels), but am unable to do so (I have read in a CSV version). I first created a transport file using the SAS code: libname ces2 'D:\CES Analysis\Data'; filename transp 'D:\CES Analysis\Data\fadata.xpt'; /* create a transport file - R cannot read file created by proc cport */ proc

Here is the website link: https://jobs.stjude.org/css_external/CSSPage_JobDetail.ASP?T=20141205131210& The position is particularly responsible for GPU (CUDA) programming to translate statistical algorithms into R CUDA package, as well as the maintenance of the R CUDA packages. The position calls for willingness to learn and adapt new methods and technologies for parallel and distributed

Since doing "yum update" this morning, an update that installed several packages, the updater shows me 3 packages to install: hpijs-1:3.12.4-4.el6_4.1 (x86_64) hplip-common-3.12.4-4.el6_4.1 (x86_64) hplip-libs-3.12.4-4.el6_4.1 (x86_64) but when I attempt to actually run the update I get: hplip-libs-3.12.4-4.el6.i686 requires libsane.so.1 hplip-libs-3.12.4-4.el6.i686 requires

Hi, I am using the rpart package to build a classification tree. I did manage to build a tree with data on a previous project. However, when attampting to build a tree on a project I am working on, I seem to be getting the error shown below: > nhg3.rp <- rpart(profitresp ~., nhg3, method="class") > plot(nhg3.rp, branch=0.4, uniform=T); text(nhg3.rp, digits=3) Error in