[R] pairs: altering pch options on upper and lower panel of pairwise scatter plots

Uwe Ligges ligges at statistik.uni-dortmund.de
Tue Jan 18 09:02:24 CET 2005


Paul Hewson wrote:

> Hello,
> 
> I can't figure out how to use the upper.panel and lower.panel options in
> pairs to alter the label options for either panel independently of the
> other.
> 
> I would like to be able to show the pairwise scatter plots for the data
> as they are (a vanilla pairs plot?) but separately to be able to label
> the points according to a factor level.   It is easy enough to do this
> independently, but I can't figure out how to combine the two settings on
> the same plot.   I have enclosed a toy example:   
> 
> Some made up data, e.g.
> 
> ## make up some multivariate normal data using mvrnorm from MASS
> ##
> require(MASS)
> mu <- c(0,0,0)
> digma <- matrix(c(0.9,0.5,0.5, 0.5,0.9,0.5, 0.5,0.5,0.9)
> dummy.data <- mvrnorm(100,mu,sigma)
> ##
> ## add a factor
> ##
> gender <- as.factor(c(rep("Male",50),rep("Female",50)))
> 
> 
> 
> 
> I can create two separate pairs plots (upper and lower panel) as
> follows:
> 
> ## pairs plot just showing the data
> 
> pairs(dummy.data, lower.panel = NULL)
> 
> ## pairs plot where the points are labelled according to a factor
> 
> pairs(dummy.data, pch = levels(gender), upper.panel = NULL)
> 
> 
> But I would like to be able to combine the two panels on the same plot.
> I would be grateful for any advice.   I've tried passing several
> combinations of functions to upper.panel and lower.panel but either
> they're the wrong function or I'm barking up the wrong tree entirely.


   pairs(dummy.data,
     lower.panel = function(x, y) points(x, y, pch = levels(gender)))

Uwe Ligges



> Thanks
> 
> Paul Hewson
> 
> -=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-
> Lecturer in Statistics
> University of Plymouth
> Drake Circus 
> Plymouth PL4 8AA
> 
> Tel ++ 44 1752 232778
> Email paul.hewson at plymouth.ac.uk
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html




More information about the R-help mailing list