semEff (version 0.6.1)

bootCI: Bootstrap Confidence Intervals

Description

Calculate confidence intervals from bootstrapped model effects.

Usage

bootCI(mod, conf = 0.95, type = "bca", digits = 3, bci.arg = NULL, ...)

Value

A summary table of the effects and bootstrapped confidence intervals (data frame), or a list or nested list of same.

Arguments

mod

A fitted model object. Alternatively, a boot object (class "boot"), containing bootstrapped model effects. Can also be a list or nested list of such objects.

conf

A numeric value specifying the confidence level for the intervals.

type

The type of confidence interval to return (defaults to "bca" – see Details). See boot.ci() for further options.

digits

The number of significant digits to return for numeric values.

bci.arg

A named list of any additional arguments to boot.ci(), excepting argument index.

...

Arguments to bootEff().

Details

bootCI() uses boot::boot.ci() to calculate confidence intervals of the specified type and level calculated from bootstrapped model effects. If a model or models is supplied, bootstrapping will first be performed via bootEff().

Nonparametric bias-corrected and accelerated confidence intervals (BCa; Efron, 1987) are calculated by default, which should provide the most accurate coverage across a range of bootstrap sampling distributions (Puth et al., 2015). They will, however, be inappropriate for parametric resampling – in which case the default will be set to the bootstrap percentile method instead ("perc").

Effects and confidence intervals are returned in a summary table, along with the bootstrap standard errors (standard deviations of the samples) and the bootstrap biases (sample means minus original estimates). Effects for which the confidence intervals do not contain zero are highlighted with a star (i.e. 'significant' at the conf level).

References

Chernick, M. R., & Labudde, R. A. (2009). Revisiting Qualms about Bootstrap Confidence Intervals. American Journal of Mathematical and Management Sciences, 29(3–4), 437–456. tools:::Rd_expr_doi("10/c8zv")

Efron, B. (1987). Better Bootstrap Confidence Intervals. Journal of the American Statistical Association, 82(397), 171–185. tools:::Rd_expr_doi("10/gfww2z")

Hesterberg, T. C. (2015). What Teachers Should Know About the Bootstrap: Resampling in the Undergraduate Statistics Curriculum. The American Statistician, 69(4), 371–386. tools:::Rd_expr_doi("10/gd85v5")

Puth, M.-T., Neuhäuser, M., & Ruxton, G. D. (2015). On the variety of methods for calculating confidence intervals by bootstrapping. Journal of Animal Ecology, 84(4), 892–897. tools:::Rd_expr_doi("10/f8n9rq")

Examples

# CIs calculated from bootstrapped SEM
(shipley.sem.ci <- bootCI(shipley.sem.boot))

# From original SEM (models)
# (not typically recommended – better to use saved boot objects)
# system.time(
#   shipley.sem.ci <- bootCI(shipley.sem, R = 1000, seed = 13,
#                            ran.eff = "site")
# )