[R] change confidence interval line length in barplot2 (plotrix package)

[David L Carlson]

Maybe a fortunate mistake. If you use the base graphics barplot(), you can
use plotCI() in plotrix to add the confidence intervals with control over
the width of the horizontal ends of the bars (if needed, the defaults are
much narrower):

out <- barplot(hh, beside = TRUE,
   col = c("lightblue", "mistyrose", "lightcyan", "lavender"),
   legend = colnames(VADeaths), ylim = c(0, 20),
   main = "Death Rates in Virginia", font.main = 4,
   sub = "Faked 95 percent error bars", col.sub = mybarcol,
   cex.names = 1.5)
plotCI(out, hh, pch="", gap=0, ui=ci.u, li=ci.l, add=TRUE)

----------------------------------------------
David L Carlson
Associate Professor of Anthropology
Texas A&M University
College Station, TX 77843-4352

> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
> project.org] On Behalf Of Martin Batholdy
> Sent: Tuesday, January 22, 2013 2:42 PM
> To: r-help at r-project.org
> Subject: Re: [R] change confidence interval line length in barplot2
> (plotrix package)
> 
> Ok, I have to apologize,
> I confused the packages.
> 
> It's the function barplot2 from the gplots package!
> 
> 
> >  It calls itself an extenstion of barplot2 and has a ci.lwd argument.
> Might save you the time of doing what I thought might be needed,
> hacking te code.
> 
> Unfortunately ci.lwd controls the thickness of the line but not the
> horizontal width.
> 
> 
> 
> On Jan 22, 2013, at 21:24 , David Winsemius <dwinsemius at comcast.net>
> wrote:
> 
> >
> > On Jan 22, 2013, at 10:28 AM, Martin Batholdy wrote:
> >
> >> Hi,
> >>
> >> is there any way to change the width of the horizontal line of
> confidence intervals
> >> in the barplot2 function in the plotrix package (independent of the
> width of the bars)?
> >>
> >>
> >> example code:
> >>
> >> library(plotrix)
> >> # Example with confidence intervals and grid
> >> hh <- t(VADeaths)[, 1]
> >> mybarcol <- "gray20"
> >> ci.l <- hh * 0.85
> >> ci.u <- hh * 1.15
> >> mp <- barplot2(hh, beside = TRUE,
> >>       col = c("lightblue", "mistyrose",
> >>               "lightcyan", "lavender"),
> >>       legend = colnames(VADeaths), ylim = c(0, 20),
> >>       main = "Death Rates in Virginia", font.main = 4,
> >>       sub = "Faked 95 percent error bars", col.sub = mybarcol,
> >>       cex.names = 1.5, plot.ci = TRUE, ci.l = ci.l, ci.u = ci.u)
> >
> > When I did an sos::findFn("barplot2") search to locate the "real"
> `barplot2` O alos noted in the same package (gplots) a function named
> `ooplot`. It calls itself an extenstion of barplot2 and has a ci.lwd
> argument. Might save you the time of doing what I thought might be
> needed, hacking te code.
> >
> > --
> > David Winsemius
> > Alameda, CA, USA
> >
> 
> ______________________________________________
> 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.