On 8/1/07, Jonathan Williams [EMAIL PROTECTED] wrote:
I thought one nice addition to a splom figure would be to have the
scatterplots in the upper triangle and a color-coordinated
correlation matrix on the bottom. So I tried my hand at customizing
panel.pairs(), and was rebuffed. Many times. Four hours of
fruitless debugging later, I turn to you for help:
panel.pairs(z=teststatfull[,6:12], pscales=0,
panel.subscripts=FALSE, subscripts=,
upper.panel=lattice.getOption(panel.splom),
lower.panel=function(x1=panel.args$x, y1=panel.args$y,
panel.args=trellis.panelArgs(),
subscripts=1:dim(teststatfull)[1],...){
panel.fill(col=brewer.pal(9,RdBu)[round(cor(x1,y1)*4 + 5)])
panel.text(mean(x1), mean(y1), round(cor(x1,y1),2),
font=2)})
This code is a bit above my level; I stole some tricks from examples
I saw elsewhere, and while it looks over-clunky, it works. Works, at
least, in creating the superpanel: you can try it yourself by
replacing z with a data frame or matrix of your choice and installing
the lattice and RColorBrewer packages. However, when I try to insert
this into the splom function, it all goes to pot. R scolds me for
either missing subscripts, improper subscripts when I try to provide
them, or missing data in the panel function, no matter how I define
one (or don't). Can anyone recommend a solution or show me how to
make my superpanel function more palatable to splom?
I haven't figured out how to run your code, even after replacing 'z',
but I think what you are making things more complicated than they have
to be:
library(lattice)
library(RColorBrewer)
foo - mtcars[c(1:6)]
splom(foo,
upper.panel = panel.splom,
lower.panel = function(x, y, ...) {
panel.fill(col = brewer.pal(9, RdBu)[ round(cor(x, y) * 4 + 5)])
panel.text(mean(x), mean(y), round(cor(x, y),2), font=2)
})
Of course, you can, if you really want to, make a custom superpanel
function to do this:
my.panel.pairs - function(..., lower.panel)
{
my.lower.panel -
function(x, y, ...) {
panel.fill(col = brewer.pal(9, RdBu)[ round(cor(x, y) * 4 + 5)])
panel.text(mean(x), mean(y), round(cor(x, y),2), font=2)
}
panel.pairs(..., lower.panel = my.lower.panel)
}
splom(foo, superpanel = my.panel.pairs)
-Deepayan
__
R-help@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
and provide commented, minimal, self-contained, reproducible code.
