Re: [R] coloring individual points in lattice xyplot

From: Gabor Grothendieck <ggrothendieck_at_gmail.com>
Date: Thu 06 Jul 2006 - 05:59:08 EST

On 7/5/06, Gabor Grothendieck <ggrothendieck@gmail.com> wrote:
> On 7/5/06, Deepayan Sarkar <deepayan.sarkar@gmail.com> wrote:
> > On 7/4/06, Gabor Grothendieck <ggrothendieck@gmail.com> wrote:
> > > OK. It looks like I need to go to the lower level llines and lpoints to
> > > do this. I wrote a panel routine, mypanel, and it seems to work (see
> > > below); however, currently it assumes types and cols are in the global
> > > environment or at least somewhere where they will be found.
> > >
> > > 1. Is it somehow possible to stick these into some structures set up by
> > > lattice already and then retrieve them from lattice from within mypanel?
> >
> > Not sure what you mean by that, but I would do something like the
> > following (this allows 'col' to be a list as in your original post).
> > Note that arguments to the panel function can be supplied directly to
> > the high level function (that's why graphical parameters can be
> > supplied to xyplot in the first place).
> >
> > mypanel <-
> > function(x, y, subscripts, groups,
> > col = 1,
> > type = "p",
> > ...)
> > {
> > col <- rep(as.list(col), length = nlevels(groups))
> > type <- rep(as.list(type), length = nlevels(groups))
> > for(g in 1:nlevels(groups)) {
> > idx <- g == groups[subscripts]
> > xx <- x[idx]; yy <- y[idx];
> > panel.xyplot(x[idx], y[idx],
> > col = col[[g]],
> > type = type[[g]],
> > ...)
> > }
> > }
> >
> > xyplot(y ~ c(x,x), groups = factor(col(y)),
> > panel = mypanel,
> > type = c("o", "p"),
> > col = list("black", 1:10))

>
>

> Thanks!!! This is quite an improvement.
>

> >
> >
> > > 2. Any other improvements to the example below?
> > >
> > > mypanel <- function(x, y, subscripts, groups, ...) {
> > > for(g in 1:nlevels(groups)) {
> > > idx <- g == groups
> >
> > This won't work for more than one panel.
>

> Could you explain this? Under what situation does it not work?
> The following, for example, appears to work:
>

> xyplot(y ~ c(x,x) | col(y), groups = factor(col(y)),
> panel = mypanel,
> type = c("o", "p"),
> col = list("black", 1:10),
> layout = 1:2)
Just to clarify mypanel refers to your version of mypanel. Were you saying that your version still does not work in certain situations with multiple panels or only mine? > > > >
> >
> > > xx <- x[idx]; yy <- y[idx]; ccols <- cols[subscripts][idx]
> > > if (any(idx)) {
> > > switch(types[g],
> > > p = lpoints(xx, yy, col = ccols),
> > > l = llines(xx, yy, col = ccols),
> > > o = { lpoints(xx, yy, col = ccols)
> > > llines(xx, yy, col = ccols) })
> > > }
> > > }
> > > }
> > >
> > > x <- 1:10
> > > y <- cbind(y1 = x, y2 = x+1)
> > > cols <- c(rep(1,10), 1:10)
> > > types <- c("o", "p")
> > > xyplot(y ~ c(x,x), groups = factor(col(y)), type = types, panel = mypanel)
> > >
> > >
> > >
> > > On 7/4/06, Deepayan Sarkar <deepayan.sarkar@gmail.com> wrote:
> > > > On 7/4/06, Gabor Grothendieck <ggrothendieck@gmail.com> wrote:
> > > > > I can get the types to work or the colors but not both:
> > > >
> > > > Sorry if I wasn't clear, but I didn't mean that you could use
> > > > panel.superpose[.2] to do what you wanted. I only meant that you could
> > > > use it as a template that may help you to write your own panel
> > > > function. What you want is not possible with tools available in
> > > > lattice.
> > > >
> > > > Deepayan
> > > >
> > > > > # this gets the types right but not the colors
> > > > > library(lattice)
> > > > > x <- 1:10
> > > > > y <- cbind(y1 = x, y2 = x+1)
> > > > > cols <- c(rep(1,10), 1:10)
> > > > > xyplot(y ~ c(x,x), groups = col(y), type = c("o", "p"),
> > > > > panel = function(x, y, subscripts, groups, ...)
> > > > > panel.superpose.2(x, y, subscripts, groups, col = cols[subscripts], ...)
> > > > > )
> > > > >
> > > > >
> > > > > # this gets the colors right but not the types
> > > > > library(lattice)
> > > > > x <- 1:10
> > > > > y <- cbind(y1 = x, y2 = x+1)
> > > > > cols <- c(rep(1,10), 1:10)
> > > > > xyplot(y ~ c(x,x), groups = col(y), type = c("o", "p"),
> > > > > panel = function(x, y, subscripts, groups, ...)
> > > > > panel.xyplot(x, y, col = cols[subscripts], ...)
> > > > > )
> > > > >
> > > > >
> > > > > On 7/4/06, Gabor Grothendieck <ggrothendieck@gmail.com> wrote:
> > > > > > On 7/4/06, Deepayan Sarkar <deepayan.sarkar@gmail.com> wrote:
> > > > > > > On 7/4/06, Gabor Grothendieck <ggrothendieck@gmail.com> wrote:
> > > > > > > > If I wish to color groups in xyplot I can do this:
> > > > > > > >
> > > > > > > > library(lattice)
> > > > > > > > x <- 1:10
> > > > > > > > y <- cbind(x, x+1)
> > > > > > > > xyplot(y ~ rep(x,2), group = col(y), col = 1:2)
> > > > > > > >
> > > > > > > > How do I color different points differently within a group.
> > > > > > > >
> > > > > > > > For example, I want to produce this plot (except that I only
> > > > > > > > want to have two groups, not 11):
> > > > > > > >
> > > > > > > > xyplot(y ~ rep(x,2), group = c(rep(1, 10), 2:11), col = 1:11)
> > > > > > > >
> > > > > > > > I am thinking of something like this (although
> > > > > > > > this does not work, its just to get the idea across):
> > > > > > > >
> > > > > > > > xyplot(y ~ rep(x,2), group = col(y), col = list(1, 2:11))
> > > > > > > >
> > > > > > > > where, in general, I have a list with one component per group
> > > > > > > > whose elements are scalars to color the whole group or
> > > > > > > > vectors one color per point in the group. I don't know
> > > > > > > > ahead of time what the list is.
> > > > > > > >
> > > > > > > > I am looking for a general approach to this within the lattice
> > > > > > > > xyplot plot framework; the above is just an example.
> > > > > > >
> > > > > > > The general approach is to write your own panel function. For a
> > > > > > > possible template, look at the functions panel.superpose and
> > > > > > > panel.superpose.2 and how they handle the 'type' argument.
> > > > > > >
> > > > > > > Deepayan
> > > > > > >
> > > > > >
> > > > > > There is no example in ?panel.superpose. Do you think you
> > > > > > could provide an example for the situation in my post?
> > > > > >
> > > > > > I have done quite a bit of RSiteSearch'ing and googling prior to
> > > > > > posting and all the examples I found had colors that depended
> > > > > > on the group, none addressed the situation in my post -- i.e.
> > > > > > coloring individual points within groups.
> > > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > http://www.stat.wisc.edu/~deepayan/
> > > >
> > >
> >
> >
> > --
> > http://www.stat.wisc.edu/~deepayan/
> >
>

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 Received on Thu Jul 06 07:22:15 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 Thu 06 Jul 2006 - 08:18:28 EST.

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