similar to: lattice code to plot columns over another variable

Hi, I have a data.frame with columns named X, D1, D2, D3 I know I can get a single plot with 3 curves by doing xyplot(D1 + D2 + D3 ~ X, data) but in some cases I might have columns D1 ... D10. Is there a way to plot all 10 columns without having to specify each individual term? (By analogy with formulae in lm, I thought, xyplot(. ~ X, data) would work, but it didn't) Thanks, --

Hi, I have a rectangular matrix and I need to check whether any columns are identical or not. Currently I'm looping over the columns and checking each column with all the others with identical(). However, as experience has shown me, getting rid of loops is a good idea :) Would anybody have any suggestions as to how I could do this job more efficiently. (It would be nice to know which

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 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

-----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, 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 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, 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'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, 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

Hi, I'm trying to take a set of rows and columns out of a matrix. I hve been using the index aray approach. My overll matrix is X and is 179 x 65. I want to take out 4 columns and 161 rows. Thus I made a 161 x 2 array I and filled it up with the row,col indices. However doing, X[ I ] gives me a vector of the extracted elements. Is there anyway I can coerce this into a 161 x 4 matrix?

Hi, I'm plotting 4 graphs on one page (2x2 matrix) but I cant seem to get the title for the whole page right. I'm doing: op <- par(mfrow = c(2,2), pty="s") hist(var$V2, breaks="FD",main="Euclidean Metric", xlab="Sum of 3NN ... hist(var$V2, breaks="FD",main="Manhattan Metric", xlab="Sum of 3NN ... hist(var$V2,

Hi, I'm using the ks.test() on two vectors. I looked up the reference and also coded up a version of the two sample Smirnov test. My question is that how can I decide from the output of R that the two vectors x & y come from the same distribution? Am I correct in assuming that smaller D values indicate that they come from the same distribution? In addition how can I use the p value that

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, I'm trying to translate some Matlab code to R and I'm trying to implement the behavior of Matlab's sort() which when applied to a matrix will sort the columns returning the column sorted matrix as well as a matrix of permuted indices. Doing: > x <- matrix(c(3,4,2,6,3,4,8,7,5), nr=3) > x [,1] [,2] [,3] [1,] 3 6 8 [2,] 4 3 7 [3,] 2 4 5

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, is there a function in R that will give me the variances of a predicted values obtained using predict.lm(). If no function is available I would need to calculate them myself - which involves taking the inverse of X'X (' indicating transpose) where X is my model matrix. I know that calculating an inverse directly is not a good idea in general - could anybody suggest a way around

Hi, I'm writing code to generate a plot, in which I draw a series of rectangles. So my code is of the form plot.new() plot.window( ... ) draw rectangle draw rectangle ... Is there a way for me to insert a PNG or PDF graphic at a specific position in the plot (ideally in plot coordinates)? I realize that this might probably be better done in a separate image editor, but if it could

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, 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 <-