mixAK (version 4.2)

cbplot: Plot a function together with its confidence/credible bands

Description

This routine typically plots a function together with its confidence or credible bands. The credible band can be indicated either by additional lines or by a shaded region or by both.

Usage

cbplot(x, y, low, upp, band.type=c("ls", "s", "l"), add=FALSE, col="darkblue", lty=1, lwd=2, cbcol=col, cblty=4, cblwd=lwd, scol=rainbow_hcl(1, start=180, end=180), slwd=5, xlim, ylim, xlab, ylab, main="", sub="", ...)

Arguments

x

a numeric vector with x coordinates corresponding to y, low, upp.

y

a numeric vector with y coordinates of the function to plot.

low

a numeric vector with y coordinates of the lower limit of the credible band.

upp

a numeric vector with y coordinates of the upper limit of the credible band.

band.type

a character which specifies the graphical way to show the credible band, “ls” stands for line and shaded region, “s” stands for shaded region only and “l” stands for line only.

add

if TRUE then everything is added to the current plot.

col, lty, lwd

graphical paramters to draw the x-y line.

cbcol, cblty, cblwd

graphical parameters to draw the x-low and x-upp lines.

scol, slwd

graphical parameters for the shaded region between the credible/confidence bounds.

xlim, ylim, xlab, ylab, main, sub

other graphical parameters.

...

additional arguments passed to the plot function.

Value

invisible(x)

Examples

### Artificial credible bands around the CDF's of N(100, 15*15)
### and N(80, 10*10)
iq <- seq(55, 145, length=100)
Fiq <- pnorm(iq, 100, 15)
low <- Fiq - 0.1
upp <- Fiq + 0.1

iq2 <- seq(35, 125, length=100)
Fiq2 <- pnorm(iq, 80, 10)
low2 <- Fiq2 - 0.1
upp2 <- Fiq2 + 0.1

cbplot(iq, Fiq, low, upp, xlim=c(35, 145))
cbplot(iq2, Fiq2, low2, upp2, add=TRUE, col="red4",
       scol=rainbow_hcl(1, start=20, end=20))