similar to: xyplot with all columns of a data.frame on a single plot

Hi, I have a data.frame that has a categorical variable, for which I would like to look at the distribution of levels of this variable, based on a grouping of two other variables. As an example: x <- data.frame(obs=sample(c('low', 'high'),100, replace=TRUE), grp1=sample(1:10, 100, replace=TRUE), grp2=runif(100)) cut.grp1 <- cut(x$grp1, 3) cut.grp2 <- cut(x$grp2, 3)

Hi, I have a series of lattice plots which I am arranging in a 2x2 grid via print: print(p.preds, split=c(1,1, 2,2), more=TRUE) print(p.comp, split=c(2,1,2,2), more=TRUE) print(p.bw, split=c(1,2,2,2), more=FALSE) What I'd like to have is a letter (A, B, ...) in the top corner of each plot. While panel.text lets me add text anywhere within a plot, I can't seem to workout how I could put

Hi, I have a lattice lot conditioned on two variables. Example code is: library(lattice) x <- data.frame(d=runif(100), f1=sample(c('yes', 'no'),100,replace=TRUE), f2=c(rep('Run1',30),rep('Run2',30),rep('Run3',40))) histogram(~d | f1 + f2, x) In the plot, for a given value of f2, there are two panels, one for 'n'

Hi, while working with decision trees and unbalanced data, I came across the use of the Hellinger distance as an alternative to information gain [1,2], when dealing with skewed data. Does anybody know of R implementations of this approach to decision trees? Thanks, [1] http://www.cse.nd.edu/Reports/2008/TR-2008-06.pdf [2] http://csmr.ca.sandia.gov/~wpk/slides/wdmda-sem.pdf -- Rajarshi Guha NIH

Hi, I have a vector of data, that I group based on two factors via tapply. For each such grouping I would like to plot a pie chart. I can layout these pie charts in a matrix layout, correpsonding to the levels of the two factors. But I am getting stuck on how to label the rows and colums. My current approach looks like this: x <- data.frame(obs=sample(c('low', 'high'),100,

Hi, I've been struggling with a lattice visualiation. I have a data.frame with 4 columns. What I'd like to have is a set of 3 panels. Ecah panel will have the first column plotted against serial number and then will superimpose the relevant column. My non-lattice version is as follows: x <- data.frame( ... ) par(mfrow=c(3,1)) for (i in 2:4) { plot(x[,1]) points(x[,i]) } Any

Hi, I am trying to use a CGI service (Pubchem PUG) via RCurl and am running into a problem where the data must be supplied via POST - but I don't know the keyword for the argument. The data to be sent is an XML fragment. I can do this via the command line using curl: I save the XML string to a file called query.xml and then do curl -d @query.xml

Hi, I'd like to use a custom color sequence (black - low values, green - high values) in am image() plot. While I can specify colors (say a sequence of grays) to the col argument, the ordering is getting messed up. I have two questions: 1. How can I get a sequence of say 256 colors starting from black and ending in green? 2. How is this specified to image() such that it uses the colors in

Hi, I am trying to see if I can use R to perform more rigorous regression analysis. I wonder if the fingerprint package is able to handle pipeline pilot fingerprints (ECFC6 etc) now. Thank you, Eric [[alternative HTML version deleted]]

Hi, I'me trying to write a function that will shuffle a vector. At the moment I'm baically making a vector of randomized indices and then making a new vector from the original one using these random indices. However, is there an alternative (more elegant) method to do this? I tried help.search('shuffle') but it does'nt return anything relevant. Thanks,

Hi, is the following possible - in a given session I make a lot of objects and save when exiting. Usually I note down seperately what each object is about. Is it possible to attach data to any object which would essentially be a short note explaining the meaning of it? Thanks, ------------------------------------------------------------------- Rajarshi Guha <rxg218 at psu.edu>

Hi, I have 3 lists, x, y, z and I'd like to perform a calculation over all the lists simultaneously. If it were a single list I could use lapply, but for more than one list I'm using a for loop. Is there an idiom that would let me use something like lapply, but the function specified to lappy would have access to an element from each list? (In Python, I would have used for a,b,c in

Hi, I have a factor and I would like to find the most frequent level. I think my current approach is a bit long winded and I was wondering if there was a more elegant way to do it: x <- factor(sample(1:0, 5,replace=TRUE)) levels(x)[ which( as.logical((table(x) == max(table(x)))) == TRUE ) ] (The length of x will always be an odd number, so I wont get a tie in max()) Thanks,

Hi, is there a way to plot distribution curves (say normal or chi sq etc) from within R? For example I looked up the *chisq family of functions but I'm not sure as to how I would use them to generate a plot of the chi sq distribution (for arbitrary d.o.f). Thanks, ------------------------------------------------------------------- Rajarshi Guha <rajarshi at presidency.com>

Hi, I am trying to build a series of rlm models. I have my data frame and the models will be built using various coulmns of the data frame. Thus a series of models would be m1 <- rlm(V1 ~ V2 + V3 + V4, data) m2 <- rlm(V1 ~ V2 + V5 + V7, data) m3 <- rlm(V1 ~ V2 + V8 + V9, data) I would like to automate this. Is it possible to use a string in place of the formula? I tried doing: fmla

Hi, this question is probably very obvious but I just cant see where I might be going wrong. I'm using the lm() function to generate a linear model and then make predictions using a different set of data. To generate the model I do (tdata & pdata are matrices of observations and parameters, tdepv, pdepv are response vectors) x <- as.data.frame(tdata) x$tdepv <- tdepv

Hi, I have some data which I was plotting using image(). I wanted to add a vertical color key to the plot and I found that heatmap.2 in gplots does let me add a color key. However, I was thinking of a vertical bar with the color range rather than the style that gplots provides. Is there any package (or code snippet) that would let me add a vertical color key to an image() or heatmap

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi, I'm a new user of R and have managed to make a plot of a histogram. Is there any way I set the title and axes labels and then save the plot as an image (png/gif)? Thanks - -- - ------------------------------------------------------------------- Rajarshi Guha <rajarshi at presidency.com> <http://jijo.cjb.net> GPG

Hi - Many organizations now make their data available as XML via a REST web service architecture. Is there any R package or facility to access this type of data directly (eg, to make the HTTP GET request and have the downloaded data put into an R data frame)? I used several R search sites to look for an answer, but came up with very little. Any help would be appreciated. Thanks very much. Gary

Hi, I'm trying to parse lines of the form: dan001.hin (0): fingerprint={256, 411, 426, 947, 973, 976} What I need is the sequence of number between {}. I'm using grep as match <- grep("{([0-9,\s]*)}",s,perl=T,value=T) where s is a character vector. But all I get is the whole string s. I tried using regexpr in an attempt to get just the sequence I wanted: match <-