Re: [R] dotplot, dropping unused levels of 'y'

From: Deepayan Sarkar <deepayan.sarkar_at_gmail.com>
Date: Tue 26 Sep 2006 - 18:28:55 GMT

On 9/26/06, Benjamin Tyner <btyner@gmail.com> wrote:
> Deepayan Sarkar wrote:
>
> > On 9/15/06, Benjamin Tyner <btyner@gmail.com> wrote:
> >
> >> In dotplot, what's the best way to suppress the unused levels of 'y' on
> >> a per-panel basis? This is useful for the case that 'y' is a factor
> >> taking perhaps thousands of levels, but for a given panel, only a
> >> handfull of these levels ever present.
> >
> >
> > It's a bit problematic. Basically, you can use
> > relation="free"/"sliced", but y behaves as as.numeric(y) would. So, if
> > the small subset in each panel are always more or less contiguous (in
> > terms of the levels being close to each other) then you would be fine.
> > Otherwise you would not. In that case, you can still write your own
> > prepanel and panel functions, e.g.:
> > -------------
> >
> > library(lattice)
> >
> > y <- factor(sample(1:100), levels = 1:100)
> > x <- 1:100
> > a <- gl(9, 1, 100)
> >
> > dotplot(y ~ x | a)
> >
> > p <-
> > dotplot(y ~ x | a,
> > scales = list(y = list(relation = "free", rot = 0)),
> >
> > prepanel = function(x, y, ...) {
> > yy <- y[, drop = TRUE]
> > list(ylim = levels(yy),
> > yat = sort(unique(as.numeric(yy))))
> > },
> >
> > panel = function(x, y, ...) {
> > yy <- y[, drop = TRUE]
> > panel.dotplot(x, yy, ...)
> > })
> >
> > ----------
> >
> > Hope that gives you what you want.
> >
> > Deepayan
>
> I've been trying to extend this to allow groups, but am running into a
> bit of trouble. For example, the following doesn't quite work: (some of
> the unused factor levels are suppressed per panel, but not all):

I don't think panel = panel.superpose is enough. Try

             panel = function(x, y, ...) {
                yy <- y[, drop = TRUE]
                yy.n <- as.numeric(yy)
                panel.superpose(x, yy.n, ...)
            },
            panel.groups =
            function(x, y, subscripts, duration, col, ...) {
                panel.abline(h = y, col = "lightgray")
                panel.xyplot(x, y, col = col, ...)
                panel.segments(x, y,
                               x + duration[subscripts], y,
                               col = col)
            },

-Deepayan

> set.seed(47905)
> temp3<-data.frame(s_port=factor(rpois(100,10)),
> POSIXtime=structure(1:100,class=c("POSIXt","POSIXct")),
> l_ipn=factor(rpois(100,10)),
> duration=runif(100),
> locality=sample(1:4,replace=TRUE,size=100),
> l_role=sample(c(-1,1),replace=TRUE,size=100))
>
> plot<-dotplot(s_port~POSIXtime|l_ipn,
> data=temp3,
> layout=c(1,1),
> pch="|",
> col=1:8,
> duration=temp3$duration,
> auto.key=list(col=1:8,points=FALSE),
> groups=locality*l_role,
> prepanel = function(x, y, ...) {
> yy <- y[, drop = TRUE]
> list(ylim = levels(yy),
> yat = sort(unique(as.numeric(yy))))
> },
> panel = panel.superpose,
> panel.groups = function(x, y, subscripts, duration, col,
> ...) {
> yy <- y[, drop = TRUE]
> yy.n <- as.numeric(yy)
> panel.abline(h=yy.n,col="lightgray")
> panel.xyplot(x=x,y=yy.n,subscripts=subscripts,col=col,...)
> panel.segments(x,
> yy.n,
> x+duration[subscripts],
> yy.n,
> col = col)
> },
> scales=list(y=list(relation="free"),
> x=list(rot=45)),
> xlab="time",
> ylab="source port")
>
>
>
> Thanks,
> Ben
>



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. Received on Wed Sep 27 04:35:07 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Tue 26 Sep 2006 - 19:30:06 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.