similar to: Help with pair plot

The following command produces red axis line in a pairs plot: pairs(iris[1:4], main = "Anderson's Iris Data -- 3 species", pch = "+", col = c("red", "green3", "blue")[unclass(iris$Species)]) Trying to fool pairs in the following way produces the same plot as above: pairs(iris[1:4], main = "Anderson's Iris Data -- 3

I have the following problem. I want to use pairs function and get a matrix of scatterplots with the correlations in the upper panel and the ordinary scatterplots in the lower panel. Moreover, I want to have points colored in five differet ways in the lower panel, because I have five subgroups. In order to do that I tried to combine examples on pairs function help. I got a colored matrix

Hi, Is it possible to create a legend in a pairs() plot? For example, suppose I have the following R codes: data(iris) pairs(iris[1:4], main = "Anderson's Iris Data -- 3 species", pch = 21, bg = c("red", "green3", "blue")[codes(iris$Species)]) How can I created a legend that shows which colour represents to which Iris species? Thanks,

Dear R users, the first 1.0 version of the rPref package is now on CRAN. rPref allows to select the Pareto-optimal tuples from a data set, also called Skylines in the database community. For example, optimal tuples from mtcars according to "high(mpg) * high(hp)" (where "*" is the Pareto operator) are those cars, for which no other dominating car exists. There, a car

> G1=barplot(A,beside=FALSE,space=0.5,col=c("red","gray","green3"),cex.axis=0.8,cex.names=0.8,ylim=c(0,140000),las=1) > legend("topright",c("Civiles","Transaccionan","Demografica"),col=c("red","gray","green3"),pch=15,bty="n",cex=1) >

Dear R-users, I would like to add a legend at the bottom of pairs plots (it's my first use of this function). With the plot function, I usually add some additional space at the bottom when I define the size of the graphical device (using mar); grid functions then allows me to draw my legend as I want. Unfortunatley, this technique does not seem to work with the pairs function as the

Hi there, I got a piece of code for the Iris data which allows to display correlation coefficients for each Iris species in the lower panel (color coded). I would now like to add e.g. a "*" to show the significance of each correlation next to the correlation coefficient. Furthermore I would like to make a t.test between the species "setosa" and "versicolor" for

Carsten, This is a common issue in technology when you have complex products that the majority of consumers of those products do not understand how they operate. If most UPS customers understood the importance of standardization they would have refused to purchase non-standard UPSes and all UPSes would have long ago standardized on a single management protocol. The normal thing is that product

Dear all, I would like to put two perspective plots into one plot. The help page for ?persp shows how one can add points and lines but not another perspective plot. data(volcano) z <- 2 * volcano # Exaggerate the relief x <- 10 * (1:nrow(z)) # 10 meter spacing (S to N) y <- 10 * (1:ncol(z)) # 10 meter spacing (E to W) ## Don't draw the grid lines : border = NA persp(x,

http://www.nabble.com/file/p14664173/at-modality.png I created the above image with R and I have one problem left: Some of the labels of the axes do not show up, probably because there's not enough space. I use the following code to create the plot: modality <- read.table("results.table", header=TRUE, sep=",") color <-

I am trying to make some plots for a poster. Each one needs to be 6 inches by inches. Right now the fonts are too small and the plots don't read well on my poster. I am pretty much an R newbie and I can't figure out how to change the fonts. It would also be nice to know how to specify the size of the plot and maybe how to save the plots as jpeg files. Any help would be much appreciated! I

Hi folks, Win7 64bit R 1.12.0 I run following command on R:- > ToothGrowth > attach(ToothGrowth) > plot(dose,len) > matrics=lm(len~dose) > abline(metrics) Error in int_abline(a = a, b = b, h = h, v = v, untf = untf, ...) : plot.new has not been called yet Only a grey diagram is displayed without content > plot(abline(metrics)) Error in int_abline(a = a, b = b, h = h, v

Hi all, Is there a function to export a dataframe to a text file? I want to store a large set of data which I have saved in a dataframe in my workspace and copy and past doesn't cut it. Thank you, Benjamin

Hi I am fairly new to R and I would appreciate some help to hopefully a trivial problem. I created a function: summary.aggregate <- function(y, ...) { temp.mean <- aggregate(y, FUN=mean, ...) temp.sd <- aggregate(y, FUN=sd, ...) temp.length <- aggregate(y, FUN=length, ...) temp <- data.frame(cbind(mean=temp.mean$x,stdev=temp.sd$x,n=temp.length$x)) } this outputs e.g.:

Hi list, not sure if this is the wanted behavior, but running the following code: > version platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 1.1 year 2005 month 06 day 20 language R > n <- 500 > d <- 4 > m <- matrix(runif(n*d, -1, 1), ncol=d) > c <- hsv(apply(m, 1, function(x) {sum(x*x)/d}),

Hi list, not sure if this is the wanted behavior, but running the following code: > version platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 1.1 year 2005 month 06 day 20 language R > n <- 500 > d <- 4 > m <- matrix(runif(n*d, -1, 1), ncol=d) > c <- hsv(apply(m, 1, function(x) {sum(x*x)/d}),

Hi, I am looking for a way to compare if every element of a vector is > 0. i.e. while(abs(U) > 0) { .. } is there a function for this or do I have to write one? I'd appreciate your help! Benjamin

Hi, is it possible to do this operation faster? I am going over 35k data entries and this takes quite some time. for(cnt in 2:length(sdata$date)) { if(sdata$value[cnt] < sdata$value[cnt - 1]) { sdata$ddtd[cnt] <- sdata$ddtd[cnt - 1] + sdata$value[cnt - 1] - sdata$value[cnt] } else sdata$ddtd[cnt] <- 0 } return(sdata) Thank you, Benjamin

Hello >From the following script I received a grafic output when I called: - xyplot.test( 'green3' ) - call.xyplot.test( 'blue3' ) I did NOT receive a grafic output when I called: - loop.xyplot.test( 'red3' ) What's the Problem? NB: I am using R 1.4.1 on Linux. --------- START OF SCRIPT ---------------- n <- 1000 x <- seq( 1, n ) y <- rnorm( n )

Dear All, [I'd already sent this message earlier (Apr 6) but as I had no response and I still see that it is the case in R1.5 I thought I'd try again] I recently encountered a problem (actually missed some data) because "white" is one of colours in the "default" palette and that does not show up to well on my transparent (on white background). Of course it is