getBRAIDbootstrap {braidrm}

Bootstrap BRAID Parameter Confidence Intervals

Description

Estimate bootstrapped confidence intervals of a BRAID surface fit's parameters, by sampling residuals and refitting.

Usage

getBRAIDbootstrap(bfit, ciLevs = c(0.025, 0.975),numBoot=NULL)

Arguments

bfit

a BRAID surface fit of class 'braidrm'

ciLevs

a two-element vector specifying the lower and upper percentiles of the desired confidence interval. The default value of c(0.025, 0.975) results in a 95% confidence interval.

numBoot

the number of bootstrapped coefficients to be calculated. If NULL (the default), number of bootstrapped coefficients is determined from the width of the interval.

Details

This function constructs bootstrapped confidence intervals for objects of class braidrm. If the input already has bootstrapped confidence intervals (or has tried and failed to generate them), this function will throw a warning and return the input model.

Value

An object of class braidrm, with the new elements:

ciPass

TRUE if a sufficient proportion of bootstrapping trials successfully converged, FALSE otherwise

ciLevs

Two-element vector specifying the lower and upper percentiles of the desired confidence interval. Defaults to 0.025 and 0.975 for a 95% confidence interval.

ciVec

Assuming 'ciPass' is true, a vector containing the lower and upper bounds on the confidence intervals of all free parameters

bCoefs

Array containing the best-fit parameter values from all bootstrapping trials. Useful for calculating confidence intervals on other values.

Author(s)

Nathaniel R. Twarog

See Also

braidrm, calcBRAIDconfint

Examples

data(es8olatmz)

brd <- braidrm(act~conc1+conc2,es8olatmz,getCIs=FALSE,fixed=c(NA,NA,NA,NA,1,NA,NA,-4,-4,-4))
summary(getBRAIDbootstrap(brd))
## Not run: summary(getBRAIDbootstrap(brd,ciLevs=c(0.05,0.95)))