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, ...)
output of bcajack
, bcajack2
, or
bcapar
The main caption (can be empty)
The x axis label (supplied if not specified)
The y axis labels (supplied if not specified)
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
y axis plot limits set automatically if not provided
x axis plot limits set automatically if not provided
add=1
adds a new plot of bca limits (in red) to
an existing plot
subtitle (can be empty)
sw=1
draws light vertical dashed lines showing the
bca intervals
further args for plot
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.