# R/plotCI.R In r4ss: R Code for Stock Synthesis

#### Documented in plotCI

```#' Plot points with confidence intervals.
#'
#' Given a set of x and y values and upper and lower bounds, this function
#' plots the points with error bars. This was Written by Venables and modified
#'
#'
#' @param x The x coordinates of points in the plot
#' @param y The y coordinates of the points in the plot.
#' @param uiw The width of the upper portion of the confidence region.
#' @param liw The width of the lower portion of the confidence region.
#' @param ylo Lower limit of y range.
#' @param yhi Upper limit of y range.
#' @param \dots Additional inputs that will be passed to the function
#' \code{plot(x,y,ylim=ylim,...)}
#' @param sfrac Fraction of width of plot to be used for bar ends.
#' @param ymax Additional input for Upper limit of y range.
#' @param col Color for the points and lines.
#' @author Bill Venables, Ian Stewart, Ian Taylor, John Wallace
#' @keywords aplot hplot
plotCI <-
function (x, y = NULL, uiw, liw = uiw, ylo = NULL, yhi = NULL,
..., sfrac = 0.01, ymax = NULL, add = FALSE, col = "black") {
# Written by Venables; modified for access to ylim, contents, and color
if (is.list(x)) {
y <- x\$y
x <- x\$x
}
if (is.null(y)) {
if (is.null(x)) {
stop("both x and y NULL")
}
y <- as.numeric(x)
x <- seq(along = x)
}
ui <- y + uiw
li <- y - liw
ylim <- range(c(y, ui, li, ylo, yhi))
if(!is.null(ymax)) ylim[2] <- ymax
if (!add) plot(x, y, type='n', ylim = ylim, col= col,...)
segments(x, li, x, ui, col= col)
smidge <- diff(par("usr")[1:2]) * sfrac
x2 <- c(x, x)
ul <- c(li, ui)
segments(x2 - smidge, ul, x2 + smidge, ul, col= col)
points(x, y, col= col, ...)
invisible(list(x = x, y = y))
}
```

## Try the r4ss package in your browser

Any scripts or data that you put into this service are public.

r4ss documentation built on May 29, 2017, 2:24 p.m.