[R] two questions re: the use of lattice

David L Carlson dcarlson at tamu.edu
Sun Jul 22 01:55:05 CEST 2012


Take a look at useOuterStrips() in package latticeExtra.

-------
David

> -----Original Message-----
> From: David L Carlson [mailto:dcarlson at tamu.edu]
> Sent: Saturday, July 21, 2012 6:51 PM
> To: 'jim holtman'; 'Ranjan Maitra'
> Cc: 'r-help at r-project.org'
> Subject: RE: [R] two questions re: the use of lattice
> 
> Run this before the bwplot() command:
> 
> z$Method <- factor(z$Method, levels = c("BIC", "ICL", "s_v", "Q_v",
>   "sig-q", "s_lsk", "s_lML", "s_mlsk", "s_mlML", "s_la8", "s_haar"))
> 
> I don't have an answer for the 2nd question. Seems like it must be
> possible.
> 
> ----------------------------------------------
> David L Carlson
> Associate Professor of Anthropology
> Texas A&M University
> College Station, TX 77843-4352
> 
> 
> 
> > -----Original Message-----
> > From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> > project.org] On Behalf Of jim holtman
> > Sent: Saturday, July 21, 2012 5:57 PM
> > To: Ranjan Maitra
> > Cc: r-help at r-project.org
> > Subject: Re: [R] two questions re: the use of lattice
> >
> > Answer to you first question, try this at the start of bwplot to
> > specify ordering:
> >
> > bwplot(Error~factor(Method, levels = unique(Method))
> >
> > On Sat, Jul 21, 2012 at 2:42 PM, Ranjan Maitra
> > <maitra.mbox.ignored at inbox.com> wrote:
> > > Dear friends,
> > >
> > > I have two questions regarding the use of lattice. First some code:
> > >
> > >
> > > ## begin code
> > >
> > > z <- cbind(rep(c("BIC", "ICL", "s_v", "Q_v", "sig-q",
> > >                       "s_lsk", "s_lML", "s_mlsk", "s_mlML",
> "s_la8",
> > > "s_haar"), each = 250), rep(c(5, 10, 20, 30, 50), each = 50))
> > > z <- rbind(cbind(z, 0), cbind(z, 20), cbind(z, 40))
> > > z <- cbind(z, rnorm(n = nrow(z)))
> > > z <- as.data.frame(z)
> > > names(z) <- c("Method", "sigma", "INU", "Error")
> > > sigma <- as.numeric(levels(z$sigma))
> > > sigmaExprList <- lapply(sigma, function(s) bquote(italic(sigma) ==
> .
> > > (s))) sigmaExpr <- as.expression(sigmaExprList)
> > > bwplot(Error~Method | sigma, data = z[z[,"INU"] == 0,],scales=list
> > > (rot=90), horiz = F, xlab = "Method", ylab = "Relative Error",
> > >        strip = function(which.given, which.panel, var.name,
> > >          strip.levels = FALSE,
> > >          strip.names = TRUE, ...) {
> > >          strip.default(which.given, which.panel,
> > >                        var.name = sigmaExpr[which.panel],
> > >                        strip.levels = FALSE,
> > >                        strip.names = TRUE, ...)
> > >         },
> > >         layout = c(5,1),  col = "red")
> > >
> > >
> > > ## end code
> > >
> > > Question 1: how do I "force" the display of the "Method" in the
> > > plotting to be in the same order (i.e.,  in the order of "BIC",
> > "ICL",
> > > "s_v", "Q_v", "sig-q", "s_lsk", "s_lML", "s_mlsk", "s_mlML",
> "s_la8",
> > > "s_haar") as the input. As you may notice, it puts them in its own
> > > merry order (I suspect in ascii alphabetical order, but that
> > conjecture
> > > is based entirely on my very few sample attempts).
> > >
> > > Question 2: I want to have 3x5 plots of the respective boxplots.
> > > Something like:
> > > Error ~ Method | sigma + INU?
> > > But I want the labels for the sigma and the INU to be only in the
> > > column and the rows (vertically here) as appropriate, in order to
> > save
> > > plotting space. How do I go about doing this?
> > >
> > > Please reply through the mailing list so that others may also
> > benefit.
> > > In any case, many thanks again for reading and for any help and
> > > pointers!
> > >
> > > Best wishes,
> > > Ranjan
> > >
> > >
> > > --
> > > Important Notice: This mailbox is ignored: e-mails are set to be
> > > deleted on receipt. For those needing to send personal or
> > professional
> > > e-mail, please use appropriate addresses.
> > >
> > > ____________________________________________________________
> > > FREE 3D MARINE AQUARIUM SCREENSAVER - Watch dolphins, sharks &
> orcas
> > on your desktop!
> > >
> > > ______________________________________________
> > > R-help at r-project.org 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.
> >
> >
> >
> > --
> > Jim Holtman
> > Data Munger Guru
> >
> > What is the problem that you are trying to solve?
> > Tell me what you want to do, not how you want to do it.
> >
> > ______________________________________________
> > R-help at r-project.org 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.



More information about the R-help mailing list