[R] forcing levelplot to use relative cuts (ie cuts for each panel)

[Mike Townsley]

Dear guRus,

I'm having trouble producing a levelplot with relative cuts for each 
panel (my data has large differences in scales, so I want to use 
quantiles for each panel).

My attempts to change the 'at'  argument in panel.levelplot function 
have not met with success.

Below is a toy example.

xy - expand.grid(x = 1:3, y = 1:3)

aaa - rbind(cbind(xy, z = 1:9, site = rep('A', 9)),
  cbind(xy, z = (1:9)/10, site = rep('B', 9)),
  cbind(xy, z = (1:9)*10, site = rep('C', 9)))

aaa

library(lattice)
levelplot(z~x+y|site, data = aaa) # using absolute cuts

# now, attempt relative cuts

levelplot(z~x+y|site, data = aaa, panel = function(...) {
   panel.levelplot(at = quantile(z),...) })

I get the following message:
Error in panel.levelplot(at = quantile(z), ...) :
 formal argument at matched by multiple actual arguments

My idea was to determine the cut points each time the panel function 
is called (ie each subset of the data), but I guess this was the 
wrong thing to do.  Can someone point out what I'm missing?

Thanks in advance,

MT



Dr Michael Townsley
Senior Research Fellow
Jill Dando Institute of Crime Science
University College London
Second Floor, Brook House
London, WC1E 7HN

Phone: 020 7679 0820
Fax: 020 7679 0828
Email: [EMAIL PROTECTED]

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

[R] monochrome mosaic plot in vcd package

[Mike Townsley]

helpeRs,

I have a nice looking mosaic plot in an article to be published 
soon.  Sadly, the published version will be in black and white and so ruin 
the advantage of the default shading scheme of tiles.

What would readers suggest as an alternative shading scheme?  If I have a 
black-and-white shading scheme graduated according to suitable cutoffs I 
won't be able to tell positive from negative residuals.  The tile borders 
can be changed of course, but I'm uncertain that is will be clear enough 
for a reader.
Another option may be to use a fill pattern of sloping lines with different 
orientations for the sign and density for the magnitude.  The problem with 
this option is I wouldn't know where to start to incorporate into a legend.

The shading_binary function is no good as I would like the cells with 
residuals less than absolute 2 to be different from other cells.  How would 
readers of this list represent a mosaic plot so that it was easily 
interpretable in monochrome?

My data can be used as an example:


library(vcd)
library(MASS)

term.1 - gl(2,1,8, labels = LETTERS[1:2])
term.2 - gl(2,2,8, labels = LETTERS[3:4])
term.3 - gl(2,4,8, labels = LETTERS[5:6])

cell.count - c(72, 19, 5, 8, 117, 115, 81, 85)

mosaic(loglm(formula = cell.count ~ term.1 + term.2 + term.3),
shade = TRUE, gp = shading_hcl, legend = TRUE,
labeling_args = list(rot_labels = rep(0,4)),
gp_args = list(lty = 1:2),legend_width = unit(0.2, npc))




Dr Michael Townsley
Senior Research Fellow
Jill Dando Institute of Crime Science
University College London
Second Floor, Brook House
London, WC1E 7HN

Phone: 020 7679 0820
Fax: 020 7679 0828
Email: [EMAIL PROTECTED]

__
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