two questions re: the use of lattice

-----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.

orcas

on your desktop!

guide.html

--
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.