getBRAIDbootstrap: Bootstrap BRAID Parameter Confidence Intervals

In nathanieltwarog/braidrm: Fitting Dose Response with the BRAID Combined Action Model

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)))

nathanieltwarog/braidrm documentation built on May 23, 2019, 12:19 p.m.