Hi:

I did this in ggplot2, which seemed easier than the approach you tried in
lattice's bwplot() - as far as I can tell, you want to plot the unique value
of v1b as a red dot in each boxplot. To that end,

ex <- data.frame(v1 = log(abs(rt(180, 3)) + 1),
                v2 = rep(c("2007", "2006", "2005"), 60),
                z  = rep(c("a", "b", "c", "d", "e", "f"), e = 30))
# the individual to be marked
ex2 <- data.frame(v1b = log(abs(rt(18, 3)) + 1),
                v2 = rep(c("2007", "2006", "2005"), 6),
                z  = rep(c("a", "b", "c", "d", "e", "f"), e = 3))
ex3 <- merge(ex, ex2, by=c("v2","z"))

library(ggplot2)
h <- ggplot(ex3)
h + geom_boxplot(aes(x = v2, y = v1)) +
    geom_point(aes(x = v2, y = v1b), shape = 17, colour = 'red', size = 3) +
    coord_flip() + facet_wrap( ~ z) +
    opts(panel.minor.grid = theme_blank()) + theme_bw()

It has a different look from bwplot(), but it appears to me that the points
are correctly located.

HTH,
Dennis

On Tue, Sep 28, 2010 at 12:29 AM, Christophe Bouffioux <
christophe.b00@gmail.com> wrote:

> Thanks for your help Peter
> but the red marks on boxplot do not correspond to ex2 dataframe
> actually, it reproduce on each panel the same marks
> that is to say the 3 first lines of ex2
> So this is not correct
>
>
> Christophe
>
> On Mon, Sep 27, 2010 at 6:18 PM, Peter Ehlers <ehlers@ucalgary.ca> wrote:
>
> > On 2010-09-27 4:54, Christophe Bouffioux wrote:
> >
> >> bwplot(v2 ~ v1 | z, data = ex3, layout=c(3,2),
> >>        pch = "|",
> >>        par.settings = list(
> >>        plot.symbol = list(alpha = 1, col = "transparent",cex = 1,pch =
> >> 20)),
> >>        panel = function(x, y){
> >>             panel.bwplot(x, y)
> >>             X<- tapply(ex3$v1b, ex3[,c(1,2)], max)
> >>             Y<- seq(length(unique(ex3[,c(1,2)])))
> >>             panel.points(X, Y, pch = 17, col = "red")
> >>             })
> >>
> >>
> > Perhaps this is what you're trying to achieve:
> >
> >
> >  bwplot(v2 ~ v1 | z, data = ex3, layout=c(3,2),
> >       panel = function(x, y){
> >            panel.bwplot(x, y, pch="|")
> >            X <- tapply(ex3$v1b, ex3[, 1:2], max)
> >            Y <- seq(nrow(unique(ex3[, 1:2])))
> >
> >            panel.points(X, Y, pch = 17, col = "red")
> >            })
> >
> > (I didn't see any need for your par.settings.)
> >
> > I'm not crazy about the way you define X,Y. I think
> > I would augment the data frame appropriately instead.
> >
> >  -Peter Ehlers
> >
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@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.
>

	[[alternative HTML version deleted]]

