similar to: Need suggestions for plotting 3D plot

Dear Rgurus, To my knowledge the best way to visualize the distribution of a discrete variable X is plot(table(X)) The problem which I have is the following. I have to discrete variables X and Y which distribution I would like to compare. To overlay the distribution of Y with lines(table(Y)) gives not satisfying results. This is the same in case of using density or histogram. Hence, I am

Dear All, I am rather rusty with my R, but I recall being able to do something like plot(x1,...);plot(x2,add=TRUE,...) to plot two series on the same axes. Any suggestions how this may be done will be appreciated.

What should I do if I need to train svm() with data having same value across all rows in some columns. These must be the important features of the class and we cant exclude these columns to build up models. The error I am getting is: Error in predict.svm(ret, xhold) : Model is empty! In addition: Warning message: In svm.default(datatrain, classtrain) : Variable(s) 'F112' and

Hi, hres <- hclust(smatr,method="single") hresd<-as.dendrogram(hres) as.dendrogram(hres) `dendrogram' with 2 branches and 380 members total, at height 2514.513 plot(hresd,leaflab="none") #<-error here. #the plotted dendrogram is incomplete. The x axis is not drawn. #The interested reader can download the save(hresd,file="hres.rda") #from the

Hi, Would like to use to the hexbin package to plot a 2D hist - because it looks really _cool_. My problem is related to drawing a pdf hexbin graphic in series of other graphics. All other figures have a par(mar=c(3.2,3.2,1,1)). So the gplot.hexbin figure in this series looks a little alienated? Was trying a to specify the _mar_ using par, viewport, hexViewport etc. a little. My question. How

Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. This makes the pdf file large. So i tried to save the file first as png and than convert it to pdf. This looks OK if printed but if viewed e.g. with acrobat as document figure the quality is bad. Anyone knows a way to reduce the size but keep the quality? /E -- Dipl. bio-chem. Witold Eryk Wolski MPI-Moleculare

This seems like a bug to me. Can someone verify this? First we define a function f that returns its second argument and lapply it to 1:2 using 9 as the second argument and all seems well. Note that "{" as a function does the same thing as f, as illustrated with f(1,9) and "{"(1,9); however, when we attempt to use "{" in the very same way we used f in lapply,

My data looks like: A,B,C,D,Class 1,2,0,2,cl1 1,5,1,9,cl1 3,2,1,2,cl2 7,2,1,2,cl2 2,2,1,2,cl2 1,2,1,5,cl2 0,2,1,2,cl2 4,2,1,2,cl2 3,5,1,2,cl2 3,2,12,3,cl2 3,2,4,2,cl2 **The steps followed are: trainfile <- read.csv("TrainFile",head=TRUE) datatrain <- subset(trainfile,select=c(-Class)) classtrain <- (subset(trainfile,select=Class)) rf <- randomForest(datatrain, classtrain)

Hi, I clustered a distance matrix and would like to draw it using plot.hclust or plot.dendrogram. The dendrogram is not informative because I have a few extremely small dissimilarities in the distance matrix (e.g. 0), but most of the other distances are in the range 1e10+-5000. I would like to show the tree only for the height of 1e10+-5000 but unfortunately their are no parameter like

I have been struggling with this problem for a while and I hope someone could help me. Or if someone could point me to a section in the manual I would be grateful. x <- "my" plot(1:10, main=expression(paste( x, Delta, "values" ))) Q : How do I get the title to say "my (triangle symbol) values" ? The following trial-and-error produced mainly errors :

Hi I'm reading a CSV file using read.table, and it keeps importing a text variable as a factor. To overcome this, I've used the as.is command referring to the variable in question (called "stim") data<-read.table(file.choose(), header=T, sep=",", as.is = "stim") However, "stim" is still imported as a factor. I notice there are other read.table

Hi, Is the ordering of the values returned something on what I can rely on, a form of a standard, that a function called unique in R (in futher versions) will return the uniq elements in order of they first occurcence. > x<-c(2,2,1,2) > unique(x) [1] 2 1 Its seems not to be the standard. E.g. matlab >> x=[2,2,1,2] x = 2 2 1 2 >> unique(x) ans =

Hi, Please consider the following trivial example. In the sample code I first define a class MVCa using composition. This class contains a slot of type list. Next I define an inheritance relation using the function setIs -> MVCa "is a" list. Then I define a function printlength which prints the list-length which is defined for the type "list". If I call this function

Hey everyone, I have been looking for a while for ways to integrate R's wonderful functions into my C++ software, but I have not found anything concrete. So finally, i post to this list to see if anyouse else knows about this, or has done it!? Is it possible? Are there C++ or C R libraries? Or is it sufficiently easy to build them? your help is much appreciated! thanks doktora

Hi, I am using R to plot the graph and the problem I am facing with my graph is that I have lots of points concentrated in one area and It is creating a visualization challenge. Is there any commands in R I could use to solve this problem. Even if there is no command, do you know how I could tackle this problem...(I want to separate these points so I could see each of them...) Any help is much

Dear List Apologies for such a simple question: I have a vector of 738 elements, coded with values between 1 and 16 (but not containing 7, 10, 11 or 13) and wish to recode value 14 to 1, 4 to 2, 1 to 3, 2 to 4 and all other values to 5. I've been trying to use the replace function (in base) and %in%, but without success. Many thanks Roy

I have a 20x3 matrix as follows: > m <- replicate(3, matrix(rnorm(20),20,1)) I need to compute, say, 95th and 99th percentiles of each column such that the resulting matrix becomes 2x3 with each row representing the respective percentile. My "best effort" is to compute one column at a time as follows: > quantile(m[,1], c(0.95, 0.99)) To do the same for columns 2 and 3, I

Hi, Executing the following code causes the prompt to enter browser() mode. >dd<-as(Npatrick1,"list") Called from: as(Npatrick1, "list") Browse[1]> R : Copyright 2004, The R Foundation for Statistical Computing Version 2.0.0 beta (2004-09-20), ISBN 3-900051-07-0 Some addtitional info about Npatrick1 > is(Npatrick1) [1] "Massvectorlist"

I have a vector with about 100.000 values representing a quite regular function (sinusoid like). I would like to find all local maxima of this function (should be about 4000). Is there a native routine for R? Thanks in advance Eryk. -- _|_ \|/ \|/ Eryk Witold Wolski tel :0049-(0)30-8413-1543 w w ?v? 'v? \'v'/ MPI Moleculare Genetik fax :0049-(0)30-8413-1139 |

I have narrowed down the problem. The error Error : 'AminoAcids' is not an exported object from 'namespace:bibliospec' Error : unable to load R code in package 'bibliospec' occurs only if I try to access the data using bibliospec::AminoAcids within the initialize method of an R reference class. It does work, as far as I tested everywhere else. In other methods of a