Plot Error Bars and Confidence Intervals
Given a set of x and y values and interval width or upper and lower bounds, plot the points with error bars. This can be a useful tool for visualizing confidence intervals.
plotCI(x, y = NULL, uiw, liw = uiw, ui, li, err='y', ylim=NULL, xlim=NULL, type="p", col=par("col"), barcol=col, pt.bg = par("bg"), sfrac = 0.01, gap=1, lwd=par("lwd"), lty=par("lty"), labels=FALSE, add=FALSE, xlab, ylab, minbar, maxbar, ... )
x,y |
coordinates for the center of error bars. |
uiw |
width of the upper or right error bar. Set to |
liw |
width of the lower or left error bar. Defaults to same value as
|
ui |
upper end of error bars. Defaults to |
li |
lower end of error bars. Defaults to |
err |
direction for error bars. Set to "y" for vertical bars. Set to "x" for horizontal bars. Defaults to "y". |
col |
color of plotting character used center marker of error bars. Default is "black". |
xlim, ylim |
range of x/y values to include in the plotting area. |
type |
point/line type; passed to |
barcol |
color of the error bars. Defaults to the same value as
|
pt.bg |
background color of points (use
|
sfrac |
width of "crossbar" at the end of error bar as a fraction of the x plotting region. Defaults to 0.01. |
gap |
space left between the center of the error bar and the lines marking the error bar in units of the height (width) of the letter "O". Defaults to 1.0 |
lwd |
width of bar lines. |
lty |
line type of bar lines. |
labels |
either a logical value indicating whether the circles
representing the x values should be replaced with text giving the
actual values or a vector containing labels to use
instead. Defaults to |
add |
logical indicating whether error bars should be added to
the current plot. If |
minbar |
minumum allowed value for bar ends. If specified,
values smaller than |
maxbar |
maximum allowed value for bar ends. If specified,
values larger than |
... |
optional plotting parameters |
xlab |
label for x axis. |
ylab |
label for y axis. |
Original version by Bill Venables wvenable@attunga.stats.adelaide.edu.au posted to r-help on Sep. 20, 1997. Enhanced version posted to r-help by Ben Bolker ben@zoo.ufl.edu on Apr. 16, 2001. This version was modified and extended by Gregory R. Warnes greg@warnes.net. Additional changes suggested by Martin Maechler maechler@stat.math.ethz.ch integrated on July 29, 2004.
plotmeans
provides an enhanced wrapper to
plotCI
.
# plot means and data(state) tmp <- split(state.area, state.region) means <- sapply(tmp, mean) stdev <- sqrt(sapply(tmp, var)) n <- sapply(tmp,length) ciw <- qt(0.975, n) * stdev / sqrt(n) # plain plotCI(x=means, uiw=ciw) # prettier plotCI(x=means, uiw=ciw, col="black", barcol="blue", lwd=1) # give mean values plotCI(x=means, uiw=ciw, col="black", barcol="blue", labels=round(means,-3), xaxt="n", xlim=c(0,5) ) axis(side=1, at=1:4, labels=names(tmp), cex=0.7) # better yet, just use plotmeans ... # plotmeans( state.area ~ state.region )
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.