[R] Use of panel.segments

[Deepayan Sarkar]

On Tue, Feb 15, 2011 at 9:44 PM, Lancaster, Vicki
<Vicki.Lancaster at fda.hhs.gov> wrote:
>
> First I have read all the previous post on error bars & segments as well as the  Lattice book by Sarkar.
>
> I am using xyplot, there are 2 panels, each panel displays the measurements from 5 matrices over time.  The matrices are identified by different colors.  My error bars show up in the 1st panel and not the second.  Any suggestions would be greatly appreciated the code follows - thanks.
>
> #Figure 2
> data2$Type<-ordered(data2$Type, c("SDM","AcSDM"))

We don't have data2, so can't run your code.

> sub2<-table(data2$Matrix, data2$Type)
> fill.palette2<-c(rep(palette2, sub2[,1]), rep(palette2, sub2[,2]))
> subscripts2<-1:dim(data2)[1]
>
> #Figure 2
> xyplot(Response~Time|Type, data=data2, upper=log10(data2$Response+data2$SD),
>        lower=log10(data2$Response-data2$SD),
>        groups=fill.palette2, cex=1.35, layout=c(2, 1), type=c("p","g"),        col="black", aspect=1,  col.line="light grey", lty=8, xlim=c(30, 120),
>        scales=list(x=list(at=seq(from=30, to=120, by=20)), y=list(log=10)),
>        ylab="Concentration in Matrices (ppb)",
>        xlab="Time after last dose (hour)",
>        main="Figure 2.  Depletion for Average Concentration by Matrix over Time
>                        (Length of the black vertical line is two standard deviations)",
>        panel=function(x, y, groups, subscripts, upper, lower, ...){
>                panel.xyplot(x, y, pch=21, fill=groups[subscripts2], ...)
>                panel.abline(v=100, col="black", lty=8)
>                panel.text(110, log10(55000), "Post-mortem", cex=0.7)
>                panel.text(110, log10(40000), "Measurements", cex=0.7)
>                switch(panel.number()
>                panel.segments(x[subscripts], lower[subscripts],
>                        x[subscripts], upper[subscripts], lty=1, col="black")},

I don't think you want to subset x (it's already subsetted).

-Deepayan


>
>        key=list(corner=c(0.02, 0.04), text=list(levels(data2$Matrix)),
>        cex=0.75, columns=1,border=TRUE,
>        points=list(pch=c(21, 21, 21, 21, 21), fill=palette2,                                           cex=c(1.5, 1.5, 1.5, 1.5, 1.5)))
>        )
>
> ______________________________________________
> R-help at 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.
>