[R] two questions re: the use of lattice

David L Carlson dcarlson at tamu.edu
Sun Jul 22 01:50:56 CEST 2012


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