similar to: title for plot contain 4 subplots

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've been generating linear models with lm(). I wanted to look at the VIF's for the coefficient. Using the vif() function from the package Design, I would get unusually high VIF's. However using vif() from the car package I get more reasonable values (ie in line with the quality of the model). What is the difference between the two vif functions? (I dont have access to the

Hi, I had two questions regarding random forests for regression. 1) I have read the original paper by Breiman as well as a paper dicussing an application of random forests and it appears that the one of the nice features of this technique is good predictive ability. However I have some data with which I have generated a linear model using lm(). I can get an RMS error of 0.43 and an R^2 of

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

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

I am trying to run R from an apache C++ module in a shell script to plot some data to display it in apache later. I get the error (reported in apache's logs): Xlib: connection to ":0.0" refused by server Xlib: No protocol specified Error in X11(paste("png::", filename, sep = ""), width, height, pointsize, : unable to start device PNG In addition: Warning

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

Dear R users, When I run the following code, R crashes: require(cclust) x <- matrix(c(0,0,0,1.5,1,-1), ncol=2, byrow=TRUE) cclust(x, centers=x[2:3,], dist="manhattan", method="kmeans") While this works: cclust(x, centers=x[2:3,], dist="euclidean", method="kmeans") I'm posting this here because I am not sure if it is a bug. I've been searching

Hello Martin Maechler and All, A simple question (I hope): How can I compute the "sum of the dissimilarities" that appears in the pam command (from the cluster package) ? Is it the "manhattan" distance (such as the one implemented by "dist") ? I am asking since I am running clustering on a dataset. I found 7 medoids with the pam command, and from it I have the

Dear All, I'm dealing with a case, where 'manhattan' distance of each of 100 vectors is calculated from 10000 other vectors. For achieving this, following 4 scenarios are tested: 1) scenario 1: > x<-read.table("query.vec") > v<-read.table("query.vec2") > d<-matrix(nrow=nrow(v),ncol=nrow(x)) > for (i in 1:nrow(v)){ + d[i,]<-

Hello all. I wish to run k-means with "manhattan" distance. Since this is not supported by the function "kmeans", I turned to the "pam" function in the "fpc" package. Yet, when I tried to have the algorithm run with different starting points, I found that pam ignores and keep on starting the algorithm from the same starting-points (medoids). For my

hi, i am using drb and it is working fine on my local pc on windows.. and after uploading the code to the server it was working for sometime but now it gives me error every time i invoke the process.. error--- uninitialized constant EmailWorker::Notifier - (NameError) /usr/lib64/ruby/gems/1.8/gems/activesupport-1.4.2/lib/active_support/ dependencies.rb:477:in `const_missing''

I use the following code to scan a (limited) parameter space of clustering strategies ... data <- read.table(... dataTranspose <- t(data) distMeth <- c("euclidean", "maximum", "manhattan", "canberra", "binary" ) clustMeth <- c("ward",

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'm making a plot in which the title takes up two lines. The title contains a subscript but when I look at the plot it does'nt seem tocome out properly. The code I'm using is: xtxt = expression(paste('Observed -log( ', IC[50], ' )')) ytxt = expression(paste('Predicted -log( ', IC[50], ' )')) mtxt = expression(paste('Plot of Observed

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>

I'm working on a project related to document clustering. I know that R has clustering algorithms such as clara, but only supports two distance metrics: euclidian and manhattan, which are not very useful for clustering documents. I was wondering how easy it would be to extend the clustering package in R to support other distance metrics, such as cosine distance, or if there was an API for

To whom it may concern: My name is Jillian Weinfeld. I am currently and undergraduate student at New York University and working at Mount Sinai School of Medicine doing research with epilepsy patients. At the moment I am creating a manhattan plot with my data set. After reading many forums and such, I have appropriately plotted my data, however, I wanted to see how I can change the colors of the