using "by" and indicies
On Mon, 18 Mar 2002, Jason Turner wrote:
I sent this to the list last week, and haven't seen it pop up. Either I deleted it when it did appear, or possibly it was destroyed as spam...? If it did appear and I somehow missed it, appologies.
[snip]
I'm using by() to call pairs(), and would like to put the factor level into the title. A toy example would be something like:
d1 <- data.frame(rnorm(100),rt(100,2),sin(1:100))
levs <- as.factor(rep(1:5,20))
levels(levs) <- c("low","med-lo","med","med-hi","high")
by(d1,levs,pairs)
Is there an elegant way to get "low"..."high" in the main title for each plot, without binding "levs" to the data sent to a custom plot function?
One way is to pass the level info with the data to be plotted:
d2 <- cbind(d1, levs)
by(d2,levs,function(x){ tmp <- levels(x[,4])[x[1,4]] ;
+ pairs( x[,1:3], main=paste("Level:",tmp) ) } )
hope this helps,
Greg Snow, PhD Office: 223A TMCB Department of Statistics Phone: (801) 378-7049 Brigham Young University Dept.: (801) 378-4505 Provo, UT 84602 email: gls at byu.edu -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._