similar to: display p-values and significance levels

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

Hello, I'd like to add custom labels to my pair() plot. These labels include math expression but they aren't correctly displayed... Further, I want that the boxes for the text.panel (diagonal) have an other background color (grey80). Is that generally possible? If yes how do I have to set it? What I've so far is: panel.cor <- function(x, y, digits=2, prefix="",

Hi there! In my case, cor(d[1:20]) makes me a good correlation matrix. Now I'd like to have it one sided, means only the left bottom side to be printed (the others are the same) and I'd like to have * where the p-value is lower than 0.05 and ** lower than 0.01. How can I do this? And another thing: Is there a way to output that table as a latex table? Thanks, Martin -- Ihr

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, 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 list. I was saving my modified pairs graphic using a custom panel from the R Graphics site, and I got an interesting difference in the final image when I save it as eps or png. This custom panel make possible to show at the left side of the pairs plot the p-value as symbol and the correlation r value with its cex proportional to the r value itself. Well, when saving as a png file

I was hoping that there might be some way to attach significance code like the ones from summary.lm to a dataframe. Anyone know how to do something like that. Here is the function i'd like to add that functionality to: add1.coef <- function(model,scope,test="F",p.value=1,order.by.p=FALSE) { num <- length(model$coefficients) add <- add1(model,scope,test=test) sub <-

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

Okay, I stumped. I'm trying to build a matrix or data.frame that consists of values by group. For example, resid, group 574 0.536196370873122 E 575 0.93191093696265 E 576 -1.27116094894117 A 577 -1.05159666055027 A 578 0.398273458301184 D 579 0.586373126102573 E 580 0.806575304513802 E 581 0.897258393118338 E 582 0.342393324256838 D 583 0.781245511773227 E I

Hi all, How do I add an asterisk representing a significant difference to a barplot? Cheers, Gil [[alternative HTML version deleted]]

Hi, I apologize if the solution is right in front of me, but I can't find anything on how to convert a class of 'noquote' to 'matrix'. I've tried as.matrix and I've tried coercing it by 'class(x)<-matrix', but of course that didn't work. I've been using the function 'symnum' which returns an object of class 'noquote'. Here is an

Hi all, Does anyone know a way to change the significant stars in mtable (package memisc)? The default is Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1, however I need it to be Signif. codes: 0 '***' 0.01 '**' 0.05 '*' 0.1 ' ' 1 Kind regards, Erich [[alternative HTML version deleted]]

Hi there, I am just starting in R and this might be a very basic question. I applied one on the examples of pairs() to my own data. The examples shows scatter plots on one side of the matrix and the correlation coefficients on the other which works well. I then modified it slightly because I want different colors in the plots. Now i am getting the error of unused arguments, because the color

R-folks, Is there a function, like aggregate, that allows users to bin values? I've got to break down a data frame into classes of 5cm (or something like it), and I only know how to do it using code like, signif <- symnum( stems$dbh, corr = FALSE, na = FALSE, cutpoints = c(0,10,20,30,40,999), symbols = c(0,10,20,30,40) )

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

I've trying to check for model bias. It was recommended to me that I should break the residuals into n groups and test to see if the mean is different from zero. What is the best way to generate the groups. I've generated a matrix that contains the residuals and groups, but the values all have quotes around them. 1) how do I get rid of the quotes? 2) is this an acceptable method for

Hello, I am trying to export the results from symnum() while maintain their readability. I tried using sink to text file and also copying and pasting but the results end up looking like this: > symnum(c5.s) bC bED bEN bLP bLS bPA bPD bPR p bbContag 1 bbED + 1 bbENN_MN + B 1 bbLPI , , , 1 bbLSI + B B , 1 bbPAFRAC , * * , * 1 bbPD , B B

Hi, I'm trying to create a plot using pair. Currently I'm doing the following pairs(~mpg + hp + wt , data=cars, labels = c("Miles per Gallon","Horsepower","Weight"),pch = c(24,25)[unclass(cars$tr + 1)], bg = c("red", "green3")[unclass(cars$tr + 1)],panel=panel.smooth) for the attached dataset. However,instead of using panel.smooth I