[R] error bars for a barchart

[ilai]

Hi,
I think the issue is not "respecting the groups" but finding the
x-location of the center of bars in panel.barchart(groups,...). Don't
know about the memisc package, but doesn't look like it provides an
easy solution. This is how I do it:
http://www.mail-archive.com/r-help@r-project.org/msg162299.html

Might need to modify some for your data.

HTH

On Tue, May 1, 2012 at 8:01 AM, Beatriz De Francisco
<Beatriz.DeFrancisco at sams.ac.uk> wrote:
> Hi
> I have the following barchart to which I want to add error bars.
>
>    library(lattice)
>    barchart(Change~fTreat,groups=Process,change,
>              auto.key=list(points=FALSE,rectangles=TRUE),
>              panel=function(x, y,...){
>                panel.barchart(x,y,origin = 0,...);
>                panel.abline(h=0,col="black",...);
>                }
>              )
>
> I have tried using the panel.errbars from the memisc package which works great for xyplots, but when I add it to my code it does not respect the groups.
>
>    library(memisc)
>    barchart(cbind(Change,lower,upper)~fTreat,groups=Process,change,
>         ylab="Pocertage change",
>         ylim=-115:50,
>         scales=list(alternating=FALSE,
>                     tick.number=7,
>                     tck=c(-1,0)),
>         panel=function(x, y,groups,...){
>           panel.barchart(x,y=change$Change,groups=change$Process,origin = 0,...);
>           panel.abline(h=0,col="black",...);
>           panel.errbars(x,y,make.grid="none",ewidth=0.2,type="n",...)
>         }
>         )
>
> Any ideas of how to add error bars to my plot either using the panel.errbars or any other function?
>
> The data:
>
>        structure(list(Treat = structure(c(3L, 4L, 1L, 2L, 3L, 4L, 1L,
>    2L), .Label = c("12-380", "12-750", "8-380", "8-750"), class = "factor"),
>        Process = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("Resp",
>        "Cal"), class = c("ordered", "factor")), Change = c(-33.05,
>        -34.74, 20.94, 18.06, 6.85, -28.57, -8.1, -78.72), upper = c(-13.22896628,
>        -28.61149669, 31.29930461, 27.30173776, 39.73271282, 9.458372948,
>        13.11035572, -47.03745704), lower = c(-52.86120694, -40.87446411,
>        10.57421563, 8.822042178, -26.03144161, -66.60447035, -29.30563327,
>        -110.3973761), fTreat = structure(c(1L, 2L, 3L, 4L, 1L, 2L,
>        3L, 4L), .Label = c("8-380", "8-750", "12-380", "12-750"), class = c("ordered",
>        "factor"))), .Names = c("Treat", "Process", "Change", "upper",
>    "lower", "fTreat"), row.names = c(NA, -8L), class = "data.frame")
>
> Regards
>
> Beatriz de Francisco Mora
> PhD Student
> The Scottish Association for Marine Science
> Scottish Marine Institute
> Oban
> PA37 1QA
> Tel: 06131 559000 (switchboard)
> Fax: 01631559001
> E. beatriz.defrancisco at sams.ac.uk<mailto:beatriz.defrancisco at sams.ac.uk>
> http://www.smi.ac.uk/beatriz-de-franciso
>
> The Scottish Association for Marine Science (SAMS) is registered in Scotland as a Company Limited by Guarantee (SC009292) and is a registered charity (9206). SAMS has an actively trading wholly owned subsidiary company: SAMS Research Services Ltd a Limited Company (SC224404). All Companies in the group are registered in Scotland and share a registered office at Scottish Marine Institute, Oban Argyll PA37 1QA. The content of this message may contain personal views which are not the views of SAMS unless specifically stated. Please note that all email traffic is monitored for purposes of security and spam filtering. As such individual emails may be examined in more detail.
>
>        [[alternative HTML version deleted]]
>
> ______________________________________________
> 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.