bcaplot uses the output of bcajack, bcajack2, or bcapar to plot bca and standard confidence limits for the parameter of interest.

bcaplot(vl, main = " ", xlab = "coverage", ylab = "limits",
  alpha = c(0.025, 0.05, 0.1, 0.16), ylim, xlim, add = 0,
  sub = "black=bca, green=standard", sw = 1, ...)

Arguments

vl

output of bcajack, bcajack2, or bcapar

main

The main caption (can be empty)

xlab

The x axis label (supplied if not specified)

ylab

The y axis labels (supplied if not specified)

alpha

Coverages are \(1-2\alpha\), e.g. alpha=c(.025,.05) plots intervals [.025,.975] and [.05,.95]. Default is alpha=c(.025,.05,.1,.16) giving coverages .95,.90,.80,.68

ylim

y axis plot limits set automatically if not provided

xlim

x axis plot limits set automatically if not provided

add

add=1 adds a new plot of bca limits (in red) to an existing plot

sub

subtitle (can be empty)

sw

sw=1 draws light vertical dashed lines showing the bca intervals

...

further args for plot

Details

confidence interval endpoints are plotted vertically versus two-sided coverages \(1-2\alpha\). Bca limits in black, Standard limits in green (dashed.). If vl$lims includes the column "jacksd" of jackknife internal standard deviations then these are indicated by vertical red bars centered at the bca limit points.